poj 2914 全局最小割 stoer wagner算法 最小割模板

题目链接

http://poj.org/problem?id=2914

题意

求最小割。
最小割：去掉一个边集，使得图分为两个部分，最小割就是这些边的权值之和最小时的边集。

/*
1. 设最小割cut=INF, 任选一个点s到集合A中, 定义W(A, p)为A中的所有点到A外一点p的权总和.
2. 对刚才选定的s, 更新W(A,p)(该值递增).
3. 选出A外一点p, 且W(A,p)最大的作为新的s, 若A!=G(V), 则继续2.
4. 把最后进入A的两点记为s和t, 用W(A,t)更新cut.
5. 合并st，即新建顶点u, 边权w(u, v)=w(s, v)+w(t, v), 删除顶点s和t, 以及与它们相连的边.
6. 若|V|!=1则继续1.
*/

#include<cstdio>
#include<queue>
#include<iostream>
#include<vector>
#include<map>
#include<cstring>
#include<string>
#include<set>
#include<stack>
#include<algorithm>
#define cle(a) memset(a,0,sizeof(a))
#define inf(a) memset(a,0x3f,sizeof(a))
#define ll long long
#define Rep(i,a,n) for(int i=a;i<=n;i++)
using namespace std;
const int INF = ( 2e9 ) + 2;
const ll maxn = 510;
int g[maxn][maxn];
int vis[maxn];// 判断i是否加入A集合中 
int w[maxn];  // w[i] 代表A集合中所有点到i点的距离 
int v[maxn];  // v[i] 代表i点所合并到的点 
int SW(int n)
{
    int ret=INF;
    for(int i=0;i<n;i++)v[i]=i;
    while(n>1)
    {
        memset(w,0,sizeof(w));
        memset(vis,0,sizeof(vis));
        int pre=0;
        for(int i=1;i<n;i++)
        {
            int k=-1;
            for(int j=1;j<n;j++)    // 寻找“距离”A集合最大的点 
            {
                if(!vis[v[j]])
                {
                    w[v[j]]+=g[v[pre]][v[j]];
                    if(k==-1||w[v[j]]>w[v[k]])k=j;
                }
            }
            vis[v[k]]=1;
            if(i==n-1)
            {
                int s=v[pre],t=v[k];
                ret=min(ret,w[t]);
                for(int j=0;j<n;j++)
                {
                    g[s][v[j]]+=g[v[j]][t];
                    g[v[j]][s]+=g[v[j]][t];
                }
                v[k]=v[--n];
            }
            pre=k;
        }
    }
    return ret;
}
int main()
{
    int n,m;
    while(~scanf("%d%d",&n,&m))
    {
        memset(g,0,sizeof(g));
        for(int i=0;i<m;i++)
        {
            int u,v,w;
            scanf("%d%d%d",&u,&v,&w);
            g[u][v]+=w;
            g[v][u]+=w;
        }
        int ans=SW(n);
        printf("%d\n",ans);
    }
}