最小费用流spfa算法模板（pascal）

以前写过，现在的码风与以前有些变化，主要是用数组模拟邻接表存图，以前是用指针存图。

以前的博文:http://www.cnblogs.com/Currier/p/6387732.html

洛谷可评测。

传送门：https://www.luogu.org/problem/show?pid=3381

program rrr(input,output);
const
  inf=123456789;
type
  etype=record
     t,c,w,next,rev:longint;
  end;
var
  e:array[0..100010]of etype;
  a,fre,frv,q,dis:array[0..5050]of longint;
  inq:array[0..5050]of boolean;
  n,m,s,t,i,j,x,y,c,w,max,ans,cnt,f:longint;
procedure add(x,y,c,w:longint);
begin
   inc(cnt);e[cnt].t:=y;e[cnt].c:=c;e[cnt].w:=w;e[cnt].next:=a[x];a[x]:=cnt;
end;
function min(a,b:longint):longint;
begin
   if a<b then exit(a) else exit(b);
end;
procedure spfa;
var
  h,t:longint;
begin
   for i:=1 to n do dis[i]:=inf;
   fillchar(inq,sizeof(inq),false);
   h:=0;t:=1;q[1]:=s;dis[s]:=0;inq[s]:=true;
   while h<>t do
      begin
         inc(h);if h>5000 then h:=1;
         i:=a[q[h]];
         while i<>0 do
            begin
               if (e[i].c>0) and (dis[q[h]]+e[i].w<dis[e[i].t]) then
                  begin
                     dis[e[i].t]:=dis[q[h]]+e[i].w;
                     frv[e[i].t]:=q[h];fre[e[i].t]:=i;
                     if not inq[e[i].t] then
                        begin
                           inc(t);if t>5000 then t:=1;
                           q[t]:=e[i].t;inq[e[i].t]:=true;
                        end;
                  end;
               i:=e[i].next;
            end;
         inq[q[h]]:=false;
      end;
end;
begin
   assign(input,'r.in');assign(output,'r.out');reset(input);rewrite(output);
   readln(n,m,s,t);
   fillchar(a,sizeof(a),0);cnt:=0;
   for i:=1 to m do
      begin
         readln(x,y,c,w);
         add(x,y,c,w);e[cnt].rev:=cnt+1;
         add(y,x,0,-w);e[cnt].rev:=cnt-1;
      end;
   max:=0;ans:=0;
   while true do
      begin
         spfa;
         if dis[t]=inf then break;
         j:=t;f:=inf;
         while j<>s do begin f:=min(f,e[fre[j]].c);j:=frv[j]; end;
         max:=max+f;
         j:=t;w:=0;
         while j<>s do begin w:=w+e[fre[j]].w;dec(e[fre[j]].c,f);inc(e[e[fre[j]].rev].c,f);j:=frv[j]; end;
         ans:=ans+w*f;
      end;
   write(max,' ',ans);
   close(input);close(output);
end.