Изменения

== Алгоритм ==Пусть дана [[Определение_сети,_потока#.D0.9E.D0.BF.D1.80.D0.B5.D0.B4.D0.B5.D0.BB.D0.B5.D0.BD.D0.B8.D0.B5_.D1.81.D0.B5.D1.82.D0.B8|сеть]] <tex> G </tex>, все рёбра которой имеют целочисленную [[Определение_сети,_потока#.D0.9E.D0.BF.D1.80.D0.B5.D0.B4.D0.B5.D0.BB.D0.B5.D0.BD.D0.B8.D0.B5_.D1.81.D0.B5.D1.82.D0.B8|пропускную способность]]. Обозначим за <tex> U </tex> максимальную пропускную способность: <tex> U = \max\limits_{{Определение(u, v) \in E} c(u, v) </tex>. Идея алгоритма заключается в нахождении путей с высокой пропускной способностью в первую очередь, чтобы сразу сильно увеличивать [[Определение_сети,_потока#.D0.9E.D0.BF.D1.80.D0.B5.D0.B4.D0.B5.D0.BB.D0.B5.D0.BD.D0.B8.D0.B5_.D0.BF.D0.BE.D1.82.D0.BE.D0.BA.D0.B0|definitionпоток]] по ним, а затем по всем остальным. Для этого воспользуемся масштабом <tex> \Delta </tex>. Изначально положим <tex> \Delta =2^{\lfloor \log_2 U \rfloor} </tex>.Алгоритм масштабирования потока — На каждой итерации в [[Дополняющая_сеть,_дополняющий_путь|дополняющей сети]] алгоритм поиска максимального потока путём регулирования находит [[Дополняющая_сеть,_дополняющий_путь|дополняющие пути]] с пропускной способности рёберспособностью не меньшей <tex> \Delta </tex> и увеличивает поток вдоль них.Этот алгоритм работает Уменьшив масштаб <tex> \Delta </tex> в предположении, что все пропускные способности рёбер целые<tex> 2 </tex> раза, так как они легко представимы в двоичном видепереходит к следующей итерации.}}

Очевидно, что при <tex> \Delta == Идея ==Идея алгоритма 1 </tex> алгоритм вырождается в нахождении путей с высокой пропускной способностью в первую очередь, чтобы сразу сильно увеличивать поток по нималгоритм [[Алоритм_Эдмондса-Карпа|Эдмондса-Карпа]], а затем по всем остальнымвследствие чего является корректным.

Пусть <tex> G </tex> — графКоличество необходимых увеличений путей, <tex> \forall(uоснованных на кратчайших путях, v) \in EG \colon c(uможет быть много больше количества увеличений,v) \in \mathbbоснованных на путях с высокой пропускной способностью.{Z_+}, U |border="0" cellpadding="5" width=30% align= \max\limits_{(u, v) \in EG} c(u, v) </tex> — максимальная пропускная способностьcenter|[[Файл:Flow_scale_1. Запишем пропускную способность каждого ребра png|550px|thumb|center|Выбор дополняющих путей в двоичном виде. Тогда каждое число будет занимать <tex> \lfloor \log_2 U \rfloor + 1 = n + 1 </tex> битпорядке длины]]|[[Файл:Flow_scale_2.png|550px|thumb|center|Выбор пути с высокой пропускной способностью в первую очередь]]|}

== Оценка времени работы =={{Лемма|about=1|statement=Максимальный поток в сети <tex> c(uG </tex> ограничен сверху значением <tex> |f_k| + 2^k E </tex>, v) = \sum\limits_где <tex> |f_k| </tex> {{i = 0---}}^n a_i(u, v) значение потока при масштабе <tex> \times Delta = 2^n, a_i(u, v) \in \{0, 1\} k </tex>.|proof=[[Файл:Flow_scale_3.png|530px|thumb|right|Разрез <tex> C_k </tex>]]

Методом Форда-Фалкерсона находим поток В конце итерации с масштабом <tex> f_0 \Delta = 2^k </tex> для графа с урезанными пропускными способностями , сеть <tex> G_{f_k} </tex> может быть разбита на два непересекающихся множества <tex> A_k </tex> и <tex> \overline{A_k} </tex> c_0(uтак, v) = a_n(uчто остаточная пропускная способность каждого ребра, идущего из <tex> A_k </tex> в <tex> \overline{A_k} </tex>, v) не превосходит масштаба <tex> \Delta </tex>.Добавим следующий бит и находим следующее приближение для графа с новыми пропускными способностями То есть образуется [[Разрез,_лемма_о_потоке_через_разрез|разрез]] <tex> c_1(u, v) C_k = 2 a_n(u\langle A_k, v) + a_\overline{n - 1A_k}(u, v) - 2 f_0(u, v) \rangle </tex>.

После При этом, количество таких рёбер не превосходит <tex> n + 1 E </tex> итерации получим ответ к задаче.Значит, так как после с каждым шагом приближение становится точнеезначение остаточного потока не может превосходить <tex> \Delta E = 2^k E </tex>.}}

{{Лемма|about=2|statement= Оценка сложности =Суммарное количество увеличивающих путей {{---}} <tex> O(E \log U) </tex>.|proof=На некоторой итерации алгоритма каждый дополняющий путь имеет пропускную способность не меньше <tex> 2^k </tex>.Дополняющий поток на предыдущем шаге ограничен значением <tex> 2^{k + 1} E </tex>. Следовательно, на каждой итерации количество дополняющих путей не превосходит <tex> 2E </tex>.}}

Время работы алгоритма — {{---}} <tex> O(E^2 \log U) </tex>.

[[ФайлВ ходе выполнения алгоритма масштаб <tex> \Delta </tex> принимает следующие значения:Scaling<tex> S = \{2^{\lfloor \log_2 U \rfloor}, \ldots, 2^k, \ldots, 2, 1, 0\} </tex>.jpgТогда <tex> |right]]S| = O(\log U) </tex> {{---}} количество итераций алгоритма.

Количество итераций — алгоритма {{---}} <tex> O(\log U) </tex>. Докажем, что сложность каждой итерации — значит, суммарное количество увеличивающих путей {{---}} <tex> O(E^2\log U) </tex>.

На каждом шаге алгоритм выполняет Алгоритм [[Обход_в_ширину|обхода в худшем случае ширину]] находит каждый дополняющий путь за время <tex>O(E)</tex> увеличений потока. Докажем это. Пусть <tex>\Delta = 2^k</tex>. В конце шага множество вершин графа можно разбить на две части: <tex>A_k</tex> и <tex>\overlineСледовательно, суммарное время работы алгоритма {A_k}</tex>. Все рёбра, выходящие из <tex>A_k</tex>, имеют остаточную пропускную способность менее <tex>2^k</tex>. Наибольшее количество рёбер между <tex>A_k</tex> и <tex>\overline{A_k---}</tex> равно <tex>E</tex>. Следовательно, остаточный поток (поток, который может быть получен на оставшихся шагах) на фазе с текущим значением <tex>k</tex> максимально составляет <tex>2^kE</tex>. Каждый увеличивающий путь при данном <tex>k</tex> имеет пропускную способность как минимум <tex>2^k</tex>. На предыдущем шаге, с масштабом <tex>k+1</tex>, остаточный поток ограничен <tex>2^{k+1}E</tex>. Значит максимальное число появившихся увеличивающих путей равно <tex>2E</tex>. Увеличивающий путь можно найти за <tex>O(E)</tex>, используя [[Обход_в_ширину | BFS]]. Количество шагов <tex>O(\log_2U)</tex>. Итоговая сложность <tex>O(E^2\log_2Ulog U)</tex>.}}

