Алгоритм нахождения максимального пути

 

При решении некоторых практических задач возникает необходимость поиска максимального пути (пути с наибольшей суммой длин дуг). Такая задача сводится к задаче нахождения минимального пути заменой знаков при длинах дуг (в матрице весов C) на противоположные. При этом необходимым является требование отсутствия в ориентированном графе контуров положительной длины.

Пример 3.16.

С помощью модифицированного алгоритма 3.1 найдем максимальный путь из верши­ны х ₁ в вершину х ₃ в графе, изображенном на рис. 3.11.

Рис. 3.11

Шаг 1. Введем число вершин графа n =5. Матрица весов этого графа после замены знаков при длинах дуг на противоположные имеет вид:

C = .

Шаг 2. Положим k = 0, l ₁(0) = 0, l ₂(0) = l ₃(0) = l ₄(0) = l ₅(0) = . Эти значения занесем в первый столбец табл. 3.2.

Шаг 3.

k = 1.

l ₁(1) = 0.

Равенство (3.1) для k = 1 имеет вид:

l_i (1) = { l_j (0) + c_ji }.

l ₂(1) = min { l ₁(0) + c ₁₂; l ₂(0) + c ₂₂; l ₃(0) + c ₃₂; l ₄(0) + c ₄₂; l ₅(0) + c ₅₂;} = min {0 – 1; ¥ + ¥; ¥ + ¥; ¥ + ¥; ¥ + ¥} = –1.

l ₃(1) = min { l ₁(0) + c ₁₃; l ₂(0) + c ₂₃; l ₃(0) + c ₃₃; l ₄(0) + c ₄₃; l ₅(0) + c ₅₃;} = min {0 + ¥; ¥ – 8; ¥ + ¥; ¥ – 2; ¥ + ¥} = ¥.

l ₄(1) = min { l ₁(0) + c ₁₄; l ₂(0) + c ₂₄; l ₃(0) + c ₃₄; l ₄(0) + c ₄₄; l ₅(0) + c ₅₄;} = min {0 + ¥; ¥ – 7; ¥ + ¥; ¥ + ¥; ¥ – 4} = ¥.

l ₅(1) = min { l ₁(0) + c ₁₅; l ₂(0) + c ₂₅; l ₃(0) + c ₃₅; l ₄(0) + c ₄₅; l ₅(0) + c ₅₅;} = min {0 – 3; ¥ – 1; ¥ + 5; ¥ + ¥; ¥ + ¥} = –3.

Полученные значения l_i (1) занесем во второй столбец табл. 3.2. Убеждаемся, что второй столбец, начиная со второго элемента, совпадает с первой строкой матрицы весов, что легко объясняется смыслом величин l_i (1), которые равны длине минимального пути из первой вершины в i -ую, содержащего не более одной дуги.

k = 2.

l ₁(2) = 0.

Равенство (3.1) для k = 2 имеет вид:

l_i (2) = { l_j (1) + c_ji }.

l ₂(2) = min {0 – 1; –1 + ¥; ¥ + ¥; ¥ + ¥; –3 + ¥} = –1.

l ₃(2) = min {0 + ¥; –1 – 8; ¥ + ¥; ¥ – 2; –3 + ¥} = –9.

l ₄(2) = min {0 + ¥; –1 – 7; ¥ + ¥; ¥ + ¥; –3 – 4} = –8.

l ₅(2) = min {0 – 3; –1 – 1; ¥ + 5; ¥ + ¥; –3 + ¥} = –3.

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

k = 3.

l ₁(3) = 0.

Равенство (3.1) для k = 3 имеет вид:

l_i (3) = { l_j (2) + c_ji }.

l ₂(3) = min {0 – 1; – 1 + ¥; – 9 + ¥; –8 + ¥; – 3 + ¥} = – 1.

l ₃(3) = min {0 + ¥; – 1 – 8; – 9 + ¥; –8 – 2; – 3 + ¥} = – 10.

l ₄(3) = min {0 + ¥; – 1 – 7; – 9 + ¥; –8 + ¥; – 3 – 4} = – 8.

l ₅(3) = min {0 – 3; – 1 – 1; – 9 + 5; –8 + ¥; – 3 + ¥} = – 4.

Полученные значения l_i (3) занесем в четвертый столбец табл. 3.2. Величины l_i (3) равны длине минимального пути из первой вершины в i -ую, содержащего не более трех дуг.

