Ford-Fulkerson算法计算如下网络中的最大流

晒个小太阳。

用Ford-Fulkerson算法计算如下网络中的最大流，每条弧上的两个数字分别表示容量和当前流量。

编写程序如下：
! h6 Y8 k9 T2 ?9 g2 f$ c8 `; g9 Vclc,clear,M=1000;
4 v; B& |6 G% u4 ?u(1,2)=1;u(1,3)=1;u(1,4)=2;
/ O- ~! w& O% |/ T' Gu(2,3)=1;u(2,5)=2;
u(3,5)=1;
/ e7 G" b- Y, Q2 Z+ qu(4,3)=3;u(4,5)=3;
4 g& `1 N# t  of(1,2)=1;f(1,3)=0;f(1,4)=1;
0 o; P$ c! k4 Yf(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
) S! {9 ?/ q  i. e3 T6 ~, O0 B   maxf=zeros(1,n);pred=zeros(1,n);
   list=1;record=list;maxf(1)=M;
4 j( U. X- y' f5 h# }9 S* {   while (~isempty(list))&(maxf(n)==0)
4 y4 c, H# j1 A* u+ \      flag=list(1);list(1)=[];
      index1=(find(u(flag,:)~=0));
      label1=index1(find(u(flag,index1)...
      -f(flag,index1)~=0));
7 t1 M1 ^) z8 e+ k! p      label1=setdiff(label1,record);
      list=union(list,label1);
      pred(label1(find(pred(label1)==0)))=flag;
& U3 L" t" E9 c      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;
$ @* i1 h# J$ o9 Y8 s      maxf(label2)=min(maxf(flag),f(label2,flag));
      record=union(record,label2);
   end
8 z: O% v+ Q! k; k; {      if maxf(n)>0
: H, B7 j' b4 h# |# X2 p. T( F         v2=n;
9 D: m3 G" M5 z8 e' x; I6 e0 H/ c         v1=pred(v2);
         while v2~=1
           if v1>0
5 o. y- d3 [- }& r              f(v1,v2)=f(v1,v2)+maxf(n);
$ @& R  x; h* R7 ]$ _! |: I           else
6 ?4 W! D8 q( c* U! N' p" y           v1=abs(v1);
           f(v2,v1)=f(v2,v1)-maxf(n);
           end
* N9 @; J8 b% }8 k* A         v2=v1;
         v1=pred(v2);
        end
      end
   end
f
7 c% E5 [4 V: z* ?- E我想知道这个程序中：
* G  e3 H+ d: e
* Z% c" f1 U. X* R/ C  u, X1.什么是当前流量？这个应该不是求最大流的通用程序吧，因为有的题目里没有当前流量这个东西。
2.还有那个输出f是什么东东？？可行流矩阵式什么意思？？矩阵的含义完全不懂。比如说这个程序运行出来的那个矩阵我就不清楚他的意思。郁闷。

* G2 |3 ^; i/ Y" f9 S9 _谢谢啦。

木兆木风

f表示当前流量，u表示容量