最高标号预流推进算法-清华大学.PPT

STEP0. (预处理) 置初始可行流x为零流; 计算精确的距离函数d；令当前节点i=s. STEP1. 若d（s） n, 继续下一步; 否则结束, 得到最优解x. STEP2. 如果存在节点i的某条出弧（i，j）为允许弧, 则转STEP3; 否则转STEP4. STEP3. 令pred(j)=i, 再令i=j; 若i=t,则找到了一条增广路, 增广并修改残量网络, 再令i=s. 转STEP1. STEP4. 修改标号(重新标号): 当{（i，j）?A(x) | uij(x) 0} ??时令d(i)=min{d(j)+1|(i，j)?A(x)且uij(x) 0}, 否则令d（i）=n; 且当i ? s时再令i=pred(i) . 转STEP1. “前进步” “回退步” 6.3.2 最短增广路算法 6.3.2 最短增广路算法 – 性质 引理6.4 在最短增广路算法中,距离标号始终保持是有效的. 归纳法： 在预处理阶段，构造的距离标号是有效的.只需证明一次增广操作和一次重新标号操作不改变距离标号的有效性. 增广操作可能引起残量网络的弧的变化只有两种情况： 允许路上的弧(i，j)可能退出残量网络:这不会影响(i，j)弧的端点上的距离标号的有效性； (j，i)弧可能进入残量网络:由于(i，j)弧是允许弧， 因此d(i)=d(j)+1，所以d(j)=d(i)-1 d(i), 这也不会影响(j，i)弧的端点上的距离标号的有效性. 重新标号将d(i)改为d’(i)=min{d(j)+1|(i,j)?A(x)且uij(x) 0}: 不会破坏任意的(i，j)弧的端点上的距离标号的有效性 ; 重新标号时节点i在残量网络中没有允许弧, 即对于任意的(i，j)?N(x)满足d(i) d(j)+1. 因此 d’(i) d(i), 即标号严格增加. 对于任意的(j，i)弧, d(j) ? d(i)+1 d’(i)+1，仍有效 6.3.2 最短增广路算法 – 性质 推论： 由于距离函数是有效的，算法结束时残量网络N（x）中从节点s到节点t没有允许路,即残量网络中不存在增广路,所以算法确实求到了最大流 例6.5 用最短增广路算法计算如下网络图6.5（a）中的最大流（s=1，t=4，弧上数字表示容量）. 1 2 3 4 1 2 3 4 5 1 2 3 4 1 2 3 4 5 1 1 2 0 1 2 3 4 1 2 3 4 1 1 2 3 4 1 2 4 1 4 0 1 1 2 1 1 2 3 4 1 1 3 4 1 4 1 0 1 2 1 2 3 4 1 2 1 4 5 1 2 3 4 1 2 1 4 5 2 0 1 2 3 1 2 3 4 1 2 1 4 5 1 2 3 4 1 2 1 4 5 2 0 1 2 4 最短增广路算法 – 首先思考：如何从一个节点出发找到从该节点出发的一条允许弧？ 数据结构：即对每个节点i，算法用链表A（i）记录该节点出发的所有弧，这些弧的顺序可以任意，但一旦确定以后在算法执行过程中就不再改变. 当前弧（Current-Arc）: 对每个节点i，用一个指针指向下一条将要被检索的弧，这条弧称为当前弧. 算法开始时，当前弧为A（i）中的第一条弧。 在查找从该节点出发的允许弧时，总是从当前弧开始判别：如果当前弧是允许弧，则返回这条弧；否则令A（i）中的下一条弧为当前弧。 如此下去，直到找到一条允许弧或判别完 A（i）的最后一条弧为止。 6.3.3 复杂度分析 1 2 4 2 2 3 1 2 3 4 5 6 0 5 4 10 0 初始当前弧:(1,2) 当前弧:(3,4) 分析：即确定“前进步”、“回退步”、增广、重标号的次数及每次的复杂度 最短增广路算法 – 采用当前弧（Current-Arc）结构时： 如果 A（i）中的一条弧在以前的迭代中不是允许弧，则在节点i的标号改变之前它仍然不是允许弧. Why? 非允许弧(i,j)的两种情况：(1) uij(x)=0; (2) d(i)d(j)+1. 这样我们得到如下结论： 复杂度分析 引理6.5 在最短增广路算法的执行过程中, 如果每个节点的重新标号操作不超过k次，则算法查找允许弧和重新标号操作的总复杂度为O（ |A（i）|）=O（km）. 如果判别完 A（i）的最后一条弧仍然没有找到允许弧，则说明 A（i）中不存在允许弧；这时，算法必须执行对该节点的重新标号操作，并重置该节点的当前弧为A（i）中的第一条弧. 由于对节点i的一次检查和一次重新标号操作分别需要对A（i）检索一遍，因此复杂度为O（| A（i）|）. 6.3.3 最短增广路算法 –复杂度分析