POJ 1815 Friendship 最小割+字典序最小割集

解题思路：

（1）求最小割。 把每个点i拆分为i和i',对于原图中每条边(u,v)建边( u' , v ),容量为1；对于每个点i，建（i , i'），容量为1。

从s'到t求一遍最小割就行了。

（2）按字典序输出 最小割集。

先记录好原图的最小割ans,从小到大枚举每各点, 如果在不建该点的( i, i' )的那条边的情况下，最小割<ans ，

说明该点属于最小割集，在以后的枚举其他点的时候建边也不建该点的（i,i'）,同时把ans更新为当前的最小割。

#include <iostream>
#include <cstdio>
#include <cstring>
#include <algorithm>
#include <queue>
#include <stack>
using namespace std;
#define FOR(i,l,r) for(int i=(l);i<=(r);++i)
#define REP(i,n) for(int i=0;i<(n);++i)
#define DSC(i,r,l) for(int i=(r);i>=(l);--i)

#define N 410
#define M 11000
#define INF 1e9
bool map[210][210];
bool flag[210];
struct
{
    int to,next;
    int c;
}edge[M];
int head[N],level[N],ip;
int que[N];
bool makelevel(int s,int t)
{
    memset(level,0,sizeof(level));
    int iq=0;
    que[iq++]=s;
    level[s]=1;
    int top;
    for(int i=0;i<iq;i++)
    {
        top=que[i];
        if(top==t)  return 1;
        for(int k=head[top];k!=-1;k=edge[k].next)
        {
            if(!level[edge[k].to]&&edge[k].c>0)
            {
                que[iq++]=edge[k].to;
                level[edge[k].to]=level[top]+1;
            }
        }
    }
    return 0;
}
int dfs(int now,int maxf,int t)
{
    if(now==t)  return maxf;
    int  ret=0,c;
    for(int k=head[now];k!=-1;k=edge[k].next)
    {
        if(edge[k].c>0&&level[edge[k].to]==(level[now]+1))
        {
            c=dfs(edge[k].to,min(maxf-ret,edge[k].c),t);
            edge[k].c-=c;
            edge[k^1].c+=c;
            ret+=c;

            if(ret==maxf)   return  ret;
        }
    }
    if(!ret) --level[now];
    return ret;
}
int dinic(int s,int t)
{
    int  ans=0;
    while(makelevel(s,t))   ans+=dfs(s,INF,t);
    return ans;
}
void add(int u,int v,int c,int f)//有向边f为0 ,否则为 c
{
    edge[ip].to=v;edge[ip].c=c;edge[ip].next=head[u];head[u]=ip++;
    edge[ip].to=u;edge[ip].c=f;edge[ip].next=head[v];head[v]=ip++;
}

int main()
{
    int n,s,t,x;
    while(cin>>n>>s>>t)
    {
        int start=s+n,end=t;
        memset(head,-1,sizeof(head)); ip=0;
        memset(map,0,sizeof(map));
        memset(flag,0,sizeof(flag));
        FOR(i,1,n)
            FOR(j,1,n)
            {
                scanf("%d",&x);
                if(x) map[i][j]=1;
            }

        if(map[s][t] || s==t )
        {
            puts("NO ANSWER!");
            continue;
        }

        FOR(i,1,n)
            FOR(j,1,n)
                if(map[i][j])
                    add(i+n,j,1,0);
        FOR(i,1,n)  add(i,i+n,1,0);
        int ans=dinic(start,end);//求最小割

        cout<<ans<<endl;
        if(ans==0) continue;
        int num=0,temp;
        for(int k=1;k<=n;++k)//从小到大进行枚举
        {
            if(k==s || k==t) continue;
            memset(head,-1,sizeof(head)); ip=0;
            FOR(i,1,n)
                FOR(j,1,n)
                    if(map[i][j])
                        add(i+n,j,1,0);
            FOR(i,1,n) if(i!=k&&!flag[i]) add(i,i+n,1,0);//之前找到的割点不建边,当前枚举的点不建边

            if((temp=dinic(start,end))<ans)//如果最小割小于ans
            {
                if(num) putchar(' ');
                printf("%d",k);
                num++;
                flag[k]=1;//把当前点更新为割点
                ans=temp;//同时更新ans
            }
        }
        cout<<endl;
    }
    return 0;
}