Algorithme de Dinic
# Algorithme de Dinic

L'algorithme de Dinic ou algorithme de Dinitz est un algorithme en temps polynomial (et même fortement polynomial) de calcul du flot maximum dans un réseau, publié en 1970 par Yefim Dinitz. Le temps de calcul est en $O(|V|^{2}|E|)$ pour un graphe dont $V$ est l'ensemble des sommets et $E$ l’ensemble des arcs. Il est semblable à l'algorithme d'Edmonds-Karp dont le temps d'exécution est en $O(|V||E|^{2})$. Comme lui, il utilise des chemins augmentants de longueur minimale. L'introduction des concepts de graphe de niveau et de flot bloquant permet d'obtenir cette meilleure performance.

## Définitions

Soit $G=((V,E),c,s,t)$ un réseau; $c(u,v)$ et $f(u,v)$ dénotent respectivement la capacité et le flot sur un arc $(u,v)$.

La capacité résiduelle est l’application $c_{f}:V\times V\to R^{+}$ définie comme suit :

* si $(u,v)\in E$,

    $c_{f}(u,v)=c(u,v)-f(u,v)$,
    $c_{f}(v,u)=f(u,v)$,
* $c_{f}(u,v)=0$ sinon.

Le graphe résiduel est le graphe $G_{f}=((V,E_{f}),c_{f}|_{E_{f}},s,t)$, où

    $E_{f}=\{(u,v)\in V\times V:c_{f}(u,v)>0\}$.

Un chemin augmentant est un chemin de la source $s$ au puits $t$ (on dit aussi un $s$-$t$-chemin) dans le graphe résiduel. On note $d(v)$ la longueur des plus courts chemins de la source $s$ au sommet $v$ dans $G_{f}$ ($d(v)$ est la distance de la source à $v$).

Le graphe de niveau de $G_{f}$ est le graphe $G_{L}=(V,E_{L},c_{f}|_{E_{L}},s,t)$, où

    $E_{L}=\{(u,v)\in E_{f}:d(v)=d(u)+1\}$.

En d'autres termes, c'est le sous-graphe contenant uniquement les arcs qui relient un sommet à un sommet de distance immédiatement supérieure.

Un flot bloquant est un flot $f$ tel que le graphe $G'=(V,E_{L}',s,t)$ avec $E_{L}'=\{(u,v):f(u,v)<c_{f}|_{E_{L}}(u,v)\}$ obtenu en ne conservant que les arcs sur lesquels le flot $f$ est strictement inférieur à la capacité $c_{f}|_{E_{L}}$ ne contient aucun $s$-$t$-chemin. Autrement dit, un flot est bloquant si tout $s$-$t$-chemin contient un arc où le flot est égal à la capacité.

## L'Algorithme

Algorithme de Dinic

    Entrée: Un réseau $G=((V,E),c,s,t)$.
    Sortie: Un flot $f$ de valeur maximale.
* Poser $f(u,v)=0$ pour tout $(u,v)\in E$.
* Construire le graphe de niveau $G_{L}$ de $G_{f}$ à partir de $G$. Si le puits n'est pas accessible, terminer et retourner le flot $f$.
* Déterminer un flot bloquant $f\;'$ dans $G_{L}$.
* Augmenter le flot $\ f$ de $f\;'$ et retourner au point 2.

## Exemple

L'exemple suivant montre l'exécution de l'algorithme de Dinic. Dans le graphe de niveau $G_{L}$ les valeurs en rouge sont les distances à la source, les arcs en bleu forment un flot bloquant.

    $G$ $G_{f}$ $G_{L}$ 1. Au départ, le flot est nul, donc $G_{f}=G$. Le flot bloquant obtenu à la fin de la première étape a pour valeur 14. Il est composé de trois chemins augmentant, tous de longueur 3. Ce sont $\{s,1,3,t\}$ avec valeur 4, $\{s,1,4,t\}$ avec valeur 6, $\{s,2,4,t\}$ avec valeur 4. 2. Le graphe résiduel $G_{f}$ a des arcs opposés à ceux du graphe de départ, notamment lorsque la capacité est égale au flot. Les distances ne sont plus les mêmes. Le flot bloquant a pour valeur 5. Il correspond au chemin seul augmentant $\{s,2,4,3,t\}$ avec valeur 5. La valeur du flot devient 14 + 5 = 19. Le chemin augmentant est composé de 4 arcs. 3. Le sommet $t$ n'est plus accessible dans $G_{f}$. L'algorithme se termine et retourne un flot dont la valeur est 19. Notons que dans chaque flot bloquant, le nombre d'arcs du chemin augmentant croît d'au moins 1.

