UVA, 563 Crimewave

来源：互联网发布：网页源文件提取数据编辑：程序博客网时间：2024/05/16 23:46

题意：团伙抢完所有银行后，撤退到S*A矩阵外就算逃出，要求一个点不能走两次，问是否可以完全逃脱

这道题是最大流问题，主要是要去构建图，然后用最大流算法得出是否银行数量和逃出的数量相等。怎么构建图呢？主要是用拆点，把一个点拆成两个点，点（i，j）可以表示为：前点(i-1)*A+j，后点(i-1)*A+j+M（M为一个较大的数，保证M大于等于S*A就行），然后连接前点和后点，方向是前到后，

相邻的点，图是无向的，用该点的后点连接相邻点的前点。最后用一个超级源点连接所有的银行点的前点，用一个超级终点连接所有的边缘点的后点，所有边的权值为1，图就构建好了。然后用你们熟悉的最大流算法。

Ford-Fulkerson最大流算法代码参考：http://blog.csdn.net/itaskyou/article/details/51331344

#include <iostream>#include<vector>#include<cstring>using namespace std;#define N 5005#define INF 1000000int q[4][2]= {1,0,0,1,-1,0,0,-1};struct edge{    int to,cap,rev;    edge(int a,int b,int c)    {        to=a;        cap=b;        rev=c;    }};vector<edge>v[N];void add_edge(int from,int to,int cap);int dfs(int a,int t,int f);int max_flow(int s,int t);int used[N];int vis[60][60];int m,n,num;int main(){    int t;    cin>>t;    while(t--)    {        int a,b,c,d;        cin>>m>>n>>num;       // memset(vis,0,sizeof(vis));        for(int i=0; i<N; i++)            v[i].clear();        for(int i=1; i<=m; i++)        {            for(int j=1; j<=n; j++)            {                add_edge((i-1)*n+j,(i-1)*n+j+m*n,1);//每个点的前点和后点相连接                for(int k=0; k<4; k++)//和临点相连                {                    int k1=i+q[k][0];                    int k2=j+q[k][1];                    if(k1>0&&k1<=m&&k2>0&&k2<=n)                    {                        add_edge((i-1)*n+j+m*n,(k1-1)*n+k2,1);                    }                }            }        }        for(int i=0; i<num; i++)//接下来的四个for循环是矩阵边缘的点和超级终点相连接。        {            cin>>c>>d;            add_edge(0,(c-1)*n+d,1);        }        for(int i=1; i<=n; i++)        {            add_edge(i+m*n,m*n*2+1,1);        }        if(m>1)        {            for(int i=1; i<=n; i++)            {                add_edge(i+(m-1)*n+m*n,m*n*2+1,1);            }        }        for(int i=2; i<m; i++)            add_edge((i-1)*n+1+m*n,m*n*2+1,1);        if(n>1)        {            for(int i=2; i<m; i++)                add_edge(i*n+m*n,m*n*2+1,1);        }        int ans=max_flow(0,m*n*2+1);     //   cout<<ans;        if(ans==num)            cout<<"possible\n";        else            cout<<"not possible\n";    }    return 0;}void add_edge(int from,int to,int cap){    v[from].push_back(edge(to,cap,v[to].size()));    v[to].push_back(edge(from,0,v[from].size()-1));}int dfs(int a,int t,int f){    if(a==t)        return f;    used[a]=1;    for(int i=0; i<v[a].size(); i++)    {        edge &e=v[a][i];        if(!used[e.to]&&e.cap>0)        {            int d=dfs(e.to,t,min(f,e.cap));            if(d>0)            {                e.cap-=d;                v[e.to][e.rev].cap+=d;                return d;            }        }    }    return 0;}int max_flow(int s,int t){    int flow=0;    while(1)    {        memset(used,0,sizeof(used));        int f=dfs(s,t,INF);        if(f==0)            return flow;        flow+=f;    }}

0 0