poj Friendship 最小割+拆点

题目链接

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

题意：

在一个关系网中，如果说A与B保持联系，那么A应该知道B的电话号码。或者存在C，A知道C的电话号码，C知道B的电话号码。可以认为如果A知道B的电话号码，那么B也知道A的电话号码。在一个关系网中可能会出现一些意外事故，比如某人更换手机号码了，导致他与其他人都失去了联系。问给定两个人S和T，问，最少多少人出现意外事故，使得S和T失去联系。

思路

如果源点和汇点是直接相连的，那么是”NO ANWER!”，然后拆点建图，对每个点u拆成入点和出点,(u’,u”,1)代表u’到u”有一条容量为1的边。在原图中，直接与u相连的点v，在新建的图中应该变成两条边(u”,v’,INF),(v”,u’,INF)。原图中源点的出点，在新图中作为源点，原图汇点的入点，作为新建图的汇点。跑一次最大流可以获得最小割，即应该删除的点数。题目要求按字典序输出应该删除的人，那么从小到大枚举I（除了s和t），将I删去，重新建图，跑最大流，看获得的最小割是否小于之前获得的，如果小于，说明删除的I点的“拆边”正是之前跑最大流时被割中的边，不用复原，将获得的最小割减一重复这个过程，知道最小割为0。否则，I不是被割中的边，那么将i复原到图中。
拆点类似这道题：http://blog.csdn.net/chen_minghui/article/details/77454041

#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 = 500;
int mp[210][210];
struct Dinic
{
    struct edge
    {
        int from,to,c;
        edge(int getu,int getv,int getc)
        {
            from=getu;
            to=getv;
            c=getc;
        }
    };
    int s,t,n;
    int d[maxn];
    int cur[maxn];
    vector<edge> e;
    vector<int> g[maxn];
    void addedge(int u,int v,int c)
    {
        e.push_back(edge(u,v,c));
        e.push_back(edge(v,u,0));
        int m=e.size();
        g[u].push_back(m-2);
        g[v].push_back(m-1);
    }
    void init()
    {
        e.clear();
        for(int i=0; i<=maxn; i++)
            g[i].clear();
    }
    bool bfs()
    {
        queue<int> q;
        q.push(s);
        memset(d,0,sizeof(d));
        d[s]=1;
        while(!q.empty())
        {
            int u=q.front();
            q.pop();
            for(int i=0; i<g[u].size(); i++)
            {
                int v=e[g[u][i]].to;
                int c=e[g[u][i]].c;
                if(!d[v]&&c)
                {
                    d[v]=d[u]+1;
                    q.push(v);
                }
            }
        }
        return d[t];
    }
    int dfs(int u,int maxf,int t)
    {
        if(u==t)return maxf;
        int ret=0;
        for(int &i=cur[u]; i<g[u].size(); i++)
        {
            int c=e[g[u][i]].c;
            int v=e[g[u][i]].to;
            int f;
            if(d[u]+1==d[v]&&c)
            {
                f=dfs(v,min(maxf-ret,c),t);
                e[g[u][i]].c-=f;
                e[g[u][i]^1].c+=f;
                ret+=f;
                if(ret==maxf)return ret;
            }
        }
        return ret;
    }
    int maxflow(int s,int t)
    {
        this->s=s;
        this->t=t;
        int flow=0;
        while(bfs())
        {
            memset(cur,0,sizeof(cur));
            int temp=dfs(s,INF,t);
            flow+=temp;
        }
        return flow;
    }
}dinic;
void build(int n)
{
    for(int i=1;i<=n;i++)
    for(int j=1;j<=n;j++)
    {
        if(mp[i][j])            // 题目中给的数据是无向图的数据，且对角线为1，所以可以这样写来拆点。 
        {
            if(i==j)
            dinic.addedge(i,i+n,1);
            else
            dinic.addedge(i+n,j,INF);
        }
    }
}
int main()
{
    int n,s,t;
    while(~scanf("%d%d%d",&n,&s,&t))
    {
        dinic.init();
        for(int i=1;i<=n;i++)
        for(int j=1;j<=n;j++)
        scanf("%d",&mp[i][j]);
        if(mp[s][t])
        {
            printf("NO ANSWER!\n");
            continue;
        }
        build(n);

        int mincut=dinic.maxflow(s+n,t);
        int ans[210];
        int cnt=0;
        int buf[210],buf2[210];
        for(int i=1;i<=n;i++)       // 去掉i点 
        {
            if(i==s||i==t)continue;
            dinic.init();

            for(int j=1;j<=n;j++)
            {
                buf[j]=mp[i][j];
                buf2[j]=mp[j][i];
                mp[i][j]=mp[j][i]=0;
            }
            build(n);

            int temp=dinic.maxflow(s+n,t);
            if(mincut>temp)
            {
                ans[cnt++]=i;
                mincut--;
            }
            else
            {
                for(int j=1;j<=n;j++)
                mp[i][j]=buf[j],mp[j][i]=buf2[j];       //把点恢复 
            }
            if(mincut==0)break;
        }
        printf("%d\n",cnt);
        for(int i=0;i<cnt;i++)
        if(i!=cnt-1)
        printf("%d ",ans[i]);
        else
        printf("%d\n",ans[i]);
    }
}