求最大流标号法的描述.docVIP

求最大流标号法的描述 ⑴数据结构 Const maxn=网络顶点数 Type node=record {可增广道路的顶点类型} l,p:integer； {第一标号，检查标号} end; arc=record {网络边类型} c,f:integer; {容量，流量} end; gtype=array[0..maxn,0..maxn] of arc； {网矩阵} ltype=array[0..maxn] of node; {可增广道路} var lt:ltype; {可增广道路} g:gtype; {网络} n,s,t:integer; {顶点数，源点，汇点} ⑵ 主要算法 初始化网络，可增广道路procedure read-graph； var I,j:integer;begin readln(n); {顶点数} fillchar(g,sizeof(g),0); {初始化网络} fillchar(lt,sizeof(lt),0); {初始化可增广道路} for i:=1 to n do for j:=1 to n do read(g[I,j].c);end; 寻找已标号而来检查的顶点序号function check:integer; var i:integer;begin i:=1; while (i=n)and not((lt[i].l0)and(lt[i].p=0)) do inc(i); if in then check:=0 {顶点不存在} else check:=;end; 标号过程，并返回是否找到可增广道路及改进量afunction ford(var a:integer):boolean; {无增广道路返回true} var I,j,m,x:integer;begin ford:=true; fillchar(lt,sizeof(lt),0); {去掉原来的标号} lt[s].l:=s; {从Vs开始} repeat i:=check; {寻找一个已标号而未检查的顶点i} if i=0 then exit; {若该顶点不存在，则退出过程，返回true} for j:=1 to n do if (lt[j].l=0)and((g[I,j].c0)or(g[j,i].c0)  {寻找与Vi相邻且标号的顶点j} then begin if (g[I,j].fg[I,j].c) then lt[j].l:=I; if (g[j,i].f0) then lt[j].l:=-I; end; lt[i].p:=1; {顶点Vi置已检查标志} until (lt[t].l0) {直到汇点Vt标号为止} m:=t; a:=maxint; {从Vt推，改进量a赋初值} repeat {求a} j:=m; m:=abs(lt[t].l); if lt[j].l0 then x:=g[j,m].f; if lt[j].l0 then x:=g[m,j].c-g[m,j].f; if ax then a:=x; until m=s; {直至例推到顶点Vs为止} ford:=false; {返回可增广道路存在标志false}end; 修正流量procedure fulkerson(a:integer); var m,j:integer;begin m:=t; {以顶点Vt出发，逆向沿可增广道路修正容量} repeat j:=m; m:=abs(lt[j].l); if lt[j].l0 then g[j,m].f:=g[j,m].f-a; {后向边p-} if lt[j].l0 then g[m,j].f:=g[m,j].f+a; {向边p+}until m=s;end; 输出所有弧的流量procedure answer; var I,j:integer;begin for i:=1 to n do for j:=1 to n do writeln(