Оптимизация точки старта в задаче последовательного обхода мегаполисов при наличии условий предшествования Текст научной статьи по специальности «Математика»

УДК 519.6

DOI: 10.14529/mmp 180207

ОПТИМИЗАЦИЯ ТОЧКИ СТАРТА В ЗАДАЧЕ ПОСЛЕДОВАТЕЛЬНОГО ОБХОДА МЕГАПОЛИСОВ ПРИ НАЛИЧИИ УСЛОВИЙ ПРЕДШЕСТВОВАНИЯ

А.Г. Чепцов1'2, П.А. Чепцов1'2

1Институт математики и механики им. H.H. Красовского УрО РАН,

г. Екатеринбург, Российская Федерация 2

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

Ключевые слова: маршрутная задача; ограничения; точка старта.

Введение

Рассматривается задача об организации системы перемещений с целью последовательного посещения системы конечных множеств, включая выбор начального соСТОЯНИЯ. Исследуется случай, когда система перемещений оценивается аддитивным критерием. Особенностью данного исследования является проблема оптимизации начального состояния без непосредственного построения начального решения, т.е. проблема нахождения экстремума. Значение последнего может оказаться полезным на этапе сравнения эвристик в задачах большой размерности. Разумеется, используемую при этом процедуру на основе динамического программирования (ДП) можно использовать и для нахождения всего оптимального процесса (включая выбор начального состояния) с любой степенью точности. В этом случае предъявляются более жесткие требования к памяти вычислителя, поскольку приходится сохранять все слои функции Беллмана, в то время как при определении только экстремума и его распределения на множестве начальных состояний достаточно использовать перезапись упомянутых слоев, что приводит к меньшей загрузке памяти.

Задачи маршрутизации находят широкие применения в инженерных исследованиях. Отметим транспортные задачи, задачи управления инструментом при листовой резке на машинах с ЧПУ, задачи о снижении дозовой нагрузки персонала АЭС при выполнении комплекса работ в условиях повышенной радиации, задачи авиапожарного патрулирования. Хотя прототипом упомянутых задач маршрутизации является известная труднорешаемая задача коммивояжера (ЗК) [1-5], здесь возникает целый ряд особенностей не только количественного, но и качественного характера: появляются ограничения, усложненные функции стоимости и др. П редставляется ? что разработка

методов решения упомянутых прикладных задач, связанных с маршрутизацией, требует разработки специальной математической теории. В связи с этим направлением сейчас отметим [6-91, где последовательно разрабатывался вариант широко понимаемого метода динамического программирования (ДП); в связи с разработкой аппарата ДП для решения ЗК отметим [10,11] (используемая в [6-9] схема решения на основе ДП допускает идейную аналогию с [10]). Среди других подходов к решению ЗК особо отметим метод ветвей и границ [12].

В настоящей работе конструкция [6-8], использующая ДП, переносится на постановку, в которой требуется еще оптимизировать начальное состояние. Последнее можно выбирать в пределах непустого и не обязательно конечного множества, на котором задана метрика. Мы полагаем получающееся метрическое пространство вполне ограниченным (если исходное множество содержится в конечномерном арифметическом пространстве, то достаточна ограниченность). При этом ориентируемся на реализацию глобального экстремума с высокой, но все же конечной степенью точности. Такой подход является естественным, поскольку исследуемая задача последовательного обхода мегаполисов является переборной, а потому дискретизация множества начальных состояний также представляется логичной и укладывающейся в общую схему решения подобных задач, связанных с дискретной оптимизацией. Сам же вариант упомянутой дискретизации непосредственно связывается со свойствами сечения функции, определяющей стоимость внешних перемещений. Все налагаемые условия выполнены в случае метрической задачи; задачи такого рода широко используются в дискретной оптимизации (имеются в виду постановки, в которых стоимости перемещений оцениваются расстояниями).

1. Общие понятия и обозначения

Используем обычную теоретико-множественную символику (кванторы, связки и

др.); через 0 обозначаем пустое множество, ёе£ заменяет фразу «по определении», = - равенство по определнию. Семейством называем множество, все элементы которого сами являются множествами. Если х и у - объекты, то {х; у} есть ёе£ множество, содержащее х у и не содержащее никаких других элементов. С учетом этого получаем, что для произвольного объекта г в гаде {г} = {г; г} реализуется синглетон, содержащий г : г £ {г}. Множества ЯВЛЯЮТСЯ объектами, а потому для произвольных двух

объектов и и ■и определена [13, с. 67] упорядоченная пара (УП) (и, у) = {{и}; {и; V}} с первым элементом и и вторым элементом у. Для произвольной УП г через рг1(г) и рг2 (г) обозначаем первый и второй элементы г соответственно ; объекты рг1(г) и рг2(г) при этом однозначно определены условием г = (рг1(г),рг2(г)). Если же а, Ь

и с - объекты, то (а,Ь, с) = ((а, Ь), с) есть [14, с. 17] триплет с первым элементом а, вторым элементом Ь и третьим элементом с. В этой связи напомним, что [14, с. 17]

для любых трех множеств А, В и С А х В х С = (А х В) х С (данное толкование существенно в последующих построениях).

Если Р и ф - непустые множества, то [13, с. 77] через фр обозначаем непустое множество всех отображений, действующих из Р в ф; тогда (Б1)[Р; ф], (Б1)[Р; ф] С фр, есть ёе£ множество всех биекций Р на ф. Если же Н - непустое множество, то элементы множества (Б1)[Н; Н] называем [15, с. 87] перестановками множества Н;

если при этом а Е (Б^[Н; И], то определена перестановка а Е (Б1)[Н; И], обратная к а:

а(а-1 (к)) = а-1(а(к)) = к Ук Е И.

Используя вышеупомянутое толкование триплетов^ заметим^ что для любых непустых множеств Л, Б, С и Б, отображения ф Е БЛхВхС, УП ^ Е Л х Б и точки V Е С определен элемент ф(^, V) Е Б, для которого используем также обозначение

ф(у>1,№, V), оде Ц1 = рг^г) и ^ = рг2(^).

Если ^ ^ ^^^^^^^тао, то через V^) обозначаем семейство всех подмножеств (п/м)

