数学建模第五章 图与网络（二）【数学建模】.docVIP

即首先给出一个初始流，这样的流是存在的，例如零流。如果存在关于它的可增广轨，那么调整该轨上每条弧上的流量，就可以得到新的流。对于新的流，如果仍存在可增广轨，则用同样的方法使流的值增大，继续这个过程，直到网络中不存在关于新得到流的可增广轨为止，则该流就是所求的最大流。 这种方法分为以下两个过程： A.标号过程：通过标号过程寻找一条可增广轨。 B.增流过程：沿着可增广轨增加网络的流量。 这两个过程的步骤分述如下。 （A）标号过程： （i）给发点标号为。 （ii）若顶点已经标号，则对的所有未标号的邻接顶点按以下规则标号： ① 若，且时，令， 则给顶点标号为，若，则不给顶点标号。 ② ，且，令，则给标号为，若，则不给标号。 （iii）不断地重复步骤（ii）直到收点被标号，或不再有顶点可以标号为止。当被标号时，表明存在一条从到的可增广轨，则转向增流过程（B）。如若点不能被标号，且不存在其它可以标号的顶点时，表明不存在从到的可增广轨，算法结束，此时所获得的流就是最大流。 （B）增流过程 （i）令。 （ii）若的标号为），则；若的标号为，则。 （iii）若，把全部标号去掉，并回到标号过程（A）。否则，令，并回到增流过程（ii）。 求网络中的最大流的算法的程序设计具体步骤如下： 对每个节点，其标号包括两部分信息 该节点在可能的增广路中的前一个节点，以及沿该可能的增广路到该节点为止可以增广的最大流量。 STEP0 置初始可行流（如零流）；对节点标号，即令=任意正值（如1）。 STEP1 若，继续下一步；否则停止，已经得到最大流，结束。 STEP2 取消所有节点的标号，即令， ；令LIST={}，对节点标号，即令充分大的正值。 STEP3 如果LIST且，继续下一步；否则：（3a）如果已经有标号（即），则找到了一条增广路，沿该增广路对流进行增广（增广的流量为，增广路可以根据pred回溯方便地得到），转STEP1。 （3b）如果没有标号（即LIST=且），转STEP1。 STEP4 从LIST中移走一个节点；寻找从节点出发的所有可能的增广弧：（4a）对非饱和前向弧，若节点没有标号（即），对进行标号，即令 ，， 并将加入LIST中。 （4b）对非空后向弧，若节点没有标号（即），对进行标号，即令 ，， 并将加入LIST中。 例14 用Ford-Fulkerson算法计算如下网络中的最大流，每条弧上的两个数字分别表示容量和当前流量。 解 编写程序如下： clc,clear,M=1000; u(1,2)=1;u(1,3)=1;u(1,4)=2; u(2,3)=1;u(2,5)=2; u(3,5)=1; u(4,3)=3;u(4,5)=3; f(1,2)=1;f(1,3)=0;f(1,4)=1; f(2,3)=0;f(2,5)=1; f(3,5)=1; f(4,3)=1;f(4,5)=0; n=length(u); list=[]; maxf=zeros(1:n);maxf(n)=1; while maxf(n)0 maxf=zeros(1,n);pred=zeros(1,n); list=1;record=list;maxf(1)=M; while (~isempty(list))(maxf(n)==0) flag=list(1);list(1)=[]; index1=(find(u(flag,:)~=0)); label1=index1(find(u(flag,index1)... -f(flag,index1)~=0)); label1=setdiff(label1,record); list=union(list,label1); pred(label1(find(pred(label1)==0)))=flag; maxf(label1)=min(maxf(flag),u(flag,label1)... -f(flag,label1)); record=union(record,label1); label2=find(f(:,flag)~=0); label2=label2; label2=setdiff(label2,record); list=union(list,label2); pred(label2(find(pred(label2)==0)))=-flag; maxf(label2)=min(maxf(flag),f(label2,flag)); record=union(record,label2);