[BeiJing2006]狼抓兔子

Time Limit: 15 Sec    Memory Limit: 162 MB

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

14

 

 

 

 

这就是一道SAP求最小割的。直接套模板。不过一开始要先bfs一遍求出标号。不然会超时。

第一次加了当前弧优化。并把流量记在边上，又省了空间。修改时用异或运算。建立链表时从数组的第二个存储单位开始，双向边放在一起。修改时就是i和i xor 1两条边一起修改。

一开始一直TLE最后检查的时候发现竟然是因为记录标号个数的数组vd在bfs后又每个都加了1。。。囧了。这样间隙优化就不一定能跑出来了。。。必然超时。。。

 

AC CODE

program hy_1001;

var g,line,next: array [ 1..8000010 ] of longint ;

     vd,d,head,hs,en,q: array [ 0..1000010 ] of longint ;

     i,tot,n,m,l,r,flow,ans,z: longint ;

//============================================================================

{procedure ins(x,y,z:longint);

begin

   inc(tot); line[tot]:=y; g[tot]:=z;

   if head[x]=0 then

   begin

     head[x]:=tot;

     en[x]:=tot;

   end else

   begin

     next[en[x]]:=tot;

     en[x]:=tot;

   end;

end;}                 //旧式链表要开head数组。这次用了zc神犇的简化链表。太简洁了~~~~

//===========================================================================

procedure ins(x,y,z: longint );       //zc神犇的简化链表！

begin

   inc(tot); line[tot]:=y; g[tot]:=z;

   next[tot]:=en[x]; en[x]:=tot;

end ;

//============================================================================

procedure init;

var i,j: longint ;

begin

   readln(n,m); tot:= 1 ;

   for i:= 1 to n do

     for j:= 1 to m- 1 do

     begin

       read(z);

       ins((i- 1 )*m+j,(i- 1 )*m+j+ 1 ,z);

       ins((i- 1 )*m+j+ 1 ,(i- 1 )*m+j,z);

     end ;

   for i:= 1 to n- 1 do

     for j:= 1 to m do

     begin

       read(z);

       ins((i- 1 )*m+j,i*m+j,z);

       ins(i*m+j,(i- 1 )*m+j,z);

     end ;

   for i:= 1 to n- 1 do

     for j:= 1 to m- 1 do

     begin

       read(z);

       ins((i- 1 )*m+j,i*m+j+ 1 ,z);

       ins(i*m+j+ 1 ,(i- 1 )*m+j,z);

     end ;

   n:=n*m;

end ;

//============================================================================

procedure bfs;

var i: longint ;

begin

   d[n]:= 0 ;

   q[ 1 ]:=n; l:= 1 ; r:= 1 ;

   while l<=r do

   begin

     inc(vd[d[q[l]]]); i:=en[q[l]];

     while i<> 0 do

     begin

       if (d[line[i]]= 0 ) and (line[i]<q[l]) then

       begin

         d[line[i]]:=d[q[l]]+ 1 ;

         inc(r); q[r]:=line[i];

       end ; i:=next[i];

     end ; inc(l);

   end ; 

  //for i:= 1 to n do inc(vd[d[i]]);     //万恶的一步、、

end ;

//============================================================================

function dfs(u,flow: longint ): longint ;

var i,v,tmp,rest: longint ;

begin

   if u=n then exit(flow);

   dfs:= 0 ;

   i:=hs[u];

   while i<> 0 do

   begin

     v:=line[i];

     if (g[i]> 0 ) and (d[u]=d[v]+ 1 ) then

     begin

       rest:=flow-dfs;

       if g[i]<rest then rest:=g[i];

       tmp:=dfs(v,rest);

       dec(g[i],tmp); inc(g[i xor 1 ],tmp);

       dfs:=dfs+tmp;

       if dfs=flow then exit(dfs);

     end ; i:=next[i]; hs[u]:=next[hs[u]];

   end ;

   if d[ 1 ]>=n then exit;

   dec(vd[d[u]]);

   if vd[d[u]]= 0 then d[ 1 ]:=n;

   d[u]:=d[u]+ 1 ; inc(vd[d[u]]);

   hs[u]:=en[u];

end ;

//============================================================================

begin

   init;

   bfs;

   for i:= 1 to n do hs[i]:=en[i];

   while d[ 1 ]<n do ans:=ans+dfs( 1 ,maxlongint);    //跑SAP。

   writeln (ans);

end .