Prim求最小生成树

一、最小生成树几种算法的对比

朴素版Prim算法主要用于稠密图，用邻接矩阵；堆优化版Prim用于稀疏图，用邻接表；
Kruskal算法是一个很优美的算法，适用于稀疏图，且时间复杂度与堆优化版Prim一致，因此最小生成树一般用朴素版Prim和Kruskal。

二、Prim算法

1.一些概念

树：连通、无圈图。
生成树：如果一棵树T是连通图G的生成子图，则称树T为图G的生成树或支撑树。
最小生成树：赋权无向连通图G的所有生成树中权重之和最小的一个。

2.基本思路

设集合S为已确认加入最小生成树的顶点，为一个连通块。
算法基本步骤：

dist[i] <- 无穷
for i in 0:n-1:
    t <- 找到集合外距离最近的点
    st[t] = true // 加入集合S
    累加生成树的权值
    用t去更新其他点到集合的距离

可以看到该算法与Dijkstra很类似，都是贪心算法。只是更新时，更新到集合的距离，而非到起点的距离。

3.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≤10⁵,
图中涉及边的边权的绝对值均不超过 10000。
输入样例：

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

输出样例：

6

思路
根据n,m的范围可知该图为稠密图，因此用朴素版Prim算法。
代码

#include <iostream>
#include <cstring>
#include <algorithm>

using namespace std;

const int N = 510, INF = 0x3f3f3f3f;

int n, m;
int g[N][N]; // 稠密图用邻接矩阵
int dist[N]; // 到集合的距离
bool st[N];  // 是否在集合S中

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; // 加入集合S
        // 用t更新其他点到集合S的距离
        for(int j = 1; j <= n; j ++ ) dist[j] = min(dist[j], g[t][j]);
        
    }
    
    return res;
}

int main()
{
    cin >> n >> m;
    
    memset(g, 0x3f, sizeof g);
    
    while(m --)
    {
        int a, b, c;
        cin >> a >> b >> c;
        g[a][b] = g[b][a] = min(g[a][b], c); // 重边取最小
    }
    
    int t = prim();
    
    if(t == INF) puts("impossible");
    else printf("%d\n", t);
    
    return 0;
}