狼抓兔子（平面图转对偶图求最短路）

一个网格图中求最小割，很特殊的是这个图是一个平面图，最小割=最大流=对偶图中最短路，把小三角换成点标号，重连边即可。

我标号用的矩阵元素标号的思想，可能看着还算清楚，乘来乘去的。

数据里还有n == 1和m == 1的数据。。。。。

插个图片帮助理解。

View Code

  1 program beijing2006(input,output);
  2 const
  3    oo = 999999999;
  4 type
  5    node    = ^link;
  6    link    = record
  7          goal,w : longint;
  8          next   : node;
  9       end;
 10 var
 11    l           : array[0..2500000] of node;
 12    d           : array[0..2500000] of int64;
 13    v           : array[0..2500000] of boolean;
 14    q           : array[0..2500000] of longint;
 15    n,m,sum     : longint;
 16    start,final : longint;
 17 procedure add(xx,yy,ww: longint );
 18 var
 19    tt : node;
 20 begin
 21    new(tt);
 22    tt^.goal:=yy;
 23    tt^.w:=ww;
 24    tt^.next:=l[xx];
 25    l[xx]:=tt;
 26 end; { add }
 27 function check():boolean;
 28 var
 29    i,x      : longint;
 30    answer : int64;
 31 begin
 32    answer:=oo;
 33    if n=1 then
 34    begin
 35       for i:=1 to m-1 do
 36       begin
 37      read(x);
 38      if x<answer then
 39         answer:=x;
 40       end;
 41       writeln(answer);
 42       exit(true);
 43    end;
 44    if m=1 then
 45    begin
 46       for i:=1 to n-1 do
 47       begin
 48      read(x);
 49      if x<answer then
 50         answer:=x;
 51       end;
 52       writeln(answer);
 53       exit(true);
 54    end;
 55    exit(false);
 56 end; { check }
 57 procedure init;
 58 var
 59    i : longint;
 60 begin
 61    readln(n,m);
 62    for i:=0 to 2*m*n+3 do
 63       l[i]:=nil;
 64    sum:=(n-1)*(m-1);
 65    Start:=0;
 66    final:=(m-1)+(m-1)*2*(n-1)+1;
 67 end; { init }
 68 procedure make_graph;
 69 var
 70    i,j : longint;
 71    ww  : longint;
 72 begin
 73    for i:=1 to n do
 74       for j:=1 to m-1 do
 75       begin
 76      read(ww);
 77      if i=1 then
 78      begin
 79         add(final,(i*(m-1)+j)+sum,ww);
 80         add((i*(m-1)+j)+sum,final,ww);
 81      end
 82      else
 83         if i=n then
 84         begin
 85            add((i-1)*(m-1)+j,start,ww);
 86            add(start,(i-1)*(m-1)+j,ww);
 87         end
 88         else
 89         begin
 90            add((i-1)*(m-1)+j,i*(m-1)+j+sum,ww);
 91            add(i*(m-1)+j+sum,(i-1)*(m-1)+j,ww);
 92         end;
 93       end;
 94    for i:=1 to n-1 do
 95       for j:=1 to m do
 96       begin
 97      read(ww);
 98      if j=1 then
 99      begin
100         add(start,i*(m-1)+j,ww);
101         add(i*(m-1)+j,start,ww);
102      end
103      else
104         if j=m then
105         begin
106            add(final,i*(m-1)+j-1+sum,ww);
107            add(i*(m-1)+j-1+sum,final,ww);
108         end
109         else
110         begin
111            add(i*(m-1)+j,i*(m-1)+j-1+sum,ww);
112            add(i*(m-1)+j-1+sum,i*(m-1)+j,ww);
113         end;
114       end;
115    for i:=1 to n-1 do
116       for j:=1 to m-1 do
117       begin
118      read(ww);
119      add(i*(m-1)+j,i*(m-1)+j+sum,ww);
120      add(i*(m-1)+j+sum,i*(m-1)+j,ww);
121       end;
122 end; { make_graph }
123 procedure spfa();
124 var
125    head,tail : longint;
126    t         : node;
127 begin
128    fillchar(v,sizeof(v),false);
129    fillchar(d,sizeof(d),63);
130    head:=0;
131    tail:=1;
132    v[start]:=true;
133    q[1]:=start;
134    d[start]:=0;
135    while head<>tail do
136    begin
137       inc(head);
138       if head=2500000 then
139      head:=0;
140       v[q[head]]:=false;
141       t:=l[q[head]];
142       while t<>nil do
143       begin
144      if d[q[head]]+t^.w<d[t^.goal] then
145      begin
146         d[t^.goal]:=d[q[head]]+t^.w;
147         if not v[t^.goal] then
148         begin
149            inc(tail);
150            if tail=2500000 then
151           tail:=0;
152            q[tail]:=t^.goal;
153            v[t^.goal]:=true;
154         end;
155      end;
156      t:=t^.next;
157       end;
158    end;
159 end; { spfa }
160 procedure print;
161 begin
162    writeln(d[final]);
163 end; { print }
164 begin
165    init;
166    if not check() then
167    begin
168       make_graph();
169       spfa();
170       print;
171    end;
172 end.

 

 

 

转载于:https://www.cnblogs.com/neverforget/archive/2012/04/16/2452190.html