^^ ^^^^гаем, что V'(Z) = V(Z) \ {0} (семейство всех непустых п/м Z) и Fin(Z) -семейство всех конечных множеств из V' ^), т.е. семейство всех непустых конечных п/м Z. Если же ^ ^ ^^^^^^^^ ^^^^таое множество, то Fin(Z) = V'^). Если Р и Q

непустые множества, к Е QP и 5 Е V(Р), то к1(5) = {к(з) : 5 Е 5} Е Vесть образ множества 5 при действии к.

В дальнейшем К - вещественная прямая, К+ = {£ Е К | 0 < £}, N = {1; 2;...} и = N и {0} = {0; 1; 2;...}; пр ир Е Н0 и д Е N

Ш = {к Е N0 | (р < к)&(к < д)} Е V(N0)

(случай р, д = 0 те исключается; в час тности, 1, 0 = 0). Есл и К - непустое конечное множество (в частности, если К Е Fin(S), оде Б - непустое множество), то КI Е N ест ь def мощность (количество эле ментов) К, а (Ы)[К ] есть def множество всех биекций «промежутка:»; 1, К | на К (bi)[K] = 0. Пр и т Е N имеем равенство

(Ы)[1,т] = (Бi)[1,m; 1,т]. Полагаем, как обычно, |0| = 0.

Если Б - непустое множество, то через обозначаем множество всех неотри-

цательных вещественнозначных (в/з) функций на т.е. = (К+)й'.

2. Основная задача

Всюду в дальнейшем фиксируем непустое множество X, а также его п/м X0 Е V(X). Мы полагаем, что на X0 задана метрика [16] р Е ^+[Х0 х X0] (отметим естественный вариант: метрикой оснащено X, а Р есть ее сужение на Xа, точнее, на X0 х X0), превращающая X0 в метрическое пространство (X°,р). Пусть

Б°р(х, е) = {у Е X0 | р(х, у) <е} Ух Е X0 Уе Е К+ \ {0}

(введены открытые шары в X , радиусы которых строго положительны). Постулируем, что (X0,р) есть вполне ограниченное пространство, т.е.

Уе Е]0, ЗК Е Fin(X0) : X0 = ^ Б0р(х,е) (1)

хек

(данное условие (1) выполняется для ограниченных множеств в конечномерных арифметических пространствах). Элементы X0 рассматриваем в качестве возможных начальных состояний для процессов, связанных с рассматриваемыми ниже задачами маршрутизации.

Пусть N Е N N > 2, заданы ^непустые конечные^ множ^вства

Mi Е Fin(X), ..,MN Е Fin(X), (2)

а также связанные с ними непустые отношения

Mi Е P'(Mi х Mi),...,Mn Е P'(Mn x Mn). (3)

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

(X0 р| Mj = 0 Vj е 1Ñ)&(Mpf| Mq = 0 Vp е 1Ñ Vq е 1N \ {p}). (4)

Условия (4) типичны для рассматриваемой ниже задачи; УП из множества Mj, где j Е

