Minimum Cut POJ

来源:互联网 发布:天书残卷披风进阶数据 编辑:程序博客网 时间:2024/06/08 15:38

无向图最小割,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;}
原创粉丝点击