【算法】最小生成树-Kruskal

引入

介绍

Kruskal算法的思想比Prin好理解一些。先把边按照权值进行排序，用贪心的思想优先选取权值较小的边，并依次连接，若出现环则跳过此边（用并查集来判断是否存在环）继续搜，直到已经使用的边的数量比总点数少一即可。

证明：如果某个连通图属于最小生成树，那么所有从外部连接到该连通图的边中的一条最短的边必然属于最小生成树。所以不难发现，当最小生成树被拆分成彼此独立的若干个连通分量的时候，所有能够连接任意两个连通分量的边中的一条最短边必然属于最小生成树

代码实现

存图

和普通的图论题不一样，最小生成树不需要链式前向星，直接输入即可。

for(int i=1;i<=m;i++)
        Read(edge[i].u,edge[i].v,edge[i].w);

判环

前置知识：并查集
我们使用并查集判环：
及若两点的祖先一样，则若再连接这两个点，图中会有环。

int find(int x)
{
    while(x!=fa[x]) 
		x=fa[x]=fa[fa[x]];
    return x;
}
 int fx=find(edge[i].u),fy=find(edge[i].v);
 	if(fx==fy)
 		continue;

细节

初始化
这个没什么好说的。

bool cmp(fy a,fy b)
{
    return a.w<b.w;
}
for(int i=1;i<=n;i++)
        fa[i]=i;
sort(edge,edge+m,cmp);

退出条件
当添加的边数=点数-1时，就可以了。
为什么？树的定义~~（QwQ）

if(cnt==n-1)
    	printf("%lld",ans);

例题

题目

给定一个 n 个点 m 条边的无向图，图中可能存在重边和自环，边权可能为负数。求最小生成树的树边权重之和，如果最小生成树不存在则输出 impossible。给定一张边带权的无向图 G=(V,E)，其中 V 表示图中点的集合，E 表示图中边的集合， n=|V|， m=|E|。由 V 中的全部 n 个顶点和 E 中 n−1 条边构成的无向连通子图被称为 G 的一棵生成树，其中边的权值之和最小的生成树被称为无向图 G 的最小生成树。
$1\le n\le 100000$,
$1\le m\le 2\ast 100000$,
图中涉及边的边权的绝对值均不超过1000。

输入格式

• 第一行包含两个整数 n 和 m。

• 接下来 m 行，每行包含三个整数 u,v,w，表示点 u 和点 v 之间存在一条权值为 w 的边。

输出格式

• 共一行，若存在最小生成树，则输出一个整数，表示最小生成树的树边权重之和，如果最小生成树不存在则输出 impossible。

输入/输出例子1
输入：

4 5
1 2 1
1 3 2
1 4 3
2 3 2
3 4 4

输出：

6

样例解释

无

Code

#include<bits/stdc++.h>
#define int long long
using namespace std;
const int N=1e6+1;
inline char gc()
{
  	static char now[1<<20],*S,*T;
  	if(T==S)
  	{
		T=(S=now)+fread(now,1,1<<20,stdin);
    	if(T==S) 
			return EOF;
  	}
  return *S++;
}
template <typename T>
inline void Read(T&x) 
{
  	x=0;
  	char c=gc();
  	while(c<'0'||c>'9') 
  		c=gc();
  	x=c-'0';
  	while((c=gc())>='0'&&c<='9') 
  		x=x*10+c-'0';
}
template <typename T, typename... Args>
inline void Read(T&x,Args&...args)
{
  	Read(x);
  	Read(args...);
}
struct fy
{
	int u,v,w;
}edge[N<<1];
int fa[N],n,m,ans,cnt;
bool cmp(fy a,fy b)
{
    return a.w<b.w;
}
int find(int x)
{
    while(x!=fa[x]) 
		x=fa[x]=fa[fa[x]];
    return x;
}
void kr()
{
    for(int i=1;i<=m;i++)
    {
        int fx=find(edge[i].u),fy=find(edge[i].v);
        if(fx==fy)
            continue;
        ans+=edge[i].w;
        fa[fx]=fy;
        cnt++;
    }
}
signed main()
{
    Read(n,m);
    for(int i=1;i<=n;i++)
        fa[i]=i;
    for(int i=1;i<=m;i++)
        Read(edge[i].u,edge[i].v,edge[i].w);
    sort(edge,edge+m,cmp);
    kr();
    if(cnt==n-1)
    	printf("%lld",ans);
    else
    	puts("impossible");	
    return 0;
}