[转帖] 网络流概念及相关算法介绍

matrixer

<TABLE height="100%" cellSpacing=0 cellPadding=0 width="100%" border=0>
4 s, c( L6 x/ \; P$ S6 d& ~  c: J0 @
<TR>
<TD  width=74><IMG src="http://www.frontfree.net/articles/pages/0000000554/title.gif" border=1></TD>
<TD vAlign=top width="100%">
<TABLE height="100%" cellSpacing=0 cellPadding=0 width="100%" border=0>
' s2 y) g3 n3 i( m3 J
<TR>
/ U3 ~' _+ h& W! z  g9 J* {<TD class=artitle vAlign=top colSpan=2>网络流概念及相关算法介绍</TD></TR>
<TR vAlign=top>
2 t& s1 n, c9 B. F<TD align=left>原创：怒火支袍 </TD>
) l4 ?7 I, \0 y/ a<TD class=text vAlign=top align=right>2003年6月17日 </TD></TR></TABLE></TD></TR>
; g( [) I1 Y( D+ f<TR>
<TD class=arcontent colSpan=3>
# z. A+ @: b' W2 y& K6 A
<STYLE type=text/css>
1 k4 P; S) v0 o  H) a4 V( B<!--
/ D. s1 X0 p0 ]" _' X.titletxt {
3 b5 e% I* S8 s        font-size: 18px;
}
% ~( d! R/ h4 `) D" H& j# W.tabletxt {
* r3 z# j8 [* c* C1 b: J* D        font-size: 14px;
        padding: 7px;
}
8 s5 f- T3 a* t! W-->
# `# N7 g, y! l, z</STYLE>
+ e  N' r( b: n0 M
<><b>一、引言</b></P>
) g, ]: B  z6 V: U/ F< align=left>如同我们可以把一个实际的道路地图抽象成一个有向图来计算两点之间的最短路径，我们也可以将一个有向图看作一个流网络来解决另一类型的问题。流网络比较适合用来模拟液体流经管道、电流在电路网络中的运动、信息网络中信息的传递等等类似的过程。</P>
( U0 v5 D, T2 L6 {% r: [4 X< align=left><b>二、网络流和最大流问题
' i  `5 C9 \3 M</b>
参看下图，给定一个有向图G=(V,E)，把图中的边看作管道，每条边上有一个权值，表示该管道的流量上限。给定源点s和汇点t，现在假设在s处有一个水源，t处有一个蓄水池，问从s到t的最大水流量是多少，类似于这类的问题都可归结为网络流问题。</P>
6 O$ k- o5 ?, [7 a1 F( J5 m: S<><IMG src="http://www.frontfree.net/articles/pages/0000000554/pic01.gif"></P>
< align=left>在流网络中，每条有向边可以被看导管。每根导管有一个固定的容量，代表物质流经这个导管的最大速率，例如一个管道每小时最多能流过200加仑液体或者一根电线最多能承载20安培的电流。流网络中的顶点可以看作是导管的连接处。除了源点和汇点之外，物质流进每个点的速率必须等于流出这个点的速率。如果我们把研究的物质特化为电流，这种“流的保持”属性就好像电路中的基尔霍夫电流定律一样。</P>
< align=left>下面我们用数学语言来进行相关概念的定义：
; e  x5 L9 x) ~2 ?. a( s* l0 @
! b8 Y9 K7 |: D5 f. ~0 [设G=(V,E)是一个流网络，设c(u, v)&gt;＝0 表示从u到v的管道的流量上限。设s为源，t为汇。G的流是一个函数f: V×V →R，且满足下面三个特征：</P>
& v) f0 s" m* Z3 b. {<TABLE cellSpacing=1 cellPadding=0 width="100%" bgColor=#000000 border=0>
* n; t# m! ]; }5 ]8 x  M
' u6 b+ Z; H) _/ j2 J<TR>
<TD class=tabletxt width="2%" bgColor=#e0e0e0>1. </TD>
<TD class=tabletxt width="98%" bgColor=#ffffff>容量限制：对于所有的 u,v ∈ V, 要求f(u, v) &lt;= c(u, v) </TD></TR>
& f: [1 p4 n( N6 |5 y6 z<TR>
<TD class=tabletxt bgColor=#e0e0e0>2. </TD>
<TD class=tabletxt bgColor=#ffffff>斜对称性：对于所有的 u,v ∈ V, 要求f(u, v) = - f(v, u)</TD></TR>
) x* y2 o/ `) E<TR>
0 i8 F5 \0 B) _! b: |" ^<TD class=tabletxt bgColor=#e0e0e0>3.</TD>
5 i1 O4 S! q9 I# |; b9 F7 ]<TD bgColor=#ffffff>
<>流的保持：对于所有的 u ∈ V - {s, t}，要求：∑ f(u, v) = 0(v∈V)
f(u,v)称为从结点u到v的网络流，它可以为正也可以为负。流 f 的值定义为：|f| = ∑ f(s, v)(v∈V)即从源出发的所有流的总和。</P></TD></TR></TABLE>
< align=left>最大流问题就是找出给定流网络的最大流。网络流问题可以归结为一类特殊的线性规划问题。</P>
< align=left><b>三、解决最大流问题常用算法一览</b></P>
9 k- }1 w) P" u7 D% \! g& B< align=left>解决最大流问题的常用到Ford-Fulkerson方法，之所以称其方法而不是算法，是因为在这种思想下包含着若干种时间复杂度不同的实现，其中较多地是使用Edmonds-Karp算法。与此相对，Push-relabel算法采用了与Ford-Fulkerson方法完全不同的思考角度，降低了渐进意义下的时间复杂度。而relabel-to-front算法则是对Push-relabel算法的改良和精炼，效率更佳。</P>
1 Q; E2 @9 M9 U$ ?5 U+ K< align=left>关于这三种常用算法的时间复杂度可见下表：（其中V表示图的顶点数，E表示边数）</P>
<TABLE cellSpacing=1 width="100%" bgColor=#000000 border=0>

<TR bgColor=#e0e0e0>
<TD class=tabletxt width="25%">
<DIV align=center>算法名称</DIV></TD>
/ s4 A8 y- ^8 w<TD class=tabletxt width="23%" bgColor=#ffffff>
) E2 w% b$ y- }6 a$ `0 Z<DIV align=center>Edmonds-Karp算法</DIV></TD>
2 x' u  A2 k8 e0 P2 F4 D<TD class=tabletxt width="26%" bgColor=#ffffff>
" f% ]2 L1 V- ]9 o7 {6 l$ I% P4 x- |<DIV align=center>一般性的push-relabel算法</DIV></TD>
<TD class=tabletxt width="26%" bgColor=#ffffff>
<DIV align=center>relabel-to-front算法</DIV></TD></TR>
' G$ G& h) c* K7 R6 |( `' f( K<TR bgColor=#e0e0e0>
<TD class=tabletxt width="25%">
<DIV align=center>时间复杂度</DIV></TD>
<TD class=tabletxt width="23%" bgColor=#ffffff>
<DIV align=center>O(V*E^2)</DIV></TD>
<TD class=tabletxt width="26%" bgColor=#ffffff>
; ~5 S  e7 e, y/ K; g! l: v6 X<DIV align=center>O(V^2*E)</DIV></TD>
* L! g% B( @( c  G/ _( |3 G<TD class=tabletxt width="26%" bgColor=#ffffff>
. W  q5 A7 M3 \0 e; }& a$ Z<DIV align=center>O(V^3)</DIV></TD></TR></TABLE>
< align=left>可以看出，当给定的有向图比较稀疏时，三种算法的效率不会相差太多，但当网络稠密时，relabel-to-front算法在效率上有着明显的优势。
<b>
4 ]5 Z' ]1 y6 a$ ]8 Y</b><b>四、基于Ford-Fulkerson方法的Edmonds-Karp实现</b></P>
% J, ~' G- }) m; H0 Z< align=left>一般的Ford-Fulkerson方法具有迭代性质，我们把顶点u和v之间的流记作f(u,v)。那么在最开始，我们对所有的u,v∈V置f(u,v)=0。在每次的迭代过程中，通过找到一条增加路径来使|f|增加。在这里，我们可以简单地认为所谓的“增加路径”就是一条可以传送比当前更多流的从源点s到汇点t的路径，一旦找到了这样的路径，我们就可以得到一个比原流数值更大的新流。重复这个过程，直到不存在增加路径为止，这就是Ford-Fulkerson方法的主要过程，可以用伪码表示如下：</P>
* C" [9 R! i$ L& ], K- t<TABLE cellSpacing=4 cellPadding=0 width="100%" border=0>
* A8 X: H5 `. d9 S. P6 |/ _4 A
* C4 n! W1 L- |/ J/ R3 ?2 T' e+ s+ @<TR>
: _" r/ B0 }1 d( K' v( O<TD bgColor=#e0e0e0>
< align=left>FORD-FULKERSON-METHOD(G,s,t)

将流f初始化为0
" H1 p8 C% x9 G0 `, I
while 存在一条增加路径p
; d, {* b" `( H& D% Y
- N5 n6 X) b( b! vdo 顺沿p增加f
1 u  t2 a# k" I. M
9 y8 g+ Y7 s$ s' B9 xreturn f</P></TD></TR></TABLE>
< align=left>实现Ford-Fulkerson的时间复杂度主要取决于如何寻找增加路径p。Edmonds-Karp实现正是通过采用了广度优先的搜索策略得以使其复杂度达到O(V*E^2)。</P>
6 ^7 T* A8 \1 `7 c, G# M4 i< align=left>由于这种算法的效率不很理想，我们在此不多着墨，而主要介绍下述push-relabel算法的思想。</P>
9 x3 h. U* M% R% z/ C& Q4 @* L# m< align=left><b>五、一般性的push-relabel算法</b></P>
5 ~# P# R  S3 w/ S9 O5 t+ V9 G9 U5 }< align=left>很多渐进意义下最优的算法都是采用了push-relabel算法的思想，而且很多其他的相关问题，比如最小费用流问题，也可以用这种方法很好的解决。首先介绍的是一般性的push-relabel算法。</P>
7 x* f# |; r+ S2 o" N< align=left>不同于Ford-Fulkerson方法在残留网络中寻找增加路径的方式，push-relabel算法在运行的过程中只关注某一个顶点以及它的相邻顶点，在这个过程中，它并不像Ford-Fulkerson方法保持着“流的保持”性质，而是以一个“先流”进行运作。这个先流同样是一个 V×V →R的函数，满足容量限制和斜对称性，同时，它对所有的u∈V-{s}满足f(V,u)&gt;=0。我们记e(u)=f(V,u)。如果e(u)&gt;0我们就说顶点u溢出。</P>
< align=left>为了步入正题，我们还需要介绍push-relabel算法引入的一个额外的高度函数。设G=(V,E)是一个流网络，源点是s，汇点是t，f是G中的一个先流。如果函数h:V→N满足h(s)=|V|，h(t)=0，而且对残留网络中所有的边(u,v)有h(u)&lt;=h(v)+1，那么称h是一个高度函数。</P>
- R+ ?% h- T$ w( w  {< align=left>正如其名称一样，push-relabel算法有两个基本操作:push和relabel。一般性的push-relabel算法就是通过往复执行这两种操作完成的：</P>
/ i* u: ^/ M) W5 m' t7 v. ]# @<TABLE cellSpacing=4 cellPadding=0 width="100%" border=0>
: q/ D- B: k1 j! a
1 m: [$ p9 o4 ?- C<TR>
<TD bgColor=#e0e0e0>
4 g7 U% H8 }$ K& y- Y+ |( ~< align=left>GENERIC-PUSH-RELABEL(G)

先流初始化

while 存在可以执行的push或relabel操作

  t& E. w8 m, A0 h　　选择一个可以执行的push或relabel操作执行</P></TD></TR></TABLE>
, i/ h* F9 l9 |: @< align=left>下面具体介绍一下这两个基本操作。</P>
% e$ M+ t& F; |* Q+ _2 C( S2 ~<TABLE cellSpacing=1 cellPadding=0 width="100%" bgColor=#000000 border=0>

<TR bgColor=#ffffff>
<TD class=tabletxt width="5%" rowSpan=3><b>USH(u,v)</b></TD>
<TD class=tabletxt width="95%">可以执行的时机：顶点u溢出，u、v之间的残留容量cf(u,v)为正，且h=h[v]+1</TD></TR>
6 z( s) i6 T5 l( ?  R# [0 l<TR>
% A% [. G5 M& I" _) C0 f<TD class=tabletxt bgColor=#ffffff>动作描述：将df(u,v)=min(e,cf(u,v))个单位的流从u压向v</TD></TR>
<TR>
<TD bgColor=#e0e0e0>
3 N- c- ~+ h# c/ n< align=left>具体步骤：
) ~& j4 {% c  K* ?( `5 x: G0 c: ~% u
* j) e; d8 j6 \# A2 |df(u,v)=min(e,cf(u,v))

6 K0 a5 o) m" M* @/ s4 p( q6 Ef[u,v]=f[u,v]+df(u,v)
1 R  N7 v4 E" O
f[v,u]=-f[u,v]

) i$ M$ b5 h& F) s0 d. N# }; |, B2 ne=e-df[u,v]
- j0 X5 n# A% R3 r9 r; o8 U& \& W
e[v]=e[v]+df[u,v]</P></TD></TR>
<TR bgColor=#ffffff>
1 a( j* l( O  `6 u) I# m9 k<TD rowSpan=3>
# q. Y' F# l4 J& d< align=left><b>RELABEL(u)</b></P></TD>
<TD class=tabletxt bgColor=#ffffff>可以执行的时机：u溢出，且对所有的残留网络中的边(u,v),有h&lt;=h[v]</TD></TR>
<TR>
* t- z6 Q6 B. Q" _" d7 }<TD bgColor=#ffffff>
: g6 o/ h& V: |3 x: x< align=left>动作描述：增加u的高度</P></TD></TR>
: Q3 P# a' ~% v  E<TR>
% B$ u7 F$ Z% V8 S- ~8 v) C<TD bgColor=#e0e0e0>
< align=left>具体步骤：

h=1+min{h[v]u,v)是残留网络中的边}</P></TD></TR></TABLE>
< align=left>通过证明，在一般性的push-relabel算法执行过程中，relabel操作的执行次数小于2|V|^2，push操作的执行次数小于2|V||E|+4|V|^3+4|E||V|^2，而每个relabel操作的耗时在O(V)级，每个push的耗时在O(1)级，选择一个可以执行的操作也可以在O(1)内完成，因此，存在具体的实现使得一般性的push-relabel算法时间复杂度达到O(V^2*E)。</P>
5 {/ D$ c! [$ `1 s( `< align=left><b>六、relabel-to-front算法</b></P>
' W7 Q  k' k* r, @. }/ P) B5 |<P align=left>通过引入邻接表和许可边的概念，relabel-to-front算法在push-relabel算法的基础上进一步提升了效率，使时间复杂度可以达到O(V^3)，但是该算法的步骤和证明的过程比较繁琐，在这里就略去了，有兴趣的读者可以参考《算法导论》。</P>
  `# t: z7 L4 l: M: [; [% x( ?6 J( _0 l& {<P align=left><b>七、二部图的最大匹配与网络最大流的关系</b></P>
+ l9 f: m: W3 u6 X0 K3 J<P align=left>有一定离散基础的读者应当对二部图的最大匹配不陌生，但它和网络流之间有什么具体的联系呢，请看下图：</P>
<P><IMG src="http://www.frontfree.net/articles/pages/0000000554/pic02_02.gif"></P>
<P><IMG src="http://www.frontfree.net/articles/pages/0000000554/pic02_01.gif"></P>
5 M5 r$ @0 _: d( c' n1 n<P align=left>是的，如果我们设二部图两部分的点集分别为L与R,现在添加源点s和汇点t,对所有的v∈L添加有向边(s,v)，再对所有的v∈R添加有向边(v,t)，再将原二部图中所有的无向边改为自L中的点指向R中的点的有向边，就构造了一个流网络。如果这个网络中每条边的容量限制都设为1，那么它的最大流数值就等于二部图的匹配数，而且这个最大流和二部图中的最大匹配是一一对应的。</P>
  s5 U: X/ ]- S1 `<P align=left>用邻接表存储的二部图可以用匈牙利算法在O(VE) 的时间内找到最大匹配，这意味着用网络流解决二部图最大匹配问题并非最具效率的选择，但是它至少向我们展示了网络流的一个侧面应用。除此之外，网络流的变形和演化还可以解决很多具有普遍意义的问题，比如最小路径覆盖等等，所以笔者建议对图论感兴趣的同学不妨多研究一下相关内容，一定会有所收获。</P></TD></TR></TABLE>

fang307

顶

blueblood

<>谢谢</P>

blueblood

<>very hard</P>

文素

这么多的？？？？？？

紫辰

恩恩，能不能给我一个二部图的最大匹配问题MATLAB程序？