$\color{green}{<–画图不易,点下这里赞一下吧}$
prim 算法干的事情是:给定一个无向图,在图中选择若干条边把图的所有节点连起来。要求边长之和最小。在图论中,叫做求最小生成树。
prim 算法采用的是一种贪心的策略。
每次将离连通部分的最近的点和点对应的边加入的连通部分,连通部分逐渐扩大,最后将整个图连通起来,并且边长之和最小。
我们将图中各个节点用数字 1 ~ n 编号。
要将所有景点连通起来,并且边长之和最小,步骤如下:
-
用一个 state 数组表示节点是否已经连通。state[i] 为真,表示已经连通,state[i] 为假,表示还没有连通。初始时,state 各个元素为假。即所有点还没有连通。
用一个 dist 数组保存各个点到连通部分的最短距离,dist[i] 表示 i 节点到连通部分的最短距离。初始时,dist 数组的各个元素为无穷大。
用一个 pre 数组保存节点的是和谁连通的。pre[i] = k 表示节点 i 和节点 k 之间需要有一条边。初始时,pre 的各个元素置为 -1。
-
从 1 号节点开始扩充连通的部分,所以 1 号节点与连通部分的最短距离为 0,即disti[1] 置为 0。
-
遍历 dist 数组,找到一个还没有连通起来,但是距离连通部分最近的点,假设该节点的编号是 i。i节点就是下一个应该加入连通部分的节点,stata[i] 置为 1。
用青色点表示还没有连通起来的点,红色点表示连通起来的点。
这里青色点中距离最小的是 dist[1],因此 state[1] 置为 1。
-
遍历所有与 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。
-
重复 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],没有可更新的。
-
此时 dist 数组中保存了各个节点需要修的路长,加起来就是。pre 数组中保存了需要选择的边。
伪代码
int dist[n],state[n],pre[n];
dist[1] = 0;
for(i : 1 ~ n)
{
t <- 没有连通起来,但是距离连通部分最近的点;
state[t] = 1;
更新 dist 和 pre;
}
代码
//2022.6.1 更新
#include <iostream>
#include <cstring>
#include <algorithm>
using namespace std;
const int N = 510;
int g[N][N];//存储图
int dt[N];//存储各个节点到生成树的距离
int st[N];//节点是否被加入到生成树中
int pre[N];//节点的前去节点
int n, m;//n 个节点,m 条边
void prim()
{
memset(dt,0x3f, sizeof(dt));//初始化距离数组为一个很大的数(10亿左右)
int res= 0;
dt[1] = 0;//从 1 号节点开始生成
for(int i = 0; i < n; i++)//每次循环选出一个点加入到生成树
{
int t = -1;
for(int j = 1; j <= n; j++)//每个节点一次判断
{
if(!st[j] && (t == -1 || dt[j] < dt[t]))//如果没有在树中,且到树的距离最短,则选择该点
t = j;
}
//2022.6.1 发现测试用例加强后,需要判断孤立点了
//如果孤立点,直返输出不能,然后退出
if(dt[t] == 0x3f3f3f3f) {
cout << "impossible";
return;
}
st[t] = 1;// 选择该点
res += dt[t];
for(int i = 1; i <= n; i++)//更新生成树外的点到生成树的距离
{
if(dt[i] > g[t][i] && !st[i])//从 t 到节点 i 的距离小于原来距离,则更新。
{
dt[i] = g[t][i];//更新距离
pre[i] = t;//从 t 到 i 的距离更短,i 的前驱变为 t.
}
}
}
cout << res;
}
void getPath()//输出各个边
{
for(int i = n; i > 1; i--)//n 个节点,所以有 n-1 条边。
{
cout << i <<" " << pre[i] << " "<< endl;// i 是节点编号,pre[i] 是 i 节点的前驱节点。他们构成一条边。
}
}
int main()
{
memset(g, 0x3f, sizeof(g));//各个点之间的距离初始化成很大的数
cin >> n >> m;//输入节点数和边数
while(m --)
{
int a, b, w;
cin >> a >> b >> w;//输出边的两个顶点和权重
g[a][b] = g[b][a] = min(g[a][b],w);//存储权重
}
prim();//求最下生成树
//getPath();//输出路径
return 0;
}
优化
上面代码的时间复杂度为 O(n^2)。
与Dijkstra类似,Prim算法也可以用堆优化,优先队列代替堆,优化的Prim算法时间复杂度O(mlogn)。适用于稀疏图,但是稀疏图的时候求最小生成树,Kruskal 算法更加实用。
画图不易,求个点赞~~
就喜欢看
Hasity
大佬写的题解,图画的太好了💕这为什么大循环套小循环变量都用i还能过啊
临时变量。就近原则。
hasity
=有能力但是这样写不易于阅读
or2or2
你这个屁股更翘了
蚌埠住了
hh
hhhhh
原来orz是这个意思hhhh,之前一直不懂hhh
什么意思啊
膜拜的意思吧大概,抽象的符号语言
中国古典象形字
nb !!!
这个图用什么软件画的?
Ppt
谢了,第一次知道😁
求个点赞~,谢谢
Orz!
g[a][b] = g[b][a] = min(g[a][b],w);//存储权重
可能有重边 而且重边权值不一样的情况。之前的代码中并没有考虑过这种情况(^///^)
这一句代码,是用最小权重更新a-b,b-a的距离
很棒!
序号6所在表格中,state[3]是不是也需要填写成为1呢
如果改成 以下,就会 WA,请问为什么?
for(int i = 1; i <= n; i++)//更新生成树外的点到生成树的距离
{
if(dt[i] > g[t][i] )//从 t 到节点 i 的距离小于原来距离,则更新。
{
dt[i] = g[t][i];//更新距离
pre[i] = t;//从 t 到 i 的距离更短,i 的前驱变为 t.
}
}
res += dt[t];
更新距离之前就要把t节点加入生成树里面,不然有自环会出问题
哦,懂了,谢谢
有个问题,时间复杂度应该是(O(nm))吧???
对于那种完全图来说,堆优化似乎会降成O(nn*log2(n))
顺带一提,正常的Prime应该是O(n*n+m)吧???
你是没学过时间复杂度吗时间复杂度不考虑常数, $n^2$ 是肯定大于 $m$ 的。说错了,是大于等于
话说memset可以用在二维数组吗?
我记得是可以的,不敢肯定,
可以
更新距离时 不加!st[i]也能过,这个影响什么吗
同问+1
从结果来说,没影响。选过的点即便更新了,也不会再选。
因为你选过点点权值放在res里面去了,即使有自环让已经放入集合的d[j]变化也不会影响res中的值
说得好
太强了QAQ
看出了动画的感觉
getpath写的有问题,应该是从i=n-1到i=1循环
下边是从 1 开始的
虽然只有十个月,但是感觉已经是远古评论了hh 你说的对
还可以在远古一点hhh,你们都说的对
再远古一点,你们说的都对
这里的i应该改成j
支持!
你一声令下,我们就拥立你为新的yxc
文化程度太高了,看不懂
orz
我悟了