1, N, определяют возможные варианты выполнения внутренних работ при посещении Mj

(ж Е X0) (x^ Е Mad) ^ x2i) Е Ma(i)) ... (xiN) Е Ma(N) - XN Е Ma(N)),

(x[i),xi1)) Е Ma(i),..., (x[N),xiN)) Е Ma(N), (5)

где а - перестановка индексов из 1, N. В (5) прямые стрелки соответствуют внешним перемещениям, а волнистые - перемещениям при выполнении внутренних работ. Вы-

а

введем множество K Е P (1,N х 1,N) (случай K = 0 не исключается и соответствует отсутствию ограничений в виде условий предшествования); УП из K именуем адресными парами. Всюду в дальнейшем предполагается, что

VKo Е P'(K) 3zo Е Ко : pr^) = pr2(z) Vz Е Ко (6)

(в [17, часть 2] указаны конкретные классы задач, удовлетворяющие (6)). Ограничение, налагаемое на выбор а, имеет следующий смысл: при z = (i,j) Е K посещение Mi должно предшествовать посещению Mj. С учетом этого в виде

A = {а Е P I a"i(pri(z)) < a_i(pr2(z)) Vz Е K},

где P = (bi)[1,N], имеем множество всех допустимых по предшествованию перестановок индексов из 1, N. При этом (см. (6); [17, (2.2.53)]) A = 0, т.е. A Е P'(P). A

та, т.е. перестановки индексов, не определяет еще конкретной реализации процесса. Следуя [17, часть 2], введем далее трассы, полагая при этом, что Vj Е 1, N

(Mj = {pri(z) : z Е Mj})&(Mj = {pr2(z) : z Е Mj}). (7)

Mj

NN

(X = (U Mj ) U X0)&(X = (U Mi^ X0), (8)

j=i i=i

получая множества из V'(X) (см. (2), (3)). Через 2 обозначаем сейчас множества всех отображений

(ъЪощ X х X.

Если х Е X0 и а Е Р, то в виде л

ад = {(ъЕ 2 I (г0 = (х,х))&& Е V Е 1, N)} Е Fin(Z) (9)

а

х. С учетом (7) - (9) получаем, что при х Е X0, а Е Р, Е [х] и т Е О, N

(рг1(гт) Е ХХ)&(рг2(гт) Е X) Введем, кроме того, X = XX и X Е V'(X); все траектории из множества (9) развиваются в X х X и соответствуют схеме (5). Пусть

Б[х] = {(а, ъ) Е А х 2 | 2 Е 2а[х]} Ух Е X0; (10)

в (10) введены множества допустимых решений (ДР) задачи последовательного обхода мегаполисов с фиксированным начальным условием (н.у.). Кроме того, пусть

Б = {(а, ъ, х) Е А х 2 х X0 | (а, ъ) Е Б[х]}

л

^введено множество ДР полной задачи). Пусть N = V'(1, N); множества из N именуем списками заданий.

Фиксируем функции с Е [X х XX х N с1 Е х X х ЭТ],... ,см Е ^[К х

X х N1, f Е R+[X]; кортеж (с,с1, ) будет использоваться для оценивания

внешних перемещений, внутренних работ и терминального состояния процесса, на содержательном уровне определяемого в (5). Если а Е Р и ъ Е 2, то полагаем, что

л N

Ca[z] = Е[ФГ2(Ф - 1)),prx(z(s)W(s,N)) + +cTS)(z(s),a1(j;N))] + f (pr2(z(N))).

Посредством (11) определяем критерии рассматриваемых ниже задач: полагаем а Е A и z Е Za[ж], где x Е X0. Так, следуя [6-8], при ж Е X0 рассматриваем задачу

Ca[z] —у min, (а, z) Е Б [ж], (12)

которой сопоставляем экстремум (значение) V [ж] Е R+ в виде наименьшего из значений Ca[z], (а, z) Е Б [ж] и множество

(SOL)[x] = {(ао, zo) Е Б [ж] | [zo] = V[ж]} Е Fin(D[x]) (13)

всех оптимальных решений задачи (12). Кроме того, имеем (полную) экстремальную задачу

Ca[z] —у inf, (а, z, ж) Е Б, (14)

для которой соответствующим значением является величина

V = inf Ca[z] = inf V[ж] Е R+; (15)

(a,z,i)£D x€X0

мнинимум в (15) может не достигаться, т.к. множество X0 не предполагалось конечным. В связи с (15) рассмотрим также задачу оптимизации н.у.:

V[ж] —► 1п£, ж € X0. (16)

Мы ставим своей целью достижение V (15) с любой наперед заданной степенью точности. В этой связи будем предполагать в дальнейшем, что функция с удовлетворяет условию равномерной непрерывности некоторых своих сечений, причем данная равномерная непрерывность будет предполагаться еще и равностепенной по набору параметров. Итак, полагаем в дальнейшем, что Уе € \ {0} € \ {0} Ух1 € X0 Чх2 € X0

(р(жьЖ2) <8) ^ (|с(хьу, 1,Ы) - с(х2,у, )| < е Уу € X). (17)

Данному предположению можно придать несколько иную форму. Для этого введем при у € X функцию

Л

c(;y, 1, N) = (c(x,y, 1, N))x£X° ER+[X0].

Тогда (17) означает, что семейство функций {с(-,у, 1, N) : у € X} является равностепенно по у € X равномерно непрерывным на (X0, р).

Пример 1. Рассмотрим следующий частный случай: множество X оснащено метрикой d € ^+[X х X], а р есть сужение ^ на X0 х X0, т.е. р(х1,х2) = d(х1,х2) Ух1 € X0 Ух2 € X0. Функция с такова, что

c(x, y, 1, N) = kd(x,y) Vx E X0 Vy E X,

где k E R+, k = 0. Тогда в силу неравенства треугольника при x' Е X0, x'' Е X0 и

y Е X ___

Ic(x',y, 1, N) - c(x'',y, 1, N)| = = kId(x',y) — d(x'',y)I < kd(x',x'') = kp(x',x'').

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

3. Динамическое программирование

Для решения задач (12), (14), (16) предполагается использовать метод ДП. В этой связи коснемся конструкций [6-8], связанных, в частности, с использованием условий предшествования в построениях, связанных с ДП. Введем прежде всего оператор вычеркивания I [17, часть 2], действующий в N то правилу: при K Е N

I(K) = K \{pr2(z): z Е ~[K]}, (19)

где S[K] = {z E K I (pri(z) E K)&(pr2(z) E K)}. Пусть, кроме того,

C = {K E N IVz E K (pri(z) E K) ^ (pr2(z) E K)} (20)

элементы (20) - существенные (для процедуры на основе ДП) списки заданий. С семейством (20) связываем подсемейства

л

Cs = {K E C j s = \Kj} Vs E 1, N.

Тогда CN = {1,N} (синглетон, содержащий 1,N) и, кроме того,

Cs_i = {K \{t} : K E Cs,t E I(K)} Vs E 2, N

(указана связь с (19)). При этом, в частности, С1 = {{Ь} : Ь € \ К1}, где К1 {рг1(г) : г € К}. Введем также

M = U Mt,Do = {(x, 0) : x EM} = M х {0}

teiN \Ki

(см. (7), (8)). Кроме того, DN = {(x, 1, N) : x E X0} = X0 х {1, N}; D0 и DN - крайние слои пространства позиций X х P (1, N).

Промежуточные слои. Если s E 1, N — 1 ж K E Cs, то последовательно конструируем

Js(K) = {j E TN \ K \{j} и K E Cs+i}, Ms[K]= и Mj,

jeJs(K)

Ds[K] = {(x,K) : x E Ms[K]} = Ms[K] х {K}.

Тогда определяем при s E 1, N — 1 елой Ds в виде объединения Ds[K], K E Cs, т.е.

Ds = U Ds[K] E P'(X х Cs). (21)

Kec

Теперь получены непустые множества П0,01, на которых последовательно

определяются в /з функции У0 € Я+[Д)], ^>1 € ^-+[^1],...,Ум € 'Я+[БМ].

