Luogu P4180 【模板】严格次小生成树[BJWC2010]
题意
题目描述
小\(C\)最近学了很多最小生成树的算法,\(Prim\)算法、\(Kurskal\)算法、消圈算法等等。正当小\(C\)洋洋得意之时,小\(P\)又来泼小\(C\)冷水了。小\(P\)说,让小\(C\)求出一个无向图的次小生成树,而且这个次小生成树还得是严格次小的,也就是说:如果最小生成树选择的边集是\(E_M\),严格次小生成树选择的边集是\(E_S\),那么需要满足:(\(value(e)\)表示边\(e\)的权值)\(\sum_{e\in E_M}value(e)<\sum_{e\in E_S}value(e)\)
这下小\(C\)蒙了,他找到了你,希望你帮他解决这个问题。
输入输出格式
输入格式:
第一行包含两个整数\(N\)和\(M\),表示无向图的点数与边数。接下来\(M\)行,每行\(3\)个数\(x\ y\ z\)表示,点\(x\)和点\(y\)之间有一条边,边的权值为\(z\)。
输出格式:
包含一行,仅一个数,表示严格次小生成树的边权和。(数据保证必定存在严格次小生成树)
输入输出样例
输入样例#1:
5 6
1 2 1
1 3 2
2 4 3
3 5 4
3 4 3
4 5 6
输出样例#1:
11
说明
数据中无向图无自环;
\(50\%\)的数据\(N\leq 2000,M\leq 3000\);
\(80\%\)的数据\(N≤50 000 M≤100 000\);
\(100\%\)的数据\(N≤100000,M≤300000\),边权值非负且不超过\(10^9\)。
思路
先考虑不严格怎么做。先建出最小生成树,然后对每一条不在树上的边进行考虑:加入这条边之后,树上就会多出来一个环,我们把环上最长的一条边删除,就可以得到另一棵生成树。这样产生的生成树里最小的,就是非严格次小生成树了。
得到环上最长的一条边,只需要在原树上求新加边的两顶点之间的\(LCA\),顺便查询最长边就好了。
之所以不是严格的,是因为加入一条边之后删除的边可能与加入的边的边权相同。这样的话,我们就在预处理倍增记录两点最长边的同时记录两点严格次长边,把它去掉,再加入新边,这样就能求得严格次小生成树了。
AC代码
#include<bits/stdc++.h>
using namespace std;
typedef long long LL;
typedef pair<LL,LL> PLL;
const LL MAXN=1e5+5,MAXM=3e5+5;
LL n,m,ans=LLONG_MAX,sum,dad[MAXN],dep[MAXN];
LL cnt,top[MAXN],to[MAXM<<1],len[MAXM<<1],nex[MAXM<<1];
LL fa[MAXN][20],mx1[MAXN][20],mx2[MAXN][20];
struct Edge
{
LL u,v,d;
bool use;
bool operator < (const Edge &sjf) const {return d<sjf.d;}
}edge[MAXM];
LL read()
{
LL re=0;char ch=getchar();
while(!isdigit(ch)) ch=getchar();
while(isdigit(ch)) re=(re<<3)+(re<<1)+ch-'0',ch=getchar();
return re;
}
void add_edge(LL x,LL y,LL z)
{
to[++cnt]=y,len[cnt]=z,nex[cnt]=top[x],top[x]=cnt;
to[++cnt]=x,len[cnt]=z,nex[cnt]=top[y],top[y]=cnt;
}
LL fd(LL x)
{
LL r=x;
while(r!=dad[r]) r=dad[r];
LL i=x,j;
while(i!=r) j=dad[i],dad[i]=r,i=j;
return r;
}
void Kruskal()
{
for(LL i=1;i<=n;i++) dad[i]=i;
for(LL i=0;i<m;i++)
{
LL fu=fd(edge[i].u),fv=fd(edge[i].v);
if(fu!=fv) dad[fu]=fv,edge[i].use=true,sum+=edge[i].d,add_edge(edge[i].u,edge[i].v,edge[i].d);
}
}
void dfs(LL now)
{
for(LL i=1;i<=18;i++)
{
fa[now][i]=fa[fa[now][i-1]][i-1];
mx1[now][i]=max(mx1[now][i-1],mx1[fa[now][i-1]][i-1]);
if(mx1[now][i-1]>mx1[fa[now][i-1]][i-1]) mx2[now][i]=max(mx2[now][i-1],mx1[fa[now][i-1]][i-1]);
else if(mx1[now][i-1]<mx1[fa[now][i-1]][i-1]) mx2[now][i]=max(mx1[now][i-1],mx2[fa[now][i-1]][i-1]);
else if(mx1[now][i-1]==mx1[fa[now][i-1]][i-1]) mx2[now][i]=max(mx2[now][i-1],mx2[fa[now][i-1]][i-1]);
}
for(LL i=top[now];i;i=nex[i])
{
if(to[i]==fa[now][0]) continue;
fa[to[i]][0]=now,mx1[to[i]][0]=len[i],dep[to[i]]=dep[now]+1;
dfs(to[i]);
}
}
PLL LCA(LL x,LL y)
{
PLL re;re.first=re.second=0;
if(dep[x]<dep[y]) swap(x,y);
for(LL i=18;i>=0;i--)
if(dep[fa[x][i]]>=dep[y])
{
if(mx1[x][i]>re.first) re.second=re.first,re.first=mx1[x][i];
else if(mx1[x][i]>re.second) re.second=mx1[x][i];
else if(mx2[x][i]>re.second) re.second=mx2[x][i];
x=fa[x][i];
}
if(x==y) return re;
for(LL i=18;i>=0;i--)
if(fa[x][i]!=fa[y][i])
{
if(mx1[x][i]>re.first) re.second=re.first,re.first=mx1[x][i];
else if(mx1[x][i]>re.second) re.second=mx1[x][i];
else if(mx2[x][i]>re.second) re.second=mx2[x][i];
x=fa[x][i];
if(mx1[y][i]>re.first) re.second=re.first,re.first=mx1[y][i];
else if(mx1[y][i]>re.second) re.second=mx1[y][i];
else if(mx2[y][i]>re.second) re.second=mx2[y][i];
y=fa[y][i];
}
if(mx1[x][0]>re.first) re.second=re.first,re.first=mx1[x][0];
else if(mx1[x][0]>re.second) re.second=mx1[x][0];
else if(mx2[x][0]>re.second) re.second=mx2[x][0];
x=fa[x][0];
if(mx1[y][0]>re.first) re.second=re.first,re.first=mx1[y][0];
else if(mx1[y][0]>re.second) re.second=mx1[y][0];
else if(mx2[y][0]>re.second) re.second=mx2[y][0];
y=fa[y][0];
return re;
}
int main()
{
n=read(),m=read();
for(LL i=0;i<m;i++) edge[i].u=read(),edge[i].v=read(),edge[i].d=read();
sort(edge,edge+m);
Kruskal();
dep[1]=1;
dfs(1);
for(LL i=0;i<m;i++)
{
if(edge[i].use) continue;
PLL hjj=LCA(edge[i].u,edge[i].v);
if(edge[i].d!=hjj.first) ans=min(ans,sum-hjj.first+edge[i].d);
else if(hjj.second) ans=min(ans,sum-hjj.second+edge[i].d);
}
printf("%lld",ans);
return 0;
}
Luogu P4180 【模板】严格次小生成树[BJWC2010]的更多相关文章
- P4180 【模板】严格次小生成树[BJWC2010]
P4180 [模板]严格次小生成树[BJWC2010] 倍增(LCA)+最小生成树 施工队挖断学校光缆导致断网1天(大雾) 考虑直接枚举不在最小生成树上的边.但是边权可能与最小生成树上的边相等,这样删 ...
- 【洛谷】4180:【模板】严格次小生成树[BJWC2010]【链剖】【线段树维护最大、严格次大值】
P4180 [模板]严格次小生成树[BJWC2010] 题目描述 小C最近学了很多最小生成树的算法,Prim算法.Kurskal算法.消圈算法等等.正当小C洋洋得意之时,小P又来泼小C冷水了.小P说, ...
- 「LuoguP4180」 【模板】严格次小生成树[BJWC2010](倍增 LCA Kruscal
题目描述 小C最近学了很多最小生成树的算法,Prim算法.Kurskal算法.消圈算法等等.正当小C洋洋得意之时,小P又来泼小C冷水了.小P说,让小C求出一个无向图的次小生成树,而且这个次小生成树还得 ...
- 【luogu P4180 严格次小生成树[BJWC2010]】 模板
题目链接:https://www.luogu.org/problemnew/show/P4180 这个题卡树剖.记得开O2. 这个题inf要到1e18. 定理:次小生成树和最小生成树差距只有在一条边上 ...
- 洛谷 P4180 【模板】严格次小生成树[BJWC2010]【次小生成树】
严格次小生成树模板 算法流程: 先用克鲁斯卡尔求最小生成树,然后给这个最小生成树树剖一下,维护边权转点权,维护最大值和严格次大值. 然后枚举没有被选入最小生成树的边,在最小生成树上查一下这条边的两端点 ...
- 【洛谷 P4180】【模板】严格次小生成树[BJWC2010](倍增)
题目链接 题意如题. 这题作为我们KS图论的T4,我直接打了个很暴力的暴力,骗了20分.. 当然,我们KS里的数据范围远不及这题. 这题我debug了整整一个晚上还没debug出来,第二天早上眼前一亮 ...
- P4180 严格次小生成树[BJWC2010] Kruskal,倍增
题目链接\(Click\) \(Here\). 题意就是要求一个图的严格次小生成树.以前被题面吓到了没敢做,写了一下发现并不难. 既然要考虑次小我们就先考虑最小.可以感性理解到一定有一种次小生成树,可 ...
- 【【模板】严格次小生成树[BJWC2010]】
树上的路径怎么能没有树剖 显然,次小生成树和最小生成树只在一条边上有差距,于是我们就可以枚举这一条边,将所有边加入最小生成树,之后再来从这些并不是那么小的生成树中找到那个最小的 我们往最小生成树里加入 ...
- 「BJWC2010」模板严格次小生成树
题目描述 小 \(C\) 最近学了很多最小生成树的算法,\(Prim\) 算法.\(Kruskal\) 算法.消圈算法等等.正当小\(C\)洋洋得意之时,小\(P\)又来泼小\(C\)冷水了.小\(P ...
随机推荐
- 从[id setValue: forKey:]了解KVC
<Objective-C基础教程> P224页有详细介绍 下边是apple官网的简单介绍 和一个应用的例子. KVC就是Key-value coding,大意是允许通过一个Key来读写一个 ...
- fedora23上安装和运行MySQL server (MySQL 已经被MariaDB取代)
[root@localhost kemin]# dnf install mysql-server Fedora 23 - x86_64 - Updates ...
- fasttext源码剖析
fasttext源码剖析 目的:记录结合多方资料以及个人理解的剖析代码: https://heleifz.github.io/14732610572844.html http://www.cnbl ...
- 5.RabbitMQ 客户端控制消息
1.生产者发送消息,消费者结束消息并回执 2.通过channel.basicConsume向服务器发送回执,删除服务上的消息 3.//不向服务器发送回执,服务器的消息一直存在 4.//消费者拒绝接受消 ...
- Java-Class-I:com.alibaba.fastjson.JSONObject
ylbtech-Java-Class-I:com.alibaba.fastjson.JSONObject 1.返回顶部 1.1.import com.alibaba.fastjson.JSON;imp ...
- Spring源码由浅入深系列四 创建BeanFactory
继上一章refresh之后,上图描述了obtainFreshBeanFactory过程.
- 重启集群的时候发现HBase的HRegionServer 服务启动失败
今天在测试环境下的集群重启了下,启动Hbase的时候报错: $ sh start-hbase.sh starting master, logging to /home/hadoop/hbase-0.9 ...
- Apsara Clouder基础技能认证:阿里巴巴编码规范 考试备考题库
考试网址: https://edu.aliyun.com/clouder/exam/intro/33 共50道题 限时90分钟 阿里云大学Apsara Clouder基础技能认证——阿里巴巴编码规范认 ...
- ionic js 加载动画 ionSpinner 提供了许多种旋转加载的动画图标。当你的界面加载时,你就可以呈现给用户相应的加载图标。 该图标采用的是SVG
ionic 加载动画 ion-spinner ionSpinner 提供了许多种旋转加载的动画图标.当你的界面加载时,你就可以呈现给用户相应的加载图标. 该图标采用的是SVG. 用法 <ion- ...
- 【学术篇】SPOJ QTREE 树链剖分
发现链剖这东西好久不写想一遍写对是有难度的.. 果然是熟能生巧吧.. WC的dalao们都回来了 然后就用WC的毒瘤题荼毒了我们一波, 本来想打个T1 44分暴力 然后好像是特判写挂了还是怎么的就只能 ...