Алгоритм Форда — Фалкерсона

Алгоритм Форда — Фалкерсона решает задачу нахождения максимального потока в транспортной сети.

Идея алгоритма заключается в следующем. Изначально величине потока присваивается значение 0: $\text{[math]}$ $\text{[math]}$ $f(u,v)=0$ $f(u,v)=0$ для всех $\text{[math]}$ $\text{[math]}$ $u,v\in V$ $u,v\in V$ . Затем величина потока итеративно увеличивается посредством поиска увеличивающего пути (путь от источника $\text{[math]}$ $\text{[math]}$ $s$ $s$ к стоку $\text{[math]}$ $\text{[math]}$ $t$ $t$ , вдоль которого можно послать больший поток). Процесс повторяется, пока можно найти увеличивающий путь.

АлгоритмПравить

Неформальное описаниеПравить

Обнуляем все потоки. Остаточная сеть изначально совпадает с исходной сетью.
В остаточной сети находим любой путь из источника в сток. Если такого пути нет, останавливаемся.
Пускаем через найденный путь (он называется увеличивающим путём или увеличивающей цепью) максимально возможный поток:
1. На найденном пути в остаточной сети ищем ребро с минимальной пропускной способностью $\text{[math]}$ $\text{[math]}$ $c_{\min }$ .
2. Для каждого ребра на найденном пути увеличиваем поток на $\text{[math]}$ $\text{[math]}$ $c_{\min }$ , а в противоположном ему — уменьшаем на $\text{[math]}$ $\text{[math]}$ $c_{\min }$ .
3. Модифицируем остаточную сеть. Для всех рёбер на найденном пути, а также для противоположных (антипараллельных) им рёбер, вычисляем новую пропускную способность. Если она стала ненулевой, добавляем ребро к остаточной сети, а если обнулилась, стираем его.
Возвращаемся на шаг 2.

Важно, что алгоритм не конкретизирует, какой именно путь мы ищем на шаге 2 или как мы это делаем. По этой причине алгоритм гарантированно сходится только для целых пропускных способностей, но даже для них при больших значениях пропускных способностей он может работать очень долго. Если пропускные способности вещественны, алгоритм может работать бесконечно долго, не сходясь к оптимальному решению (см. пример ниже).

Если искать не любой путь, а кратчайший, получится алгоритм Эдмондса — Карпа или алгоритм Диница. Эти алгоритмы сходятся для любых вещественных весов за время $\text{[math]}$ $\text{[math]}$ $O(|V||E|^{2})$ и $\text{[math]}$ $\text{[math]}$ $O(|V|^{2}|E|)$ соответственно.

Формальное описаниеПравить

Дан граф $\text{[math]}$ $\text{[math]}$ $G(V,E)$ с пропускной способностью $\text{[math]}$ $\text{[math]}$ $c(u,v)$ и потоком $\text{[math]}$ $\text{[math]}$ $f(u,v)=0$ для рёбер из $\text{[math]}$ $\text{[math]}$ $u$ в $\text{[math]}$ $\text{[math]}$ $v$ . Необходимо найти максимальный поток из источника $\text{[math]}$ $\text{[math]}$ $s$ в сток $\text{[math]}$ $\text{[math]}$ $t$ . На каждом шаге алгоритма действуют те же условия, что и для всех потоков:

$\text{[math]}$ $\text{[math]}$ $\ 0\leqslant f(u,v)\leqslant c(u,v)$ . Поток из $\text{[math]}$ $\text{[math]}$ $u$ в $\text{[math]}$ $\text{[math]}$ $v$ не превосходит пропускной способности.
$\text{[math]}$ $\text{[math]}$ $\ f(u,v)=-f(v,u)$ .
$\text{[math]}$ $\text{[math]}$ $\ \sum _{v}f(u,v)=0\Longleftrightarrow f_{in}(u)=f_{out}(u)$ для всех узлов $\text{[math]}$ $\text{[math]}$ $u$ , кроме $\text{[math]}$ $\text{[math]}$ $s$ и $\text{[math]}$ $\text{[math]}$ $t$ . Поток не изменяется при прохождении через узел.

Остаточная сеть $\text{[math]}$ $\text{[math]}$ $G_{f}(V,E_{f})$ — сеть с пропускной способностью $\text{[math]}$ $\text{[math]}$ $c_{f}(u,v)=c(u,v)-f(u,v)+f(v,u)$ и без потока.

Вход Граф $\text{[math]}$ $\text{[math]}$ $G$ с пропускной способностью $\text{[math]}$ $\text{[math]}$ $c$ , источник $\text{[math]}$ $\text{[math]}$ $s$ и сток $\text{[math]}$ $\text{[math]}$ $t$
Выход Максимальный поток $\text{[math]}$ $\text{[math]}$ $f$ из $\text{[math]}$ $\text{[math]}$ $s$ в $\text{[math]}$ $\text{[math]}$ $t$

