Kruscal 算法是求最小生成树的基础算法,很容易求得最小生成树,但是怎么利用这个算法求非严格次小生成树与严格次小生成树呢?
Kruskal 算法
做法
这个算法很基础,也很好理解,首先对所有边的权值从小到大排序(既然是要求最小生成树),从边权最小的边开始枚举,对于每条边,如果其两个端点目前不联通,则加入这条边,否则说明加入这条边会成环,就不加。“加入这条边”的操作(也就是合并联通块)用并查集解决。
证明
(其实这个算法 yy 一下就大概知道了……)
Wikipedia 上说:
可以用反证法:如果这条边不在最小生成树中,它连接的两个连通分量最终还是要连起来的,通过其他的连法,那么另一种连法与这条边一定构成了环,而环中一定有一条权值大于这条边的边,用这条边将其替换掉,图仍旧保持连通,但总权值减小了。也就是说,如果不选取这条边,最后构成的生成树的总权值一定不会是最小的。
讲得详细点:
假设按照 Kruskal 算法加入的边集为 ,形成了树 ;
(用反证法)假设这时候求得的 不是该图的最小生成树,则存在最小生成树 ;假设其边集为 ;
那么这两个集合有一个前缀是相同的(可能前缀长度为0),即 ,直到 。
假设现在我们强制把 也加入 ,那么现在一定会形成一个环,这个环里 ,也就是说把 去掉会形成一棵更小的生成树。
说明前面的假设错误,证得这时 是最小生成树。
次小生成树
先介绍一种超级暴力的做法,枚举最小生成树(MST)上的边,把这条边去掉重新跑 MST……
现在怎么用 Kruskal 求次小生成树呢?首先说说非严格次小生成树的做法(所谓非严格次小生成树,就是可能有多个最小生成树的情况下次小生成树可以被认为是最小生成树)。
在 Kruskal 算法里我们可以标记哪些边在 MST 里,那么现在我们可以枚举每条不在 MST 上的边,考虑把这条边加入的情况。
显然,加入这条边会形成环。我们要求的可是次小生成树,可仍然是一棵树。那么在加入这条边后形成的环里,我们肯定要去掉一条边。去掉那条边呢?
假设去掉边 ,加入 ,为了确保我们得到的是次小生成树,我们要使 最小(显然 )。所以对于每个 ,让 尽量大就可以了,换句话说我们要找到加入 所在的环里的最大边,替换之,看看答案。最后找到最接近最小生成树的答案就可以了。
如何实现“加入 所在的环里的最大边”这一神奇的操作呢?用树上倍增(类似 LCA 的做法)就可以了。
再考虑严格次小生成树的做法。我们需要确保 “加入 所在的环里的最大边” 不与 相同,多写一个特殊判断的 DFS 就可以实现了。
参考代码
洛谷模板题:P4180 【模板】严格次小生成树[BJWC2010]
#include<cstdio>
#include<cstring>
#include<iostream>
#include<algorithm>
using namespace std;
const int maxn=100005,maxe=300005;
int n,m,fa[maxn],f[maxn][18],to[maxn][18],deep[maxn];
int tot=0,lnk[maxn],nxt[maxe*2],son[maxe*2],w[maxe*2];
bool vis[maxn];
long long ans=(long long)1<<60,sum=0;
struct SideData{
int x,y,w;
bool used;
}a[maxe];
inline int read(){
int ret=0,f=1;char ch=getchar();
while (ch<'0'||ch>'9') {if (ch=='-') f=-1;ch=getchar();}
while (ch>='0'&&ch<='9') ret=ret*+ch-,ch=();
ret*f;
}
{
aa.w<bb.w;
}
{
tot++;son[tot]=y;w[tot]=z;nxt[tot]=lnk[x];lnk[x]=tot;
}
{
(fa[x]==x) x;
fa[x]=(fa[x]);
fa[x];
}
{
( i=;i<=n;i++) fa[i]=i;
cnt=;
( i=;i<=m&&cnt<n;i++){
fx=(a[i].x),fy=(a[i].y);
(fx==fy) ;
cnt++;a[i].used=;sum+=( )a[i].w;
fa[fx]=fy;
(a[i].x,a[i].y,a[i].w);
(a[i].y,a[i].x,a[i].w);
}
}
{
vis[x]=;
f[x][]=w[id];to[x][]=lst;
( i=;i<;i++) f[x][i]=(f[x][i],f[to[x][i]][i]),to[x][i]=to[to[x][i]][i];
( i=lnk[x];i;i=nxt[i]) (!vis[son[i]]) deep[son[i]]=deep[x],(son[i],i,x);
}
{
ret=;
(deep[x]<deep[y]) (x,y);
( i=;i>=;i--) (deep[to[x][i]]>=deep[y]) ret=(ret,f[x][i]),x=to[x][i];
( i=;i>=;i--) (to[x][i]!=to[y][i]) ret=(ret,(f[x][i],f[y][i])),x=to[x][i],y=to[y][i];
ret=(ret,(f[x][],f[y][]));
ret;
}
{
ret=;
(deep[x]<deep[y]) (x,y);
( i=;i>=;i--) (deep[to[x][i]]>=deep[y]){
(f[x][i]!=z) ret=(ret,f[x][i]);
x=to[x][i];
}
( i=;i>=;i--) (to[x][i]!=to[y][i]){
(f[x][i]!=z) ret=(ret,f[x][i]);
(f[y][i]!=z) ret=(ret,f[y][i]);
x=to[x][i];y=to[y][i];
}
(f[x][]!=z) ret=(ret,f[x][]);
(f[y][]!=z) ret=(ret,f[y][]);
ret;
}
{
n=();m=();
( i=;i<=m;i++) a[i].x=(),a[i].y=(),a[i].w=();
(a,a+m,cmp);
();
deep[]=;
(,,);
( i=;i<=m;i++) (!a[i].used){
qry=(a[i].x,a[i].y);
(qry==a[i].w) qry=(a[i].x,a[i].y,qry);
now=sum-( )qry+( )a[i].w;
(now>sum) ans=(ans,now);
}
(,ans);
;
}
参考
Kruskal's algorithm - Wikipedia
(维基百科英文版解禁了吗?!不用魔法居然可以访问……)
克鲁斯克尔算法 - 维基百科,自由的百科全书