При этом у0(х, 0) = f (х) Ух € Л4. Далее учитываем свойство (рг2(г),К \ {]}) € В3-1 Уз € Т^М У(х,К) € В3 У] € 1(К) У г € М,. Это свойство позволяет определить преобразование ь3-1 € /К+[03-1] ъ ь3 € ^+[03\. Итак, если з € и функция ь3-1 € /К+{03-1\ уже построена, то ь3 € ^+[03\ определяем

условиями: при (х, К) € 03

vs(x,K) = min min[c(x,pr1(z),K)+ „ ч

sV ; jei(K) zeM^ F 1W' ; (22)

+cj(z,K)+ Vs_i(pr2(z),K \{j})].

Получили рекуррентную процедуру v0 —> v1 —> ... —> vN. При этом vN E R+[DN] определяется согласно (22) следующим образом:

vN (x, 1,N )= min min[c(x, pr1(z), 1,N) +

__jei(i,N) zeM__(23)

Cj(z, 1,N) + vn-i(pr2(z), 1, N \ {j})] Vx E X0.

Из результатов [7] вытекает, что

V[x] = vn(x, 1, N) Vx E X0. (24)

Таким образом, решение задачи (16) сводится (см. (24)) к минимизации значений (23). Реально, однако, такое решение может быть найдено с использованием дискретизации X 0 на основе (1). Рассмотрим соответствующую процедуру.

Итак, фиксируем в дальнейшем произвольное (положительное) число £0 Е К+\{0} в качестве параметра точности. Тогда с учетом (17) подберем 60 Е \ {0} такое, что Ухг Е X0 4x2 Е X0

(Р(Х1,Х2) < 6о) ^ (|с(хьу, ) - с(х2,у, )| < £0 4у Е X). (25)

Теперь, используя (1), подберем К Е Пп(Х0) так, что при этом

X0 = и В0(хЛ). (26)

жеж

Мы зафиксируем (непустое конечное) множество Ж и рассмотрим зависимость

ж — vN(x, 1,N): K — R+, (27)

имея в виду (24).

Предложение 1. Справедливо неравенство

| min vn(x, 1, N) - V| < eg. (28)

жеж

Доказательство практически очевидно (см. (23) - (26)), но все же рассмотрим его схему. Если x Е X0, то для некото poro x* Е K имеем в силу (26), что x Е Bj°(x*,$0).

Тогда в силу (25) ___

\c(x,y, 1, N) - c(x*,y, 1, N)| <£о Уу Е X. (29)

С учетом (8), (19) и (23) получаем

VN(x„ 1, N) - eg < c(X, pri(z), 1, N) + Cj(z, 1, N) + vn-i(pr2(z), 1,N \ {j})

Vj E 1(1, N) yz E Mj.

Вновь используя (23) получаем, что

vn (ж*, 1, N) — eg < vn (X, 1, N). Как следствие получаем, что

min vn (ж, 1, N) <vn (X, 1,N)+ eo. (30)

жеж

Поскольку X выбиралось произвольно, установлено (см. (30))

—e0 + min vN(ж, 1, N) <vN(ж, 1,N) Vx E Xo; (31)

жеж

в свою очередь, из (15), (24) и (31) вытекает, что

—e0 + min vN(ж, 1, N) < V.

жеж

Иными словами, установлено неравенство

min vn (ж, 1, N) — V < eo. (32)

жеж

С другой стороны, по выбору K имеем из (15), что

V = inf vN(x, 1,N) < minvN(x, 1,N),

x&x° xeK

а потому с учетом (32) получаем (28). □

Итак, ограничиваясь рассмотрением н.у. x Е K, мы получаем искомый экстремум V с точностью до £0. Данное значение может использоваться в двух направлениях: 1) тестирование эвристик; 2) построение £-оптим^ьного решения при £ = £о- В обоих случаях мы выбираем из соображений минимизации функции (27) точку x0 Е K, для которой

vN (x0,1,N) = min vN (x, 1,N). (33)

В силу (28) получаем тогда неравенство 1ум(х0,1, N) — VI < £0.

1) Располагая такой точкой х0 и значением УМ (х0,1, N), мы рассматрпваем «быстрые:»; эвристики, для которых регистрируем значения (11): имеются в виду УП (а*, ъ*) Е Б[х0], определяемые тем или иным эвристическим алгоритмом, и получающиеся в ходе вычислительного эксперимента величины Са* [ъ*] Е определяемые посредством (11). Мы выбираем эвристическое решение из соображений наибольшей близости к УМ(х0,1, N), чем обеспечивается (см. (27)) и должное приближение к V по

