bzoj1001 [BeiJing2006]狼抓兔子

1001: [BeiJing2006]狼抓兔子

Time Limit: 15 Sec Memory Limit: 162 MB
Submit: 19687 Solved: 4870
[Submit][Status][Discuss]

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

题解

平面图最大流……周冬论文里说的很清楚了……贴一下代码吧……

  1 /**************************************************************
  2     Problem: 1001
  3     User: 1090900715
  4     Language: Pascal
  5     Result: Accepted
  6     Time:7744 ms
  7     Memory:103760 kb
  8 ****************************************************************/
  9  
 10 program j01;
 11 const maxn=2000086;
 12 var b:array[0..1000,0..1000,0..1]of longint;
 13     q,next,data:array[0..6000086]of longint;
 14     dis,head,l:array[0..maxn]of longint;
 15     inl:array[0..maxn]of boolean;
 16     n,m,x,i,j,tt,tot:longint;
 17  
 18 procedure add(u,v,w:longint);
 19 begin
 20   inc(tt);
 21   q[tt]:=v;
 22   next[tt]:=head[u];
 23   head[u]:=tt;
 24   data[tt]:=w;
 25 end;
 26  
 27 procedure spfa;
 28 var h,tail,i,j:longint;
 29 begin
 30   fillchar(dis,sizeof(dis),$3f);
 31   fillchar(inl,sizeof(inl),0);
 32   dis[0]:=0;h:=0;tail:=1;l[1]:=0;inl[0]:=true;
 33   while h<>tail do
 34   begin
 35     inc(h);if h>maxn then h:=1;
 36     i:=l[h];
 37     j:=head[i];
 38     while j>0 do
 39     begin
 40       if dis[i]+data[j]<dis[q[j]] then
 41       begin
 42         dis[q[j]]:=dis[i]+data[j];
 43         if inl[q[j]]=false then
 44         begin
 45           inc(tail);if tail>maxn then tail:=1;
 46           l[tail]:=q[j];
 47           inl[q[j]]:=true;
 48         end;
 49       end;
 50       j:=next[j];
 51     end;
 52     inl[i]:=false;
 53   end;
 54 end;
 55  
 56 begin
 57   readln(n,m);
 58   tot:=0;
 59   for i:=1 to n-1 do
 60     for j:=1 to m-1 do
 61     begin
 62       inc(tot);
 63       b[i,j,0]:=tot;
 64       inc(tot);
 65       b[i,j,1]:=tot;
 66     end;
 67   inc(tot);
 68   fillchar(head,sizeof(head),0);
 69   tt:=0;
 70   for i:=1 to n do
 71     for j:=1 to m-1 do
 72     begin
 73       read(x);
 74       if i=1 then add(b[i,j,1],tot,x);
 75       if i=n then
 76         add(0,b[i-1,j,0],x);
 77       if (i<>1)and(i<>n) then
 78       begin
 79         add(b[i,j,1],b[i-1,j,0],x);
 80         add(b[i-1,j,0],b[i,j,1],x);
 81       end;
 82     end;
 83   for i:=1 to n-1 do
 84     for j:=1 to m do
 85     begin
 86       read(x);
 87       if j=1 then
 88         add(0,b[i,j,0],x);
 89       if j=m then add(b[i,j-1,1],tot,x);
 90       if (j<>1)and(j<>m) then
 91       begin
 92         add(b[i,j,0],b[i,j-1,1],x);
 93         add(b[i,j-1,1],b[i,j,0],x);
 94       end;
 95     end;
 96   for i:=1 to n-1 do
 97     for j:=1 to m-1 do
 98     begin
 99       read(x);
100       add(b[i,j,0],b[i,j,1],x);
101       add(b[i,j,1],b[i,j,0],x);
102     end;
103   spfa;
104   writeln(dis[tot]);
105 end.
106

View Code