求网络最大流及最小费用最大流问题的 Ford-Fulkerson标号算法.docVIP

最优化方法上机实验3 求网络最大流及最小费用最大流Ford-Fulkerson标号算法 上机时间：2014.01.07  实验 问题 描述 1例 求下图所示网络的最大流 2例求下图所示网络的最小费用最大流，弧旁中的分别表示顶点i到顶点j间弧的费用和容量。 b=[0 10 6 3 0 0; 0 0 4 0 0 3; 0 0 0 0 12 4; 0 0 5 0 2 0; 0 0 0 0 0 5; 0 0 0 0 0 0]; c=[0 12 7 8 0 0; 0 0 8 0 0 5; 0 0 0 0 10 7; 0 0 9 0 7 0; 0 0 0 0 0 8; 0 0 0 0 0 0]; [f wf zwf]=BGf(c,b) 原理 及 算法 算法原理 步骤： STEP1 若maxf(t) 0，继续下一步；否则停止，已经得到最大流，结束。 STEP2 取消所有节点 j∈V 的标号，即令maxf(j) = 0， pred(j) = 0；令 LIST={ s }，对节点s标号，即令maxf(s) =充分大的正值。 STEP3 如果 LIST≠ Φ且maxf(t) = 0，继续下一步；否则： （3a）如果t已经有标号（即maxf(t) 0），则找到了一条增广路，沿该增广路对流x进行增广（增广的流量为maxf(t)，增广路可以根据 pred 回溯方便地得到），转 STEP1。 （3b）如果t没有标号（即 LIST=Φ且maxf(t) = 0），则停止，已得到最大流。 STEP4 从 LIST 中移走一个节点i；寻找从节点i出发的所有可能的增广弧： （4a）对非饱和前向弧(i, j)，若节点 j没有标号（即pred(j) = 0），对 j进行标号，并将 j加入 LIST 中。 （4b）对非空后向弧(j,i)，若节点 j没有标号（即pred(j) = 0），对 j进行标号，并将 j 加入 LIST 中。 完 整 程 序 清 单 （ 含 注 释 ） 3.1Ford-Fulkerson标号算法 clc,clear u=[0 5 4 3 0 0 0 0; 0 0 0 0 5 3 0 0; 0 0 0 0 0 3 2 0; 0 0 0 0 0 0 2 0; 0 0 0 0 0 0 0 4; 0 0 0 0 0 0 0 3; 0 0 0 0 0 0 0 5; 0 0 0 0 0 0 0 0]; n=length(u); f=zeros(n,n); list=[];maxf(n)=1; while maxf(n)0 maxf=zeros(1,n);pred=zeros(1,n); list=1;record=list;maxf(1)=inf; %list是未检查邻接点的标号点，record是已标号点 while (~isempty(list))(maxf(n)==0) flag=list(1);list(1)=[]; label1= find(u(flag,:)-f(flag,:)); label1=setdiff(label1,record); list=union(list,label1); pred(label1)=flag; maxf(label1)=min(maxf(flag),u(flag,label1)... -f(flag,label1)); record=union(record,label1); label2=find(f(:,flag)); label2=label2; label2=setdiff(label2,record); list=union(list,label2); pred(label2)=-flag; maxf(label2)=min(maxf(flag),f(label2,flag)); record=union(record,label2); end if maxf(n)0 v2=n; v1=pred(v2); while v2~=1 if v10 f(v1,v2)=f(v1,v2)+maxf(n); else v1=abs(v1); f(v2,v1)=f(v2,v1)-maxf(n); end v2=v1; v1=pred(v2); end end end f 3.2最小费用最大流问题程序： function[f wf zwf]=BGf(C,b) %C表示弧容量矩阵 %b表示弧上单位流量的费用 %f表示最小费用最大流矩阵 %wf表示最大流量 %zwf表示最小费用 n=size(C,2); wf=0;wf0=Inf; %wf表示