k = 4.

l ₁(4) = 0.

Равенство (3.1) для k = 4 имеет вид:

l_i (4) = { l_j (3) + c_ji }.

l ₂(4) = min {0 – 1; – 1 + ¥; – 10 + ¥; – 8 + ¥; – 4 + ¥} = – 1.

l ₃(4) = min {0 + ¥; – 1 – 8; – 10 + ¥; – 8 – 2; – 4 + ¥} = – 10.

l ₄(4) = min {0 + ¥; – 1 – 7; – 10 + ¥; – 8 + ¥; – 4 – 4} = – 8.

l ₅(4) = min {0 – 3; – 1 – 1; – 10 + 5; – 8 + ¥; – 4 + ¥} = – 5.

Полученные значения l_i (4) занесем в пятый столбец табл. 3.2. Величины l_i (4) равны длине минимального пути из первой вершины в i -ую, содержащего не более четырех дуг.

Таблица 3.2

i (номер вершины)	li (0) li (1) li (2) li (3) li (4)
	0 0 0 0 0 ¥ – 1 – 1 – 1 1 ¥ ¥ – 9 – 10 – 10 ¥ ¥ – 8 – 8 – 8 ¥ – 3 –3 – 4 – 5

Заменив в табл. 3.2 отрицательные числа положительными, получим таблицу индексов максимальных путей (табл. 3.3). При этом l_i (k) определяет длину максимального пути из первой вершины в i -ую, содержащего не более k дуг.

Таблица 3.3

i (номер вершины)	li (0) li (1) li (2) li (3) li (4)
	0 0 0 0 0 ¥ 1 1 1 1 ¥ ¥ 9 10 10 ¥ ¥ 8 8 8 ¥ 3 3 4 5

Шаг 5. Восстановление максимального пути производится по тому же правилу, что и для минимального пути.

Длина максимального пути равна 10. Этот путь состоит из трех дуг, т. к. l_i (3) = l_i (4) = 10. Поэтому в соотношении (3.2) будет выполнено, начиная с n – 1.

Учитывая это замечание, для последней вершины x ₃предшествующую ей вершину x_r определим из соотношения (3.2) полученного при s =3:

l_r (2) + c_{r 3} = l ₃(3), (3.7)

x_r Î G^{- 1}(x ₃), где G^{- 1}(x ₃) - прообраз вершины x ₃.

G^{- 1}(x ₃)= { x ₂, x ₄}.

Подставим в (3.7) последовательно r = 2 и r = 4, чтобы определить, для какого r это равенство выполняется:

l ₂(2) + c ₂₃ = 1 + 8 ¹ l ₃(4) = 10,

l ₄(2) + c ₄₃ = 8 + 2 = l ₃(4) = 10.

Таким образом, вершиной, предшествующей вершине x ₃, является вершина x ₄.

Для вершины x ₄предшествующая ей вершина x_r определяется из соотношения (3.2) полученного при s =4:

l_r (1) + c_{r 4} = l ₄(2), x_r Î G^{- 1}(x ₄), (3.8)

где G^{- 1}(x ₄) - прообраз вершины x ₄.

G^{- 1}(x ₄)= { x ₂, x ₅}.

Подставим в (3.8) последовательно r = 2, r = 3 и r = 5, чтобы определить, для какого r это равенство выполняется:

l ₂(1) + c ₂₄ = 1 + 7 = l ₄(3) = 8,

l ₅(1) + c ₅₄ = 3 + 4 ¹ l ₄(3) = 8,

Таким образом, вершиной, предшествующей вершине x ₄, является вершина x ₂.

Для вершины x ₂предшествующая ей вершина x_r определяется из соотношения (3.2) полученного при s =2:

l_r (0) + c_{r 2} = l ₂(1), x_r G^{- 1}(x ₂), (3.9)

где G^{- 1}(x ₂) - прообраз вершины x ₂.

G^{- 1}(x ₂)= { x ₁}.

Подставим в (3.9) r = 1, чтобы определить, выполняется ли это равенство:

l ₁(1) + c ₁₂ = 0 + 1 = l ₂(1) = 1.

Таким образом, вершиной, предшествующей вершине x ₂, является вершина x ₁.

Итак, найден максимальный путь – x ₁, x ₂, x ₄, x ₃, его длина равна 10.