大流问题的标号.pptVIP

最大流问题的标号法Ford-Fulkerson增广路径算法;一、标号法的基本思路;1、标号过程;（2）若在弧（vj,vi）上，fji0,则给vj标号（-vi,L(vj)）, L(vj)=min{L(vi),fji}。这时顶点vj成为标号而未检查的顶点。 在vi的全部相邻顶点都已标号后，vi成为标号而已检查过的顶点。重复上述步骤，一旦vt被标上号，表明得到一条从vs到vt的可改进路P，转入调整过程；若所有标号都已检查过致使标号过程无法继续时，则算法结束。这时的可行流即最大流。 ;2、调整过程;并去掉所有的标号，得到新的可行流F’={fij’},重新进入标号过程，直至检查完所有标号为止。;二、标号法的算法流程;可改进路的数据类型：将当前已标号点的可改进路存储在数组path中，数据元素为可改进路上的结点，结点信息包括标号（前驱结点指针p，当前可改进量d）和弧指针w: Ptype=record p,d:integer w:link; end; Var path:array[1..maxn] of ptype; 队列：采用宽度优先搜索的方法求最大流。 Var q:array[1..maxn] of integer; op,cl:integer; {队列指针};2、构造网络D的邻接表。 每读入一条弧（u,v）的两个端点u,v及其容量c后，则在D[u]的邻接表中插入流量为0，容量为c的前向弧（u，v）；在D[v]对应的邻接表中插入流量为0，容量为0的后向弧（v,u）。这个过程可以用子过程insert(u,v,c)完成。 构造过程： 读定点数n，源点序号s和汇点序号t； For i:=1 to n do d[i]:=nil; While D网未读完 do begin 读入当前弧的两个端点u,v及其容量c； insert(u,v,c)； end;;Procedure insert(u,v,c:integer); Var x:link; Begin x:=d[u]; While (xnil) and(x^.kv) do x:=x^.next; If xnil then x^.c:=c Else begin new(x); x^.k:=v; x^.c:=c; x^.f:=0; x^.next:=d[u]; d[u]:=x; new(x^.g); x^.g^.k:=u; x^.g^.c:=0; x^.g^.f:=0; x^.g^.g:=x; x^.g^.next:=d[v]; d[v]:=x^.g; end;;3、采用宽度优先搜索求最大流。 Procedure ford; Var u,a:integer; x:link; Begin repeat fillchar(path,sizeof(path),0); path[s].d:=maxint; path[s].p=s; cl:=0; op:=1; q[op]:=s; while (clop) and(path[t].p=0) do begin inc(cl); u:=q[cl]; x:=d[u]; while xnil do begin a:=x^.c-x^.f; if (path[x^.k].p=0) and (a0) then begin inc(op); q[op]:=x^.k; path[x^.k].p:=u; path[x^.k].w:=x; if path[u].da then path[x^.k].d:=path[u].d else path[x^.k].d:=a; end; x:=x^.next; end; end;;If path[t].p0 then Begin u:=t; a:=path[u].d; repeat path[u].w^.f:=path[u].w^.f+a; path[u].w^.g^.f:=-path[u].w^.f; u:=path[u].p; until u=s; End Until path[t].p=0; End; ;4、输出