【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[6*N+10];
int cur[N+10],deep[N+10],head[2*N+10],q[N+10];
int tn,n,m,ans; 
 
inline int get()
{
    char c;int anss=0;
    while ((c=getchar())==' ' || c=='
' || c=='
');
    anss=c-'0';
    while (isdigit(c=getchar()))    anss=anss*10+c-'0';
    return anss;
}
 
inline bool bfs(int s,int t)
{
    memset(deep,-1,sizeof(deep));
//  for (int i=0;i<=n;i++)   cur[i]=head[i];
    deep[s]=0;
    int tt=0,w=1;
    q[tt]=s;
    while (tt<w)
    {
        int now=q[tt]; tt++;
        for (int tn=head[now];tn!=-1;tn=e[tn].next)
            if (deep[e[tn].point]==-1 && e[tn].v)
                deep[e[tn].point]=deep[now]+1,q[w++]=e[tn].point;   
    }
    if (deep[t]==-1)    return 0;
     return 1;
}
 
int dfs(int now,int t,int limit)
{
    if (!limit || now==t)   return limit;
    int flow=0,f,used=0;
    for (int tn=head[now];tn!=-1;tn=e[tn].next)
    {
        //cur[now]=tn;
        if (deep[e[tn].point]==deep[now]+1 && 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^1].v+=f;
            if (used==limit)    return limit;
        }
    }
    if (!used)  deep[now]=-1;
    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=-1;
    memset(e,-1,sizeof(e));
    memset(head,-1,sizeof(head));
    n=get(); m=get();
    for (int i=1;i<=n;i++)
        for(int j=1;j<m;j++)
            w=get(),se(m*(i-1)+j,m*(i-1)+j+1,w);
    for (int i=1;i<n;i++)
        for (int j=1;j<=m;j++)
            w=get(),se(m*(i-1)+j,m*i+j,w);
    for (int i=1;i<n;i++)
        for (int j=1;j<m;j++)
            w=get(),se(m*(i-1)+j,m*i+j+1,w);
    dinic(1,n*m);
    printf("%d
",ans);
    return 0;
} 

14

HINT

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

Source