POJ1815 Friendship

题意：有n个人，人之间可以相互联系，如果A可以联系B，B可以联系C，那么A也可以联系C，现在给出一个人s和t，问最少要几个人消失，他们两就不能再联系了，他们自己本身不能消失，如果无论消失多少人（除了他们两），他们都能联系就输出NO ANSWER！。。如果有多组情况，按最小字典序输出。

拆点建图，题目要求消失多少个点，拆点之后每个点之间的边的流量就是1，所以就是最少消失多少条边，也就成了求最小割，然后从小到大枚举每一个点，把这个点删去，再重新建图，再跑一遍最大流，看流会不会减小，减小了就保持删去，没有减小就把这条边回复，最后输出删去的点就行了。。

/*************************************************************************
> File Name: poj1815.cpp
> Author: tjw
> Mail:
> Created Time: 2014年11月01日 星期六 12时35分54秒
************************************************************************/

#include<cstdio>
#include<cstring>
#include<algorithm>
#include<cmath>
#include<queue>
#include<vector>
#include<stack>
#include<map>
#define ll long long
#define ls k<<1
#define rs k<<1|1
using namespace std;
const int MAXN=500;
const int MAXE=100010;
const int INF=1<<30;
struct EDGE
{
    int v,next;
    int cap;
}edge[MAXE];
int head[MAXN],size;
void init()
{
    memset(head,-1,sizeof(head));
    size=0;
}
void add_edge(int u,int v,int cap)
{
    edge[size].v=v;
    edge[size].cap=cap;
    edge[size].next=head[u];
    head[u]=size++;
    edge[size].v=u;
    edge[size].cap=0;
    edge[size].next=head[v];
    head[v]=size++;
}
int gap[MAXN],dist[MAXN],pre[MAXN],pe[MAXN];
int nv,n,s,t;
int ISAP()
{
    int cur_flow,temp,u,i,neck;
    int max_flow;
    memset(gap,0,sizeof(gap));
    memset(dist,0,sizeof(dist));
    memset(pre,-1,sizeof(pre));
    for(i=0;i<=nv;i++)
    pe[i]=head[i];
    max_flow=0;
    u=s;
    while(dist[s]<nv)
    {
        if(u==t)
        {
            cur_flow=INF;
            for(i=s;i!=t;i=edge[pe[i]].v)
            {
                if(cur_flow>edge[pe[i]].cap)
                {
                    cur_flow=edge[pe[i]].cap;
                    neck=i;
                }
            }
            for(i=s;i!=t;i=edge[pe[i]].v)
            {
                temp=pe[i];
                edge[temp].cap-=cur_flow;
                edge[temp^1].cap+=cur_flow;
            }
            max_flow+=cur_flow;
            u=neck;
        }
        for(i=pe[u];i!=-1;i=edge[i].next)
        {
            if(edge[i].cap&&dist[u]==dist[edge[i].v]+1)
            break;
        }
        if(i!=-1)
        {
            pe[u]=i;
            pre[edge[i].v]=u;
            u=edge[i].v;
        }
        else
        {
            if(0==--gap[dist[u]])
            break;
            pe[u]=head[u];
            temp=nv;
            for(i=head[u];i!=-1;i=edge[i].next)
            {
                if(edge[i].cap)
                temp=min(temp,dist[edge[i].v]);
            }
            dist[u]=temp+1;
            ++gap[dist[u]];
            if(u!=s)
            u=pre[u];
        }
    }
    return max_flow;
}
int mp[210][210];
bool del[210];
void build()
{
    int i,j;
    init();
    for(i=1;i<=n;i++)
    {
        if(del[i])
            continue;
        for(j=1;j<=n;j++)
        {
            if(del[j])
                continue;
            if(i==j)
                add_edge(i,i+n,1);
            else if(mp[i][j]==1)
                add_edge(i+n,j,INF);
        }
    }
}
int main()
{
    int i,j;
    while(scanf("%d%d%d",&n,&s,&t)==3)
    {
        for(i=1;i<=n;i++)
        {
            for(j=1;j<=n;j++)
            {
                scanf("%d",&mp[i][j]);
            }
        }
        if(mp[s][t])
        {
            printf("NO ANSWER!\n");
            continue;
        }
        s+=n;
        build();
        nv=2*n+1;
        int pre=ISAP();
        int ans=0;
        memset(del,0,sizeof(del));
        for(i=1;i<=n;i++)
        {
            if(pre==0)
                break;
            if(i==s-n||i==t)
                continue;
            del[i]=1;
            build();
            int temp=ISAP();
            if(temp<pre)
            {
                ans++;
                pre=temp;
            }
            else
                del[i]=0;
        }
        printf("%d\n",ans);
        for(i=1;i<=n;i++)
        {
            if(del[i])
                printf("%d ",i);
        }
        printf("\n");
    }
    return 0;
}