$\text{[math]}$ $\text{[math]}$ $f(u,v):=0$ для всех рёбер $\text{[math]}$ $\text{[math]}$ $(u,v)$
Пока есть путь p из s в t в G f , такой что c f ( u , v ) > 0 для всех рёбер ( u , v ) ∈ p :
1. Найти $\text{[math]}$ $\text{[math]}$ $c_{f}(p)=\min\{c_{f}(u,v)\;|\;(u,v)\in p\}$
2. Для каждого ребра ( u , v ) ∈ p
  1. $\text{[math]}$ $\text{[math]}$ $f(u,v):=f(u,v)+c_{f}(p)$
  2. $\text{[math]}$ $\text{[math]}$ $f(v,u):=f(v,u)-c_{f}(p)$
  3. $\text{[math]}$ $\text{[math]}$ $c_{f}(u,v)=c(u,v)-f(u,v)$
  4. $\text{[math]}$ $\text{[math]}$ $c_{f}(v,u)=f(u,v)$

Путь может быть найден, например, поиском в ширину (алгоритм Эдмондса — Карпа) или поиском в глубину в $\text{[math]}$ $\text{[math]}$ $G_{f}(V,E_{f})$ .

СложностьПравить

На каждом шаге алгоритм добавляет поток увеличивающего пути к уже имеющемуся потоку. Если пропускные способности всех рёбер — целые числа, легко доказать по индукции, что и потоки через все рёбра всегда будут целыми. Следовательно, на каждом шаге алгоритм увеличивает поток по крайней мере на единицу, следовательно, он сойдётся не более чем за $\text{[math]}$ $\text{[math]}$ $O(|f|)$ шагов, где f — максимальный поток в графе. Можно выполнить каждый шаг за время $\text{[math]}$ $\text{[math]}$ $O(|E|)$ , где $\text{[math]}$ $\text{[math]}$ $|E|$ — число рёбер в графе, тогда общее время работы алгоритма ограничено $\text{[math]}$ $\text{[math]}$ $O(|E|\max |f|)$ .

Если величина пропускной способности хотя бы одного из рёбер — иррациональное число, то алгоритм может работать бесконечно, даже не обязательно сходясь к правильному решению. Пример дан ниже.

Пример несходящегося алгоритмаПравить

Рассмотрим приведённую справа сеть, с источником $\text{[math]}$ $\text{[math]}$ $s$ , стоком $\text{[math]}$ $\text{[math]}$ $t$ , пропускными способностями рёбер $\text{[math]}$ $\text{[math]}$ $e_{1}=1$ , $\text{[math]}$ $\text{[math]}$ $e_{2}=r=({\sqrt {5}}-1)/2$ , $\text{[math]}$ $\text{[math]}$ $e_{3}=1$ и пропускной способностью всех остальных рёбер, равной целому числу $\text{[math]}$ $\text{[math]}$ $M\geq 2$ . Константа $\text{[math]}$ $\text{[math]}$ $r$ выбрана так, что $\text{[math]}$ $\text{[math]}$ $r^{2}=1-r$ . Мы используем пути из остаточного графа, приведённые в таблице, причём $\text{[math]}$ $\text{[math]}$ $p_{1}=\{s,v_{4},v_{3},v_{2},v_{1},t\}$ , $\text{[math]}$ $\text{[math]}$ $p_{2}=\{s,v_{2},v_{3},v_{4},t\}$ и $\text{[math]}$ $\text{[math]}$ $p_{3}=\{s,v_{1},v_{2},v_{3},t\}$ .

Шаг	Найденный путь	Добавленный поток	Остаточные пропускные способности
$\text{[math]}$ $\text{[math]}$ $e_{1}$	$\text{[math]}$ $\text{[math]}$ $e_{2}$	$\text{[math]}$ $\text{[math]}$ $e_{3}$
0			$\text{[math]}$ $\text{[math]}$ $r^{0}=1$	$\text{[math]}$ $\text{[math]}$ $r$	$\text{[math]}$ $\text{[math]}$ $1$
1	$\text{[math]}$ $\text{[math]}$ $\{s,v_{2},v_{3},t\}$	$\text{[math]}$ $\text{[math]}$ $1$	$\text{[math]}$ $\text{[math]}$ $r^{0}$	$\text{[math]}$ $\text{[math]}$ $r^{1}$	$\text{[math]}$ $\text{[math]}$ $0$
2	$\text{[math]}$ $\text{[math]}$ $p_{1}$	$\text{[math]}$ $\text{[math]}$ $r^{1}$	$\text{[math]}$ $\text{[math]}$ $r^{2}$	$\text{[math]}$ $\text{[math]}$ $0$	$\text{[math]}$ $\text{[math]}$ $r^{1}$
3	$\text{[math]}$ $\text{[math]}$ $p_{2}$	$\text{[math]}$ $\text{[math]}$ $r^{1}$	$\text{[math]}$ $\text{[math]}$ $r^{2}$	$\text{[math]}$ $\text{[math]}$ $r^{1}$	$\text{[math]}$ $\text{[math]}$ $0$
4	$\text{[math]}$ $\text{[math]}$ $p_{1}$	$\text{[math]}$ $\text{[math]}$ $r^{2}$	$\text{[math]}$ $\text{[math]}$ $0$	$\text{[math]}$ $\text{[math]}$ $r^{3}$	$\text{[math]}$ $\text{[math]}$ $r^{2}$
5	$\text{[math]}$ $\text{[math]}$ $p_{3}$	$\text{[math]}$ $\text{[math]}$ $r^{2}$	$\text{[math]}$ $\text{[math]}$ $r^{2}$	$\text{[math]}$ $\text{[math]}$ $r^{3}$	$\text{[math]}$ $\text{[math]}$ $0$

