正题

题目链接:https://www.luogu.com.cn/problem/P4180


题目大意

nnn个点mmm条边的一张无向图,求它的严格次小生成树。

1≤n≤105,1≤m≤3×1051\leq n\leq 10^5,1\leq m\leq 3\times 10^51≤n≤105,1≤m≤3×105


解题思路

一定存在一种严格次小生成树和最小生成树只差一条边,感性理解的话大概就是如果有两条不同那么肯定有一条可以替换成另一条要么更优要么不变。

所以我们可以枚举一条不选的边(u,v,w)(u,v,w)(u,v,w)然后找到最小生成树上u,vu,vu,v路径最大的权值kkk替换。

但是发现有可能恰好w=kw=kw=k,所以我们不只需要统计最大的权值,还要统计一个严格次大的,如果等于就选择严格次大的。

时间复杂度O(nlog⁡n)O(n\log n)O(nlogn)


code

#include<cstdio>
#include<cstring>
#include<algorithm>
#define ll long long
using namespace std;
const ll N=1e5+10,T=18;
struct enode{ll x,y,w,v;
}e[3*N];
struct node{ll to,next,w;
}a[N<<1];
ll n,m,tot,ls[N],fa[N],dep[N];
ll ans,sum,f[N][T],g[N][T][2];
bool cmp(enode x,enode y)
{return x.w<y.w;}
ll find(ll x)
{return (fa[x]==x)?x:(fa[x]=find(fa[x]));}
void addl(ll x,ll y,ll w){a[++tot].to=y;a[tot].next=ls[x];ls[x]=tot;a[tot].w=w;return;
}
void dfs(ll x,ll fa){dep[x]=dep[fa]+1;for(ll i=ls[x];i;i=a[i].next){ll y=a[i].to;if(y==fa)continue;f[y][0]=x;g[y][0][0]=a[i].w;dfs(y,x);}return;
}
void calc(ll &mx,ll &mi,ll x){if(x>mx)mi=mx,mx=x;else if(x>mi&&x!=mx)mi=x;return;
}
void calccc(ll x,ll y,ll w){if(x==y)return;if(dep[x]>dep[y])swap(x,y);ll mx=-1,mi=-1;for(ll i=T-1;i>=0;i--)if(dep[f[y][i]]>=dep[x]){calc(mx,mi,g[y][i][0]);calc(mx,mi,g[y][i][1]);y=f[y][i];}if(x!=y){for(ll i=T-1;i>=0;i--)if(f[x][i]!=f[y][i]){calc(mx,mi,g[x][i][0]);calc(mx,mi,g[x][i][1]);calc(mx,mi,g[y][i][0]);calc(mx,mi,g[y][i][1]);x=f[x][i];y=f[y][i];}calc(mx,mi,g[x][0][0]);calc(mx,mi,g[x][0][1]);calc(mx,mi,g[y][0][0]);calc(mx,mi,g[y][0][1]);}if(w!=mx)ans=min(ans,sum+w-mx);else if(mi>=0)ans=min(ans,sum+w-mi);return;
}
signed main()
{scanf("%lld%lld",&n,&m);for(ll i=1;i<=m;i++)scanf("%lld%lld%lld",&e[i].x,&e[i].y,&e[i].w);sort(e+1,e+1+m,cmp);for(ll i=1;i<=n;i++)fa[i]=i;for(ll i=1;i<=m;i++){ll x=find(e[i].x),y=find(e[i].y);if(x==y)continue;addl(e[i].x,e[i].y,e[i].w);addl(e[i].y,e[i].x,e[i].w);fa[x]=y;e[i].v=1;sum+=e[i].w;}memset(g,-1,sizeof(g));dfs(1,1);for(ll j=1;j<T;j++)for(ll i=1;i<=n;i++){f[i][j]=f[f[i][j-1]][j-1];calc(g[i][j][0],g[i][j][1],g[i][j-1][0]);calc(g[i][j][0],g[i][j][1],g[i][j-1][1]);calc(g[i][j][0],g[i][j][1],g[f[i][j-1]][j-1][0]);calc(g[i][j][0],g[i][j][1],g[f[i][j-1]][j-1][1]);}ans=1e18;for(ll i=1;i<=m;i++)if(!e[i].v)calccc(e[i].x,e[i].y,e[i].w);printf("%lld\n",ans);return 0;
}

