Программная реализация метода ветвей и границ

     Решение задачи

     Для удобства изложения везде ниже в  платежной матрице заменим имена городов (A, B, …, F) номерами соответствующих строк и столбцов (1, 2, …, 6).

     Найдем  нижнюю границу длин множества всех маршрутов. Вычтем из каждой строки число, равное минимальному элементу этой строки, далее вычтем из каждого столбца  число, равное минимальному элементу этого столбца, и таким образом приведем матрицу по строкам и столбцам. Минимумы по строкам: r₁=15, r₂=1, r₃=0, r₄=16, r₅=5, r₆=5.

     После их вычитания по строкам получим: 
 

 

     Минимумы  по столбцам: h₁=5, h₂=h₃=h₄=h₅=h₆.

     После их вычитания по столбцам получим  приведенную матрицу: 

     Найдем  нижнюю границу φ(Z) = 15+1+0+16+5+5+5 = 47.

     Для выделения претендентов на включение  во множество дуг, по которым производится ветвление, найдем степени Θ_ij нулевых элементов этой матрицы (суммы минимумов по строке и столбцу). Θ₁₄ = 10 + 0,  
Θ₂₄ = 1 + 0, Θ₃₆ = 5+0, Θ₄₁ = 0 + 1, Θ₄₂ = 0 + 0, Θ₅₆ = 2 + 0, Θ₆₂ = 0 + 0,  
Θ₆₃ = 0 + 9, Θ₆₅ = 0 + 2. Наибольшая степень Θ₁₄ = 10. Ветвление проводим по дуге (1, 4).

     Нижняя  граница для множества  остается равной 47. Для всех маршрутов множества из города A мы не перемещаемся в город D. В матрице это обозначается выставлением в ячейку (1, 4) знака ∞. В этом случае выход из города A добавляет к оценке нижней границы по крайней мере наименьший элемент первой строки. φ ( ) = 47 + 10.

     В матрице, соответствующей  полагаем c₁₄= ∞. 

 

     После проведения процедуры приведения с r₁=10 получим новую нижнюю границу 57 + 10 = 67.

     В матрице, соответствующей  , вычеркиваем первую строку и четвертый столбец и положим c₄₁= ∞, чтобы предотвратить появления цикла 1→ 4 → 1. Получим новую платежную матрицу {c¹_ij}: 

 

     Для приведения надо вычесть минимум  по первому столбцу: h₁=1. При этом нижняя граница станет равной 47+1 = 48. Сравнивая нижние границы 
φ ( ) = 67 и φ ( ) = 48 < 67 выделяем подмножество маршрутов , которое с большей вероятностью содержит маршрут минимальной длины.

     

     Рис. 1. Ветвление на первом шаге

 

      Приведенная платежная матрица для 

 

     Далее продолжаем процесс ветвления. Найдем степени Θ_ij нулевых элементов этой матрицы Θ₂₁ =16, Θ₃₆ = 5, Θ₄₂ = 2, Θ₅₆ = 2, Θ₆₂ = 0, Θ₆₃ =9, Θ₆₅ = 2. Наибольшая степень Θ₂₁. Затем множество разбивается дуге (2, 1) на два новых и .

     В матрице для  вычеркиваем строку 2 и столбец 1. дуги (1, 4) и (2, 1) образуют связный путь (2, 1, 4), положим c₄₂= ∞, чтобы предотвратить появления цикла 2→1→ 4 → 2. 

 

     Для приведения надо вычесть минимум  по строке 4: r₄=2. При этом нижняя граница станет равной 48+2 = 50.

     Нижняя  граница для , полученная как на предыдущем шаге ветвления, равна 48 + 16 = 64. Сравнивая нижние границы φ ( ) = 64 и φ ( ) = 50 < 64 выбираем для дальнейшего разбиения подмножество маршрутов .

     

     Рис. 2. Ветвление на втором шаге 

     Приведенная платежная матрица для   

 

     Степени Θ_ij нулевых элементов этой матрицы Θ₃₆ = 5, Θ₄₅ = 0, Θ₅₆ = 22, Θ₆₂ = 13, Θ₆₃ =7, Θ₆₅ = 0. Наибольшая степень Θ_56. Затем множество разбивается дуге (2, 1) на два новых и .

     Нижняя  граница для равна 50 + 22 = 72. В матрице для вычеркиваем строку 5 и столбец 6 и полагаем c₆₅= ∞. Получим матрицу: 

 

     Для приведения надо вычесть минимум  по строке 3: r₃=5. При этом нижняя граница станет равной 50+5 = 55. Выбираем для дальнейшего разбиения подмножество маршрутов .