poj 1815 最小割+枚举

题意：给一个无向图，求最少删除多少个点，使得从S到T不连通，并输出字典序最小的方案数。

思路：拆点，从S到T跑一遍最小割，求出最少删除的点数。

关键是如何求字典序最小的方案，可以从1到n枚举每个点（拆成的边）是否可以在最小割中，即把这个点所连的所有边断掉，再跑一遍最小割。如果最小割减小，就把这个点从割集中删除；否则，把这个点恢复。

代码：

#include<cstdio>
#include<algorithm>

using namespace std;
typedef long long ll;

const int N=210;
const int M=1e5+10;
const int inf=0x3f3f3f3f;

int head[N<<1],to[M],nxt[M],w[M];
int cur[N<<1],d[N<<1],gap[N<<1];
int tot;

bool vis[N];

int g[N][N];

int n,s,t;
int ans[N];

inline void add(int u,int v,int x){
    to[++tot]=v,nxt[tot]=head[u],w[tot]=x,head[u]=tot;
    to[++tot]=u,nxt[tot]=head[v],w[tot]=0,head[v]=tot;
}

int dfs(int u,int in){
    if(u==t) return in;
    int out=0;
    for(int &i=cur[u];i;i=nxt[i]){
        int v=to[i];
        if(d[v]+1==d[u]){
            int res=dfs(v,min(in,w[i]));
            w[i]-=res,w[i^1]+=res;
            in-=res,out+=res;
            if(!in) return out;
        }
    }
    if(!(--gap[d[u]])) d[s]=(n<<1)+1;
    ++gap[++d[u]],cur[u]=head[u];
    return out;
}

int main()
{
    scanf("%d%d%d",&n,&s,&t);
    for(int i=1;i<=n;i++){
        for(int j=1;j<=n;j++){
            scanf("%d",&g[i][j]);
        }
    }
    if(g[s][t]){
        printf("NO ANSWER!
");
        return 0;
    }
    tot=1;
    for(int i=1;i<=n;i++){
        if(i==s||i==t) add(i,i+n,inf);
        else add(i,i+n,1);
    }
    for(int i=1;i<=n;i++){
        for(int j=i+1;j<=n;j++){
            if(g[i][j]){
                add(i+n,j,inf);
                add(j+n,i,inf);
            }
        }
    }
    int res=0;
    while(d[s]<=(n<<1)) res+=dfs(s,inf);
    int cnt=0;
    for(int i=1;i<=n;i++){
        if(i==s||i==t) continue;
        vis[i]=1;
        tot=1;
        for(int j=1;j<=(n<<1);j++){
            head[j]=0;
            cur[j]=gap[j]=d[j]=0;
        }
        for(int j=1;j<=n;j++){
            if(!vis[j]){
                if(j==s||j==t) add(j,j+n,inf);
                else add(j,j+n,1);
            }
        }
        for(int j=1;j<=n;j++){
            if(!vis[j])
            for(int k=j+1;k<=n;k++){
                if(g[j][k]&&!vis[k]){
                    add(j+n,k,inf);
                    add(k+n,j,inf);
                }
            }
        }
        int tmp=0;
        while(d[s]<=(n<<1)) tmp+=dfs(s,inf);
        if(tmp<res){
            res=tmp;
            ans[++cnt]=i;
        }
        else vis[i]=0;
    }
    printf("%d
",cnt);
    for(int i=1;i<=cnt;i++){
        printf("%d%c",ans[i],i==cnt? '
':' ');
    }
    return 0;
}