£0

мым из набора эвристик можно выбрать £-наилучшую в смысле приближения к V при £ > £0- Заметим, ч[то для достаточно больших значений N при решении задачи (14) приходится использовать эвристики. В [18] отмечено, однако, что размерность задач исследуемого типа, допускающих точное построение экстремума, все-таки несколько выше той размерности, при которой удается, помимо экстремума, найти оптимальное решение, т.е. элемент соответствующего множества (13). Это связано с тем, что при построении функции (27) можно использовать (при реализации (22)) процедуру перезаписи слоев и, как следствие, более рационально использовать память компьютера (в памяти при этом находится только один слой значений функции Беллмана; подробнее см. в [18]). Данная процедура допускает естественную аналогию с [19].

2) Построение £-оптимш1ьного (£ = £0) решения соответствует [7] при выборе х0 Е Ж со свойством (33) в качестве н.у. Наметим данную процедуру совсем кратко

(см. [7, раздел 4]), фиксируя х0 (см. (33)). Полагаем ъ(0) = (х0,х0) и с учетом (23) выбираем п1 Е 1(1, N) И ъ(1) е м так! • ч[то

Ум (х0, ^) = с(х0, рГ1(ъ(1)), Т^) + +с^ (ъ(1), 1, N) + Ум-1(рг2(ъ(1)), 1, N \ {щ}). ^

При этом (рг2(ъ(1)), \ {гц}) Е -1,/ что позволяет применить для нахождения УМ-1(рг2(ъ(1)), \ {щ}) надлежащий вариант (22). Далее, выбираем п2 Е 1(1, N \ {п1}) и ъ(2) Е МП2 из условия

УМ-1(рг2(ъ(1)), Т^ \ {щ}) = с(рг2(ъ(1)), рг1(ъ(2)), ^ \Ы) +

+ст(ъ(2), 1, N \ {щ}) + Ум-2(рг2(ъ(2)), 1, N \ {щ; т}) ^

(см. (22)), получая также (рг2(ъ(2)), \ {п1; п2}) Е ОМ-2. Дальнейшие построения

подобны (34), (35) и доставляют в итоге ДР (п, (ъ(з))зе0м) Е 1Э[х0], где п = (Пз^еТМ) со свойством [(2(з))з€0м] = V[х0] (см. [20, раздел 7]). С учетом (13) получаем, что (П, Зесм) Е (ЯОЬ)[х0].

4. Вычислительный эксперимент

В ходе вычислительного эксперимента были произведены несколько расчетов, демонстрирующих работу алгоритма выбора начальной точки. Вычисления производились на компьютере с процессором Intel i5-2400 и 8 Гб оперативной памяти, работающем под управлением Linux Ivubuntu 16.04 (64-bit). Программа разработана на языке С • • с использованием библиотеки Qt. Рассматривался модельный вариант задачи управления инструментом при листовой резке на машинах с ЧПУ. соответствующий схеме [7, р аздел ы 5, 6] (процедура раскроя предполагается выполненной). Предполагается, что мегаполисы получены дискретизацией контуров, подлежащих резке. С каждым «городом> мегаполиса (т.е. с точкой на эквидистанте) связываются точка врезки и точка выключения инструмента; получаемый триплет именуем системой врезки. Время резки по эквидистантам не учитываем, поскольку оно одно и то жб для всех ДР. Конкретный выбор ДР влияет на время внешних перемещений, время перемещений для каждой системы врезки и время перемещения к заданной терминальной точке.

Предполагаем, что множество X раздела 2, содержащее мегаполисы, есть прямоугольник на плоскости (т.е. в R х К), определяемый точками ( — 25 мм; —25 мм), (—25 мм; 1025 мм), (1600 мм; 1025 мм) и (1600 мм; —25 мм); X0 - граница упомянутого прямоугольника. Используемые функции стоимости соответствуют [7, разделы 5, 6] и характеризуют время движения инструмента. При этом стоимость внешних перемещений, определяемая функцией с, соответствовала скорости 500 мм/с. Аналогично задавалась стоимость перемещения в терминальную точку, характеризуемая посредством f. Стоимость внутренних работ для каждой системы врезки определялась в виде суммы двух слагаемых, одно из которых определялось временем перемещений в системе врезки (от точки врезки до соответствующей ей точки на эквидистанте и от этой последней точки до точки выключения инструмента) со скоростью 10 мм/с, а второе определялось штрафом, связанным с использованием тепловых допусков. Соответствующий вариант оценивания внутренних работ приведен в [7, разделы 5, 6] (в частности, см. [7, (6.1)]). Различие в скоростях соответствует движениям в режиме

с

соответствует примеру 1 раздела 2.

Предполагалось, что N = 28, множество К адресных пар имело мощность |K| = 21. Найденный экстремум 62,44 секунды, время счета 8 минут 7 секунд. Оптимальная начальная точка (см. (33)) есть (—25; 675). Маршрут и трасса приведены на рисунке.

Финии

Старт,

Маршрут и трасса обхода множеств

Работа проводилась при финансовой поддержке комплексной программы УрО

РАН, проект, 18-1-1-9 «Оценивание динамики нелинейных управляемых систем, и

маршрутная оптимизациям.

Литература

1. Gutin, G. The Traveling Salesman Problem and Its Variations / G. Gutin, A.P. Punnen. -N.Y.: Springer, 2002.

2. Cook, W.J. In Pursuit of the Traveling Salesman. Mathematics at the Limits of Computation / W.J. Cook. - New Jersey: Princeton University Press, 2012.

3. Меламед, И.И. Задача коммивояжера. Вопросы теории / И.И. Меламед, С.И. Сергеев, И.Х. Сигал // Автоматика и телемеханика. - 1989. - Т. 50, № 9. - С. 3-34.

4. Меламед, И.И. Задача коммивояжера. Точные алгоритмы / И.И. Меламед, С.И. Сергеев, И.Х. Сигал // Автоматика и телемеханика. - 1989. - Т. 50, № 10. - С. 3-29.

5. Меламед, И.И. Задача коммивояжера. Приближенные алгоритмы / И.И. Меламед, С.И. Сергеев, И.Х Сигал // Автоматика и телемеханика. - 1989. - Т. 50, № 11. - С. 3-26.

6. Ченцов, А.Г. Задача маршрутизации с ограничениями, зависящими от списка заданий / А.Г. Ченцов, A.A. Ченцов // Доклады Академии наук. - 2015. - Т. 465, № 2. - С. 154-158.

7. Ченцов, А.Г. Маршрутизация в условиях ограничений: задача о посещении мегаполисов / А.Г. Ченцов, П.А. Ченцов // Автоматика и телемеханика. - 2016. - № 11. -С. 96-117.

8. Ченцов, А.Г. Одна параллельная процедура построения функции Беллмана в обобщенной задаче курьера с внутренними работами / А.Г. Ченцов // Автоматика и телемеханика. - 2012. - № 3. - С. 134-149.

9. Коробкин, В.В. Методы маршрутизации и их приложения в задачах повышения безопасности и эффективности эксплуатации атомных станций / В.В. Коробкин, А.И. Сесекин, О.Л. Ташлыков, А.Г. Ченцов. - М.: Новые технологии, 2012.

10. Беллман, Р. Применение динамического программирования к задаче о коммивояжере / Р. Беллман // Кибернетический сборник. - 1964. - Т. 9. - С. 219-228.

11. Хелд, М. Применение динамического программирования к задачам упорядочения / М. Хелд, P.M. Карп // Кибернетический сборник. - 1964. - Т. 9. - С. 202-218.

12. Литл, Дж. Алгоритм для решения задачи о коммивояжере / Дж. Литл, К. Мурти, Д. Суини, К. Кэрел // Экономика и математические методы. - 1965. - Т. 1, вып. 1. -С. 94-107.

13. Куратовский, К. Теория множеств / К. Куратовский, А. Мостовский. - М.: Мир. - 1970.

14. Дьедонне, Ж. Основы современного анализа / Ж. Дьедонне. - М.: Мир. - 1964.

15. Кормен, Т. Алгоритмы: Построение И Э.НЭЛИЗ / Т. Кормен, Ч. Лейзерсон, Р. Ривест. -М: МЦНМО, 2002.

16. Энгелькинг, Р. Общая топология / Р. Энгелькинг. - М.: Мир, 1986.

17. Ченцов, А.Г. Экстремальные задачи маршрутизации и распределения заданий: вопросы теории / А.Г. Чепцов. - Ижевск: НИЦ «Регулярная и хаотическая динамика:», 2008.

18. Ченцов, A.A. К вопросу о нахождении значения маршрутной задачи с ограничениями / A.A. Ченцов, А.Г. Ченцов // Проблемы управления и информатики. - 2016. - № 1. -С. 41-54.

19. Lawler, E.L. Efficient Implementation of Dynamic Programming Algorithms for Sequencing Problems / E.L. Lawler // CWI. Technical Reports. Stichting Mathematish Centrum. Mathematische Besliskunde. - 1979. - BW 106/79 - P. 1-16.

20. Ченцов, AT. К вопросу о маршрутизации комплексов работ / А.Г. Ченцов // Вестник Удмуртского университета. Математика. Механика. Компьютерные науки. - 2013. -№ 1. - С. 59-82.

Александр Георгиевич Ченцов, член-корреспондент РАН, главный научный сотрудник, отдел «Управляемые системы:», Институт математики и механики им.

Н.Н. Красовского УрО РАН; профессор, кафедра «Прикладная математика и механи-

:

[email protected].

Павел Александрович Ченцов, К&НДИДсХТ физико-математических наук, старший

«:

«

тимальный раскрой промышленных материалов и оптимальные маршрутные техно:

ция), [email protected].