'''Capacity-Scalingfunction''' maxFlowByScaling(G: '''graph''', s: '''int''', t: '''int'''): '''int''' '''int''' flow = 0 <texfont color=darkgreen> f \leftarrow 0 // поток в сети </texfont> '''int''' scale = <tex> \Delta \leftarrow 2^{\lfloor \log_2 U log_2U\rfloor}</tex> <font color=darkgreen> // текущий минимальный размер потока, который пытаемся пустить </font> '''while''' scale <tex> \Delta >0geqslant </tex>1 '''do''' '''while''' в <tex>G_f</tex> существует увеличивающий путь <tex>s-tp </tex> путь с пропускной способностью большей <tex>\Delta</tex>не меньше, чем scale '''doint''' <tex>P\leftarrow</tex> путь с пропускной способностью большей minCapacity = <tex>\Delta</tex> <tex>\delta\leftarrow\min\{c_{ij}c(u, v) \colon(iu,jv)\in Pp\}</tex> <font color=darkgreen> // минимальная пропускная способность в увеличивающем пути </font> увеличить поток по рёбрам <tex>Pp </tex> на <tex>\delta</tex>minCapacity обновить <tex>G_f</tex> <tex>f\leftarrow f flow = flow +\delta</tex>minCapacity <tex>\Delta\leftarrow\Delta scale = scale /2</tex> '''return''' flow

== Литература См. также ==* [[Определение_сети,_потока|Определение сети, потока]]* [[Алоритм_Эдмондса-Карпа|Алоритм Эдмондса-Карпа]]* [[Алгоритм_Форда-Фалкерсона,_реализация_с_помощью_поиска_в_глубину|Алгоритм Форда-Фалкерсона]] == Источники информации ==

* [http://wwwlogic.cs-seminarpdmi.spbras.ru/reportsics/34talks/21stream.pdf ''Андрей Станкевич'': Задача о максимальном потоке]* [https://youtu.be/sEwp5ZAJJps?t=18m9s ''Андрей Станкевич'': Лекториум, дополнительные главы алгоритмов, лекция 12] [[Категория: Алгоритмы и структуры данных]][[Категория: Задача о максимальном потоке]]