网络流算法专题优秀PPT.ppt
《网络流算法专题优秀PPT.ppt》由会员分享,可在线阅读,更多相关《网络流算法专题优秀PPT.ppt(52页珍藏版)》请在淘文阁 - 分享文档赚钱的网站上搜索。
1、图论算法图论算法 -最大流问题最大流问题长沙市雅礼中学朱 全 民 运输网络现在想将一些物资从S运抵T,必需经过一些中转站。连接中转站的是马路,每条马路都有最大运载量。每条弧代表一条马路,弧上的数表示该马路的最大运载量。最多能将多少货物从S运抵T?4248473 621STV1V2V3V4公路运输图公路运输图基本概念这是一个典型的网络流模型。为了解答此题,我们先了解网络流的有关定义和概念。若有向图G=(V,E)满足下列条件:有且仅有一个顶点S,它的入度为零,即d-(S)=0,这个顶点S便称为源点,或称为发点。有且仅有一个顶点T,它的出度为零,即d+(T)=0,这个顶点T便称为汇点,或称为收点。每
2、一条弧都有非负数,叫做该边的容量。边(vi,vj)的容量用cij表示。则称之为网络流图,记为G=(V,E,C)可行流可行流对于网络流图G,每一条弧(i,j)都给定一个非负数fij,这一组数满足下列三条件时称为这网络的可行流,用f表示它。1.每一条弧(i,j)有fijCij2.流量平衡除源点S和汇点T以外的全部的点vi,恒有:j(fij)=k(fjk)该等式说明中间点vi的流量守恒,输入与输出量相等。3.对于源点S和汇点T有,i(fSi)=j(fjT)=V(f)可增广路 给定一个可行流f=fij。若fij=Cij,称为饱和弧;否则称为非饱和弧。若fij=0,称为零流弧;否则称为非零流弧。定义一条
3、道路P,起点是S、终点是T。把P上全部与P方向一样的弧定义为正向弧,正向弧的全体记为P+;把P上全部与P方向相悖的弧定义为反向弧,反向弧的全体记为P-。譬如在图中,P=(S,V1,V2,V3,V4,T),那么P+=,P-=给定一个可行流f,P是从S到T的一条道路,假如满足:fij是非饱和流,并且 P+,fij是非零流,并且 P-那么就称P是f的一条可增广路。之所以称作“可增广”,是因为可改进路上弧的流量通过确定的规则修改,可以令整个流量放大。剩余图(残余网络)剩余图G=(V,E)流量网络G=(V,E)中,对于随意一条边(a,b),若flow(a,b)0则(a,b)E可以沿着a-b方向增广l剩余
4、图中,从源点到汇点的每一条路径都对应一条增广路Capacity=5Capacity=6Capacity=2Flow=2Flow=2Flow=2有向图32224剩余图l剩余图中,每条边都可以沿其方向增广剩余图的权值代表能沿边增广的大小G=(V,E,C)是已知的网络流图,设U是V的一个子集,W=VU,满足S U,TW。即U、W把V分成两个不相交的集合,且源点和汇点分属不同的集合。对于弧尾在U,弧头在W的弧所构成的集合称之为割切,用(U,W)表示。把割切(U,W)中全部弧的容量之和叫做此割切的容量,记为C(U,W),即:割切 割切示例v上例中,令U=S,V1,则W=V2,V3,V4,T,那么,vC(
5、U,W)=+=8+4+4+1=17流量算法的基本理论定理定理1:对于已知的网络流图,设随意一可行流为:对于已知的网络流图,设随意一可行流为f,随意一割切为随意一割切为(U,W),必有:,必有:V(f)C(U,W)。定理定理2:可行流:可行流f是最大流的充分必要条件是:是最大流的充分必要条件是:f中不中不存在可改进路。存在可改进路。定理定理3:整流定理。:整流定理。假如网络中全部的弧的容量是整数,则存在整数假如网络中全部的弧的容量是整数,则存在整数值的最大流。值的最大流。定理定理4:最大流最小割定理。:最大流最小割定理。最大流等于最小割,即最大流等于最小割,即max V(f)=min C(U,W
6、)。最大流算法第1步,令x=(xij)是随意整数可行流,可能是零流,给s一个永久标号(-,)。第2步(找增广路),假如全部标号都已经被检查,转到第4步。找到一个标号但未检查的点i,并做如下检查,对每一个弧(i,j),假如xij0,且j未标号,则给j一个标号(-i,(j),其中,(j)=minxji,(i)第3步(增广),由点t起先,运用指示标号构造一个增广路,指示标号的正负则表示通过增加还是削减弧流量来增加还是削减弧流量来增大流量,抹去s点以外的全部标号,转其次步接着找增广轨。第4步(构造最小割),这时现行流是最大的,若把全部标号的集合记为S,全部未标号点的集合记为T,便得到最小割(S,T)。
7、实例困难度分析设图中弧数为m,每找一条增广轨最多须要进行2m次弧的检查。假如全部弧的容量为整数,则最多须要v(其中v为最大流)次增广,因此总的计算量为O(mv)。procedure maxflow;最大流最大流var i,j,delta,x:integer;last:tline;可改进路中的前趋可改进路中的前趋 check:array0.maxn of boolean;检查数组检查数组begin repeat fillchar(last,sizeof(last),0);fillchar(check,sizeof(check),false);last1:=maxint;repeat i:=0;r
8、epeat inc(i)until(i n)or(lasti 0)and not checki;找到一个已检查而未标号的点找到一个已检查而未标号的点 if i n then break;for j:=1 to n do if lastj=0 then if flowi,j 0 then lastj:=-i;反向弧反向弧 checki:=true;until lastn 0;if lastn=0 then break;delta:=maxint;i:=n;repeat j:=i;i:=abs(lastj);if lastj 0 then x:=limiti,j-flowi,j else x:=f
9、lowj,i;if x 0 then inc(flowi,j,delta)else dec(flowj,i,delta);until i=1;放大网络流放大网络流 until false;end;利用找增广路的其他流量算法增广路的思想在于每次从源点搜寻出一条前往汇点的增广路,并变更路上的边权,直到无法再进行增广:一般增广路方法:在剩余图中,每次随意找一条增广路径增广。O(nmU)容量缩放增广路方法:在剩余图中,每次随意找一条最大可增广容量和的增广路径增广。O(nm*logU)最短增广路方法(MPLA):在剩余图中,每次随意找一条含结点数最少的增广路径增广。O(nm2)连续最短增广路方法(DIN
10、IC):在剩余图中,每次BFS找增广路径增广路径时,记录每个点的距离标号。在距离标号最短路图上,不断dfs找增广路,即一次标号,多次增广。O(n2m)DINIC算法演示:算法演示:源点汇点422532汇点32对增广路进行增广,增广后退回到源点1汇点232汇点1找到增广路路线,(红色路线)找到增广路路线,(红色路线)对增广路进行增广,增广后退回到源点,再无增广路线3用预流推动方法求网络流预流推动算法给每一个顶点一个标号h(v),表示该点到t的最短路(在残量网络中)。第一步hights()过程,就是BFS出初始最短路,计算出每一个顶点的h(v)。预流推动算法的特征是运用了预流来加快运算。预流说明图
11、中的节点(除s,t),仅须要满足流入量=流出量。其中流入量流出量的结点,我们称之为活动节点。我们的算法就是不断地将活动结点,变为非活动结点,使得预流成为可行流。预流推动算法流程算法过程prepare(),即首先将与s相连的边设为满流,并将这时产生的活动结点加入队列Q。这是算法的起先。以后便重复以下过程直到Q为空:(1).选出Q的一个活动顶点u。并依次推断残量网络G中每条边(u,v),若h(u)=h(v)+1,则顺着这里条边推流,直到Q变成非活动结点(不存在多余流量)。(Push推流过程)(2).假如u还是活动结点。则须要对u进行重新标号:h(u)=minh(v)+1,其中边(u,v)存在于G
12、中。然后再将u加入队列。(relable过程)可以证明,通过以上算法得到的结果就是最大流。预流推动算法示例顶点u的通过量g(u):剩余图中,找入边权和与出边权和的较小值l增广时,每次找一个通过量最小的点v,从点v向源点“推”大小为g(v)的流量向汇点“拉”大小为g(v)的流量尽量使剩余图中的边饱和34578g(u)=12用预流推动方法的一些网络流算法预流推动的算法核心思想是以边为单元进行推流操作:一般的预流推动算法:在剩余图中,维护一个预流,不断对活跃点执行push操作,或者relable操作来重新调整这个预流,直到不能操作。O(nm2)先进先出预流推动算法:在剩余图中,以先进先出队列维护活跃
13、点。O(n3)最高标号预流推动算法:在剩余图中,每次检查最高标号的活跃点,须要用到优先队列。O(n2m1/2)费用流流最重要的应用是尽可能多的分流物资,这也就是我们已经探讨过的最大流问题。然而实际生活中,最大配置方案确定不止一种,一旦有了选择的余地,费用的因素就自然参与到决策中来。右图是一个最简洁的例子:弧上标的两个数字第一个是容量,其次个是费用。这里的费用是单位流量的花费,譬如fs1=4,所需花费为3*4=12。简洁看出,此图的最大流(流量是8)为:fs1=f1t=5,fs2=f2t=3。所以它的费用是:3*5+4*5+7*3+2*3=62。(6,3)(5,4)(3,7)(8,2)STV1V
14、2费用流问题费用流定义设有带费用的网络流图G=(V,E,C,W),每条弧对应两个非负整数Cij、Wij,表示该弧的容量和费用。若流f满足:流量V(f)最大。满足a的前提下,流的费用Cost(f)=E(fij*Wij)最小。就称f是网络流图G的最小费用最大流。最小费用可改进路设P是流f的可改进路,定义P+Wij-P-Wij 为P的费用(为什么如此定义?)假如P是关于f的可改进路中费用最小的,就称P是f的最小费用可改进路。费用流算法求最小费用最大流的基本思想是贪心法。即:对于流f,每次选择最小费用可改进路进行改进,直到不存在可改进路为止。这样的得到的最大流必定是费用最小的。算法可描述为:第1步.令
15、f为零流。第2步.若无最小费用可改进路,转第5步;否则找到最小费用可改进路,设为P。第3步.依据P求delta(改进量)。第4步.放大f。转第2步。第5步.算法结束。此时的f即最小费用最大流。如何求最小费用可改进路 设带费用的网络流图G=(V,E,C,W),它的一个可行流是f。我们构造带权有向图B=(V,E),其中:V=V。若E,fijCij,那么E,权为Wij。若E,fij0,那么E,权为-Wij。明显,B中从S到T的每一条道路都对应关于f的一条可改进路;反之,关于f的每条可改进路也能对应B中从S到T的一条路径。即两者存在一一映射的逻辑关系。故若B中不存在从S到T的路径,则f必定没有可改进路
16、;不然,B中从S到T的最短路径即为f的最小费用可改进路。现在的问题变成:给定带权有向图B=(V,E),求从S到T的一条最短路径。迭代法求最短路经考虑到图中存在权值为负数的弧,不能接受Dijkstra算法;Floyd算法的效率又不尽如人意所以,这里接受一种折衷的算法:迭代法(bellman算法)。设Shorti表示从S到i顶点的最短路径长度;从S到顶点i的最短路径中,顶点i的前趋记为Lasti。那么迭代算法描述如下:(为了便于描述,令n=|V|,S的编号为0,T的编号为n+1)step 1.令Shorti +(1in+1),Short0 0。step 2.遍历每一条弧。若Shorti+Short
17、j,则令Shortj Shorti+,同时Lastj i。重复做step 2直到不存在任何任何弧满足此条件为止。step 3.算法结束。若Shortn+1=+,则不存在从S到T的路径;否则可以依据Last记录的有关信息得到最短路径。一次迭代算法的时间困难度为O(kn2),其中k是一个不大于n的变量。在费用流的求解过程中,k大部分状况下都远小于n。procedure costflow;求最小费用最大流求最小费用最大流var i,j,x,delta:integer;best,last:tline;best:最短路长度;:最短路长度;last:可改进路中的前趋顶点:可改进路中的前趋顶点 more:b
18、oolean;begin repeat fillword(best,sizeof(best)shr 1,maxint);fillchar(last,sizeof(last),0);last1:=maxint;best1:=0;赋初值赋初值 repeat more:=false;for i:=1 to n do if besti maxint then for j:=1 to n do begin if(flowi,j limiti,j)and (besti+costi,j 0)and (besti-costj,i 0 then x:=limiti,j-flowi,j else x:=flowj
19、,i;if x 0 then inc(flowi,j,delta)else dec(flowj,i,delta);until i=1;until false;根据改进量放大流根据改进量放大流end;思维发散与探究 可改进路费用:“递增!递增?”设f从零流到最大流共被改进了k次,每i次选择的可改进路的费用为pi,那么会不会有p1p2p3pk呢?迭代法:“当心死循环!嘿嘿”迭代法会出现死循环吗?也就是说,构造的带权有向图B中会存在负回路吗?费用:“你在乎我是负数吗?”容量:“我管的可不仅是弧。”网络流图中的“容量”都是对弧而言的,但若是给每个顶点也加上一个容量限制:即通过此顶点的流量的上限;任务仍
20、旧是求从S到T的最小费用最大流。你能解决吗?有上下界的最大流 上面探讨的网络流都只对每条弧都限定了上界(其实其下界可以看成0),现在给每条弧加上一个下界限制Aij(即必需满足Aijfij)。弧上数字对第一个是上界,其次个是下界。若是撇开下界不看,此图的最大流如图所示,流量是6;但若是加入了下界的限制,它的最大流量就只有5了。(3,0)(3,0)(3,0)(3,0)(10,1)原问题原问题33330(a)32231(b)增广增广怎样找可行流一种自然的想法是去掉下界,将其转化为只含上界的网络流图。这种奇妙的愿望是可以实现的。具体方法如下:设原网络流图为G=(V,E,C,A),构造不含下界的网络流图
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 网络 算法 专题 优秀 PPT
限制150内