首页 > 编程语言 >最小生成树(Kruskal和Prim算法)

最小生成树(Kruskal和Prim算法)

时间:2024-07-24 15:28:50浏览次数:19  
标签:Prim 权重 int Kruskal 最小 生成 算法

最小生成树(Kruskal和Prim算法)

部分资料来源于:最小生成树(Kruskal算法)_kruskal算法求最小生成树-CSDN博客【算法】最小生成树——Prim和Kruskal算法-CSDN博客

关于图的几个概念定义:

连通图:在无向图中,若任意两个顶点vi与vj都有路径相通,则称该无向图为连通图。

强连通图:在有向图中,若任意两个顶点vi与vj都有路径相通,则称该有向图为强连通图。

连通网:在连通图中,若图的边具有一定的意义,每一条边都对应着一个数,称为权;权代表着连接连个顶点的代价,称这种连通图叫做连通网。

生成树:一个连通图的生成树是指一个连通子图,它含有图中全部n个顶点,但只有足以构成一棵树的n-1条边。一颗有n个顶点的生成树有且仅有n-1条边,如果生成树中再添加一条边,则必定成环。

最小生成树:在连通网的所有生成树中,所有边的代价和最小的生成树,称为最小生成树。

img

1.Kruskal算法

此算法可以称为“加边法”,初始最小生成树边数为0,每迭代一次就选择一条满足条件的最小代价边,加入到最小生成树的边集合里。

把图中的所有边按代价从小到大排序;

把图中的n个顶点看成独立的n棵树组成的森林;

按权值从小到大选择边,所选的边连接的两个顶点ui,vi,应属于两颗不同的树,则成为最小生成树的一条边,并将这两颗树合并作为一颗树。

重复(3),直到所有顶点都在一颗树内或者有n-1条边为止。

img

img

parent158778706

例子:

输入

第1行:2个数N,M中间用空格分隔,N为点的数量,M为边的数量。(2 <= N <= 1000, 1 <= M <= 50000) 第2 - M + 1行:每行3个数S E W,分别表示M条边的2个顶点及权值。(1 <= S, E <= N,1 <= W <= 10000)

输出

输出最小生成树的所有边的权值之和。

输入示例

9 14
1 2 4
2 3 8
3 4 7
4 5 9
5 6 10
6 7 2
7 8 1
8 9 7
2 8 11
3 9 2
7 9 6
3 6 4
4 6 14
1 8 8

输出示例

7->8 权重:1
6->7 权重:2
3->9 权重:2
1->2 权重:4
3->6 权重:4
3->4 权重:7
1->3 权重:8
1->5 权重:9
总权重:37
代码示例:
#include <iostream>
#include <algorithm>

using namespace std;
const int MAXN = 5000 + 5;
const int MAXM = 200000 + 5;
const int INF = 0x3fffffff;

struct edge {
	int u;
	int v;
	int w;
} e[MAXM];
int f[MAXN],cnt,m,n,ans;

bool cmp(edge x,edge y)//根据权重排序
{
	return x.w<y.w;
}

int find(int x)//查找并查集
{
	if(f[x]==x) {
		return x;
	} else {
		f[x]=find(f[x]);//路径压缩
		return f[x];
	}
}

void Kruskal()
{
	sort(e+1,e+m+1,cmp);
	for(int i=1; i<=m; i++) {
		int u = find(e[i].u);
		int v = find(e[i].v);
		if(u==v)continue;//判断两个点是否在同一颗树,同一棵树则成环跳过
		ans+=e[i].w;
		printf("%d->%d 权重:%d\n",u,v,e[i].w);
		f[v]=u;//v点的父亲为u,意思为(u,v)这条边加入
		cnt++;
		if(cnt==n-1)break;//所有的点构成构成一棵树
	}
}

int main()
{
	printf("输入点的数目:");
	cin>>n;
	printf("输入边的数目:");
	cin>>m;
	for(int i=1; i<=n; i++) {//初始化,开始时每一个点都在各自的集合
		f[i]=i;
	}
	printf("输入边:\n");
	for(int i=1; i<=m; i++) {
		cin>>e[i].u>>e[i].v>>e[i].w;
	}
	Kruskal();
	printf("总权重:%d",ans);
	return 0;
}

