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

来源:互联网 发布:怎样才能约到炮 知乎 编辑:程序博客网 时间:2024/06/05 18:45

题目链接

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);    }}
阅读全文
0 0
原创粉丝点击