P4180-[BJWC2010]严格次小生成树【Kruskal,倍增】相关推荐

  1. 解题报告:luogu P4180 [BJWC2010]严格次小生成树(次小生成树、倍增LCA优化、O(mlogn) )

    P4180 [BJWC2010]严格次小生成树 次小生成树有两种,一种是不严格次小生成树,也就是可以数值上等于最小生成树,一种是严格次小生成树,是权值严格大于最小生成树,两种求法大同小异. 方法2在严 ...

  2. 【kruskal】【倍增】严格次小生成树(P4180)

    正题 P4180 题目大意 求严格次小生成树 解题思路 一定存在一种严格次小生成树,和最小生成树只差一条边,不然可以用一条最小生成树上的边代替,从而使边权和更小 那么可以先构造出最小生成树,然后枚举每 ...

  3. BZOJ 1977: [BeiJing2010组队]次小生成树(Kruskal+树上倍增)

    1977: [BeiJing2010组队] 次小生成树 Tree Time Limit: 10 Sec Memory Limit: 512 MB Description 小 C 最近学了很多最小生成树 ...

  4. 模板 - LCA最近公共祖先(倍增法、Tarjan、树上差分、LCA优化的次小生成树)

    整理的算法模板合集: ACM模板 注意x和y的LCA可以是x或者y本身 一.LCA的在线倍增算法 /*给定一棵包含 n个节点的有根无向树,有 m个询问,每个询问 给出了一对节点的编号 x和 y,询问 ...

  5. 最小生成树(kruskal、prim、最小生成森林问题、严格次小生成树)

    整理的算法模板合集: ACM模板 目录 一.kruskal算法 二.prim算法 三.Boruvka算法 四.生成森林问题(K颗树) 五.最小生成树的唯一性 六.严格次小生成树 LCA优化的次小生成树 ...

  6. `Computer-Algorithm` 最小生成树MST,Prim,Kruskal,次小生成树

    Contents 最小生成树 Algorithm Prim Code Kruskal Prim&KruskalPrim \& KruskalPrim&Kruskal算法的性质 ...

  7. 一棵树的生成树有几颗_次小生成树(树剖,生成树)

    生成树的概念: 在一个无向图中,设顶点数为\(n\),取其中\(n-1\)条边并使所有点相连,所得到的一棵树即为生成树. 最小生成树: 如果还没有接触过生成树的同学,欢迎戳->最小生成树详解 次 ...

  8. BZOJ1977: [BeiJing2010组队]次小生成树 Tree

    1977: [BeiJing2010组队]次小生成树 Tree 题意:求严格次小生成树 我为什么要单独发这篇呢 因为愚蠢的我不停换写法最后发现是因为没开long long所以wa掉的 很简单,次小生成 ...

  9. BZOJ 1977 [BeiJing2010组队]次小生成树 Tree

    严格次小生成树.一开始没有特批一圈都相等的情况,一直WA,十分难受. 先生成最小生成树,枚举每条非树边,连上它构成一个环,拆掉环上树边中最大的一条(若和该边相等则次大的一条)换上这条. 用倍增维护一条 ...

最新文章

  1. Datawhale组队学习周报(第047周)
  2. 对于注册中心,ZooKeeper、Eureka 哪个更合适?
  3. python2.7输入函数_Python2.7的用户输入函数有问题,无法让这些输入与程序一起工作...
  4. 1040 有几个PAT
  5. 注:以前我的博客,因为丢了用户名和口令,无法使用,声明作废;现转于此。...
  6. Kaggle入门——使用scikit-learn解决DigitRecognition问题
  7. 【电路仿真】基于matlab BP神经网络三相逆变器故障诊断【含Matlab源码 1655期】
  8. boost电路输出电流公式_boost升压电路参数计算
  9. Google Guice依赖注入框架使用
  10. TL-WAN726N_V2_免驱版Win10开启AP网络
  11. 【NeoVim Coc.nvim】禁用从文档内提取字段的补全选项(禁用带有“yank”的补全选项)
  12. 免手机验证码网站备案
  13. emacs java_Emacs配置与使用
  14. 第一P2P收购中国典当联盟发力供应链金融模式
  15. oracle时分秒修改值_oracle优化(一) oracle数据库使用 TIMESTAMP(6)类型保存年月日时分秒...
  16. mysql级联更新_Mysql实现级联操作(级联更新、级联删除)(转)
  17. 80端口、443端口、8080端口、8000端口的区别
  18. 微信机器人-定制消息
  19. QuartusII-引脚分配与下载
  20. 视频:插入U盘后Win7系统蓝屏,报iusb3xhc.sys

热门文章

  1. java当前时间推前三个月_获取当前时间的前三个月 java
  2. java中string 和stringbuffer的区别_Java中的String,StringBuilder,StringBuffer三者的区别...
  3. linux docker导入镜像,Docker镜像的导入和导出
  4. long 对应oracle,【转】oracle number与java中long、int的对应
  5. java 进程同步代码_java 实现进程间的同步(源代码)
  6. 经纬度 c代码中定义_如何将TXT文本格式的批量经纬度值导入到奥维成为标签
  7. mysql脚本的制作_制作脚本实现mysql自动备份
  8. [MyBatisPlus]常用注解_@TableName_@TableId_@TableField_@TableLogic通过全局配置配置主键生成策略
  9. 数据结构与算法--字符串的排列组合问题
  10. WordList02