Kruskal、Prim

引用了算法导论--最小生成树(Kruskal和Prim算法)

稠密用Prim,稀疏用Kruskal

1.Kruskal

此算法可以称为“加边法”,初始最小生成树边数为0,每迭代一次就选择一条满足条件的最小代价边,加入到最小生成树的边集合里。 
1. 把图中的所有边按代价从小到大排序; 
2. 把图中的n个顶点看成独立的n棵树组成的森林; 
3. 按权值从小到大选择边,所选的边连接的两个顶点ui,vi,应属于两颗不同的树,则成为最小生成树的一条边,并将这两颗树合并作为一颗树。 
4. 重复(3),直到所有顶点都在一颗树内或者有n-1条边为止。

struct Edge{
    int from, to, w;
    Edge(int x, int y, int z): from(x), to(y), w(z){}
    bool operator < (const Edge &rhs) const{
        return w > rhs.w;
    }
};
int parent[1005];
priority_queue<Edge> q;
queue<Edge> q_res;
void make_set(int x){parent[x]=x;}
int findset(int x){return parent[x] != x ? parent[x] = findset(parent[x]) : x;}
void unionset(int x, int y){
    if(x < y)
        parent[y] = x;
    else
        parent[x] = y;
}
void kruskal(){
    for(int i = 1; i <= N; i++)
        make_set(i);
    while(!q.empty()){
        Edge tmp = q.top(); q.pop();
        int x = findset(tmp.from);
        int y = findset(tmp.to);
        if(x != y){
            if(res < tmp.w)
                res = tmp.w;
            unionset(x, y);
            q_res.push(tmp);
        }
    }
}
void init(){
    while(!q.empty())
        q.pop();
    while(!q_res.empty())
        q_res.pop();
}
int main(){
    while(scanf("%d%d", &N, &M) != EOF){
        init();
        res = -1;
        int a, b, wt;
        for(i = 0; i < M; i++){
            scanf("%d%d%d", &a, &b, &wt);
            q.push(Edge(a, b, wt));
        }
        kruskal();
        printf("%d\n", res);
        printf("%d\n", q_res.size());
        while(!q_res.empty()){
            Edge tmp = q_res.front(); q_res.pop();
            printf("%d %d\n", tmp.from, tmp.to);
        }
    }
    return 0;
}

2.Prim

此算法可以称为“加点法”,每次迭代选择代价最小的边对应的点,加入到最小生成树中。算法从某一个顶点s开始,逐渐长大覆盖整个连通网的所有顶点。

  1. 图的所有顶点集合为V;初始令集合u={s},v=V−u;
  2. 在两个集合u,v能够组成的边中,选择一条代价最小的边(u0,v0),加入到最小生成树中,并把v0并入到集合u中。
  3. 重复上述步骤,直到最小生成树有n-1条边或者n个顶点为止。

由于不断向集合u中加点,所以最小代价边必须同步更新;需要建立一个辅助数组closedge,用来维护集合v中每个顶点与集合u中最小代价边信息,:

#define INF 0x3f3f3f3f
#define VNUM  1e5 + 5
 
struct Node{
    int v,next,w;
    bool operator < (const Node &a) const{
        return w > a.w;
    }
} p[N],t1,t2;
 
int dis[N], vis[N], head[N], cnt;
int res;
 
void addedge(int u,int v,int w){
    p[cnt].v = v;
    p[cnt].next = head[u];
    p[cnt].w = w;
    head[u] = cnt++;
}
 
void prim(){
    priority_queue<Node> q;
    for(int i = head[0] ; i != -1 ; i = p[i].next){
        int v = p[i].v;
        if(p[i].w < dis[v]){
            dis[v] = p[i].w;
            t1.w = dis[v];
            t1.v =  v;
            q.push(t1);
        }
    }
    vis[0] = 1;
    while(!q.empty()){
        t1 = q.top();
        q.pop();
        int u = t1.v;
        if(vis[u]) continue;
        vis[u] = 1;
        res += dis[u];
        for(int i = head[u]; i != -1; i = p[i].next){
            int v = p[i].v;
            if(!vis[v] && dis[v] > p[i].w){
                dis[v] = p[i].w;
                t2.v = v;
                t2.w = dis[v];
                q.push(t2);
            }
        }
    }
}
 
int main(){
    int n, m, w;
    while(scanf("%d", &n), n){
        memset(p, 0, sizeof(p));
        memset(head, -1, sizeof(head));
        memset(vis, 0, sizeof(vis));
        char u, v;
        for(int i = 0; i < n - 1; i++){
            cin >> u >> m;
            for(int j = 0; j < m; j++){
                cin >> v >> w;
                addedge(u, v, w);
                addedge(v, u, w);
            }
        }
        for(int i = 0; i < n; i++) dis[i] = INF;
        res = 0;
        prim();
        printf("%d\n", res);
    }
    return 0;
}

猜你喜欢

转载自blog.csdn.net/aaakkk_1996/article/details/82085539