function?[fwfNo]=MaxFlowMinCut_Me（nC）?
%?利用Ford--Fulkerson?標(biāo)號法求最大流算法的MATLAB?程序代碼?
%?f?%顯示最大流?
%?wf?%顯示最大流量?
%?No?%顯示標(biāo)號?由此可得最小割?
%?n?節(jié)點(diǎn)個(gè)數(shù)?
%?C?%弧容量?
%?Example:?
%?????n=8;?
????C=[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];??
%??????[fwfNo]=MaxFlowMinCut_Me（nC）?
??????
for（i=1:n）for（j=1:n）f（ij）=0;
????end;
end?%取初始可行流f?為零流?
for（i=1:n）No（i）=0;
????d（i）=0;
end?%Nod?記錄標(biāo)號?
while（1）?
No（1）=n+1;
d（1）=Inf;?%給發(fā)點(diǎn)vs?標(biāo)號?
while（1）pd=1;?%標(biāo)號過程?
for（i=1:n）if（No（i））?%選擇一個(gè)已標(biāo)號的點(diǎn)vi?
for（j=1:n）
????if（No（j）==0&f（ij）No（j）=i;
d（j）=C（ij）-f（ij）;
pd=0;?
if（d（j）>d（i））d（j）=d（i）;
end?
elseif（No（j）==0&f（ji）>0）?%對于未給標(biāo)號的點(diǎn)vj?當(dāng)vjvi?為非零流弧時(shí)?
No（j）=-i;
d（j）=f（ji）;
pd=0;?
if（d（j）>d（i））d（j）=d（i）;
end;
???end;
end;
????end;
end?
if（No（n）|pd）
????break;
end;
end?%若收點(diǎn)vt?得到標(biāo)號或者無法標(biāo)號?終止標(biāo)號過程?
if（pd）
????break;
end?%vt?未得到標(biāo)號?f?已是最大流?算法終止?
dvt=d（n）;
t=n;?%進(jìn)入調(diào)整過程?dvt?表示調(diào)整量?
while（1）?
if（No（t）>0）f（No（t）t）=f（No（t）t）+dvt;?%前向弧調(diào)整?
elseif（No（t）<0）f（No（t）t）=f（No（t）t）-dvt;end?%后向弧調(diào)整?
if（No（t）==1）
????for（i=1:n）No（i）=0;
????????d（i）=0;?
????end;
????break;
end?%當(dāng)t?的標(biāo)號為vs?時(shí)?終止調(diào)整過程?
t=No（t）;
end;
end;?%繼續(xù)調(diào)整前一段弧上的流f?
wf=0;
for（j=1:n）wf=wf+f（1j）;
end??
?
end

?屬性????????????大小?????日期????時(shí)間???名稱
-----------?---------??----------?-----??----
?????文件????????1573??2013-08-01?10:46??MaxFlowMinCut_Me.m