次小生成树

最小生成树

---

设G=(V, E,ω)是连通的无向图，T是图G的一个最小生成树。如果有另一棵树T1，满足不存在树T’,T’≠T, ω（T’）<ω（T1），则称T1是图G的次小生成树。

---

定理

定理:设T是图G的最小生成树，如果T1满足ω(T1)=min{ω(T’)| T’∈Not(T)},则T1是G的次小生成树。

也就是说，最小生成树邻集中权值和最小的一棵生成树即为该图的次小生成树。

证明：
可以证明下面一个强一些的结论：
T是某一棵最小生成树，T0是任一棵异于T的树，通过变换 T₀ –> T₁ –> T₂ –> … –> T_n (T) 变成最小生成树。
所谓的变换是，每次把Ti中的某条边换成T中的一条边, 而且树T_(i+1)的权小于等于Ti的权。

操作

在Ti中任取一条不在T中的边(u, v)。

把边(u, v)去掉，就剩下两个连通分量A和B，在T中，必有唯一的边(u‘, v’) 连结A和B。

显然(u‘, v’)的权比(u, v)小 (否则,(u, v)就应该在T中)。把(u‘, v’)替换(u, v)即得树T(i+1)。

特别地：取T0为任一棵次小生成树，则T(n-1)满足定义且跟T差一条边。上述定理得证。

算法

Step 1.
先用prim求出最小生成树T。在prim的同时，用一个矩阵max_val[u][v] 记录在T中连结任意两点u,v的唯一的路中权值最大的那条边的权值。

这是很容易做到的，因为prim是每次增加一个结点s，而设已经标号了的结点集合为W，则易求W中所有点到s的路中的最大边权值

设 u 属于W,且 s是被连接到W中的v点的，

则

Max_val[v][s] = 边(v,s)的权
Max_val[u][s] = Max( Max_val[v][s], Max_val[u][v])

用时O(V^2)。

const int MAXN=110;
const int INF=0x3ffffff;
bool vis[MAXN];
int lowc[MAXN];
int pre[MAXN];
int Max[MAXN][MAXN];

int n;
int cost[MAXN][MAXN];
bool used[MAXN][MAXN];
int prim(){
    int mst = 0;
    memset(vis, false, sizeof(vis));
    memset(Max, 0, sizeof(Max));
    memset(used, 0, sizeof(used));
    vis[0] = true;
    pre[0] = -1;
    for (int i = 1; i < n; i++){
        lowc[i] = cost[0][i];
        pre[i] = 0;
    }
    lowc[0] = 0;
    for (int i = 1; i < n; i++){
        int minc = INF;
        int p = -1;
        for (int j = 0; j < n; j++){
            if (!vis[j] && minc>lowc[j]){
                p = j;
                minc = lowc[j];
            }
        }
        if (p == -1){
            return -1;
        }
        mst += minc;
        vis[p] = true;
        used[p][pre[p]] = used[pre[p]][p] = true;
        for (int j = 0; j < n; j++){
            if (vis[j]){
                Max[j][p] = Max[p][j] = max(Max[j][pre[p]], lowc[p]);
            }else if (lowc[j] > cost[p][j]){
                lowc[j] = cost[p][j];
                pre[j] = p;
            }
        }
    }
    return mst;
}
int secondTree(){
    int res = INF;
    int mst = prim();
    for (int i = 0; i < n; i++){
        for (int j = 0; j < n; j++){
            if (i != j && !used[i][j]){
                res = min(res, mst-Max[i][j]+cost[i][j]);
            }
        }
    }
    return res;
}

---

内容参考北大讲义