算法模板——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[..] of point;
    c,d:array[..] 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:=;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),);
               d[]:=s;f:=;r:=;c[s]:=;
               while f<r do
                     begin
                          p:=a[d[f]];
                          while p<>nil do
                                begin
                                     if (p^.w<>) and (c[p^.g]=-) then
                                        begin
                                             c[p^.g]:=c[d[f]]+;
                                             e[p^.g]:=p;
                                             d[r]:=p^.g;inc(r);
                                        end;
                                     p:=p^.next;
                                end;
                          inc(f);
                     end;
               exit(c[t]<>-);
          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:= to n do a[i]:=nil;
      for i:= to m do
          begin
               readln(j,k,l);
               add(j,k,l);
          end;
      flow:=;while spfa do calc;
      writeln(flow);
      readln;
 end.