Заметим, что после шага 1, как и после шага 5, остаточные способности рёбер $\text{[math]}$ $\text{[math]}$ $e_{1}$ , $\text{[math]}$ $\text{[math]}$ $e_{2}$ и $\text{[math]}$ $\text{[math]}$ $e_{3}$ имеют форму $\text{[math]}$ $\text{[math]}$ $r^{n}$ , $\text{[math]}$ $\text{[math]}$ $r^{n+1}$ и $\text{[math]}$ $\text{[math]}$ $0$ , соответственно, для какого-то натурального $\text{[math]}$ $\text{[math]}$ $n$ . Это значит, что мы можем использовать увеличивающие пути $\text{[math]}$ $\text{[math]}$ $p_{1}$ , $\text{[math]}$ $\text{[math]}$ $p_{2}$ , $\text{[math]}$ $\text{[math]}$ $p_{1}$ и $\text{[math]}$ $\text{[math]}$ $p_{3}$ бесконечно много раз, и остаточные пропускные способности этих рёбер всегда будут в той же форме. Полный поток после шага 5 равен $\text{[math]}$ $\text{[math]}$ $1+2(r^{1}+r^{2})$ . За бесконечное время полный поток сойдётся к $\text{[math]}$ $\text{[math]}$ $\textstyle 1+2\sum _{i=1}^{\infty }r^{i}=3+2r$ , тогда как максимальный поток равен $\text{[math]}$ $\text{[math]}$ $2M+1$ . Таким образом, алгоритм не только работает бесконечно долго, но даже и не сходится к оптимальному решению.^[1]

ПримерПравить

Следующий пример показывает первые шаги алгоритма Форда — Фалкерсона в транспортной сети с четырьмя узлами, источником A и стоком D. Этот пример показывает худший случай. При использовании поиска в ширину алгоритму потребуется всего лишь 2 шага.

Путь	Пропускная способность	Результат
	Начальная транспортная сеть
$\text{[math]}$ $\text{[math]}$ $A, B, C, D$	$\text{[math]}$ $\text{[math]}$ $\min(c_{f}(A,B),c_{f}(B,C),c_{f}(C,D))=$ $\text{[math]}$ $\text{[math]}$ $\min(c(A,B)-f(A,B),c(B,C)-f(B,C),c(C,D)-f(C,D))=$ $\text{[math]}$ $\text{[math]}$ $\min(1000-0,1-0,1000-0)=1$
$\text{[math]}$ $\text{[math]}$ $A, C, B, D$	$\text{[math]}$ $\text{[math]}$ $\min(c_{f}(A,C),c_{f}(C,B),c_{f}(B,D))=$ $\text{[math]}$ $\text{[math]}$ $\min(c(A,C)-f(A,C),c(C,B)-f(C,B),c(B,D)-f(B,D))=$ $\text{[math]}$ $\text{[math]}$ $\min(1000-0,0-(-1),1000-0)=1$
	После 1998 шагов …
	Конечная сеть

См. такжеПравить

Теорема Форда — Фалкерсона

СсылкиПравить

Пример решения задачи на YouTube
Визуализатор алгоритма
Реализация поиска максимального потока методом Форда-Фалкерсона на Java
На Викискладе есть медиафайлы по теме Алгоритм Форда — Фалкерсона

ЛитератураПравить

Томас Х. Кормен и др. Алгоритмы: построение и анализ = INTRODUCTION TO ALGORITHMS. — 2-е изд. — М.: «Вильямс», 2006. — С. 1296. — ISBN 0-07-013151-1.

ПримечанияПравить

↑ Zwick, Uri. The smallest networks on which the Ford-Fulkerson maximum flow procedure may fail to terminate (англ.) // Theoretical Computer Science (англ.) (рус. : journal. — 1995. — 21 August (vol. 148, no. 1). — P. 165—170. — doi:10.1016/0304-3975(95)00022-O.

[1] Zwick, Uri. The smallest networks on which the Ford-Fulkerson maximum flow procedure may fail to terminate (англ.) // Theoretical Computer Science (англ.) (рус. : journal. — 1995. — 21 August (vol. 148, no. 1). — P. 165—170. — doi:10.1016/0304-3975(95)00022-O.

[1]