算法模板——Dinic网络最大流 2

实现功能：同Dinic网络最大流 1

这个新的想法源于Dinic费用流算法。。。

在费用流算法里面，每次处理一条最短路，是通过spfa的过程中就记录下来，然后顺藤摸瓜处理一路

于是在这个里面我的最大流也采用这种模式，这样子有效避免的递归，防止了爆栈么么哒

type
point=^node;
node=record
g,w:longint;
next,anti:point;
end;
var
i,j,k,l,m,n,s,t,flow:longint;
a,e:array[0..10000] of point;
c,d:array[0..10000] of longint;
function min(x,y:longint):longint;
begin
if x<y then min:=x else min:=y;
end;
procedure add(x,y,z:longint);
var p:point;
begin
new(p);p^.g:=y;p^.w:=z;p^.next:=a[x];a[x]:=p;
new(p);p^.g:=x;p^.w:=0;p^.next:=a[y];a[y]:=p;
a[x]^.anti:=a[y];a[y]^.anti:=a[x];
end;
function spfa:boolean;
var p:point;f,r:longint;
begin
fillchar(c,sizeof(c),255);
d[1]:=s;f:=1;r:=2;c[s]:=0;
while f<r do
begin
p:=a[d[f]];
while p<>nil do
begin
if (p^.w<>0) and (c[p^.g]=-1) then
begin
c[p^.g]:=c[d[f]]+1;
e[p^.g]:=p;
d[r]:=p^.g;inc(r);
end;
p:=p^.next;
end;
inc(f);
end;
exit(c[t]<>-1);
end;
procedure calc;   //“顺藤摸瓜”模式有效避免了递归带来的爆栈隐患
begin
i:=t;l:=maxlongint;
while i<>s do
begin
l:=min(l,e[i]^.w);
i:=e[i]^.anti^.g;
end;
i:=t;inc(flow,l);
while i<>s do
begin
if e[i]^.w<>maxlongint then dec(e[i]^.w,l);
if e[i]^.anti^.w<>maxlongint then inc(e[i]^.anti^.w,l);
i:=e[i]^.anti^.g;
end;
end;
begin
readln(m,n,s,t);
for i:=1 to n do a[i]:=nil;
for i:=1 to m do
begin
readln(j,k,l);
add(j,k,l);
end;
flow:=0;while spfa do calc;
writeln(flow);
readln;
end.