Поступила в редакцию 10 апреля 2018 г.

MSC 93CXX DOI: 10.14529/mmpl80207

OPTIMIZATION OF THE START POINT IN THE GTSP

WITH THE PRECEDENCE CONDITIONS

A.G. Chentsov1'2, P.A. Chentsov1'2

1 Krasovskii Institute of Mathematics and Mechanics UrB RAS, Ekaterinburg,

Russian Federation

2

E-mail: [email protected], [email protected]

The paper is devoted to the routing problem with constraints and cost functions that can depend on the list of tasks. It is assumed that the initial condition for the process with discrete time can be selected within a metric space that satisfies the condition of complete boundedness. It is supposed that the problem includes a visiting of a finite system of megalopolises (non-empty finite sets) with the fulfillment of some works. The cost of these works each time depend on the point of arrival and the point of departure. The costs of movement and work are aggregated additively. For the problem solution widely understood dynamic programming method providing e-optimal solution for any e > 0 is used.

Keywords: route task; restrictions; start point.

References

1. Gutin G., Punnen A.P. The Traveling Salesman Problem and Its Variations, N.Y., Springer, 2002.

2. Cook W.J. In Pursuit of the Traveling Salesman. Mathematics at the Limits of Computation. New Jersey, Princeton University Press, 2012.

