Задача о составлении маршрута коммивояжера. Метод ветвей и границ

Актуальность данной темы заключается в следующем, Для решения оптимизационных и других задач строительства нередко прибегают к формулировке поставленной задачи в виде каких-то хорошо известных математических задач: транспортная задача, задача поиска оптимального пути (задача коммивояжера) и другие. Сформулированную таким образом задачу решают, используя такие математические методы, как метод ветвей и границ, симплексный метод,метод Фогеля (приближенного решения), метод дефектов и другие методы.

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

Знаменитая задача коммивояжера, поставленная еще в 1934 г., является одной из самых важнейших задач в теории графов. В своей области (оптимизации дискретных задач) задача коммивояжера служит своеобразным полигоном, на котором испытываются все новые методы.

Целью данной работы будет:

1. Познакомить читателя с основными понятиями теории графов.

2. Дать представление о задаче коммивояжера.

3. Описать метод ветвей и границ.

4. Привести пример использования метода ветвей и границ для решения задачи коммивояжера.

1. Постановка задачи

Коммивояжер (бродячий торговец) должен выйти из первого города, посетить по разу в неизвестном порядке города 2,3,4тАжn и вернуться в первый город. Расстояния между всеми городами известны. В каком порядке следует обходить города, чтобы замкнутый путь коммивояжера был кратчайшим? В терминах теории графов: найти гамильтонов цикл в графе минимальной длины.

Задача коммивояжера является так называемой NP-трудной задачей, т.е. задачей, точное решение которой в общем случае может быть получено только за экспоненциальное время. Следовательно, решать ее переборным алгоритмом неэффективно при большом количестве вершин графа.

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

Впервые метод ветвей и границ был предложен Лендом и Дойгом в 1960 для решения общей задачи целочисленного линейного программирования. Интерес к этому методу и фактически его Влвторое рождениеВ» связано с работой Литтла, Мурти, Суини и Кэрела, посвященной задаче коммивояжера. Начиная с этого момента, появилось большое число работ, посвященных методу ветвей и границ и различным его модификациям. Столь большой успех объясняется тем, что авторы первыми обратили внимание на широту возможностей метода, отметили важность использования специфики задачи и сами воспользовались спецификой задачи коммивояжера. В основе метода ветвей и границ лежит идея последовательного разбиения множества допустимых решений на подмножества (стратегия Влразделяй и властвуйВ»). На каждом шаге метода элементы разбиения подвергаются проверке для выяснения, содержит данное подмножество оптимальное решение или нет. В качестве примера конкретного применения метода может быть предложена прикладная задача, связанная с проблемой размещения и обслуживания оборудования, в которой требуется определить оптимальную траекторию циклического маршрута движения робота-транспортера по траектории цеха с целью периодического обслуживания оборудования.

2. Обзор других методов решения задачи коммивояжера

