蓝桥杯备考：图论之Prim算法

嗯。通过我们前面的学习，我们知道了，一个具有n个顶点的连通图，它的生成树包括n-1个边，如果边多一条就会变成图，少一条就不连通了

接下来我们来学一下把图变成生成树的一个算法

Prim算法，我们从任意一个结点开始构造最小生成树，先把离生成树最近的点拉进来，然后更新所有结点距离生成树的距离，然后再拉进来最近的点，然后再更新所有结点到这颗生成树的距离，一直到所有结点全部拉进来的时候，我们的算法就结束了

这里我们需要开两个数组，一个是dist数组来判断某个点距离生成树的距离，另一个是st数组，来判断某个点在不在生成树里

如图，我们把起始点加入进去了，接下来我们更新一下所有点离最小生成树的距离

我们还是找距离生成树最近的结点，如图，是5，我们把5拉进去

接下来把5标记上，然后更新所有点距离最小生成树的距离

然后继续找最小的，应该是2结点了

把2结点拉进去之后，我们再次更新dist数组，....如此往复，直到所有点都遍历完之后，我们的最小生成树也就构建完了；

话不多说我们来实现一下代码

我们先来实现一下邻接矩阵的代码；

这道题我们还要考虑一下连通不连通的情况

比如这张图，当我们用prim算法把2，1，5，4都加到生成树里的时候，3这个结点距离生成树的值还是无穷，就说明他不连通，没法找最小生成树

下面就是我们的代码

#include <iostream>
#include <cstring>
using namespace std;
const int N = 5010;
const int INF = 0x3f3f3f3f;
int dist[N];
bool st[N];
int edges[N][N];
int n, m;
int ret = 0;
int prim()
{
	dist[1] = 0;
	for (int i = 1; i <= n; i++)//循环把n个结点加入到生成树里面 
	{
		//1.找最近点
		int t = 0;
		for (int j = 1; j <= n; j++)
		{
			if (dist[j] < dist[t] && !st[j])
			{
				t = j;
			}
		}
		if (dist[t] == INF) return INF;//不连通
		//把最小结点加入到生成树里
		st[t] = true;
		ret += dist[t];
		//更新所有结点离最小生成树的距离
		for (int i = 1; i <= n; i++)
		{
			dist[i] = min(dist[i], edges[t][i]);
		}
	}
	return ret;
}
int main()
{
	cin >> n >> m;
	memset(dist, 0x3f, sizeof(dist));
	memset(edges, 0x3f, sizeof(edges));
	for (int i = 1; i <= m; i++)
	{
		int a, b, c; cin >> a >> b >> c;
		edges[a][b] = edges[b][a] = min(edges[a][b],c);
	}
	ret = prim();
	if (ret == INF) cout << "orz" << endl;
	else cout << ret << endl;





	return 0;
}

接下来，我们用邻接表实现一下这道题

#include <iostream>
#include <cstring>
#include <vector>
using namespace std;
typedef pair<int, int> PII;
const int N = 5010, INF = 0x3f3f3f3f;
vector <PII> path[N];
int dist[N];
bool st[N];
int ret = 0;
int n, m;
int prim()
{
	dist[1] = 0;
	for (int i = 1; i <= n; i++)
	{
		int t = 0;
		for (int j = 1; j <= n; j++)
		{
			if (dist[t] > dist[j] && !st[j])
			{
				t = j;
			}
		}
		if (dist[t] == INF) return INF;
		st[t] = true;
		ret += dist[t];
		for (auto &e : path[t])
		{
			int x = e.first, y = e.second;
			dist[x] = min(dist[x], y);
		}
	}
	return ret;
}
int main()
{
	memset(dist, 0x3f, sizeof(dist));
	cin >> n >> m;
	for (int i = 1; i <= m; i++)
	{
		int x, y, z; cin >> x >> y >> z;
		path[x].push_back({ y,z });
		path[y].push_back({ x,z });
	}
	int r = prim();
	if (r == INF) cout << "orz" << endl;
	else cout << r << endl;






	return 0;
}

根据我们的代码可以看到，我们的prim算法时间复杂度是N²，如果结点很少边很多的话，我们用prim算法最合适