查看: 2531|回复: 0

[问题求助] 最大流通用函数问题。

字体大小: 正常放大

6 主题	7 听众	140 积分

升级 20%

TA的每日心情

	郁闷 2014-2-7 13:28

签到天数: 47 天

[LV.5]常住居民I

自我介绍: 好好学习，天天向上。

电梯直达

1^#

发表于 2013-1-20 17:40 |只看该作者 |倒序浏览

|招呼Ta 关注Ta

最大流通用程序
%function [f,s]=maxflow(startp,endp,c)
%c为容量网络
%对容量网络的填写做一下说明
%容量具有方向性，比如弧(i,j)的容量为10，弧(j,i)为0
%即矩阵无须有对称性
function [f,s]=maxflow(startp,endp,c)
n=length(c);
f=zeros(size(c));
l=zeros(1,n);d=zeros(1,n);examine=zeros(1,n);
l(startp)=0.5;d(startp)=inf;
while 1
ifexam=0;ifl=0;
for i=1:n
      if l(i)~=0
         ifl=ifl+1;
         if examine(i)==1
            ifexam=ifexam+1;
         end
      end
end
if ifl==ifexam
      break;
end
for i=1:n
      if l(i)~=0&examine(i)==0
         break;
      end
end
for j=1:n
      if c(i,j)~=0
         if f(i,j)<c(i,j)&l(j)==0
            l(j)=i;
            d(j)=min(d(i),c(i,j)-f(i,j));
         end
      end
      if c(j,i)~=0
         if f(j,i)>0&l(j)==0

            l(j)=-i;
            d(j)=min(d(i),f(i,j));
         end
      end
end
examine(i)=1;
if l(endp)~=0
      j=endp;
      while 1
         if l(j)~=0.5
            if l(j)>0
                  i=l(j);
                  f(i,j)=f(i,j)+d(endp);
                  j=i;
            end
            if l(j)<0
                  i=-l(j);
                  f(j,i)=f(j,i)-d(endp);
                  j=i;
            end
         else
            l=zeros(1,n);break;
         end
      end
      l(startp)=0.5;d(startp)=inf;examine=zeros(1,n);
end
end
s=[];ns=0;
for i=1:n
if l(i)~=0
      ns=ns+1;
      s(ns)=i;
end
end
fprintf('f为最大可行流\n');
fprintf('图的最小截划分得到的一个子集s为：\n');
disp(s);

我想知道这个程序中function [f,s]=maxflow(startp,endp,c)，f,s,startp,endp,c都代表什么意思？？
最好能举个例子。麻烦啦。

zan

function, examine