【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了，传—送—门

 1  
 2 #include <iostream>
 3 #include <cstdio>
 4 #include <cstring>
 5 #include <cmath>
 6 #include <cstdlib>
 7 #include <queue>
 8 #define N 1000000
 9 #define inf 0x7fffffff 
10 using namespace std;
11  
12 struct data{int next,point,v;}e[6*N+10];
13 int cur[N+10],deep[N+10],head[2*N+10],q[N+10];
14 int tn,n,m,ans; 
15  
16 inline int get()
17 {
18     char c;int anss=0;
19     while ((c=getchar())==' ' || c=='\n' || c=='\r');
20     anss=c-'0';
21     while (isdigit(c=getchar()))    anss=anss*10+c-'0';
22     return anss;
23 }
24  
25 inline bool bfs(int s,int t)
26 {
27     memset(deep,-1,sizeof(deep));
28 //  for (int i=0;i<=n;i++)   cur[i]=head[i];
29     deep[s]=0;
30     int tt=0,w=1;
31     q[tt]=s;
32     while (tt<w)
33     {
34         int now=q[tt]; tt++;
35         for (int tn=head[now];tn!=-1;tn=e[tn].next)
36             if (deep[e[tn].point]==-1 && e[tn].v)
37                 deep[e[tn].point]=deep[now]+1,q[w++]=e[tn].point;   
38     }
39     if (deep[t]==-1)    return 0;
40      return 1;
41 }
42  
43 int dfs(int now,int t,int limit)
44 {
45     if (!limit || now==t)   return limit;
46     int flow=0,f,used=0;
47     for (int tn=head[now];tn!=-1;tn=e[tn].next)
48     {
49         //cur[now]=tn;
50         if (deep[e[tn].point]==deep[now]+1 && e[tn].v)
51         {
52             f=limit-used;
53             f=dfs(e[tn].point,t,min(e[tn].v,f));
54             used+=f;
55             e[tn].v-=f;
56             e[tn^1].v+=f;
57             if (used==limit)    return limit;
58         }
59     }
60     if (!used)  deep[now]=-1;
61     return used;
62 }
63  
64 void se(int x,int y,int w)
65 {
66     tn++; e[tn].next=head[x]; head[x]=tn; e[tn].point=y; e[tn].v=w;
67     tn++; e[tn].next=head[y]; head[y]=tn; e[tn].point=x; e[tn].v=w;
68 }
69  
70 void dinic(int s,int t)
71 {
72     while (bfs(s,t))    ans+=dfs(s,t,inf);
73 }
74  
75 int main()
76 {
77     //freopen("bjrabbit.in","r",stdin);
78     //freopen("bjrabbit.out","w",stdout);
79     int w;
80     tn=-1;
81     memset(e,-1,sizeof(e));
82     memset(head,-1,sizeof(head));
83     n=get(); m=get();
84     for (int i=1;i<=n;i++)
85         for(int j=1;j<m;j++)
86             w=get(),se(m*(i-1)+j,m*(i-1)+j+1,w);
87     for (int i=1;i<n;i++)
88         for (int j=1;j<=m;j++)
89             w=get(),se(m*(i-1)+j,m*i+j,w);
90     for (int i=1;i<n;i++)
91         for (int j=1;j<m;j++)
92             w=get(),se(m*(i-1)+j,m*i+j+1,w);
93     dinic(1,n*m);
94     printf("%d\n",ans);
95     return 0;
96 }
View Code

 

14

HINT

 

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

 

Source

posted @ 2016-02-28 19:00  DMoon  阅读(209)  评论(0)  编辑  收藏  举报