## Analyse

### Cas général

On peut montrer que le nombre d'arcs dans chaque flot bloquant croit d'au moins 1 à chaque étape; il y a donc au plus $n-1$ étapes, et au plus $n-1$flots bloquants à calculer, où $n$ est le nombre de sommets du réseau.

Le graphe de niveau $G_{L}$ peut être calculé par un parcours en largeur; un tel parcours prend un temps en $O(m)$, où $m$ est le nombre d'arcs. Un flot bloquant peut être calculé en temps $O(nm)$ . Le temps total de l'algorithme de Dinic est donc en $O(n^{2}m)$.

On peut diminuer le temps de calcul d'un flot bloquant en utilisant une structure de données de la famille des arbres dynamiques, soit les arbres adaptatifs (self-adjusting trees) soit les link/cut tree (en), inventés par Sleator et Tarjan. La détermination du flot bloquant de chaque étape peut être ramenée à $O(m\log n)$, et le temps total à $O(nm\log n)$.

### Cas particuliers

Capacités 0 ou 1. Dans un réseau où les capacités sont 0 ou 1, la complexité en temps peut être nettement mieux bornée. Chaque flot bloquant peut être trouvé en temps $O(m)$, et on peut montrer que, de plus, le nombre d'étapes est majoré par $O({\sqrt {m}})$ et $O(n^{2/3})$. L'algorithme est donc en temps $O(\min(n^{2/3},m^{1/2})m)$.

Réseaux de couplages. Dans les réseaux qui interviennent dans la résolution du couplage biparti, le nombre de phases est borné par $O({\sqrt {n}})$, ce qui donne un temps borné par $O({\sqrt {n}}m)$. L'algorithme est aussi connu sous le nom d'algorithme de Hopcroft-Karp.

Réseaux unitaires. Plus généralement, cette borne reste valable pour des réseaux dits unitaires, qui sont des réseaux où chaque sommet autre que la source et le puits possèdent soit un arc entrant de capacité un, soit un arc sortant de capacité un, les autres arcs ayant des capacités qui sont des entiers arbitraires.

## Historique

### Dinic ou Dinitz ?

L'algorithme de Dinic a été publié en 1970 par Yefim (Chaim) A. Dinitz (Ефим Диниц), alors qu'il était en Russie. Comme le dit Dinitz dans son article historique en hommage au mathématicien et informaticien Shimon Even, le rideau de fer séparait alors efficacement les chercheurs soviétiques des autres; son algorithme a été traduit dans la même année, et dans cette traduction, la lettre finale de son nom est remplacée par un « c ». C'est Shimon Even et Alon Itai, alors élève d'Even, qui ont découvert la traduction assez obscure de l'article, lui ont apporté des modifications et compléments, et en ont répandu la connaissance dans le monde scientifique occidental. Dinitz dit que l'algorithme ainsi reformulé - et qui circule encore aujourd'hui sous le nom d'algorithme de Dinic - doit son succès à Even et Atai, et que son algorithme original, algorithme de Dinitz, n'aurait peut-être pas eu ce succès. L'algorithme continue à être appelé algorithme de Dinic dans les ouvrages de référence, et par Dinitz lui-même, sur un ton amusé. Dinitz est maintenant membre du département d'informatique de l’université Ben Gourion du Néguev en Israël. Ses publications scientifiques paraissent toujours sous son vrai nom.

### Dinic et Edmonds-Karp

Deux ans plus tard, en 1972, paraît l'algorithme d'Edmonds-Karp mais qui avait déjà été présenté plus tôt. Ils ont montré indépendamment que si l'on choisit comme chemin augmentant le plus court à chaque pas, alors la longueur des chemins augmentant est croissante au sens large. L'algorithme est moins efficace que celui de Dinic parce qu'il ne fait pas usage de la notion de distance. Il avait déjà été présenté à une conférence aux États-Unis en 1968 (donc avant la publication de l'article de Dinic), mais pas publié à l'époque.

