25. Задача коммивояжера. Решение задачи методом ветвей и границ.

задача коммивояжера.

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

В общем случае предполагается, что c_ij  c_ji, то есть граф может быть и ориентированным, причем веса дуг (u,v) и (v,u) могут отличаться.

Решение задачи «в лоб» предполагает перебор всех возможных вариантов, то есть p! операций. Но если граф большой, это слишком дорогое удовольствие, поэтому разработаны приемы сокращения перебора. Наиболее распространенный из таких методов ограниченного перебора носит название «метод ветвей и границ».

Решение задачи коммивояжера методом ветвей и границ.

Рассмотрим дерево перебора (дерево решений) и назовем его T, как мы называли упорядоченные корневые деревья. Общий принцип решения методом ветвей и границ такой: на каждом шаге дерево решений разбивается на подмножества ветвей Затем вычисляется некоторая оценка, которая служит критерием отсечения заведомо бесперспективных ветвей.

	1	3	1	2
2		5	4	6
2	4		6	1
3	1	5		4
3	2	4	1	

Рассмотрим, как это выглядит в данной задаче. Будем иллюстрировать рассуждения на полном пятивершинном графе (точнее, ориентированном мультиграфе, где c_ij  c_ji) с матрицей весов C. Очевидно, что этот граф гамильтонов. Матрица весов у нас будет следующая:

Основной операцией, на основе которой выполняется необходимая оценка, является приведение матрицы весов. Что это означает? Рассмотрим по шагам.

Введем обозначения: z – искомый гамильтонов цикл, l(z) – вес такого гамильтонова цикла.

Шаг 1. Получение нижней l₀(z) и верхней L(z) оценок веса искомого гамильтонова цикла (построение границ решения).

Начнем с получения нижней оценки, так как именно она наиболее важна для решения задачи.

1.1. Выбираем min по каждой строке матрицы весов. Строка i у нас, как всегда, соответствует выходам из вершины v_i, а столбец i - входам. То есть мы выбираем сначала минимальный по весу выход из каждой вершины.

Теперь из каждой строки матрицы вычтем соответствующее этой строке минимальное значение.

Получаем матрицу:

	0	2	0	1
0		3	2	4
1	3		5	0
2	0	4		3
2	1	3	0	

Если бы в каждом столбце после вычитания оказалось по нулю, это означало бы, что из каждой вершины можно выйти по дуге минимального веса и в каждую вершину можно войти по дуге минимального веса. В этом случае вес гамильтонова цикла был бы равен l₀(z)=6, то есть той сумме, которую мы получили. И задача была бы решена уже на первом шаге. Но у нас в третьем столбце нулей нет. Поэтому придется принять, что l(z)> l₀(z), и продолжать поиск l(z).

1.2. Выбираем min по каждому столбцу, то есть рассматриваем все дуги, по которым можно войти в каждую вершину v_j. Повторяем то, что делали на первом шаге, только теперь для столбцов.

	0	0	0	1
0		1	2	4
1	3		5	0
2	0	2		3
2	1	1	0	

Теперь вычитаем из столбцов минимальные значения.

Полученный результат показывает, что минимальная дуга, по которой мы можем войти в вершину 3, - это c₁₃=2. Теперь уже ясно, что l(z) l₀(z)+ c₁₃=8.

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

1.3. Для получения верхней оценки L(z) построим произвольный цикл методом поиска в глубину и посчитаем его вес по исходной матрице весов.

L(z)= c₁₂+ c₂₃+ c₃₄+ c₄₅+ c₅₁=1+5+6+4+3=19.

Если бы мы получили, что l(z)=L(z), то это и был бы искомый гамильтонов цикл минимального веса. Но в нашем примере это не так. Мы, всего-навсего, имеем 8l(z)19, и нам придется продолжать поиск.

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

Рассмотрим два случая(рис.3.27).: 1) гамильтонов цикл z проходит через некоторую минимальную по весу дугу (v_i,v_j); 2) гамильтонов цикл z не проходит через эту дугу. Назовем эти два варианта z(i,j) и соответственно. Они образуют две ветви решений.

Рассмотрим оценки для каждой из этих ветвей. Для ветви z(i,j) остается нижняя оценка l(z)=8. Рассмотрим теперь оценку для ветви . Если путь не проходит через (i,j), то мы должны выйти из вершины i по какой-то дуге (i,k) и войти в вершину j по какой-то дуге (m,j). В наиболее удачном случае это будут дуги минимальных весов. Тогда

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

Вычисляем такие оценки для всех дуг, обнуленных в процессе приведения матрицы, и в качестве альтернативы перехода выбираем вариант с максимальным значением _ij. То есть, включаем в цикл то ребро, у которого альтернативный путь заведомо будет самым плохим.

Посмотрим, что получается в нашем примере. Для всех нулевых элементов нашей приведенной матрицы c_ij берем ребро с минимальным значением (по i-й строке), по которому можно выйти из вершины i в любую вершину, кроме j, и ребро с минимальным значением (по j-му столбцу), по которому можно войти в вершину j из любой вершины, кроме i.

₁₂=0+0=0; ₁₃=0+1=1; ₁₄=0+0=0; ₂₁=1+1=2; ₃₅=1+1=2; ₄₂=2+0=2; ₅₄=1+0=1.

У нас получилось 3 одинаковых значения, берем первое: ₂₁. После первого шага z={(2,1)}.

Теперь «склеим» вершины 2 и 1, чтобы вывести дугу (2,1) из рассмотрения. Для этого удалим из матрицы 2-ю строку и 1-й

	0	0	1
3		5	0
0	2		3
1	1	0	

столбец, учитывая направление дуги в цикле. Получим матрицу весов меньшей размерности. В этой матрице заменим элемент c₂₁ на символ , чтобы про него забыть, поскольку назад от 1 к 2 мы уже переходить не имеем права.

Теперь опять приводим полученную матрицу по строке и столбцу и вычисляем оценки _ij.

Итак, нижняя оценка l(z) остается прежней: l¹₀(z)=8.

Вычисляем оценки: ₁₃=0+1=1; ₁₄=0+0=0; ₃₅=3+1=4; ₄₂=2+1=3; ₅₄=1+0=1.

Имеем: max=4, то есть, выбираем c₃₅.

	0	0
0	2	
1		0

Склеиваем вершины 3 и 5, убирая из матрицы третью строку и пятый столбец, а обратный путь из 5 в 3 заменяем на .

После приведения этой матрицы нижняя оценка веса цикла остается прежней. Вычисляем оценки _ij: ₁₃=0+2=2; ₁₄=0+0=0; ₄₂=0+0=0.

Выбираем c₁₃. Склеиваем вершины 1 и 3, удаляя первую строку и третий столбец. Получаем

0	
1	0

Здесь уже без приведения ясно, что нужно брать (4,2) и (5,4).

Итак, мы получили следующий гамильтонов цикл: (3,5),(5,4),(4,2),(2,1),(1,3). Вес этого цикла равен l(z)=1+1+1+2+3=8, что удовлетворяет ранее найденной нижней оценке.

Дерево решений для рассмотренного примера выглядит так, как показано на рис.3.28.

В рассмотренном нами случае мы имели полный набор вариантов перехода, поэтому использовали только нижнюю оценку границы решений. Но если граф не полный, то, действуя по рассмотренной схеме, на каком-то шаге после приведения матрицы мы можем получить превышение нижней оценки веса цикла над L(z). В этом случае цикл, построенный для верхней оценки L(z), и будет искомым решением задачи.