3. Melamed I.I., Sergeev S.I., Sigal I. The Traveling Salesman Problem. Issues in the Theory. Automation and Remote Control, 1989, vol. 50, no. 9, pp. 1147-1173.

4. Melamed I.I., Sergeev S.I., Sigal I. The Traveling Salesman Problem. Exact Methods. Automation and Remote Control, 1989, vol. 50, no. 10, pp. 1303-1324.

5. Melamed I.I., Sergeev S.I., Sigal I. The Traveling Salesman Problem. Approximate Algorithms. Automation and Remote Control, 1989, vol. 50, no. 11, pp. 1459-1479.

6. Chentsov A.G., Chentsov A.A. Route Problem with Constraints Depending on a List of Tasks. Doklady Mathematics, 2015, vol. 92, no. 3, pp. 685-688. DOI: 10.1134/S1064562415060083

7. Chentsov A.G., Chentsov P.A. Routing Under Constraints: Problem of Visit to Megalopolises. Automation and Remote Control, 2016, vol. 77, no. 11, pp. 1957-1974. DOI: 10.1134/S0005117916110060

8. Chentsov A.G. One Parallel Procedure for the Construction of the Bellman Function in the Generalized Problem of the Courier with the Inner Workings. Automation and Remote Control, 2012, vol. 3, pp. 134-149.

