Minimum Cut POJ - 2914 无向图最小割,stoer_wagner算法

19 篇文章 0 订阅
19 篇文章 0 订阅

无向图最小割,stoer_wagner算法

题意
给一副无向图,求一组最小割集,将图分成两份
题解
某算法刚好可解决 stoer_wagner算法
算法流程:
1.min=MAXINT,固定一个顶点P
2.从点P用类似prim的算法扩展出“最大生成树”,记录最后扩展的顶点和最后扩展的边
3.计算最后扩展到的顶点的切割值(即与此顶点相连的所有边权和),若比min小更新min
4.合并最后扩展的那条边的两个端点为一个顶点(当然他们的边也要合并)
5.转到2,合并N-1次后结束
6.min即为所求,输出min
理解?
分割后会有两个集合S和T,
每次将两个点合并
据说按照以上步骤起初只会 合并S或T中的两个点,
这样总会在某次合并后刚好得到答案
复杂度O(n^3)
代码

/*
 * Author       :  Echo
 * Email        :  1666424499@qq.com  
 * Description  :   
 * Created Time :  2017/10/16 17:43:16
 * Last Modify  :  2017/10/17 17:43:23
 * File Name    :  write.cpp
 */
#include <iostream>
#include <cstdio>
#include <cstring>
#include <cmath>
#include <algorithm>
#include <string>
using namespace std;

const int maxn=505;
const int INF=1e9;
int mp[maxn][maxn];
int vis[maxn];
int dis[maxn],dot[maxn];
int N;
int stoer_wagner(int now){
    int ret=INF;
    for(int i=0;i<now;i++) dot[i]=i;
    while(now>1){//重复合并n-1个点
        int pre=0;
        for(int i=0;i<=N;i++) vis[i]=dis[i]=0;
        for(int i=1;i<now;i++){//向集合A中添加元素
            int k= -1;
            for(int j=1;j<now;j++){ //修改跨界边权
                if(vis[dot[j]]) continue;
                dis[dot[j]]+=mp[dot[pre]][dot[j]];
                if(k==-1||dis[dot[j]]>dis[dot[k]]) k=j;//找到最大边权点
            }
            vis[dot[k]]=true;
            if(i==now-1){//只剩两个的时候 合并点
                int s=dot[pre],t=dot[k];
                if(ret>dis[t]) ret=dis[t];
                for(int j=0;j<now;j++){
                    mp[s][dot[j]]+=mp[t][dot[j]];
                    mp[dot[j]][s]+=mp[t][dot[j]];
                }
                dot[k]=dot[--now];
            }
            pre=k;
        }
    }
    return ret;
}
int main(){
    int n,m;
    while(~scanf("%d%d",&n,&m)){
        N=n;
        for(int i=0;i<=n;i++)
            for(int j=0;j<=n;j++)
                mp[i][j]=0;
        for(int i=1;i<=m;i++){
            int u,v,w;
            scanf("%d%d%d",&u,&v,&w);
            mp[u][v]+=w;
            mp[v][u]+=w;
        }
        printf("%d\n",stoer_wagner(n));
    }
    return 0;
}

评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值