题目链接 : 点击查看
题目描述 :
给定一个 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。
输入输出样例 :
输入
4 5
1 2 1
1 3 2
1 4 3
2 3 2
3 4 4
输出
6
题目分析 :
Prim算法用于处理稠密图(m ≈ n ^ 2)的最小生成树的问题,与Dijkstra算法类似,Prim算法主要通过贪心策略从某个点每一次找到边权最小的边,并将这条边的另一个顶点加入最小生成树点集之中,重复n次,直到点集中加入所有顶点。从上述解释可知,Prim算法中要维护一个最小生成树的顶点集,我们用st[N]数组来表示当前点是否在点集之中,而Dijsktra算法中st[N]数组表示某点是否已经确定到源点的最短路径。同样不同的是还有dist[N]数组,在Dijkstra算法中表示的某点到源点的最短距离,而在Prim算法中却表示到点集的最短距离。也就是,相比于Dijkstra算法来说,我们每一次找到点集的最短路径点t,我们要更新的dist数组不再是d[j]与dist[t] + g[t][j]相比较,而是d[j]与g[t][j]相比较,这就是点j到源点的距离与点j到点集的距离的不同。详见如下代码。
代码 :
#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;
const int N = 510, INF = 0x3f3f3f3f;
int n, m;
int g[N][N];
int dist[N];//dist维护的是距离点集的最短距离
bool st[N];//st[i]用来记录点i是否在最小生成树的点集中
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 (!st[j] && (t == -1 || dist[t] > dist[j]))
t = j;
}
if (i && dist[t] == INF) return INF;//当前找到的最短距离为INF说明图不连通无最小生成树
if (i) res += dist[t];//计算最小生成树的权重和
st[t] = true;//将t点加入最小生成树点集中
for (int j = 1; j <= n; j ++ ) dist[j] = min(dist[j], g[t][j]);//更新 因为dist表示是距离点集的距离所以只需更新点集外的边权最小的那条边即可
}
return res;
}
int main() {
cin >> n >> m;
memset(g, 0x3f, sizeof(g));
while (m -- ) {
int v1, v2, w;
cin >> v1 >> v2 >> w;
g[v1][v2] = g[v2][v1] = min(g[v1][v2], w);
}
int t = prim();
if (t == INF) cout << "impossible" << endl;
else cout << t << endl;
return 0;
}
下面我们给出Prim算法的相关模板
时间复杂度是 O(n^2+m), n 表示点数, m 表示边数
int n; // n表示点数
int g[N][N]; // 邻接矩阵,存储所有边
int dist[N]; // 存储其他点到当前最小生成树的距离
bool st[N]; // 存储每个点是否已经在生成树中
// 如果图不连通,则返回INF(值是0x3f3f3f3f), 否则返回最小生成树的树边权重之和
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 (!st[j] && (t == -1 || dist[t] > dist[j]))
t = j;
if (i && dist[t] == INF) return INF;
if (i) res += dist[t];
st[t] = true;
for (int j = 1; j <= n; j ++ ) dist[j] = min(dist[j], g[t][j]);
}
return res;
}