再来个例子:

原题链接:P3366 【模板】最小生成树 - 洛谷 | 计算机科学教育新生态 (luogu.com.cn)

题目描述

如题,给出一个无向图,求出最小生成树,如果该图不连通,则输出 orz

输入格式

第一行包含两个整数 N , M N, M N,M ,表示该图共有 N N N 个结点和 M M M 条无向边。
接下来 M M M 行每行包含三个整数 X i , Y i , Z i X_{i}, Y_{i}, Z_{i} Xi​,Yi​,Zi​ ,表示有一条长度为 Z i Z_{i} Zi​ 的无向边连接结点 X i , Y i X_{i}, Y_{i} Xi​,Yi​ 。

输出格式

如果该图连通,则输出一个整数表示最小生成树的各边的长度之和。如果该图不连通则输出 orz。

输入输出样例

输入 #1

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

输出 #1

7

说明/提示

数据规模:

对于 20 % 20 \% 20% 的数据, N ≤ 5 , M ≤ 20 N \leq 5 , M \leq 20 N≤5,M≤20 。
对于 40 % 40 \% 40% 的数据, N ≤ 50 , M ≤ 2500 N \leq 50 , M \leq 2500 N≤50,M≤2500 。
对于 70 % 70 \% 70% 的数据, N ≤ 500 , M ≤ 1 0 4 N \leq 500, M \leq 10^{4} N≤500,M≤104 。
对于 100 % 100 \% 100% 的数据: 1 ≤ N ≤ 5000 , 1 ≤ M ≤ 2 × 1 0 5 , 1 ≤ Z i ≤ 1 0 4 1 \leq N \leq 5000 , 1 \leq M \leq 2 \times 10^{5} , 1 \leq Z_{i} \leq 10^{4} 1≤N≤5000,1≤M≤2×105,1≤Zi​≤104 。

样例解释:

img

所以最小生成树的总边权为 2+2+3=72+2+3=7。

