最小生成树专题总结

最小生成树的概念　

　　　给定无向图G = (V, E)，连接G中所有点，且边集是E的子集的树称为G的生成树，而权值和最小的生成树称为最小生成树，即MST。

　　　构造MST的方法有很多种。常用的有Kruskal算法和Prim算法，前者好写，时间复杂度为O(m)，后者稍微难写，时间复杂度O(n*n)。(n为树的节点数，m为边数)。

　　　

Kruskal算法(摘自刘汝佳白书P199)：

　　　算法的第一步是给所有边按照从小到大的顺序排序，然后从小到大考查所有边。考查到边(u, v)的时候有两种情况：

　　　情况1：u, v此时属于同一个连通分量中，则加入(u, v)会形成环，不能加入该边。

　　　情况2：u, v属于不同的连通分量。那么加入(u, v)一定是最优情况。为什么呢？用反证法证。如果不加这条边得到最优解T，则T+(u,v)一定有且只有一个环，而且环中至少有一条边(u', v')权值大于或等于(u, v)。删除该边后，得到新树T' = T + (u,v) - (u', v')，权值和T <= T，所以加入(u, v)不会比不加入差。

　　　算法中，判断是否属于同一个连通分量用并查集即可。

//Kruskal算法求MST，返回MST的权值和，如果不联通返回-1
//时间复杂度O(m)

struct Pat{                 //表示边
    int s, e, w;            //s, e表示边的两个点，w为权值
};

int f[N];                   //并查集

int kruskal(int n, int m)       //n为点的数量，m为边的数量
{
    sort (p, p+m, cmp);
    for (int i = 0; i < n; ++ i)
        f[i] = i;
    
    int cost = 0;
    for (int i = 0; i < m; ++ i){
        int t1 = find(p[i].s), t2 = find(p[i].e);
        if (t1 != t2){
            cost += p[i].w;
            f[t1] = t2;
        }
    }

    int tmp;
    for (int i = 0; i < n; ++ i){
        if (!i) tmp = find(i);
        else if (tmp != find(i)) return -1;
    }
    return cost;
}

 

Prim算法：

　　　详见我的另一篇文章，http://www.cnblogs.com/plumrain/p/Prim_Dijskra.html。比较讨论Dijskra和Prim算法的。

//Prim算法求MST，返回MST的权值和，如果图不联通返回-1
//时间复杂度O(n*n)
//这代码是用优先队列priority_queue写的，其实用set来写会快一些。

typedef pair<int, int> pii;
const int N = ;                         //N表示节点数

int c[N];                               //c[i]表示i点目前到已选出的点中最短边的权值
bool v[N];                              //v[i] = 0表示i点目前还没有被选出
vector<pii> p[N];                       //p[i]表示一点为i的边，另一点为p[i][j].first，该边权值为p[i][j].second

int prim(int n, int m)                  //n为节点数，m为边数
{
    int cost = 0;
    CLR (v);
    priority_queue<pii, vector<pii>, cmp> q;
    while (!q.empty()) q.pop();

    for (int i = 1; i < n; ++ i)
        c[i] = maxint - 1000;
    c[0] = 0;
    v[0] = 1;
    for (int i = 0; i < (int)p[0].size(); ++ i){
        pii tmp = p[0][i];
        c[tmp.first] = tmp.second;
        q.push (make_pair(tmp.first, tmp.second));
    }

    for (int i = 0; i < n-1; ++ i){
        if (q.empty()) break;
        pii tmp = q.top(); q.pop();
        while (v[tmp.first] && !q.empty()){
            tmp = q.top(); q.pop();
        }
        if (v[tmp.first]) break;

        int t1 = tmp.first, t2 = tmp.second;
        cost += t2;
        v[t1] = 1;
        c[t1] = 0;
        for (int j = 0; j < (int)p[t1].size(); ++ j) if (!v[p[t1][j].first]){
            pii t = p[t1][j];
            if (c[t.first] > t.second){
                q.push (make_pair(t.first, t.second));
                c[t.first] = t.second;
            }
        }
    }

    for (int i = 0; i < n; ++ i)
        if (!v[i]) return -1;
    return cost;
}

 

题目训练：

虽然在网上挂了一个专题，但是做了以后才发现里面几乎都是裸题。所以就挑其中一些题来写题解，以后碰到好题再添加。

1、HDU 1863，HDU 1875，HDU 1879，HDU 1233

　　　畅通工程系列，都是很裸的MST的题。贴一篇题解，http://www.cnblogs.com/plumrain/p/HDU_changtong.html。

 

2、POJ 2349 Arctic Network

　　　对于给定图，求其生成树的第k大的边最小为多少。其实很容易想到用MST，但是关键是为什么。要证明也不容易。http://www.cnblogs.com/plumrain/p/POJ_2349.html

 

3、POJ 1679 The Unique MST

　　　问最小生成树是否唯一。http://www.cnblogs.com/plumrain/p/POJ_1679.html

 

4、这是我目前做过的一些水题：ZOJ 1586，HDU 1102(POJ 2421)，POJ 1251(HDU 1301)，POJ 1287，POJ 2031，POJ 1789，POJ 1751，POJ 1258，POJ 3026