Dinic's algorithm explained

Dinic's algorithm or Dinitz's algorithm is a strongly polynomial algorithm for computing the maximum flow in a flow network, conceived in 1970 by Israeli (formerly Soviet) computer scientist Yefim Dinitz.^[1] The algorithm runs in

O(|V|^2|E|)

time and is similar to the Edmonds–Karp algorithm, which runs in

O(|V||E|²⁾

time, in that it uses shortest augmenting paths. The introduction of the concepts of the level graph and blocking flow enable Dinic's algorithm to achieve its performance.

History

Dinitz invented the algorithm in January 1969, as a master's student in Georgy Adelson-Velsky's group. A few decades later, he would recall:

In 1970, Dinitz published a description of the algorithm in Doklady Akademii Nauk SSSR. In 1974, Shimon Even and (his then Ph.D. student) Alon Itai at the Technion in Haifa were very curious and intrigued by Dinitz's algorithm as well as Alexander V. Karzanov's related idea of blocking flow. However it was hard for them to decipher these two papers, each being limited to four pages to meet the restrictions of journal Doklady Akademii Nauk SSSR. Even did not give up, and after three days of effort managed to understand both papers except for the layered network maintenance issue. Over the next couple of years, Even gave lectures on "Dinic's algorithm", mispronouncing the name of the author while popularizing it. Even and Itai also contributed to this algorithm by combining BFS and DFS, which is how the algorithm is now commonly presented.

For about 10 years of time after the Ford–Fulkerson algorithm was invented, it was unknown if it could be made to terminate in polynomial time in the general case of irrational edge capacities. This caused a lack of any known polynomial-time algorithm to solve the max flow problem in generic cases. Dinitz's algorithm and the Edmonds–Karp algorithm (published in 1972) both independently showed that in the Ford–Fulkerson algorithm, if each augmenting path is the shortest one, then the length of the augmenting paths is non-decreasing and the algorithm always terminates.

Definition

Let

G=((V,E),c,f,s,t)

be a network with

c(u,v)

and

f(u,v)

the capacity and the flow of the edge

(u,v)

, respectively.

The residual capacity is a mapping

c_f\colonV x V\toR⁺

defined as,

(u,v)\inE

c_f(u,v)=c(u,v)-f(u,v)

(v,u)\inE

c_f(u,v)=f(v,u)

c_f(u,v)=0

otherwise.

The residual graph is an unweighted graph

G_f=((V,E_f),c_f|


	E_f

,s,t)

, where

E_f=\{(u,v)\inV x V\colon c_f(u,v)>0\}

An augmenting path is an

–

path in the residual graph

G_f

Define

\operatorname{dist}(v)

to be the length of the shortest path from

G_f

. Then the level graph of

G_f

is the graph

G_L=((V,E_L),c_f|


	E_L

,s,t)

, where

E_L=\{(u,v)\inE_f\colon \operatorname{dist}(v)=\operatorname{dist}(u)+1\}

A blocking flow is an

–

flow

such that the graph

G'=((V,E_L'),s,t)

with

E_L'=\{(u,v)\colon f'(u,v)<c_f|


	E_L

(u,v)\}

contains no

–

path. ^[2]

Algorithm

Dinic's Algorithm

Input: A network

G=((V,E),c,s,t)

Output: An

–

flow

of maximum value.

f(e)=0

for each

e\inE

Construct

G_L

from

G_f

. If

\operatorname{dist}(t)=infty

, stop and output

Find a blocking flow

G_L

Add augment flow

and go back to step 2.

Analysis

It can be shown that the number of layers in each blocking flow increases by at least 1 each time and thus there are at most

|V|-1

blocking flows in the algorithm. For each of them:

the level graph

G_L

can be constructed by breadth-first search in

O(E)

time

a blocking flow in the level graph

G_L

can be found in

O(VE)

time^[3]

with total running time

O(E+VE)=O(VE)

for each layer. As a consequence, the running time of Dinic's algorithm is

O(V²E)

.^[4]

Using a data structure called dynamic trees, the running time of finding a blocking flow in each phase can be reduced to

O(ElogV)

and therefore the running time of Dinic's algorithm can be improved to

O(VElogV)

Special cases

In networks with unit capacities, a much stronger time bound holds. Each blocking flow can be found in

O(E)

time, and it can be shown that the number of phases does not exceed

O(\sqrt{E})

and

O(V^2/3)

.^[5] Thus the algorithm runs in

O(min\{V^2/3,E^1/2\}E)

time.^[6]

In networks that arise from the bipartite matching problem, the number of phases is bounded by

O(\sqrt{V})

, therefore leading to the

O(\sqrt{V}E)

time bound. The resulting algorithm is also known as Hopcroft–Karp algorithm. More generally, this bound holds for any unit network — a network in which each vertex, except for source and sink, either has a single entering edge of capacity one, or a single outgoing edge of capacity one, and all other capacities are arbitrary integers.

Example

The following is a simulation of Dinic's algorithm. In the level graph

G_L

, the vertices with labels in red are the values