【BZOJ1001】【BeiJing2006】狼抓兔子最大流

Description

现在小朋友们最喜欢的"喜羊羊与灰太狼",话说灰太狼抓羊不到，但抓兔子还是比较在行的，而且现在的兔子还比较笨，它们只有两个窝，现在你做为狼王，面对下面这样一个网格的地形：

左上角点为(1,1),右下角点为(N,M)(上图中N=4,M=5).有以下三种类型的道路 1:(x,y)<==>(x+1,y) 2:(x,y)<==>(x,y+1) 3:(x,y)<==>(x+1,y+1) 道路上的权值表示这条路上最多能够通过的兔子数，道路是无向的. 左上角和右下角为兔子的两个窝，开始时所有的兔子都聚集在左上角(1,1)的窝里，现在它们要跑到右下解(N,M)的窝中去，狼王开始伏击这些兔子.当然为了保险起见，如果一条道路上最多通过的兔子数为K，狼王需要安排同样数量的K只狼，才能完全封锁这条道路，你需要帮助狼王安排一个伏击方案，使得在将兔子一网打尽的前提下，参与的狼的数量要最小。因为狼还要去找喜羊羊麻烦.

Input

第一行为N,M.表示网格的大小，N,M均小于等于1000.接下来分三部分第一部分共N行，每行M-1个数，表示横向道路的权值. 第二部分共N-1行，每行M个数，表示纵向道路的权值. 第三部分共N-1行，每行M-1个数，表示斜向道路的权值. 输入文件保证不超过10M

Output

输出一个整数，表示参与伏击的狼的最小数量.

Sample Input

3 4
5 6 4
4 3 1
7 5 3
5 6 7 8
8 7 6 5
5 5 5
6 6 6

Sample Output

Solution:记得第一次看时觉得炒鸡高大上，其实现在也是，因为我并不会平面图转对偶图，一开始感觉就是裸的最大流，然而加当前弧优化后跑不粗来，估计写惨了，看了黄学长的blog用了某种优化后才A了，传—送—门

 #include <iostream>

 #include <cstdio>

 #include <cstring>

 #include <cmath>

 #include <cstdlib>

 #include <queue>

 #define N 1000000

 #define inf 0x7fffffff

 using namespace std;

 struct data{int next,point,v;}e[*N+];

 int cur[N+],deep[N+],head[*N+],q[N+];

 int tn,n,m,ans; 

 inline int get()

 {

     char c;int anss=;

     while ((c=getchar())==' ' || c=='\n' || c=='\r');

     anss=c-'';

     while (isdigit(c=getchar()))    anss=anss*+c-'';

     return anss;

 }

 inline bool bfs(int s,int t)

 {

     memset(deep,-,sizeof(deep));

 //  for (int i=0;i<=n;i++)   cur[i]=head[i];

     deep[s]=;

     int tt=,w=;

     q[tt]=s;

     while (tt<w)

     {

         int now=q[tt]; tt++;

         for (int tn=head[now];tn!=-;tn=e[tn].next)

             if (deep[e[tn].point]==- && e[tn].v)

                 deep[e[tn].point]=deep[now]+,q[w++]=e[tn].point;

     }

     if (deep[t]==-)    return ;

      return ;

 }

 int dfs(int now,int t,int limit)

 {

     if (!limit || now==t)   return limit;

     int flow=,f,used=;

     for (int tn=head[now];tn!=-;tn=e[tn].next)

     {

         //cur[now]=tn;

         if (deep[e[tn].point]==deep[now]+ && e[tn].v)

         {

             f=limit-used;

             f=dfs(e[tn].point,t,min(e[tn].v,f));

             used+=f;

             e[tn].v-=f;

             e[tn^].v+=f;

             if (used==limit)    return limit;

         }

     }

     if (!used)  deep[now]=-;

     return used;

 }

 void se(int x,int y,int w)

 {

     tn++; e[tn].next=head[x]; head[x]=tn; e[tn].point=y; e[tn].v=w;

     tn++; e[tn].next=head[y]; head[y]=tn; e[tn].point=x; e[tn].v=w;

 }

 void dinic(int s,int t)

 {

     while (bfs(s,t))    ans+=dfs(s,t,inf);

 }

 int main()

 {

     //freopen("bjrabbit.in","r",stdin);

     //freopen("bjrabbit.out","w",stdout);

     int w;

     tn=-;

     memset(e,-,sizeof(e));

     memset(head,-,sizeof(head));

     n=get(); m=get();

     for (int i=;i<=n;i++)

         for(int j=;j<m;j++)

             w=get(),se(m*(i-)+j,m*(i-)+j+,w);

     for (int i=;i<n;i++)

         for (int j=;j<=m;j++)

             w=get(),se(m*(i-)+j,m*i+j,w);

     for (int i=;i<n;i++)

         for (int j=;j<m;j++)

             w=get(),se(m*(i-)+j,m*i+j+,w);

     dinic(,n*m);

     printf("%d\n",ans);

     return ;

 }

HINT

2015.4.16新加数据一组，可能会卡掉从前可以过的程序。

Source