网络流——增广路算法（dinic）模板 [BeiJing2006]狼抓兔子

 #include<iostream>
 #include<cstring>
 #include<algorithm>
 #include<cmath>
 #include<cstdio>
 #include<queue>
 using namespace std;
 struct data
 {
     int from,to,next,cup;
     data(){from=-,to=-,next=-,cup=-;}
 }e[];
 int cnt=,head[];
 int dis[];
 bool vis[];
 void add(int u,int v,int w)
 {e[cnt].from=u,e[cnt].to=v,e[cnt].next=head[u],head[u]=cnt,e[cnt].cup=w;cnt++;}
 int s,t;
 bool bfs()
 {
     memset(dis,-,sizeof(dis));
     int f=;
     queue<int> q;
     q.push(s);
     dis[s]=;
     while(!q.empty())
     {
         int n=q.front();
         q.pop();
         for(int i=head[n];i>=;i=e[i].next)
         {
             if(dis[e[i].to]==-&&e[i].cup>)
             {
                 dis[e[i].to]=dis[n]+;
                 q.push(e[i].to);
             }
         }
     }
     return dis[t]!=-;
 }
 int dfs(int n,int a)
 {

     if(n==t||a==) return a;
     int flow=,f;
     for(int i=head[n];i>=;i=e[i].next)
     {
         if(dis[e[i].to]==dis[n]+&&e[i].cup>)
         {
             f=dfs(e[i].to,min(a,e[i].cup));
             e[i].cup-=f;
             e[i^].cup+=f;
             flow+=f;
             a-=f;
             if(a==) break;
         }
     }
     if(!flow) dis[n]=-;
     return flow;
 }
 int dinic()
 {
     int sum=;
     while(bfs()) sum+=dfs(s,);
     return sum;
 }
 int main()
 {
     memset(head,-,sizeof(head));
     int n,m;
     scanf("%d%d",&n,&m);
     int x;
     for(int i=;i<=n;i++)
         for(int j=;j<m;j++)
         {
             scanf("%d",&x);
             add(m*(i-)+j,m*(i-)+j+,x);
             add(m*(i-)+j+,m*(i-)+j,x);
         }
     for(int i=;i<n;i++)
         for(int j=;j<=m;j++)
         {
             scanf("%d",&x);
             add(m*(i-)+j,m*(i)+j,x);
             add(m*(i)+j,m*(i-)+j,x);
         }
     for(int i=;i<n;i++)
         for(int j=;j<m;j++)
         {
             scanf("%d",&x);
             add(m*(i-)+j,m*(i)+j+,x);
             add(m*(i)+j+,m*(i-)+j,x);
         }
     s=,t=m*n;
     printf("%d\n",dinic());
 }