搜索与图论：Prim

Prim算法求最小生成树

给定一个 n 个点 m 条边的无向图，图中可能存在重边和自环，边权可能为负数。

求最小生成树的树边权重之和，如果最小生成树不存在则输出 impossible。

给定一张边带权的无向图 G = (V,E)，其中 V 表示图中点的集合，

E 表示图中边的集合，n = |V|，m = |E|。

由 V 中的全部 n 个顶点和 E 中 n − 1 条边构成的无向连通子图被称为 G 的一棵生成树，

其中边的权值之和最小的生成树被称为无向图 G 的最小生成树。

输入格式
第一行包含两个整数 n 和 m。

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

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

数据范围
1 ≤ n ≤ 500,
1 ≤ m ≤ 105,
图中涉及边的边权的绝对值均不超过 10000。

输入样例：

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

输出样例：

6

prim 算法干的事情是：给定一个无向图，在图中选择若干条边把图的所有节点连起来。要求边长之和最小。在图论中，叫做求最小生成树。

prim 算法采用的是一种贪心的策略。

每次将离连通部分的最近的点和点对应的边加入的连通部分，连通部分逐渐扩大，最后将整个图连通起来，并且边长之和最小。

我们将图中各个节点用数字 1 ~ n 编号。

要将所有景点连通起来，并且边长之和最小，步骤如下：

1. 用一个 state 数组表示节点是否已经连通。state[i] 为真，表示已经连通，state[i] 为假，表示还没有连通。初始时，state 各个元素为假。即所有点还没有连通。
   用一个 dist 数组保存各个点到连通部分的最短距离，dist[i] 表示 i 节点到连通部分的最短距离。初始时，dist 数组的各个元素为无穷大。
   用一个 pre 数组保存节点的是和谁连通的。pre[i] = k 表示节点 i 和节点 k 之间需要有一条边。初始时，pre 的各个元素置为 -1。
   
2. 从 1 号节点开始扩充连通的部分，所以 1 号节点与连通部分的最短距离为 0，即disti[1] 置为 0。
   
3. 遍历 dist 数组，找到一个还没有连通起来，但是距离连通部分最近的点，假设该节点的编号是 i。i节点就是下一个应该加入连通部分的节点，stata[i] 置为 1。
   用青色点表示还没有连通起来的点，红色点表示连通起来的点。
   这里青色点中距离最小的是 dist[1]，因此 state[1] 置为 1。
   
4. 遍历所有与 i 相连但没有加入到连通部分的点 j，如果 j 距离连通部分的距离大于 i j 之间的距离，即 dist[j] > w[i][j]（w[i][j] 为 i j 节点之间的距离），则更新 dist[j] 为 w[i][j]。这时候表示，j 到连通部分的最短方式是和 i 相连，因此，更新pre[j] = i。
   与节点 1 相连的有 2， 3， 4 号节点。1->2 的距离为 100，小于 dist[2]，dist[2] 更新为 100，pre[2] 更新为1。1->4 的距离为 140，小于 dist[4]，dist[4] 更新为 140，pre[2] 更新为1。1->3 的距离为 150，小于 dist[3]，dist[3] 更新为 150，pre[3] 更新为1。
   
5. 重复 3 4步骤，直到所有节点的状态都被置为 1.
   这里青色点中距离最小的是 dist[2]，因此 state[2] 置为 1。
   
   与节点 2 相连的有 5， 4号节点。2->5 的距离为 80，小于 dist[5]，dist[5] 更新为 80，pre[5] 更新为 2。2->4 的距离为 80，小于 dist[4]，dist[4] 更新为 80，pre[4] 更新为2。
   
   选dist[4]，更新dist[3]，dist[5]，pre[3]，pre[5]。
   
   
   选dist[5]，没有可更新的。
   
   选dist[3]，没有可更新的。
   
6. 此时 dist 数组中保存了各个节点需要修的路长，加起来就是。pre 数组中保存了需要选择的边。
   
   伪代码

int dist[n],state[n],pre[n];
dist[1] = 0;
for(i : 1 ~ n)
{
    t <- 没有连通起来，但是距离连通部分最近的点;
    state[t] = 1;
    更新 dist 和 pre;
}

#include <iostream>
#include <cstring>
#include <algorithm>
using namespace std;
const int N = 510, INF = 0x3f3f3f3f;
int n, m;
int num[N][N];
int dist[N];
bool status[N];

int prim()
{
	memset(dist, 0x3f, sizeof dist);
	int res = 0; // 最小生成树的长度和 
	for (int i = 0; i < n; i ++ ) {
		int t = -1; // 记录当前找到点的下标 
		for (int j = 1; j <= n; j ++ ) 
			// 如果当前点不在集合内，且是第一个点或者距离更小，则记录下标 
			if (!status[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], num[t][j]); // 更新每个点到集合的最短距离 
		/*
		Dijkstra 更新最短距离为点到起点的距离：dist[j] = min(dist[j], dist[t] + num[t][j]); 
		Prim 更新最短距离为点到集合的距离：dist[j] = min(dist[j], num[t][j]); 
		*/
		status[t] = true;
	}
	return res;
}

int main()
{
	cin >> n >> m;
	memset(num, 0x3f, sizeof num);
	while (m -- ) {
		int a, b, c;
		cin >> a >> b >> c;
		num[a][b] = num[b][a] = min(num[a][b], c);
	}
	int t = prim();
	if (t == INF) puts("impossible");
	else cout << t << endl;
	return 0;
}

笔记学习：
作者：Hasity
链接：https://www.acwing.com/solution/content/38312/
来源：AcWing
著作权归作者所有。商业转载请联系作者获得授权，非商业转载请注明出处。

---