代码示例:
#include<iostream>
#include<algorithm>
using namespace std;
const int N = 1e6 + 10;
int fa[N];
int sum, u, v, w,n,m;
int tot;
struct node```cpp
#include<iostream>
#include<algorithm>
using namespace std;
const int N = 1e6 + 10;
int fa[N];
int sum, u, v, w,n,m;
int tot;
struct node
{
	int x, y, z;
}e[N];
bool cmp(node a, node b)
{
	return a.z < b.z;
}
void init(int n)
{
	for (int i = 1; i <= n; i++)
	{
		fa[i] = i; // 初始化并查集,每个节点的父节点为自己
	}
		
}
int find(int u)
{
	if (u == fa[u]) // 如果节点的父节点就是自己,则返回该节点
		return u;
	else
		return fa[u] = find(fa[u]); // 否则递归查找父节点,并路径压缩
}
int main()
{
	cin >> n >> m; // 输入节点数和边数
	init(n); // 初始化并查集
	for (int i = 1; i <= m; i++)
	{
		cin >> e[i].x >> e[i].y >> e[i].z; // 输入边的起点、终点和权重
	}
	sort(e + 1, e + m + 1, cmp); // 对边按权重从小到大进行排序
	for (int i = 1; i <= m; i++)
	{
		int q = find(e[i].x); // 查找起点所在集合的代表元素
		int p = find(e[i].y); // 查找终点所在集合的代表元素
		if (q != p) // 如果起点和终点不在同一个集合
		{
			fa[q] = p; // 将起点所在集合合并到终点所在集合
			sum += e[i].z; // 累加权重
			tot++; // 记录边的数量
		}
		else
			continue;
	}
	if (tot == n - 1) // 如果最小生成树的边数为n-1,则输出最小生成树的权重
		cout << sum << endl;
	else
		cout << "orz" << endl; // 否则输出"orz"
	return 0;
}

Prim算法:

Prim算法解释

是采用从点方面考虑来构建MST的一种算法,Prim 算法在稠密图中比Kruskal优,通常步骤如下

1.从源点出发,将所有与源点连接的点加入一个待处理的集合中
2.从集合中找出与源点的边中权重最小的点,从待处理的集合中移除标记为确定的点
3.将找到的点按照步骤1的方式处理
4.重复2,3步直到所有的点都被标记

Prim算法的正确性证明

  1. 源点单独一个点时,可以作为MST
  2. 在与源点相连的点中找出权重最小的加入,显然加入MST后仍然成立
  3. 将所有的点一个一个加入,如果成环,如果新的边权重更小,则替代权重大的,得到的仍然是MST
  4. 综上Prim算法中,MST总不被破坏,所以Prim算法总是正确的

图解

作为坚信图比文字更容易理解,接下来我们将通过例题,更直观的了解Prim算法,给定下图,要我们求出从V1点出发的最小生成树

img

第一步,我们首先将源点V1点加入,并把所有相连的点都加入带处理的集合
惯例的我们将处理完的点标记为红点,待处理的点几位蓝点,如下图

在这里插入图片描述

接下来从待处理的点中找出权重最小的边,将所连的点标记为红色,并且重复第一步得到,在找到(V4,V3)边的时候,因为V4->V3的权重小于V1->V3的权重,所以这里V1->V3这个边可以扔了

在这里插入图片描述

同理对其他的点都进行上诉操作,最后得到总权重和为16

在这里插入图片描述

Prim算法和Dijkstra算法极其相似,特别注意的是:

  • Dijkstra算法的dist数组存的是该点到起点的距离

  • Prim算法的dist数组存的是该点到"树"的集合的距离

就造成了细微的差别,其中该语句:

Prim算法中是这样的:

for (int j = 1; j <= n; j++) // 更新距离数组
       dist[j] = min(dist[j], g[t][j]);

而Dijkstra算法中是这样的:

for (int j = 1; j <= n; j++) // 更新距离数组
      dist[j] = min(dist[j], dis[t]+g[t][j]);

代码示例:

#include<iostream>
#include<algorithm>
#include<cstring>
using namespace std;
const int N = 5100;
int n, m;
int a, b, c;
int dist[N];
const int INF = 0x3f3f3f3f; // 定义一个极大值
bool st[N];

int g[N][N]; // 存储图的邻接矩阵
int prim()
{
    memset(dist, INF, sizeof dist); // 将距离数组初始化为极大值
    int res = 0; // 记录最小生成树的权值和
    for (int i = 0; i < n; i++) // 进行 n 次循环,每次找到一个新节点加入最小生成树
    {
        int t = -1; // t 表示当前距离最小的节点编号
        for (int j = 1; j <= n; j++) // 在所有未加入最小生成树的节点中,找到距离最小的节点
        {
            if(!st[j]&&(t==-1||dist[t]>dist[j]))
                t=j;
        }
        if(i&&dist[t]==INF) // 如果当前加入的节点距离为极大值,说明无法到达该节点,返回极大值表示无解
            return INF;
        if (i)
            res += dist[t]; // 将当前加入节点的距离加入最小生成树的权值和
        for (int j = 1; j <= n; j++) // 更新距离数组
            dist[j] = min(dist[j], g[t][j]);
        st[t] = true; // 标记该节点已加入最小生成树
    }
    return res; // 返回最小生成树的权值和
}
int main()
{
    cin >> n >> m;
    memset(g, INF, sizeof g); // 初始化邻接矩阵为极大值
    while (m--)
    {
        cin >> a >> b >> c;
        g[a][b] = g[b][a] = min(g[a][b], c); // 如果有重边,取最小值
    }
    int t = prim(); // 求最小生成树的权值和
    if (t == INF)
        cout << "orz" << endl;
    else
        cout << t << endl;
    return 0;
}

四.总结

以上就是最小生成树的其中的两个最经常使用的算法了,两种各有优劣,Kruskal算法容易理解,编写容易,但是在稠密图时间效率比Prim算法差,但Prim算法又相对于Kruskal难理解,代码量也要多些。

标签:Prim,权重,int,Kruskal,最小,生成,算法
From: https://blog.csdn.net/Mr____Wangk/article/details/140664441

相关文章