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 ...
随机推荐
- pytorch处理模型过拟合
演示代码如下 import torch from torch.autograd import Variable import torch.nn.functional as F import matpl ...
- c++ 枚举简单举例
#include <iostream> enum Enumeration{ VAL1, VAL2, VAL3=100, VAL4 }; int main() { using namespa ...
- 关于合并pdf文件出现的问题
输出端出现以下问题: PdfReadWarning: Xref table not zero-indexed. ID numbers for objects will be 解决方案: import ...
- jmeter接口测试(基础)
一.jmeter创建请求 1.运行jmeter:jmeter.bat 2.右键测试计划,添加一个:线程组 3.右键前面添加的线程组,添加一个:http请求 4.填写请求信息: 5.右键线程组, ...
- Python分布式爬虫必学框架scrapy打造搜索引擎✍✍✍
Python分布式爬虫必学框架scrapy打造搜索引擎 整个课程都看完了,这个课程的分享可以往下看,下面有链接,之前做java开发也做了一些年头,也分享下自己看这个视频的感受,单论单个知识点课程本身 ...
- D3.js(v3)+react 制作 一个带坐标与比例尺的柱形图 (V3版本)
现在用D3.js + react做一个带坐标轴和比例尺的柱形图.我已经尽力把代码全部注释上了,最后我也会把完整柱形图代码奉上.如果还有疑惑的,可以去翻看一下我之前介绍的方法,以下方法都有介绍到. 还有 ...
- 关于UIPageViewController去除边缘点击手势
如果page上方还有一层UI控件的话,不去除边缘点击手势会造成手势的冲突干扰. 首先我做的处理是设置pageView的手势代理 for (UIGestureRecognizer *gr in _pag ...
- 2019-8-31-C#-通过编程的方法在桌面创建回收站快捷方式
title author date CreateTime categories C# 通过编程的方法在桌面创建回收站快捷方式 lindexi 2019-08-31 16:55:58 +0800 201 ...
- Hadoop Pig组件
- nginx代理配置备份
补充: 代理对文件大小的限制, server {client_max_body_size 100M;listen 9096;server_name gate.chaohuoyy.com; locati ...