Русская Википедия:Алгоритм Беллмана — Форда
Шаблон:Значения Шаблон:Алгоритм
Алгоритм Беллмана — Форда — алгоритм поиска кратчайшего пути во взвешенном графе. За время <math>O(|V|\cdot |E|)</math> алгоритм находит кратчайшие пути от одной вершины графа до всех остальных. В отличие от алгоритма Дейкстры, алгоритм Беллмана — Форда допускает рёбра с отрицательным весом. Предложен независимо Ричардом Беллманом и Лестером Фордом.
Алгоритм маршрутизации RIP (алгоритм Беллмана — Форда) был впервые разработан в 1969 году, как основной для сети ARPANET.
Формулировка задачи
Дан ориентированный или неориентированный граф <math>G</math> со взвешенными рёбрами. Длиной пути назовём сумму весов рёбер, входящих в этот путь. Требуется найти кратчайшие пути от выделенной вершины <math>s</math> до всех вершин графа.
Заметим, что кратчайших путей может не существовать. Так, в графе, содержащем цикл с отрицательным суммарным весом, существует сколь угодно короткий путь от одной вершины этого цикла до другой (каждый обход цикла уменьшает длину пути). Цикл, сумма весов рёбер которого отрицательна, называется отрицательным циклом.
Решение задачи на графе без отрицательных циклов
Решим поставленную задачу на графе, в котором заведомо нет отрицательных циклов.
Для нахождения кратчайших путей от одной вершины до всех остальных, воспользуемся методом динамического программирования. Построим матрицу <math>A_{ij}</math>, элементы которой будут обозначать следующее: <math>A_{ij}</math> — это длина кратчайшего пути из <math>s</math> в <math>i</math>, содержащего не более <math>j</math> рёбер.
Путь, содержащий 0 рёбер, существует только до вершины <math>s</math>. Таким образом, <math>A_{i0}</math> равно 0 при <math>i = s</math>, и <math>+\infty</math> в противном случае.
Теперь рассмотрим все пути из <math>s</math> в <math>i</math>, содержащие ровно <math>j</math> рёбер. Каждый такой путь есть путь из <math>j-1</math> ребра, к которому добавлено последнее ребро. Если про пути длины <math>j-1</math> все данные уже подсчитаны, то определить <math>j</math>-й столбец матрицы не составляет труда.
Так выглядит алгоритм поиска длин кратчайших путей в графе без отрицательных циклов:
for <math>v \in V</math>
do <math>d[v] \gets +\infty</math>
<math>d[s] \gets 0</math>
for <math>i \gets 1</math> to <math>|V| - 1</math>
do for <math>(u, v) \in E</math>
if <math>d[v] > d[u] + w(u, v)</math>
then <math>d[v] \gets d[u] + w(u, v)</math>
return <math>d</math>
Здесь <math>V</math> — множество вершин графа <math>G</math>, <math>E</math> — множество его рёбер, а <math>w</math> — весовая функция, заданная на рёбрах графа (возвращает длину дуги, ведущей из вершины <math>u</math> в <math>v</math>), <math>d</math> — массив, содержащий расстояния от вершины <math>s</math> до любой другой вершины.
Внешний цикл выполняется <math>|V| - 1</math> раз, поскольку кратчайший путь не может содержать большее число рёбер, иначе он будет содержать цикл, который точно можно выкинуть.
Вместо массива <math>d</math> можно хранить всю матрицу <math>A</math>, но это требует <math>O(V^2)</math> памяти. Зато при этом можно вычислить и сами кратчайшие пути, а не только их длины. Для этого заведем матрицу <math>P_{ij}</math>.
Если элемент <math>A_{ij}</math> содержит длину кратчайшего пути из <math>s</math> в <math>i</math>, содержащего <math>j</math> рёбер, то <math>P_{ij}</math> содержит предыдущую вершину до <math>i</math> в одном из таких кратчайших путей (ведь их может быть несколько).
Теперь алгоритм Беллмана-Форда выглядит так:
for <math>v \in V</math>
for <math>i \gets 0</math> to <math>|V| - 1</math>
do <math>A_{vi} \gets +\infty</math>
<math>A_{s0} \gets 0</math>
for <math>i \gets 1</math> to <math>|V| - 1</math>
do for <math>(u, v) \in E</math>
if <math>A_{vi} > A_{u, i-1} + w(u, v)</math>
then <math>A_{vi} \gets A_{u, i-1} + w(u, v)</math>
<math>P_{vi} \gets u</math>
После выполнения этого алгоритма элементы <math>A_{i, j}</math> содержат длины кратчайших путей от <math>s</math> до <math>i</math> с количеством рёбер <math>j</math>, и из всех таких путей следует выбрать самый короткий. А сам кратчайший путь до вершины <math>i</math> с <math>j</math> рёбрами восстанавливается так:
while <math>j > 0</math>
<math>p[j] \gets i</math>
<math>i \gets P_{ij}</math>
<math>j \gets j - 1</math>
return p
Граф с отрицательными циклами
Алгоритм Беллмана-Форда позволяет очень просто определить, существует ли в графе <math>G</math> отрицательный цикл, достижимый из вершины <math>s</math>. Достаточно произвести внешнюю итерацию цикла не <math>|V| - 1</math>, a ровно <math>|V|</math> раз. Если при исполнении последней итерации длина кратчайшего пути до какой-либо вершины строго уменьшилась, то в графе есть отрицательный цикл, достижимый из <math>s</math>. На основе этого можно предложить следующую оптимизацию: отслеживать изменения в графе и, как только они закончатся, сделать выход из цикла (дальнейшие итерации будут бессмысленны).
Литература
- R. Bellman: On a Routing Problem // Quarterly of Applied Mathematics. 1958. Vol 16, No. 1. C. 87-90, 1958.
- L. R. Ford, Jr., D. R. Fulkerson. Flows in Networks, Princeton University Press, 1962.
См. также
Ссылки
- Реализация алгоритма Форда-Беллмана на e-maxx.ru
- Реализация алгоритма Поиска отрицательного цикла на e-maxx.ru
- Псевдокод из книги «Introduction To Algorithms» (Thomas H Cormen,Charles E Leiserson,Ronald L Rivest,Clifford Stein)
Шаблон:Алгоритмы поиска на графах
