无向图最小割,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;
}