网络流之isap.docVIP

program isap; const maxm 1020;maxn 500; var n,m,i,j,k,x,y,z,s,t,sum,tot:longint; next,other,u,v,w:array[0..maxm]of longint; point,last,d,num:array[0..maxn]of longint; d为距离标号，表示其到t的距离；num统计d函数的每个值是多少，如果某次修改后出现了num[i] 0，则说明出现了断层，可以直接退出 function min a,b:longint :longint; begin if a b then exit a ;exit b ; end; function dfs now,flow:longint :longint; //now节点可提供的流量为flow，返回值为最从now中流出的流。 var j,p,ans:longint; begin if now t then exit flow ; ans: 0;j: last[now]; while j 0 do begin if w[j] 0 and d[v[j]]+1 d[u[j]] then begin last[now]: j;//当前弧优化，下次直接从此开始寻找 p: dfs v[j],min flow-ans,w[j] ; dec w[j],p ; inc w[other[j]],p ; inc ans,p ; if ans flow then exit ans ;//已经满流就退出。 if d[s] tot then exit ans ; end; j: next[j]; end; dec num[d[now]] ; if num[d[now]] 0 then d[s]: tot;//GAP优化。 inc d[now] ;//寻找新的可增广路 inc num[d[now]] ; last[now]: point[now];//重置last[now] exit ans ; end; begin readln n,m ;s: 1;t: n;tot: n; for i: 1 to m do begin readln x,y,z ; inc k ;u[k]: x;v[k]: y;w[k]: z; next[k]: point[x];point[x]: k;other[k]: k+1; inc k ;u[k]: y;v[k]: x;w[k]: 0; next[k]: point[y];point[y]: k;other[k]: k-1; //u,v,w分别记录起点、终点和残余容量，other记录它的反向边。 end; fillchar d,sizeof d ,0 ; fillchar num,sizeof num ,0 ; num[0]: tot; sum: 0; for i: 1 to tot do last[i]: point[i];//一开始要从第一条边开始 while d[s] tot do inc sum,dfs s,maxlongint ; //d[s] tot就一定不存在增广路了 writeln sum ; end.