网络流+拆点+ISAP+uva563

题意：有一个矩形城市，坐标为整数的的点为路口，之间的连线为街道，现在有b个银行要被抢劫，能不能在经过每个路口至多一次的情况下，劫匪能逃到城市边缘

思路：网络流。建立一个源点s和一个汇点t。将n*m个点的每个点A拆成两个点A1和A2且A1向A2连一条边。如果点A处是银行那么s与A1连一条边，如果B与A相邻，那么A2与B1连一条边，对于所有的边界点C，C2与汇点t连一条边，所有的边的容量都是1，最后做网络流，如果最大流量等于k就可行，否则不可行。

下面是代码：

#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;
const int maxn=50010;
const int INF=1000000000;
int n,m,b;
struct node
{
    int v,next,f;
}edge[500010];
int head[maxn],pre[maxn],gap[maxn],dis[maxn],cur[maxn];
int num,nn;
void init()
{
    num=0;
    memset(head,-1,sizeof(head));
}
void add_edge(int x,int y,int f)
{
    edge[num].v=y;
    edge[num].next=head[x];
    edge[num].f=f;
    head[x]=num++;
    edge[num].v=x;
    edge[num].next=head[y];
    edge[num].f=0;
    head[y]=num++;
}
void Add_adjcant(int x,int y)
{
    int p=(x-1)*m+y+n*m;
    if(x>1){int q=(x-2)*m+y;add_edge(p,q,1);}
    if(x<n){int q=x*m+y;add_edge(p,q,1);}
    if(y>1){int q=(x-1)*m+y-1;add_edge(p,q,1);}
    if(y<m){int q=(x-1)*m+y+1;add_edge(p,q,1);}
}
int sap(int s,int t)
{
    for(int i=0;i<=nn;i++)
    {
        cur[i]=head[i];
        dis[i]=gap[i]=0;
    }
    int flow=0,aug=INF;
    int u=pre[s]=s;
    gap[s]=nn;
    bool flag;
    while(dis[s]<nn)
    {
        flag=false;
        for(int &i=cur[u];i!=-1;i=edge[i].next)
        {
            int v=edge[i].v;
            if(edge[i].f>0&&dis[u]==dis[v]+1)
            {
                flag=true;
                if(aug>edge[i].f)
                    aug=edge[i].f;
                pre[v]=u;
                u=v;
                if(u==t)
                {
                    flow+=aug;
                    while(u!=s)
                    {
                        u=pre[u];
                        edge[cur[u]].f-=aug;
                        edge[cur[u]^1].f+=aug;
                    }
                    aug=INF;
                }
                break;
            }
        }
        if(flag)continue;
        int mindis=nn;
        for(int j=head[u];j!=-1;j=edge[j].next)
        {
            int v=edge[j].v;
            if(dis[v]<mindis&&edge[j].f>0)
            {
                mindis=dis[v];
                cur[u]=j;
            }
        }
        if((--gap[dis[u]])==0)break;
        gap[dis[u]=mindis+1]++;
        u=pre[u];
    }
    return flow;
}
int main()
{
    #ifndef ONLINE_JUDGE
        freopen("in.txt","r",stdin);
    #endif

    int T,x,y,s,t;
    scanf("%d",&T);
    while(T--)
    {
        scanf("%d%d%d",&n,&m,&b);
        s=0,t=2*n*m+1;
        init();
        for(int i=1;i<=n;i++)
        {
            int p;
            for(int j=1;j<=m;j++)
            {
                p=(i-1)*m+j;
                add_edge(p,p+n*m,1);
                Add_adjcant(i,j);
            }
        }
        for(int i=1;i<=m;i++)
        {
            add_edge(i+n*m,t,1);
            add_edge((n-1)*m+i+n*m,t,1);
        }
        for(int i=2;i<n;i++)
        {
            add_edge((i-1)*m+1+n*m,t,1);
            add_edge(i*m+n*m,t,1);
        }
        for(int i=0;i<b;i++)
        {
            scanf("%d%d",&x,&y);
            int p=(x-1)*m+y;
            add_edge(s,p,1);
        }
        nn=t+1;
        if(sap(s,t)==b)cout<<"possible"<<endl;
        else cout<<"not possible"<<endl;
    }
    return 0;
}