Successive shortest paths

Abstract view

Definition:

1. For a node [math]\displaystyle{ v\in V }[/math], let [math]\displaystyle{ \Delta f(v):=\sum_{w:(v,w)\in A}f(v,w)-\sum_{w:(w,v)\in A}f(w,v) }[/math].
2. The imbalance of a node [math]\displaystyle{ v\in V }[/math] is defined as [math]\displaystyle{ I_f(v):=\Delta f(v)-b(v) }[/math].
3. The imbalance of a node [math]\displaystyle{ v\in V }[/math] is underestimating if [math]\displaystyle{ 0\leq \Delta f(v)\leq b(v) }[/math] or [math]\displaystyle{ 0\geq\Delta f(v)\geq b(v) }[/math].
4. The total imbalance of [math]\displaystyle{ f }[/math] is the defined as [math]\displaystyle{ \sum_{v\in V}|I_f(v)| }[/math].

Invariant:

1. The capacity constraints are fulfilled, that is, [math]\displaystyle{ 0\leq f(a)\leq u(a) }[/math] for all [math]\displaystyle{ a\in A }[/math].
2. There is no negative cycle in the residual network of [math]\displaystyle{ f }[/math].
3. The imbalance of every node is underestimating.

Variant: The total imbalance strictly decreases.

Break condition: The imbalances of all nodes are zero.

Induction basis

Abstract view: Start with the zero flow.

Proof: Obvious.

Induction step

1. In the residual network of [math]\displaystyle{ f }[/math], find a shortest path [math]\displaystyle{ p }[/math] from the set of nodes [math]\displaystyle{ v\in V }[/math] with [math]\displaystyle{ I_f(v)\lt 0 }[/math] to the set of nodes [math]\displaystyle{ w\in V }[/math] with [math]\displaystyle{ I_f(v)\gt 0 }[/math].
2. Let [math]\displaystyle{ v_0 }[/math] be the node where [math]\displaystyle{ p }[/math] actually starts and [math]\displaystyle{ w_0 }[/math] the node where [math]\displaystyle{ p }[/math] actually ends.
3. Increase the flow on all arcs of [math]\displaystyle{ p }[/math] by the minimum of [math]\displaystyle{ |I_f(v_0)| }[/math], of [math]\displaystyle{ I_f(w_0)\gt 0 }[/math], and of the residual capacities of all arcs on [math]\displaystyle{ p }[/math].

Correctness

The variant follows from the fact that the amount by which the flow is increased in step 3 is strictly positive. The first and third points of the invariant result from the choice of this amount as the minimum of imbalances and residual capacities. Thus, it remains to show point 2 of the invariant.

By induction hypothesis, there were nonegative cycles immediately before the iteration. So, any negative cycle [math]\displaystyle{ C }[/math] must have been introduced by the flow augmentation along [math]\displaystyle{ p }[/math]. Such a cycle must contain at least one arc [math]\displaystyle{ (v,w) }[/math] on [math]\displaystyle{ p }[/math] such that the residual capacity [math]\displaystyle{ w\rightarrow v }[/math] was zero immediately before the current iteration (and is positive afterwards, of course). Suppose for a contradiction that a negative cycle [math]\displaystyle{ C }[/math] has been introduced. For some [math]\displaystyle{ k\gt 0 }[/math], let [math]\displaystyle{ (v_1,w_1),\ldots,(v_k,w_k) }[/math] denote, in this (cyclic) order, the arcs on [math]\displaystyle{ p }[/math] such that, for each [math]\displaystyle{ k\in\{1,\ldots,,k\} }[/math], [math]\displaystyle{ (w_i,v_i) }[/math] belongs to [math]\displaystyle{ C }[/math] and the residual capacity of [math]\displaystyle{ (w_i,v_i) }[/math] was zero immediately before the current iteration.

Let [math]\displaystyle{ p_1,\ldots,p_k }[/math] denote the maximal subpaths of [math]\displaystyle{ p }[/math] whose internal nodes do not belong to [math]\displaystyle{ p }[/math]. In other words, each [math]\displaystyle{ p_i }[/math] starts and ends on [math]\displaystyle{ p }[/math] but does not share any further node with [math]\displaystyle{ p }[/math]. All arcs no all of these subpaths are in the residual network immediately before and immediately after the current iteration. For [math]\displaystyle{ i\in\{1,\ldots,k\} }[/math], let [math]\displaystyle{ p' }[/math] denote

1. the corresponding subpath of [math]\displaystyle{ p }[/math] if the start node of [math]\displaystyle{ p_i }[/math] appears on [math]\displaystyle{ p }[/math] before its end node;
2. the corresponding subpath of the transpose of [math]\displaystyle{ p\lt /math\lt if the end node of \lt math\gt p_i }[/math] appears on [math]\displaystyle{ p }[/math] before its start node.

Let [math]\displaystyle{ C' }[/math] be the cycle that results from [math]\displaystyle{ C }[/math] when each subpath [math]\displaystyle{ p_i }[/math] of the first type is replaced by the corresponding subpath [math]\displaystyle{ p'_i }[/math] of [math]\displaystyle{ p }[/math]. Obviously, [math]\displaystyle{ C' }[/math] is augmenting immediately before the current iteration.

Since [math]\displaystyle{ p }[/math] is shortest, each [math]\displaystyle{ p_i }[/math] of the first type is at least as long as the corresponding [math]\displaystyle{ p'_i }[/math] (recall that each [math]\displaystyle{ p_i }[/math] was augmenting when [math]\displaystyle{ p }[/math] was found).