Другие методы, предложенные для поиска кратчайших гамильтоновых циклов тАУ алгебраический метод, основанный на работе Йоу, Даниэльсона и ДхаванаВаВаВаВаВаВа (включает в себя построение всех простых цепей с помощью последовательного перемножения матриц), метод перебора Робертса и Флореса (метод перебора имеет дело с одной цепью, непрерывно продлеваемой вплоть до момента, когда либо становится ясно, что эта цепь не может привести к гамильтонову циклу. Тогда цепь модифицируется некоторым систематическим способом, после чего продолжается поиск гамильтонова цикла. Другой подход тАУ мультицепной метод, предложенный Селби.

3. Теория графов

3.1 Основные понятия теории графов. Определения

Граф G задается множеством точек или вершин x
₁, x₂, тАжx_n(которое обозначается через X) и множеством линий или ребер a₁, a₂, тАжa_m(которое обозначается символом А), соединяющих между собой все или часть этих точек. Таким образом, граф G полностью задается парой (X, A).

Если ребра из множества А ориентированны, что обычно показывается стрелкой, то они называются дугами, и граф с такими ребрами называют ориентированным графом. Другое, употребляемое чаще описание ориентированного графа G состоит в задании множества вершин Х и соответствия Г, которое показывает, как между собой связаны вершины. Граф в этом случае обозначается парой G=(X, Г).

Путем (или ориентированным маршрутом) ориентированного графа называется последовательность дуг, в которой конечная вершина всякой дуги, отличной от последней, является начальной вершиной следующей.

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

Иногда дугам графа G сопоставляются (приписываются) числа тАУ дуге (x_i, x_j) ставится в соответствие некоторое число c_i_j, называемое весом, или длинной. Тогда граф G называется графом со взвешенными дугами. Иногда веса приписываются вершинам x_i графа, и тогда получается граф со взвешенными вершинами. Если в графе веса приписаны и дугам, и вершинам, то он называется просто взвешенным. При рассмотрении пути , представленного последовательностью дуг, за его вес (или длину) принимается число , равное сумме весов всех дуг, входящих в , т.е. .

Гамильтонов цикл в орграфе тАУ это ориентированный цикл (контур), проходящий ровно один раз через каждую вершину графа (т.е. простая орцепь).

коммивояжер граф задача метод

Если G = (X, A) тАУ неориентированный граф с n вершинами, то остовным деревом (или, короче, остовом) графа G называется всякий остовный подграф G, являющийся деревом (т.е. графом не имеющим циклов, в котором каждая пара вершин соединена одной и только одной простой цепью). Например, если G тАУ граф, показанный на рис. 1а, то графы на рис. 1.б, в являются остовом этого графа.

3.2 Задача коммивояжера

Постановка задачи. Коммивояжер (бродячий торговец) должен выйти из первого города, посетить по разу в неизвестном порядке города 2,3,4тАжn и вернуться в первый город. Расстояния между всеми городами известны. В каком порядке следует обходить города, чтобы замкнутый путь коммивояжера был кратчайшим? В терминах теории графов: найти гамильтонов цикл в графе минимальной длины.

Решение обобщенной задачи коммивояжера

Пусть каждой дуге (i, j) графа (I, U) поставлено в соответствие число Ваназываемое длиной дуги.

Рассмотрим задачу: определить кратчайший путь из множества А в множество В, пересекающий каждое множество разбиения один и только один раз. Очевидно, что если каждое множество разбиения состоит из одного элемента и , то имеем обычную задачу коммивояжера.

Определим функцию : положим для произвольного пути . Итак, значениями функции Вабудут множества номеров подмножеств разбиения, которые пересекает путь . Пусть каждое множество Ф_i, , состоит из всевозможных подмножеств множества {1, 2,тАж, p}, a . Применим для решения этой задачи следующий алгоритм.

Достаточной системой функций в данном случае будут функции

Обозначим через Вачисло элементов произвольного конечного множества А.

Шаг 0. Положим . Пометим вершины Вапризнаками . Для помеченных вершин увеличим на 1. Рассмотрим одну из пометок и перейдем к шагу 1.

Шаг 1. Пусть ВатАУ рассматриваемая пометка. Пометим признаками Вавсе те вершины, для которых . Для вновь помеченных вершин увеличим Вана 1.

Рассмотрим следующую помеченную вершину, и будем повторять помечивания до тех пор, пока не пометим некоторую вершину Ватак, чтобы для признака пометки этой вершины выполнялось условие Ваили пока нельзя будет сделать дальнейших пометок. В первом случае перейдем к шагу 2. а во втором к шагу 3.

Шаг 2. Строим кратчайший допустимый путь от вершины . Пусть ВатАУ пометка вершины, для которой . Перейдём к вершине Ваи рассмотрим пометку Вавершины,

для которой . Далее перейдем к вершине , с пометкой , для которой . Последовательность и является кратчайшим допустимым путем.

Шаг 3. Пусть ВатАУ множество помеченных вершин. Рассмотрим все возможные числа при . Определим среди этих чисел наименьшее и возьмем его за новое приближение к длине искомого пути. Затем перейдем к шагу 1.

Этот алгоритм можно изменить. Если для некоторой пометки при всех j, для которых Ваили , то путь соответствующий этой пометке, уже продленво все смежные с вершины. Следовательно, для таких пометок признаки Ваможно стирать.

3.3 Метод ветвей и границ (МВГ)

Представим, что необходимо обойти все города страны. Так как их много, то определить кратчайший путь затруднительно. Тогда можно выбрать некоторое разбиение множества городов, например, рассматривать республики, области или районы, и определить кратчайший путь, пересекающий каждое из выбранных подмножеств разбиения только один раз. Затем уже в пределах выбранных подмножеств достроить полученный путь до требуемого. Такой подход используется в методе ветвей и границ. Этот метод позволяет опознать бесперспективные частичные решения, в результате чего от дерева поиска на одном шаге отсекается целая ветвь. Тем не менее, удовлетворительных оценок быстродействию алгоритма Литтла, основанного на этом методе, и родственных алгоритмов нет, хотя практика показывает, что на современных ЭВМ они иногда позволяют решить задачу коммивояжера для графов с количеством вершин, меньшим 100.

Описание метода ветвей и границ

Пусть х¹ тАУ центр куба Х. Вычисляем Ваи присваиваем это значение рекорду . Разбиваем куб Х ¹ⁱ со стороной ½ и вычисляем значения целевой функции в их центрах: , i= 1,тАж 2ⁿ, обновляя по ходу вычислений значение рекорда . Проверяем выполнение условия для i=1,тАж, 2ⁿ и отбрасываем соответствующие подкубы. Каждый из оставшихся разбиваем на 2ⁿ одинаковых подкубов Х²^ij со стороной ¼ и поступаем как прежде. На любом шаге у нас формируется множество К ВлкубиковВ» со сторонами 2^-^l, l>=2, целое. Правило выбора очередного кубика для разбиения называется правилом ветвления тАУ возможные варианты приводятся ниже. Кубики со стороной не больше Ваисключаются из множества К тАУ дробление кубика заканчивается. Также исключаются кубики, попавшие в множество Ва(с индексом k тАУ номером кубика) для текущего значения рекорда, тАУ правило отсечения ветвей. Рекорд обновляется при получении меньшего значения целевой функции (правило получения границ, т.е. оценок). Значения целевой функции вычисляются в центре каждого нового подкубика, включаемого в К после разбиения выбранного для этого кубика. Алгоритм останавливается, когда К пусто.

Рис. 2. Граф-дерево

Указанная терминология и название метода определяются тем, что визуально данная схема перебора представляется в виде графа-дерева, корневая вершина которого соответствует кубу Х, вершины первого яруса тАУ подкубам X^li, вершины второго яруса тАУ кубикам X²^li, подсоединенным к своим порождающим вершинам X^li-го яруса, и т.д. (см. рис. 2). Если кубик исключается из К, его вершина закрывается тАУ из нее не будут идти ветви на следующий ярус. Порядок закрытия вершины определяется правилом отсечения (своим для каждой массовой задачи), порядок раскрытия тАУ правилом ветвления (своим для каждой индивидуальной задачи). Различают два вида правил ветвления по типу построения дерева решений (выбора вершин для раскрытия): Влв ширинуВ», когда сначала раскрываются все вершины одного яруса до перехода к следующему, и в ВлглубинуВ» тАУ всякий раз раскрывается лишь одна (обычно с лучшим значением рекорда) вершина на ярусе до конца ветви. На практике реализуют некоторую смесь, например, первое правило пока хватает машинной памяти (в К не слишком много элементов), затем переключаются на второе. Предпочтительность той или иной стратегии ветвления оценивается каждым вычислителем по-своему, исходя из главной задачи метода ветвей и границ тАУ быстрее получить лучший рекорд, чтобы отсечь больше ветвей. Удачный выбор стратегии ветвления в МВГ (например, на основе имеющейся у вычислителя дополнительной информации или эвристических соображений об объекте) позволяет (хотя и не гарантированно) решать задачи большой размерности.

Отметим, что в худшем случае ВатАУ не удается отбросить ни одной точки х тАУ и приходим к полному перебору; т.е. указанная экспоненциальная оценка точна на классе всех липшицевых функций.

4. Выбор объекта управления. анализ аспектов его работы

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

Для рассмотрения можно представить себе, что робот развозит по цеху некоторый материал, требуемый для работы агрегатов. Предположим, что уже имеется некоторый образец робота, т.е. нет необходимости покупать новый тАУ капитальные затраты отсутствуют. Пусть робот имеет электрический привод и в качестве источника питания тАУ переносные аккумуляторы.

5. Постановка и решение оптимизационной задачи

5.1 Выбор критерия оптимальности

Основным будем считать экономический критерий, т.е. будем стремиться снизить все статьи денежных расходов, связанных с работой данного робота. Отметим, что схема работы робота должна соответствовать технологической схеме работы цеха, т.е. поставка требуемого материала должна осуществляться в поставленные сроки, определяемые видом работ цеха. Условно будем считать это требование удовлетворенным.

Перечислим основные экономические показатели, связанные с работой робота:

1. Прямые затраты

1.1 Единовременные тАУ по доставке машин.

1.2 Амортизационные суммы.

2. Накладные расходы. Заработная плата рабочих.

3. Условно-постоянные расходы.

Содержание служебных помещений.

Пожарной и сторожевой охраны.

Бытовое благоустройство цеха и подсобных помещений.

Расходы, связанные с техникой безопасности.

4. Косвенные расходы.

Расходы на содержание обслуживающего персонала.

Ремонтных мастерских. Расходы на ремонт.

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

Таким образом, следует по возможности сократить путь, проходимый роботом.

Итак, задачу уменьшения денежных затрат мы свели к задаче поиска пути минимальной длинны. Имеем задачу коммивояжера.

5.2 Выявление основных особенностей рассматриваемого объекта

Будем считать, что у нас имеются собранные статистические данные, показывающие время движения робота между агрегатами цеха (См. табл. 1). Здесь ВатАУ номера агрегатов. ВатАУ соответствует времени движения выраженном в некоторых условных единицах. Таблица симметрична. Незаполненные поля говорят о невозможности данного маршрута по каким-то причинам.

Таблица 1.

	1	2	3	4	5
1	*	4	2		5
2		*		1	9
3			*	3	4
4				*	11
5					*

Вместе с этим смотрят:

РЖнварiантнi пiдпростори. Власнi вектори i власнi значення лiнiйного оператора

Актуальные проблемы квантовой механики

Алгебра и алгебраические системы

Волоконно-оптические датчики температуры на основе решеток показателя преломления

Время и пространство - идеалистические понятия