9. Korobkin V.V., Sesekin A.N., Tashlykov O.L., Chentsov A.G. Routing Methods and Their Applications to the Enhancement of Safety and Efficiency of Nuclear Plant Operation. Moscow, Novye tekhnologii, 2012. (in Russian)

10. Bellman R. Dynamic Programming Treatment of the Travelling Salesman Problem. Journal of the Association for Computing Machinery, 1962, vol. 9, pp. 61-63. DOI: 10.1145/321105.321111

11. Held M., Karp R.M. A Dynamic Programming Approach to Sequencing Problems. Journal of the Society for Industrial and Applied Mathematics, 1962, vol. 10, no. 1, pp. 196-210. DOI: 10.1137/0110015

12. Little J., Murty K., Sweeney D., Karel C. An Algorithm for the Traveling Salesman Problem, Operations Research, 1963, vol. 11, no. 6, pp. 972-989. DOI: 10.1287/opre.ll.6.972

13. Kuratowski K., Mostowski A. Set Theory. Amsterdam, North-Holland Publishing Company, 1967.

14. Dieudonne J. Foundations of Modern Analysis. New York, Academic Press, 1960.

15. Cormen T., Leiserson C., Rivest R. Introduction to Algorithms. MIT Press and McGraw-Hill, 1990.

16. Engelking R. General Topology. Polish Scientific Publishers, 1977.

17. Chentsov A.G. Ekstremalnye zadachi marshrutizacii i raspredeleniya zadaniy voprosy teorii [Extreme Routing and Task Distribution Problems: Theory Questions]. Izhevsk, NITs "Regulyarnaya i Khaoticheskaya Dinamika", 2008. (in Russian)

18. Chentsov A.G., Chentsov A.A. On the Problem of Obtaining the Value of Routing Problem with Constraints. Journal of Automation and Information Sciences, 2016, vol. 6, pp. 41-54.

19. Lawler E.L. Efficient Implementation of Dynamic Programming Algorithms for Sequencing Problems. CWI Technical report. Stichting Mathematisch Centrum. Mathematische Besliskunde-BW, 1979, vol. 106, no. 79, pp. 1-16.

20. Chentsov A.G. To Question of Routing of Works Complexes. Vestnik Udmurtskogo universiteta. Matematika. Mekhanika. Kompyuternye nauki, 2013, no. 1, pp. 59-82. (in Russian)

Received April 10, 2018