POJ1815 Friendship(求最小割点 dinic)

时间:2022-03-15 04:27:31

传送门:http://poj.org/problem?id=1815
题意就是求s点到t点,最少去掉几个点使得他们不连通。如果无解输出NO ANSWER!
       因为最小割只能求割掉几条边的解,我们要求的是割掉几个点。那么我们可以这样考虑:把每个点拆成入点和出点。入点->出点权值为1。那么割掉这条边就相当于割掉这个点了,就能把这题转化成最小割。那么原来的边,我们是不希望去割它的,所以我们将原来的边的权值设置为INF。例如,原来边为u->v那么就变成out(u)->in(v)。这样跑一次最大流,跑出来的答案就是最小割掉的点数。
       这题还要我们求出一个点集,就是割掉的点,按字典序排列。因为割点可能会有多种组合,所以最简单的思路就是枚举。如果这个点割掉之后,跑出来的答案能少割一个点,那么这个点就是要求的点之一,割掉它,继续枚举,一直到枚举出所有答案为止。
       应该还有不枚举的更好方法,求教学。。。
代码:

#include<cstdio>
#include<cstring>
#include<string>
#include<queue>
#include<map>
#include<iostream>
#include<vector>
#include<cstdlib>
#include<cmath>
#include<stack>
#include<cctype>
#include<set>
#include<ctime>
#include<cassert>
#include<algorithm>
using namespace std;
typedef long long ll;
const int INF=1e9+7;
typedef pair<int,int> pii;
#define MP make_pair
#define PB push_back
#define in(x) (x)
#define out(x) (x+n)
int n,st,ed;
int start[210][210],pic[555][555],d[555];
bool del[555];
vector<int> pr;
void build(){
    for(int i=1;i<=n;i++){
        if(!del[i])
            pic[in(i)][out(i)]=1;//拆点,后面枚举了割掉的点,所以用del数组标记
        //容量是1,希望在割边的时候割它,割掉点中的边相当于割点,就把最小割转化成割点了
    }
    for(int i=1;i<=n;i++){
        for(int j=1;j<=n;j++)if(start[i][j]){
            pic[out(i)][in(j)]=INF;
            //割这条边的流量是INF,意味着不希望去割它
        }
    }
}
bool BFS(){
    queue<int> Q;
    memset(d,-1,sizeof d);
    Q.push(st);d[st]=0;
    while(!Q.empty()){
        int s=Q.front();Q.pop();
        for(int i=1;i<=n+n;i++){
            if(pic[s][i]>0&&d[i]<0){
                d[i]=d[s]+1;if(i==ed)return true;
                Q.push(i);
            }
        }
    }
    return false;
}
int DFS(int s,int t,int flow){
    if(s==t||flow==0)return flow;
    int ans=0;
    for(int i=1;i<=n+n;i++){
        if(d[i]==d[s]+1&&pic[s][i]>0){
            int ff=DFS(i,t,min(pic[s][i],flow));
            if(ff>0){
                pic[s][i]-=ff;
                pic[i][s]+=ff;
                ans+=ff;
                flow-=ff;
                if(!flow)break;
            }
        }
    }
    if(!ans)d[s]=-1;
    return ans;
}
int dinic(){
    int ans=0;
    while(BFS()){
        ans+=DFS(st,ed,INF);
    }
    return ans;
}
int main()
{
// freopen("D://input.txt","r",stdin);
    scanf("%d%d%d",&n,&st,&ed);
    for(int i=1;i<=n;i++){
        for(int j=1;j<=n;j++){
            scanf("%d",&start[i][j]);
        }
    }
    if(start[st][ed]){printf("NO ANSWER!\n");return 0;}
    //只有直接相连的时候会无解
    st=out(st);ed=in(ed);
    build();
    int ans=dinic();
    printf("%d\n",ans);
    if(!ans)return 0;
    for(int i=1;i<=n;i++){
    //因为割点可能会有很多组合,而且要求字典序最小,所以只能割一个就标记上,一直到割完为止
        if(i==st-n||i==ed)continue;
        del[i]=true;
        build();
        int t=dinic();
        if(t<ans){
            ans--;pr.PB(i);
            if(ans==0)break;
        }
        else del[i]=false;;
    }
    for(int i=0;i<(int)pr.size();i++){
        if(i)printf(" ");
        printf("%d",pr[i]);
    }
    printf("\n");
    return 0;
}