CC BY
12
Известия Института математики и информатики Удмуртского государственного университета
2021. Том 58. С. 94-126
УДК 517.958, 530.145.6
© А. Г. Ченцов, А. А. Ченцов, А. Н. Сесекин
ОДНА ЗАДАЧА МАРШРУТИЗАЦИИ РАБОТ В УСЛОВИЯХ ПОВЫШЕННОЙ РАДИАЦИИ
Исследуется задача последовательного обхода мегаполисов, ориентированная на проблему демонтажа системы радиационно опасных объектов при ограничениях в виде условий предшествования. Радиационное воздействие на исполнителей оценивается дозами, получаемыми при перемещениях и при выполнении работ по демонтажу. Рассматривается маршрутная задача минимизации дозовой нагрузки работников, осуществляющих демонтаж в той или иной последовательности операций. Исследуется процедура построения оптимального решения с использованием варианта динамического программирования. На этой основе построен алгоритм, реализованный на ПЭВМ. Приведены примеры численного решения модельной задачи на минимум дозовой нагрузки.
Ключевые слова: динамическое программирование, демонтаж, маршрут. 001: 10.35634/2226-3594-2021-58-06
Введение
Задачи, связанные с маршрутизацией перемещений возникают в самых различных сферах человеческой деятельности. Можно отметить задачу управления инструментом при фигурной листовой резке на машинах с ЧПУ (см. [1-3]), задачу авиапожарного патрулирования лесов, транспортные задачи. В настоящей работе мы ориентируемся на применения, связанные с атомной энергетикой; имеется в виду задача минимизации дозовой нагрузки при демонтаже системы радиационно опасных элементов. Эта задача может возникать в связи с авариями на АЭС, подобных Чернобылю и Фукусиме, а также в связи с необходимостью демонтажа энергоблока, выведенного из эксплуатации. В своих построениях мы ставим своей целью оптимизацию процессов, включающих выбор точки старта, варианта очередности выполнения заданий по демонтажу и конкретной траектории. Возникает очень сложная задача маршрутизации, в которой даже насчитывание функций стоимости (здесь — совокупных доз радиации) составляет серьезную проблему; см. в этой связи [4,5]. В этой связи в статье уделяется данной проблеме серьезное внимание; предполагается, что построение функций затрат осуществляется в значительной степени по мере их использования для целей решения основной задачи. Используемые методы восходят к [6]. Монография [4] содержит необходимую детализацию в части применения в атомной энергетике; имеется в виду задача минимизации дозовой нагрузки исполнителей при выполнении комплекса работ по демонтажу радиационно опасных элементов.
Отметим, что рассматриваемая задача маршрутизации, ориентированная на инженерные приложения, имеет своим прототипом известную труднорешаемую задачу коммивояжера (см. [7-9] и др.), однако содержит целый ряд существенных особенностей. Среди методов решения задачи коммивояжера (ТБР в англоязычной литературе) отметим сейчас метод ветвей и границ [10] и динамическое программирование (ДП) [11,12]. Отметим также обстоятельный обзор [13-15]. Существенный момент настоящего исследования связан с условиями предшествования, когда для некоторых пар объектов посещения требуется, чтобы один из этих объектов посещался раньше другого. Данное ограничение удается использовать для снижения вычислительной сложности (см. [6, § 4.9]). В исследуемой задаче о минимизации дозовой нагрузки появление условий предшествования может быть связано,
например, с тем, что один из демонтируемых объектов размещается на другом и поэтому верхний объект должен демонтироваться раньше нижнего (имеются и другие условия содержательного характера, приводящие к требованиям предшествования).
В настоящей работе при исследовании задачи демонтажа, включая часть, связанную с конкретными вычислениями, используется достаточно подробная детализация функций стоимости, значениями которых являются дозы радиации. Как уже отмечалось, само построение функций стоимости является достаточно сложным. В статье предполагается, что критерий отвечает аддитивному агрегированию затрат (в связи с «неаддитивными» задачами такого типа см. [3]), при котором естественным образом выделяются компоненты, оценивающие внешние перемещения, и работы, проводимые в пределах объектов посещения (здесь мегаполисов, то есть непустых конечных множеств) и именуемые внутренними, а также компонента, оценивающая терминальное состояние. Дозы, получаемые при перемещениях (здесь — прямолинейных), находятся всякий раз интегрированием зависимости, обратно пропорциональной квадрату расстояния до источника. Дозы, получаемые при выполнении внутренних работ получаются суммированием динамической компоненты, оценивающей перемещения (при этом раздельно оценивается перемещение к источнику с целью демонтажа и перемещение после демонтажа), и статической, оценивающей дозу, получаемую вблизи источника при непосредственном проведении демонтажа. Терминальная компонента может учитывать остаточный фон, но, возможно, и воздействие оставшихся после выполнения всех заданий источников, осуществляющееся на этапе перемещения работника к пункту эвакуации (в примерах, однако, терминальная компонента предполагалась тождественно равной нулю).
Представляется, что решение упомянутой задачи в ее оптимизационной постановке может быть полезным на этапе планирования работ по демонтажу. В частности, это касается вопросов замены исполнителей в связи с получением достаточно высоких доз радиации, когда «одним составом» удается демонтировать лишь часть радиационно опасных объектов. Эта проблема, однако, требует дополнительного исследования.
§ 1. Общие понятия
В дальнейшем используется стандартная теоретико-множественная символика (кванторы, связки и т. п.); 0 — пустое множество, = — равенство по определению. Семейством называем множество, все элементы которого сами являются множествами. Если х и у — объекты, то {х; у} есть их неупорядоченная пара, то есть множество, содержащее х, у и не
содержащее никаких других элементов. Для всякого объекта г в виде {г} = {г; г} имеем
синглетон, содержащий г: г € {г}. Если и и V — объекты, то [16, с. 67] (и, V) = {{и}; {и; V}} есть упорядоченная пара (УП) с первым элементом и и вторым элементом V. Если к есть какая-либо УП, то через рг1(к) и рг2(к) обозначаем соответственно первый и второй элементы к; к = (рг1(к),рг2(к)). Каждым трем объектам а, Ь и с сопоставляем их упорядоченный триплет (а,Ь,с) = ((а, Ь),с). Соответственно, полагаем А х В х С = (А х В) х С для трех любых множеств А, В и С .В связи с последними определениями см. [16, гл. II] и [17, с. 17]. Если Н — множество, то через Р(Н) обозначаем семейство всех п/м Н и полагаем, что Р'(Н) = Р(Н) \ {0} (семейство всех непустых п/м Н); через Пп(Н) обозначаем семейство всех непустых множеств из Р'(Н), то есть семейство всех непустых конечных п/м Н. Если А и В — непустые множества, то через ВА обозначаем (см. [16, гл. II, § 6]) множество всех функций, действующих из А в В; если £ € ВА (а это означает, что £: А ^ В)
и С € Р(А), то £1 (С) = {£(х): х € С} есть образ множества С при действии £. В случае, когда С и Н — непустые множества,
(Б1)[С,Н] = £ € Нс | (/ 1(С) = Н)&(У01 € С У£2 € С £Ы = £(52)) ^ (£1 = 52))}
есть множество всех биекций множества С на Н; см. [18, с. 86]. Тогда для всякого непустого множества Б в виде (В1)[Б; имеем [18, с. 87] множество всех перестановок Б. Как
обычно, К+ = {£ е К | 0 ^ £} = [0, то [ (К — вещественная прямая), N = {1;2;...}
и N0 = {0} и N = {0; 1; 2;...}; при р е N0 и д е N0
рТд = {к е N0 | (р ^ к)&(к ^ д)}
(ясно, что р, д = 0 при д < р). Для всякого непустого множества Б через "&+[£] обозначаем множество всех неотрицательных вещественнозначных (в/з) функций на Б, то есть
^+[Б] = (К+)5. Непустому конечному множеству К сопоставляем его мощность |К| е N и непустое множество (Ы)[К] = (В1)[1, |К|; К]; как обычно, полагаем, что |0| = 0. Ясно, что (Ы)[1,т] = (В1) [1, т; 1, т] при т е N есть множество всех перестановок 1,т; если а е (Ы)[1,т], то а-1 е (Ы)[1,т] есть по определению перестановка, обратная к а:
а(а-1(к)) = а-1(а(к)) = к Ук е 1,т.
§ 2. Постановка задачи
Мы введем общую постановку и кратко обсудим некоторые детализации, связанные с задачей о демонтаже радиационно опасных элементов. Фиксируем непустое множество X и (непустое конечное) его п/м X0 е Пп(Х). Элементы X0 играют роль стартовых точек; в пределах X могут осуществляться рассматриваемые ниже перемещения. Пусть N е N таково, что 2 ^ N, заданы множества
м1 е Пп^),...,е Пп^), (2.1)
то есть мегаполисы, подлежащие посещению, а также
М1 е Р'(М1 х М1),..., е Р'(Мм х ) (2.2)
((2.2) — система непустых отношений: 0 = М, С М, х М, при ] е 1, N). Полагаем, что
(X0 п М, = 0 У? е ^)&(Мр п Мя = 0 Ур е Уд е Т^\ {р}). (2.3)
Условия (2.3) типичны для задач маршрутизации. Отношения (2.2) будут определять условия выполнения работ, связанных с посещением мегаполисов (см. (2.1)) и именуемых внутренними. При ] е 1, N полагаем, что
(М, = {рг1(г): г е М,})&(М, = {рг2(г): г е М,}); (2.4)
ясно, что М, е Р'(М,) и М, е Р'(М,). Кроме того, имеем, что
N N
(х = X0 и (ие Пп^))& (X = X0 и (и М^ е Пп^)). (2.5)
г=1 г=1
Пусть Р = (Ы)[1,Ж|. Рассматриваемые ниже системы перемещений содержательно характеризуются схемой
(ж° е X0) ^ (^1,1 е М а(1) ^ Х1,2 е Ма(1)) ^ ... ^ (^N,1 е ^ ^N,2 е Ma(N)), (2.6)
где а е Р, (ж1;1,ж1;2) е Ма(1), ..., ) е Ма(^. Выбор а е Р может быть
стеснен условиями предшествования, для введения которых зафиксируем множество К е е Р (1,N х 1,N), элементы последнего называем адресными парами (напомним, что К С С 1, N х 1, N). Полагаем далее, что
УК0 е Р'(К) 3^0 е К0: рг^) = рг2(г) Уг е К0. (2.7)
Данное условие типично для практических задач (см. обсуждение в [6, ч. 2]). Называя перестановки из P маршрутами, полагаем, что маршрут допустим, если при всяком выборе адресной пары z £ K мегаполис с номером pr1(z) £ 1, N посещается раньше чем мегаполис с номером pr2(z) £ 1, N. Тогда (см. [6, ч. 2])
A = {а £ P | Vti £ VN Vi2 £ VN ((a(ti), a(t2)) £ K) ^ (ti < ¿2)} (2.8)
= {а £ P | a-1(pr1(z)) < a-1(pr2(z)) Vz £ K} £ P'(P)
есть множество всех маршрутов, допустимых по предшествованию. Возвращаясь к (2.6), введем траектории, согласованные с маршрутами. В силу условий, отмеченных после (2.6), уместно рассматривать перемещения в X х X. Итак, траекториями называем перемещения в пространстве УП. В силу (2.5), (2.6) в качестве фазового пространства процесса логично теперь рассматривать X х X.
Отметим, что траектории, используемые ниже, являются кортежами, а точнее, функциями, определенными на конечных п/м No. В этой связи введем множество Z всех отображений из 0, N в X х X. Тогда при x £ X0 и а £ P
= {z £ Z | (z(0) = (x,x))&(z(t) £ Ma(t) Vt £ VN)} £ Fin(Z) (2.9)
есть пучок траекторий со стартом в x, согласованных с маршрутом а (следуя традиции задачи коммивояжера, называем перестановку индексов маршрутом). Тогда при x £ X0
D[x] = {(а, z) £ A х Z | z £ Z«[x]} £ Fin(A х Z) (2.10)
есть множество всех допустимых решений (ДР) задачи со стартом в x. Наконец,
D = {(а, z,x) £ A х Z х X0 | (а, z) £ D[x]} £ Fin(A х Z х X0). (2.11)
Триплеты, являющиеся элементами D (2.11), будем называть (допустимыми) маршрутными процессами.
Функции стоимости. Пусть N = P'(1, N). Множества — элементы семейства N — называем далее списками (заданий). Используем (2.5) и фиксируем следующие N + 2 функции:
c £R+[X х X х N], С1 £ R+M-1 х N], ..., cn £ R+[Mn х N], f £ R+[M],(2.12)
где M есть по определению объединение всех множеств Mj, i £ 1,N. Используем c для оценивания внешних перемещений, с1,..., cn — для оценивания внутренних работ, связанных с посещением мегаполисов, а f — для оценивания терминального состояния (точки xN 2 в (2.6)). Если x £ X0, а £ P и z £ Za[x], то полагаем (см. (2.4), (2.5), (2.9)), что
N
Ca[z] = [c(pr2(z(t - 1)), pr1(z(t)), + Са^^а1^))] (2.13)
t=1
+f (pr2(z(N))).
Полагаем, что (2.13) определяет аддитивный критерий. Тогда при x £ X0 в виде
Ca[z] ^ min, (а,z) £ D[x], (2.14)
получаем x-задачу, характеризуемую экстремумом
V[x]= min Ca [z] £ R+, (2.15)
(a,z)eD[x]
а также (непустым конечным) множеством
(SOL)[x] = {(a0, z0) G D[x] | C«o[z0] = V[x]} G P'(dD[x]) (2.16) оптимальных решений. Соответственно, полная задача
Ca[z] ^ min, (a, z,x) G D, (2.17) характеризуется глобальным экстремумом
V = min Ca[z] G R+, (2.18)
а также множеством
SOL = {(ao, zo,Xo) G D | [zo] = V} G Fin(D) (2.19)
всех оптимальных маршрутных процессов. Отметим, наконец, задачу
V[x] ^ min, x G X0, (2.20)
оптимизации точки старта, для которой V (2.18) является экстремумом, а
Xo0pt = {x0 G X0 | V[x0] = V} G P'(X0) (2.21)
есть множество всех оптимальных точек старта. Нашей основной целью является решение задачи (2.17); задачи (2.14) и (2.20) являются по смыслу вспомогательными, но их роль в последующих построениях весьма существенна.
§ 3. Динамическое программирование, 1
В этом параграфе мы конструируем систему частичных задач, которую рассматриваем как расширение основной задачи маршрутизации. Для построения этого расширения напомним прежде всего определение специального оператора вычеркивания (заданий из списка), связанного с реализацией условий предшествования; см. [2-6] и др. Итак, следуя [6, ч. 2] определяем оператор I, действующий в N посредством правила: если K G N, то
I(K) = K \{pr2(z): z G S[K]}, (3.1)
где S[K] = {z G K | (pr1(z) G K)&(pr2(z) G K)}. Итак, I сопоставляет каждому списку заданий некоторый подсписок. Отметим, что (см. [19, замечание 3.2]) при t G 1,N одноэлементный список {t} является неподвижной точкой I, то есть I({t}) = {t}; подчеркнем, что данное свойство следует из того, что (при условии (2.7)) pr1(z) = pr2(z) Vz G K (см. [19, (3.13)]).
С (3.1) связана схема построения частичных маршрутов со свойством допустимости по вычеркиванию. Итак, при K G N мы имеем значение |K| G 1,N и получаем (см. [19, (3.16)]) множество
(I - bi)[K] = {a G (bi)[K] | a(m) G I(a1(^7^)) Vm G lHK]} G P'((bi)[K]). (3.2) Напомним важное свойство [19, (3.17)]: справедливо равенство
A = (I - bi)[l7Ñ]. (3.3)
В (3.2) определены частичные маршруты, допустимые по вычеркиванию. Далее, рассмотрим траектории, согласованные с этими маршрутами. Пусть X = X U X.
Если K £ N, то через ZK обозначим множество всех отображений из 0, |K | в X х X; итак, ZK есть множество всех кортежей (zi)ie0"fK|, где Zj £ X х X Vj £ 0, |K|. Если x £ X, K £ N и а £ (bi)[K], то полагаем, что
Z(x, K, а) = {z £ Zk | (z(0) = (x, x))&(z(t) £ MaW Vt £ 1, |K|)} £ Fin(ZK); (3.4)
отметим, что при z £ Z(x,K, а) и t £ 1, |K| непременно z(t) £ X х X, z(0) £ X х X и, в частности, z(0) £ XX х X. В связи с (3.3) представляется важным случай K = 1,N. Кроме того, при x £ X0 и а £ (bi) [1, N] имеем, что Z(x, 1, N, а) С Z. Иными словами, как видно из (3.3) и (3.4), при x £ X0 и а £ A определено (см. (2.5)) множество
Z(x, 17Жа) = {z £ Ъщ | (z(0) = (x,x))&(z(t) £ M«(t) Vt £ 1,N)} (3.5)
= {z £ Z | (z(0) = (x,x))&(z(t) £ Ma(t) Vt £ 1,N)}
(мы учитываем, что |1, N| = N и Z С Z^n); тогда, как видно из (2.9) и (3.5),
Za[x] = Z (x, 17Жа). (3.6)
Полезно учитывать (3.3) и (3.6) в их естественной совокупности.
Введем в рассмотрение варианты аддитивного критерия частичных задач. Итак, при x £ X, K £ N, а £ (bi)[K] и z £ Z(x, K, а) полагаем, что
|K|
Ca [z | K] = J>(pr2(z(t - 1)), pr1(z(t)), (3.7)
t=1
а^^Л^})) + с^та1^))] + f (pr2(z(|K |))).
Заметим, что (3.7) определено, в частности, при x £ X0, K = 1, N, а £ P и z £ Za[x] (здесь мы учитываем, конечно, что P = (bi) [1, N], и (3.6)); получаем тогда, что (см.(2.13))
N
Ca [z | iTN] = J>(pr2(z(t - 1)), pr1(z(t)), (3.8)
t=1
а1^)) + Ca(t)(z(t), а1^))] + f (pr2(z(N))) = Ca[z].
Введем теперь в рассмотрение частичные (укороченные) задачи, учитывая (3.2), (3.5) и (3.7). Итак, при x £ X и K £ N рассматриваем задачу
Ca[z | K] ^ min, а £ (I - bi)[K], z £ Z(x,K^). (3.9)
В связи с (3.9) удобно ввести в рассмотрение ДР частичных задач, полагая при x £ X и K £ N
Dk[x] = {(а, z) £ (I - bi)[K] х Zk | z £ Z(x,K,a)} £ Fin((I - bi)[K] х Zk). (3.10) Тогда при x £ X и K £ N задача (3.9) может быть переписана в виде
Ca[z | K] ^ min, (а, z) £ Dк[x]; (3.11)
сопоставляем данной задаче ее экстремум
v(x,K) = min Ca[z | K] = min min Ca[z | K] £ R+, (3.12)
(a,z )sDk [x] ae(I-bi)[K] z eZ(x,K,a)
а также (непустое) множество оптимальных решений
(sol)[x | K] = {(ао, zo) G Dк[x] | [zq | K] = v(x,K)} G P'(Dк[x]). (3.13)
Задачу (3.11) можно рассматривать при ж С X0 и К = ; тогда согласно (3.3) и (3.6) получаем после простых преобразований
[ж] = {(а, г) с А х Z | 2 С ЗД} = (3.14)
Далее, учтем (3.8). Тогда при ж € X0 и К = 1, N имеем, что (см.(2.15), (3.8), (3.14))
v(x, 1, N)= min Ca[z] = V[x] (3.15)
(a,z)€-D[x]
и, кроме того,
(sol)[x | 1, N] = {(ао, zq) G Din[x] | C«o[*o | 1, N] = v(x, 1,N)} (3.16)
= {(ао, zq) G D[x] | C«0[zq] = V[x]} = (SOL)[x].
Свойства (3.15), (3.16) показывают, что система задач (3.11) действительно может рассматриваться как расширение системы задач (2.14); иными словами, при x G X0 задачи (2.14) и (3.11), где K = 1,N, совпадают: см. (3.8), (3.14). Введем, наконец, значения v(x, 0), x G M. А именно, полагаем, что
v(x, 0) = f (x) Vx G M. (3.17)
Теперь определена функция (см. (3.12), (3.17))
v: (X х N) U (M x {0}) ^ R+ (3.18)
(заметим, что (X x N) П (M x {0}) = 0)). Если x G X, K G N, j G I(K) и z G M^, то
((pr2(z), K \ {j}) G X x N) V ((pr2(z), K \ {j}) G M x {0}). (3.19)
Замечание 1. В целях полноты изложения проверим (3.19), фиксируя x G X, K G N,
j G I(K) и z G Mj. Пусть n = |KТогда (n = 1) V (n G 2, N). Рассмотрим два этих случая отдельно.
1) Пусть п = 1. Тогда К = {т}, где т € и, как уже отмечалось, 1(К) = 1({т}) = {т}. Поэтому ^ = т и К \ {_?'} = {т} \ {т} = 0. При этом
рг2(г) € Мт,
где Мт С М. Тогда рг2(г) € М. Кроме того К \ {_?'} = {т} \ {т} = 0. Поэтому в рассматриваемом сейчас случае
(рг2(г),К \{;}) € М х {0}. Получили следующую импликацию
(п = 1) ^ ((рг2(г), К \ {;}) € М х {0}). (3.20)
2) Пусть п € 2, N. Тогда, поскольку 1(К) С К, имеем по выбору ], что |К \ {^}| = 0, а потому К \ {_?' } € N. Далее, поскольку г € М^-, то в силу (2.4) рг2 (г) € М^- и, в частности, согласно (2.5) рг2(г) € X. В итоге
(рг2(г),К \{;}) € X х N
в рассматриваемом случае 2). Итак, имеем импликацию
(n £ 2, N) ^ ((pr2(z), K \ {j}) £ X х N).
C учетом (3.20) получаем теперь требуемое свойство (3.19). □
Отметим, что согласно (3.18) при x £ X, K £ N, j £ I(K) и z £ Mj определено значение
v(pr2(z),K \{j}) £ R+.
Теорема3.1. Если x £ X и K £ N, то справедливо равенство
v(x,K)= min min[c(x,pr1(z),K) + Cj(z,K)+ v(pr2(z),K \ {j})]. (3.21)
jei(K) zeMj
Доказательство. Фиксируем x £ X и K £ N. Пусть
(n = |K|)&(ш = min mMi[c(x,pr1(z),K) + Cj(z,K)+ v(pVz),K \ {j})]). (3.22)
j€i(K) zeMj
Тогда n £ 1, N и ш £ R+. Получаем, что
(n = 1) V (n £ 2, N). (3.23)
Оба случая в (3.23) рассмотрим отдельно.
1) Пусть n =1. Здесь рассуждение соответствует [19, c. 66]. Тем не менее напомним его в краткой форме. Имеем с очевидностью, что K = {r}, где r £ 1, N. Далее, как уже отмечалось,
I(K) = I({r}) = {r} = K.
Из (3.22) следует поэтому, что справедливо равенство ш = min[c(x, pr1(z), {r})+cr(z, {r}) +
+ v(pr2(z), 0)], где (pr2(z), 0) £ Mr х {0} при z £ Mr и Mr С M. Тогда в силу (3.17) имеем, что
ш = min[c(x pr1(z^ {r}) + сr(z, {r}) + f (pr2(z))]. (3.24)
Рассмотрим представление v(x,K) = v(x, {r}). Прежде всего заметим, что (bi)[K] = = (bi)[{r}] есть одноэлементное множество, содержащее биекцию a: 1,1 ^ {r}, для которой a(1) = r. При этом (I - bi)[K] £ P'((bi)[K]), где (bi)[K] = {a}. Тогда (I - bi)[K] = 0 и (I - bi)[K] С {a}, а потому a £ (I - bi)[K] и, более того, (I - bi)[K] = {a}. Поэтому в силу (3.12) имеем равенство
v(x,K)= min Ca[z | K]. (3.25)
z€Z (x,K,a)
Согласно (3.4) и определению a имеем, что
Z(x, K, a) = {z £ Zk | (z(0) = (x, x))&(z(1) £ Ma(1))}. (3.26)
В то же время согласно (3.7) и (3.26) при z £ Z(x, K, a)
Ca[z | K] = c(x, pr1(z(1)), a1({1}) + Ca(1)(z(1), a1({1})) + f (pr2(z(1))) (3.27)
= c(x, pr1(z(1)), {r}) + Cr(z(1), {r}) + f (pr2(z(1))).
Из (3.26) и (3.27) получаем, что (см. (3.25))
v(x, K) = min [c(x, pr1 (z(1)), {r}) + ca(1)(z(1), {r}) + f (pVz(1)))] (3.28)
z€Z (x,K,a)
= min [c(x, prl(z), {r}) + Cr(z {r}) + f (pr2 (z))],
z€Z(1)(x,K,a)
где множество Z(1)(x,K, a) = {z(1): z £ Z(x,K, a)}. С учетом (3.26) получаем, однако, что Z (1)( x, K, a) = Ma(1) (легкопроверяемое равенство). Поэтому (см. (3.28))
v(x,K)= min [c(x, pr1(z), {r}) + Cr(z, {r}) + f (pr2(z))].
zeMa(i)
С учетом (3.24) получаем, что справедливо равенство ^(ж, К) = ш в случае 1). Итак, истинна импликация
(п = 1) ^ (ф,К)= ш). (3.29)
2) Рассмотрим теперь случай
n £ 2, N. (3.30)
Тогда получаем (см. (3.30)) очевидное свойство
n - 1 £ 1, N - 1. (3.31)
При этом определены множества
0, n - 1 = {k £ N0 | k ^ n - 1} и 1, n - 1 = {k £ N | k ^ n - 1},
причем
0, п - 1 = {0}и 1, п - 1. (3.32)
Поскольку |К| = п ^ 2, то (см. (3.1)) имеем с очевидностью, что
|К }| = п - 1 V? е 1(К). (3.33)
C учетом (3.13) выберем (а0,2°) е (во1)[ж | К], получая при этом, что (а0, 2°) е [ж] и, кроме того,
Са0 [2° | К]= ^(ж, К). (3.34)
Из (3.10) вытекает, что справедливы свойства
а° е (I - Ы)[К]: 2° е 2(ж, К, а°). (3.35)
Тогда, в частности (см. (3.2)), а° е (Ь1)[К] и с учетом (3.22)
а°: Т7п ^ К, причем а° биективно. Кроме того, в силу (3.2)
а°(т) е 1((а°)1(тТп)) Ут е ^п. (3.36)
Далее, из (3.4) и (3.35) имеем с очевидностью,что е , а потому (см. (3.22))
2°: 0Тп ^ XX х X;
при этом согласно (3.4) и (3.35) реализуются свойства
(zq(0) = (x,x))&(z°(t) G Ma0(t) Vt G Нуга). (3.37)
По выбору а0 и z0 имеем, что (а0)1 (1,n) = K, а потому (см. (3.36)) а0(1) G I(K). Кроме того, из (3.37) вытекает, что z0(1) G Ma0(1). Из (3.22) получаем, как следствие, что
w ^ c(x,pri(z0(1)),K) + c«0(i)(z°(1),K)+ v(pr2(z0(1)),К), (3.38)
где К = K \ {а°(1)}. Заметим, что согласно (3.33) имеем равенство
|К| = n — 1. (3.39)
Введем в рассмотрение следующее отображение:
а° = (а°(г + 1))гещ—1 (3.40)
а°: 1, n — 1 ^ K и при этом а°(j) = а°(? + 1) Vj G 1, n — 1. Тогда по выбору а° (см. (3.35)) имеем из (3.40) свойство
а° G (I — bi)[K]. (3.41)
Данное свойство установлено в [20]. Итак, мы получили новый частичный маршрут. Заметим теперь, что pr2(z°(1)) G X, а потому определено (см. (3.4), (3.39), (3.41)) множество
Z(pr2(z°(1)), K, а°) = {z G Zk | (z(0) = (pr2(z°(1)), pr2(z°(1)))) (3.42)
&(z(t) G M«0(t) Vt G 1, n — 1)}.
Отметим, что z°(t + 1) G X x X при t G 0,n — 1. Из (3.37) вытекает, что z°(t + 1) G G MQ0(t+1) Vt G 0, n — 1. С учетом (3.40) получаем теперь, что
z°(t + 1) G M«0(t) Vt G 1,n — 1. (3.43)
Введем в рассмотрение картеж z* : 0,n — 1 ^ X x X по следующему правилу
(z*(0) = (pr2(z°(1)),pr2(z°(1))))&(z^(t) = (z°(t + 1) Vt G 1) (3.44)
(учитываем, что X С XX). Тогда согласно (3.43), (3.44) z*(t) G Ma0(t) Vt G 1, n — 1. Заметим, что (см. (3.39), (3.42), (3.44))
z* GZ(pr2(z°(1)), K, а°), (3.45)
а тогда (см. (3.41), (3.45)) имеем из (3.10), что (а°, z*) G ¿k[pr2(z0(1))]. В силу (3.12)
v(pr2(z°(1)),K) ^ ¿«0[z* | К]. (3.46)
В связи с (3.46) напомним, что (pr2(z°(1)),К) G X x N. Напомним (3.7). Тогда
n— 1
¿«0[z* | К] = J>(pr2(z*(t — 1)),pri(z*(t)), (а°)1(^,¥—Г)) (3.47)
t=i
+c«0(t)(z*(t), (ао)1^—1))] + f (pr2(z*(n — 1))). 103
Из (3.38), (3.46) вытекает, что справедливо неравенство
w ^ c(x, рГ1 (z0(1)),K) + сао(1)(z0(1), K) + Cao[z" | K]. (3.48)
Рассмотрим (3.47), учитывая (3.31) и (3.44). Прежде всего отметим, что f (pr2(z"(n — 1))) = = f (pr2(z°(n)). Кроме того, из (3.44) вытекает, что
cao(t)(z"(t), (ao)1(t;n—Г)) = cao(t+i)(z°(í + 1), (ao)1(t7¥—Г)) Vt G Г. (3.49)
Пусть т G 1, n — 1. Тогда т + 1 G 2, n. Имеем с очевидностью (см. (3.40))
(ao)1(r,n — 1) = {ao(j): j G т, n — 1} = {a0(j + 1): j G т, n — 1} (3.50)
= {a0(k): k G т + 1,n} = (a0)1(т + 1,n),
где учитывается то, что по выбору т в силу (3.35) a0(j) = a0(j + 1) при j G т, n — 1. Поскольку выбор т был произвольным, установлено, что
(a0)1(t,n — 1) = (a0)1(t + 1,n) Vt G 1,n — 1. (3.51)
Поэтому из (3.49) вытекает свойство
c«o(t)(z*(t), (a0)1(t7n—Г)) = c«0(t+1)(z°(t + 1), (a0)1(íT^)) Vt G 1,n — 1. (3.52)
Пусть f G 1, n — 1. Рассмотрим представление значения
c(pVz"(f — 1)),pr!(z*(f)), (a0)1 (fn—!))
При этом в силу (3.44) z"(f) = z0(f + 1) и согласно (3.51)
(a0)1(f,n — 1) = (a0 )1(f + 1,n). (3.53)
Отметим, что (f — 1 = 0) V (f — 1 G 1,n — 2). При этом (см. (3.44)) (f — 1 = 0) ^ ^ (pr2(z"(f — 1)) = pr2(z0(f))); кроме того, из (3.44) имеем также импликацию
(f — 1 G 1,n — 2) ^ (z"(f — 1) = z0(f)).
В итоге во всех возможных случаях (pr2(z"(f — 1)) = pr2(z0(f)). Итак, получили, что (во всех возможных случаях)
c(pr2(z*(f — 1),pr^f)), (a0)1 (fn—1)) = c(pr2(z0(f)),pr1(z0(f + 1)), (a0)1(f + ^)).
Поскольку f выбиралось произвольно, установлено, что Vt G 1, n — 1
c(pr2(z*(t — 1)),pr^t)), Ы1^—I))= (3.54)
= c(pr2(z0(t)),pr1(z0(t + 1)), (a0)1(t + 1,n)). Из (3.47), (3.51), (3.52) и (3.54) получаем цепочку равенств
n— 1
Cao [z" | K] = Дф^ (z0(t)), pr1(z0(t + 1)), (a0)1 (¿ТЩ)) + t=1
+ c«o(t+1)(z0(t + 1), (a0)1 (T+M))] + f (pr2(z0(n))) =
n
^[c(pr2(z0(t — 1)), pr1(z0(t)), (a0)1^)) + Cao(t)(z0(t), (a0)1^))] + f (pr2(z0(n))). t=2
Тогда в силу (3.37) и (3.48) получаем неравенство
и ^ ф^ъ0(0)),рГ1(2°(1)), (а0)1^)) + с«0(1)(2°(1), (а0)1(Т,п)) +
п
+ Дф^ОЛ* - 1)),рг^*)), (а0)1^)) + с«0и(ъ°(*), (а0)1 (1^))] + /(рг2(ъ0(п))) =
4=2
п
= Дф^ъ0^ - 1)),РГ1(2°(*)), (а0)1 (¿7П)) + с«о№(ъ°(*), (а0)1^))] + /(р^п))).
4=1
Из (3.7), (3.22) и последнего неравенства получаем, что и ^ Сао[ъ° | К], и, как следствие (см. (3.34)), имеем неравенство
и ^ г>(ж,К). (3.55)
Вернемся к (3.22). Используя (3.22), выберем и зафиксируем д € 1(К) и и € Мд, для которых при
Я = К \ {д} (3.56)
справедливо следующее равенство
и = с(ж, рг1(и), К) + сд(и, К) + фг2(и),Я). (3.57)
Тогда (см. (3.33), (3.56)) имеем, что справедливо равенство
|Я| = п - 1 € 1, N - 1. (3.58)
Ясно, что Я € N. При этом по выбору и имеем, что
(рг1(и) € М) Л (рг2(и) € Мд). (3.59)
Тогда в силу (2.5) и (3.59) рг2(и) € X. Поэтому (рг2(и),Я) € X х N. Тогда согласно (3.4) имеем при а € (Ы) [Я] равенство
2(рг2(и),Я,а) = {ъ € ^^^ | (г(0) = (рг2(и),рг2(и))) (3.60)
&(ъ(*) € Ма(() V* € 1, п - 1)} € Fin(ZQ),
где ^^^ есть множество всех кортежей (г^)^е0)П-1: 0,п - 1 ^ XX х X. Кроме того, согласно (3.2) имеем, что
(I - Ы)[Я] = {а € (Ы)[Я] | а(т) € 1(а1(т,п - 1)) Vm € 1,п - 1} € Р((Ы)[Я]). (3.61) Согласно (3.7) имеем при а € (Ы)[Я] и ъ € 2(рг2(и), Я, а), что
п- 1
£а[ъ | Я] = Дф^ - 1)),рг1(ъ(*)),а1(*7п-Г)) + (3.62)
4=1
+ с«(4)(ъ(*), а1(*;п-Т))] + /(рг2(ъ(п - 1))).
Из (3.10) следует в рассматриваемом случае, что
ЗДг2(и)] = {(а, ъ) € (I - Ы)[Я] х ZQ | ъ € 2(р^и), Я, а)} (3.63)
€ Кп((1 - Ы)[Я] х ZQ).
Согласно (3.12) имеем следующее равенство
^(р^(и),Я) = шт ¿а[ъ | Я].
(а,г)€-0 ^[рг2(и)]
С учетом этого выберем (в°,№ € -Од[рг2(и)] со свойством
^(рг2(и),Я) = ^во| Я]. (3.64)
Тогда в силу (3.63) имеем, что
в° € (I - Ы)[Я] (3.65) и, кроме того, справедливо свойство
№ €2(рг2(и),Я,в°). (3.66) Из (3.61) и (3.65) следует, что в° С (Ы)[Я] и при этом
в°(т) € 1((в°)1(т,п - 1)) Vm € 1,п - 1. (3.67) С учетом (3.58) получаем, что
в° : 1, п - 1 ^ Я
таково, что справедливы свойства
(Я = (в°)1(1, п - 1))&(^1 € 1, п - 1 ^2 € 1, п - 1 (в°(¿1) = в°Ы) ^ (¿1 = ¿2)). (3.68) Далее, из (3.60) и (3.66) вытекает, что № € ZQ, то есть
№: 0, п - 1 ^ XX х X. При этом согласно (3.60) и (3.66) имеем также, что
И0) = (рг2(и), рг2(и)))&(№(*) € Мво(4) V* € 1,п - 1). (3.69)
Заметим, что * - 1 € 1, п - 1 при * € 2, п. Поэтому определены значения
в°(* - 1) € Я V* € 2, п.
Из (3.56) получаем, как следствие, что в°(* - 1) С К V* € 2, п. С учетом этого введем в рассмотрение отображение
в°: 1, п ^ К
посредством следующего правила:
(в°(1) = д)&(в°(*) = в°(* - 1) V* € 2, п). (3.70)
Тогда, как нетрудно проверить, имеет место
в° € (I - Ы)[К] (3.71)
(в связи с проверкой (3.71) отметим подробное рассуждение в [6, с. 67-69]; имеется в виду обоснование свойства [6, (3.2.88)]). Отметим теперь, что согласно (3.69) при * € 2,п имеет
место * - 1 е 1,п - 1, а потому (см. (3.70)) - 1) е Мво(4). Введем в рассмотрение вспомогательный кортеж
г: 1,п ^ XX х X,
определяемый по следующему правилу
А ч п / ~ /А
(гг(1) = и)&(гг(£) = - 1) V* е 2, п). (3.72)
Заметим, что по выбору и имеем в силу (3.70), что гг(1) е М^0(1); кроме того, согласно (3.69) и (3.70) имеем при * е 2, п свойство гг(£) = - 1) е М^0(4). Следовательно, у нас
г(£) е Мво(4) V* е 1, п. (3.73)
С учетом (3.73) дополним кортеж г начальной УП (ж,ж). Итак, введем кортеж гг е , полагая (см. (3.22)), что
А / ч ч „ / Л /.ч А
(г(0) = (ж, ж))&(гг(£) = гг(£) V* е 1, п). Из (3.4), (3.22) и (3.73) вытекает следующее очевидное свойство
г е2 (ж, К, в °), (3.74)
а тогда (см. (3.10), (3.71)) получаем включение
(в °,г) е Б к [ж]. Тогда в силу (3.12) реализуется следующее неравенство:
г(ж, К) ^ Сво [г | К]. (3.75)
С учетом (3.7) и (3.22) имеем (см. (3.71), (3.74)), что
п
С во [г I К] = 5>(рг2(г(* - 1)), рг1(гг(£)), (в °)1(£7п)) +
4=1
+ св>(4)И*), (вУ^))] + /(рг2(гг(п))). Поэтому с учетом (3.75) получаем неравенство
п
г(ж,К) ^[с(рг2(гг(* - 1)),рг^)), (вУ^)) + (3.76)
4=1
+ Сво(4)(гг(*), (в°)1 (*, п))] + /(рг2(гг(п))). Отметим, что в силу биективности в° (см. (3.71)) справедливо равенство
(в°)1(1,п) = К. (3.77)
При этом по построению гг имеем, что рг2(гг(0)) = ж и рг1(гг(1)) = рг1(гг(1)); с учетом (3.72) и (3.77) получаем теперь, что
с(рг2(гг(0)),рг1 (гг(1)), (вУ^п)) = с(ж,рг1(гу(1)),К) = с(ж,рг1(и),К),
а тогда в силу (3.76) реализуется неравенство
п
г(ж,К) ^ с(ж,рг1(и),К) + Св0(1)(№(1), (в°)1(1;п)) ^[с(рг2(^Й(* - 1)), (3.78)
4=2
рг1 (№(*)), (вУ^)) + ов0(4)(№(*), (в0)1^))] + /(рг2(гЙ(п))). Заметим, что согласно (3.70), (3.77) имеем по построению №, что
^(№(1), (в°)1(1;п) = сд(№(1),К). Поэтому из (3.78) вытекает неравенство
п
г(ж,К) ^ с(ж, рг1 (и), К) + сд(№(1),К) + ^[с(рг2(№(* - 1)), р^(№(*)), (вУ^)) +
4=2
+ Ов0(4)(гЙ(*), (в0)1^))] + /(рг2(№(п))). Учитывая (3.72) и то, что п € 1, п, получаем теперь, что
п
г(ж,К) ^ с(ж, рг1(и),К) + сд(и, К) + ^[с(рг2(№(* - 1)), рг1(гу(*)), (в0)1^)) +
4=2
+ св0(4)(гЙ(*), (в0)1^))] + /(рг2(№(п))).
Но в нашем случае п € 2, п (см. (3.30)), а потому (см. (3.72)) №(п) = №(п - 1) и, следовательно,
п
г(ж, К) ^ с(ж, рг1(и),К) + сд(и, К) + Дф^И* - 1)), рг1 (№(*)), (вУИ) + (3.79)
4=2
+ Св0(4)(гЙ(*), (вУ^))] + /(рг2(№(п - 1))).
Заметим, что при * € 2, п имеем п С 2, п, а потому согласно (3.70)
(в°)1(*7п)) = {в°(т): т € = {в°(т - 1): т € ¿7^} = (3.80)
= {в°(£): £ € * - 1,п - 1} = (в°)1(* - 1,п - 1).
Далее по построению № имеем, что (см. (3.72))
(№¿(1) = №(1) = и) & (№¿(2) = №(2) = №(1)).
Поэтому имеем следующее равенство
с(рг2(№(1)), рг1(№(2)), (в°)1(2;п)) = с(рг2(и), рг^^)), (в°)1(17п-1)),
где учтено также (3.80). Вместе с тем, используя (3.69) получаем теперь, что
с(рг2(№(1)),рг^)), (в°)1(2Тп)) = с(рг2(№(0)), рг^^)), (в°)1(^Г-Г)). (3.81)
С другой стороны, при * € 2, п \ {2} имеем, что * € 3, п, а потому (см. определение №)
ф^Н* - 1)), р^ (№(*)), (в0)^)) = - 1)), р^ (№(*)), (в° - 1,п - 1)) =
= с(рг2(№(* - 2)), рг1(№(* - 1)), (в°)1(* - 1,п - 1)).
С учетом (3.81) получаем теперь, что
с(рг2(гг(* - 1)), рГ1(гг(*)), (в°)1(*,п)|= (3.82)
с(рг2(г(* - 2)), рг1(г(* - 1)), (в°)1(* - 1,п - 1)) V* е 2, п.
Заметим теперь, что при * е 2,п согласно (3.70), (3.72) и определению гг
С0О(4)(г(*), (в°)1(*7п)) = сво(4_1)(гг(*), (в°)1(* - 1,п - 1)) = = сво(4_1)(г(* - 1), (в°)1(* - 1,п - 1)).
С учетом (3.79) и (3.82) получаем теперь, что
г(ж, К) ^ с(ж, рг1(и), К) + сд(и, К) + ^[с(рг2(г(* - 2)), рг1(г(* - 1)), (3.83)
4=2
(в°)1(* - 1,п - 1)) + сво(4_1)(г(* - 1), (в°)1(* - 1,п - 1))] + /(рг2(г(п - 1))) =
п_ 1
= с(ж, рг1(и), К) + сд(и, К) + ^][с(рг2(г(г - 1)), рг^Цт)), (в°)1(т,п - 1)) +
+ сво(т)(г(т), (в°)1 (т, п - 1))] + /(рг2(г(п - 1))). Вместе с тем согласно (3.62), (3.65) и (3.66) имеем равенство
п_ 1
| VI = ^[с(рг2(г(т - 1)), рг (г(
т=1
+ сво(т)(г(т), (в°)1(Т;п-Т))] + /(рг2(г(п - 1))).
Сво [г | V] = ^>(рг2(г(т - 1)), рг1 (г(т)), (в°)1(т,п - 1)) +
Тогда из (3.83) вытекает следующее неравенство:
г(ж, К) ^ с (ж, рг1(и),К) + сд (и, К) + ¿во [г | V]. В силу (3.64) получаем неравенство
г(ж, К) ^ с(ж, рг1(и), К) + сд(и, К) + г(рг2(и), V). C учетом (3.57) получаем следующее неравенство
г(ж, К) ^ ш,
откуда в силу (3.55) имеем, что г(ж, К) = ш и при условии (3.30). Итак, установлена импликация
(п е 2,Ж) ^ (г(ж, К) = ш). (3.84)
Тогда согласно (3.23), (3.29) и (3.84) имеем во всех возможных случаях равенство г (ж, К) = = ш, а потому (см. (3.22)) получаем требуемое равенство (3.21). □
С учетом (3.15) и теоремы 3.1 получаем, что
V[ж] = Ш1п шт[с(ж, рг1(г), ) + с,(г, )+ г(рг2(г), \ {?})] Уж е Х°. (3.85) ,е1(1,м) ¿ем,
§ 4. Динамическое программирование, 2
В этом параграфе конструируется оптимальный алгоритм на базе ДП. Oсновная идея состоит в том, чтобы при наличии условий предшествования подменить построение всего массива значений функции Беллмана построением системы некоторых ее слоев, которые получаются насчитыванием ее значений на некоторых специальных п/м области определения (данной функции). Эти п/м будем называть слоями пространства позиций. Упомянутая процедура восходит к [6, §4.9]. Ее расмотрение начинаем с введения существенных (для дальнейшего) списков заданий. Итак, следуя [6, (4.9.1)] полагаем, что
д = {К е N IV* е К (рг1(г) е К) ^ (рг2(*) е К)}; (4.1)
множества — элементы семейства (4.1) — называем существенными списками (заданий); см. также [19, (6.1)]. Данные списки ранжируем по мощности, полагая, что
д* = {к е д | в = |к|} V« е туж (4.2)
Семейства (4.2) составляют разбиение д (4.1). Легко видеть, что
(д* = {т, N})&(д1 = {{г}: г е т, N \ К1}), (4.3)
. ш
где К1 = {рг1(*): * € К}. Кроме того (см. [21]),
д*-1 = {к \{г}: к ед*е 1(к)} V* е 2,N. (4.4) Посредством (4.3), (4.4) определена рекурентная процедура
д* ^д*-1 ^... (4.5) Перейдем к построению слоев пространства позиций. Пусть (см. [3, (3.4.18)])
М = у М; е Р'(М). (4.6)
Введем в рассмотрение множества Д0, Д1;..., в пространстве позиций; эти множества
именуем слоями пространства позиций. Проще всего определяются Д0 и — крайние слои:
(До = {(х, 0): х е М})&(Д* = {(х, Т^): х е X0}). (4.7)
Если в е 1, N — 1 и к е дв, то последовательно определяем
Л(к) = О е т^ \ к I О} и к е д^},
млк]= и м^,
¿€Л[К]
] = {(х,к): х е М*[к]}.
Тогда при в е 1, N — 1 слой Д* определяем правилом
Д* О*[к]. (4.8)
к
Данное построение слоев Д0, Дь..., Ду соответствует [22, раздел 2.9]. Все упомянутые слои — суть непустые множества (см. [22, (2.73)]). При этом Д0 С М х {0}, Ду = X0 х х } и Да С М х при ^ € — 1. Отметим важное свойство [22, (2.74)]: если
5 € 17Ж (ж, К) € Да, 3 € 1(К) и г € М,, то
(рг2(г),К \{з}) € Да-1- (4.9)
Слои функции Беллмана: построение и анализ. Рассмотрим сначала алгоритм построения функций
зд €Я+[До], V! €Я+[Д1], ..., зд €^+[Ду]. (4.10)
Затем будет показано, что построенные функции являются сужениями функции Беллмана (3.18). Наконец, совсем кратко обсудим вариант построения только одной функции зд. Итак, возвращаясь к (4.10), полагаем, что функция г>0 € ^+[Д0] определяется правилом
1>о (ж, 0) = /(ж) Уж € М (4.11)
(учитываем (4.7)). Дальнейшее построение осуществляется посредством рекуррентной процедуры, использующей (4.9).
Итак, пусть 5 € и функция г>а-1 € ^+[Да-1] уже построена. Тогда согласно (4.9) определены значения
^а-1(рг2(г),К \{3}) € Е+ У(ж, К) € Да У? € 1(К) Уг € М,.
С учетом этого при в € 1, N полагаем, что функция € [Да] определяется правилом
^а(ж,К)= шт ш1п[е(ж, рг1(г),К) + с,(г, К)+ ^а-1(рг2(г), К \ {3})] (4.12)
У(ж, К) € Да.
Таким образом (см. (4.11), (4.12)) реализуется рекуррентная процедура
^0 ^ г>1 ^ ... ^ гу. (4.13)
В связи с (4.13) отметим, что при в € 0^ и (ж, К) € Да определено согласно (3.18) значение г>(ж, К) € Е+.
Предложение 4.1. Если в € 0, N и (ж, К) € Да, то г>а(ж, К) = г>(ж, К).
Доказательство. Из (3.17), (4.7) и (4.11) следует, что г>0 (ж, К) = г>(ж, К) У(ж, К) € € Д0. Пусть теперь
NN = {в € 0, N | иа(ж, К) = и(ж, К) У(ж, К) € Да}. (4.14)
Тогда 0 € N. Покажем, что NN = 0, N .В самом деле, допустим противное: пусть NN = 0, N.
Тогда в силу (4.14) 0, N \ NN = 0. В силу конечности множества 0, N \ NN имеем, что
△ .
п = т£ (0, N \ N € 0, N \ NN. (4.15)
Тогда п € 0, N \ NN и при этом п ^ к Ук € 0, N \ NN. Поскольку 0 € 5:, то п = 0. Это означает, что п € 1, N. Далее, п — 1/ 0, N \ NN по выбору п. Следовательно,
п — 1 € 0, N \ (0, N \ ЭТ),
то есть п — 1 е N. Поэтому п — 1 е 0, N и согласно (4.14)
ига_1 (х, к) = ^(х, к) V(x, к) е Д„-1. (4.16)
При этом (см. (4.9)) V(x, к) е Д„ V; е 1(к) V* е М,
(рг2(*),к \{;}) е Д„_1. (4.17)
Поэтому согласно (4.16), (4.17) имеем, что
^_1(рг2(*),к \ {;}) = v(рг2(*), к \ {;}) ^х,к) е Д„ V; е 1(к) V* е М,. (4.18)
При этом согласно (4.12) и (4.18) получаем представление
^(х,к)= .тш тт[с(х, рг1(*),к) + с, (*,к)+ (4.19)
+ ^_1(рг2(*),к})] = гап тт[с(х, рг1(*),к) +
,е1(к) ¿ем^
+ с,(*, к)+ v(pг2(z), к \ {;})] V(x, к) е Дга.
При этом Дп С X х дп С X х N. Поэтому согласно теореме 3.1
тт тт[с(х, рг1(*), к) + с,(*, к) + v(pг2(z), к \ {;})] = v(x, к) V(x, к) е Дп. ,е1(к) ¿ем^
С учетом (4.19) получаем, что ^(х,к) = v(x,K) V(x,K) е Дп. В силу (4.14) получаем, что п е N что противоречит выбору п. Полученное противоречие показывает, что свойство N = 0^ невозможно и, следовательно, N = 0^. Поэтому (см. (4.14)) ^(х, к) = v(x, к) V* е 0, N V(x, к) е Д*. Предложение доказано. □
Итак, при в е 0^ функция vs е есть сужение функции v на множество Д*.
В частности,
^(х, к) = v(x, к) V(x, к) е Д*. (4.20)
Из (3.15) и (4.20) вытекает следующее свойство:
V[х] = ^(х, 1, N) Vx е X0. (4.21)
С учетом (4.7), (4.12) и (4.21) получаем, что
V[х] = т]п тт[с(х, рг1(*), 1,N) + с,(*, 1,N)+ зд_1(рг2(*), \ {;})] (4.22)
Vx е X0;
при этом, конечно, в силу (4.7) и (4.9) при ; е 1(1, N) и * е М,
(рг2(*), 1, N \{;}) е Д*_1 (4.23)
(здесь полезно учитывать свойство X0 = 0). Итак, рекуррентная процедура (4.13) приводит, в частности, к функции
VН = ^[х])хехо 0], (4.24)
что позволяет достаточно просто решить задачу (2.20) и определить глобальный экстремум V (2.18) и множество X(0pt (2.21). В этой связи, однако, полезно следующее замечание.
Замечание 2. Отметим естественную аналогию с конструкцией [23], применяемой для решения другой задачи. Для этого заметим, что согласно (4.12) для построения функции vs, где s € 1, N, требуется (помимо c и ci,..., cn) только функция vs-i, то есть предыдущий слой функции Беллма-на. В этой связи возникает следующий алгоритм решения задачи (2.20).
Итак, используя (4.11), определяем функцию vo € R+[Do]. Пусть s € 1,N таково, что функция vs-i € R+[Ds-i] уже построена. Тогда определяем функцию vs € R+[Ds] по правилу (4.12). Если при этом s = N, то построение vn завершено. Если же s < N, то массив значений функции vs-i уничтожается и заменяется массивом значений построенной функции vs, который используется для построения vs+i по правилу
Vs+i(x,K)= min min [c(x, pri(z),K) + Cj (z, K) + vs(pr2(z),K \ {j})] V(x, K) € Ds+i, jei(K) zeMj
после чего данный массив (значений функции vs) уничтожается. Как видно из построения, в памяти вычислителя при данном варианте находится всякий раз только один слой функции Беллмана, что доставляет определенную экономию ресурсов памяти (см. в этой связи [24]).
Результатом данной процедуры является уже функция vn (а не кортеж v0, vi,..., vn), что, однако, вполне достаточно для определения функции V[ ] (4.24) и последующего решения задачи (2.20), включая определение V и построения множества (2.21). □
Вернемся к основной задаче (2.17), предполагая завершенной процедуру (4.13) с нахождением всех функций v0, vi,..., vN. Рассмотрим схему построения решения из множества SOL, то есть схему построения оптимального решения полной задачи (2.17).
Прежде всего мы решаем задачу (2.20), располагая функцией vN (см. (4.21)). Итак, (2.20) сводится к задаче
Vn(x, ^ min, x € X0. (4.25)
Поскольку функция vN известна, определяем V и (оптимальную) точку старта x0 € X(°¡lpt:
vn(x0,17Ñ) = min vn (x, VÑ) = V. (4.26)
xei0
Итак, задача (4.25) доставляет решение задачи (2.20) и мы фиксируем, наряду с V, точку x0 € X0 со свойством (4.26) (это и означает, что x0 € X(pt). Мы полагаем, что
z(0) = (x0,x0). (4.27)
В силу (4.22), (4.26) имеем следующую цепочку равенств
V = V [x0 ] = VN (x0,17Ñ) = min min[c(x0, pri(z), 1,Ñ) +
jei(i,N) zeMj + Cj (z, 1,Ñ) + VN-i(pr2(z), \ {j})].
С учетом этого находим € I(1, N) и z(i) € МП1, для которых
V = c(x0, pri(z(i)), 1,N) + cn1 (z(i), 17N) + VN-i(pr2(z(i)), 1,N \{ni}) (4.28)
(решаем задачу на минимум); при этом согласно (4.23)
(pr2(z(i)), 17N \{ni}) € Dn-i.
Поэтому согласно (4.12) имеем равенство
vN-i(pr2(z(i)), 1,N \ {ni}) = min min[c(pr2(z(i)), pri(z), 1,N \ {^i}) +
jei(i,N\{ni}) zeM
+ Cj (z, i;N \{ni}) + Vn-2(pr2(z), ITN \{W j})].
С учетом этого определяем (из решения задачи на минимум) щ2 £ I(1, N\ {ni}) и z(2) £ Mn2 со свойством
VN-1 (pr2(z(1)), тух\{ni}) = c(pi2(z(1)), pri(z(2)),^ \ {ni}) + (4.29)
+ СП2 (z(2), VÑ \{ni})+ VN-2(pr2(z(2)), ^ \ {ni; П2 }) •
Теперь из (4.28) и (4.29) получаем следующее равенство
V = c(x0, pri(z(i)), VÑ) + c(pr2(z(i)), pri(z(2)),тух \ {ni}) + C41 (z(i),T,Ñ) + (4.30) + СП2 (z(2), \ {ni} ) + VN-2 (pr2 (z(2)), t;ñ\{ ni; n2} ) = = c(pr2(z(0)), pri(z(i)), 1, N) + c(pr2(z(i)), pri(z(2)), 1, N \ {ni}) + + Cni(z(i), T^Ñ) + СП2(z(2),T;Ñ \ {ni}) + VN-2(pr2(z(2)),I^N \ {щ; n2})•
Замечание 3. При N = 2 имеем свойство
((n)¿gi,2, (z^) £ (SOL)[x0].
Проверим данное простое следствие (4.30). Прежде всего, по выбору n2 имеем (см. (3.1)), что ni = n2, а тогда в силу цепочки равенств 1,N = 1, 2 = {1;2} имеем, что {ni; n2} = = 1, 2 = 1,N, а (nOjgi^ есть перестановка 1,N = 1, 2, то есть (n^gi^ £ (bi) [1, N]. Далее, отметим, что ni £ I({n: i £ 1, 2}) и 1, N \ {ni} = 1, 2 \ {ni} = {n2}, а потому n2 £ I({n2}). Поэтому в нашем случае nm £ I({n : i £ m, N}) Vm £ 1, N. Тогда согласно (3.2) (ni)ie^ £ £ (I-bi)[17N], а потому (см. (3.3)) (щ)iei2 £ A. Далее, z(0) = (x0,x0) £ X0 x X0, z(i) £ M41 и z(2) £ Mn2, а потому согласно (2.9) (z(i))ie02 £ Z(ni).eY2 [x0]. В итоге (см. (2.10))
(Ыгда (z^W £ dd[x0].
Вернемся к (4.30) (при N = 2). При этом 1, N \ {щ1; щ2} = 0, а тогда
VN-2(pr2(z(2)), 1TN \ {ni; n2}) = V0(pr2(z(2)), 0) = f (p^(z(2))).
Поэтому в нашем случае из (4.30) следует, что
V = c(pr2(z(0)),pr1(z(1)), H^W) + c(pr2(z(1)),pr1(z(2)), \ {ni}) + + Cni 1z(1), l^) z(2), ^ \ {ni}) + f (pr2(z(2))) = = c(pr2(z(0)), pri (z(1)), {ni: i £ 1,N}) + c(pr2(z(1)), pri(z(2)), {^}) + + Cni(z(1), {ni: i £ ^}) + Cn2 (z(2), {%}) + f (pVz(2))).
C учетом (2.13) получаем (при N = 2) равенство V = C(ni) _ [(z(i))ie0Ñ]. В силу (2.19) это означает, что
((ni^iÑ, (z(iW,x0) = (Ы^, (z^c^x0) £ SOL,
где учитывается, конечно, равенство (2.11). Итак, в нашем простейшем случае N = 2 построения (4.28), (4.29) непосредственно приводят к оптимальному решению. □ Отметим, что в общем случае N, N ^ 2, процедуры решения локальных задач, подобные (4.28), (4.29), следует продолжать вплоть до исчерпывания полного списка 1,N. В результате будут построены
n = (n)jeiÑ £ A £ Zn[x0]
со свойством [(z(j))je0Ñ] = V = V[ж0]. Тогда
(П, G (SOL)[x0]. (4.31)
При этом, в частности, (n, (z(j))jeo;Ñ) G DD[x0] согласно (2.10), а потому (n, (z(j))jeo;Ñ,x0) G G D. Более того, поскольку [(z^^go^] = V, имеем из (2.19), что (см. (4.31))
(П, (z(j))jeo;Ñ,x0) G SOL. (4.32)
Итак, определены глобальный экстремум V и оптимальное решение (4.32) нашей основной (полной)задачи.
§ 5. К вопросу о планировании системы работ в условиях повышенной радиации
В этом параграфе на содержательном уровне обсудим одну особенность задачи о последовательном демонтаже системы радиационно опасных элементов. Это обсуждение позволяет и несколько дополнить весьма общие построения предыдущего параграфа. Речь идет о том, что в условиях упомянутой конкретной задачи демонтажа исполнитель по достижении некоторого уровня d суммарной радиации должен сниматься с выполнения работ и выводиться из зоны повышенной радиации. В этой связи существенным является соотношение V и d. Если V ^ d, то ранее упомянутая система работ по демонтажу может быть завершена при использовании оптимального ДР (4.32). В случае, когда d < V, осуществить все намеченные работы не предоставляется возможным, так как их проведение в любом случае (то есть при использовании любых ДР) будет сопряжено с недопустимо высокой дозовой нагрузкой. В этой связи представляется важным уже на этапе постановки осуществление прогноза в вопросе соотношения значений V и d. Представляется, что построения предыдущего раздела и на этом этапе могут быть полезны; в частности, здесь представляется полезным вариант процедуры, отмеченный в замечании 2.
Будем предполагать, что после завершения всех работ, связанных с посещением мегаполисов, или по достижении дозы радиации, определяемой значением d, исполнитель перемещается в заданный пункт эвакуации ж* G X .В случае демонтажа всех источников в процессе этого перемещения он уже не испытывает ощутимого радиационного воздействия; имеется лишь некоторый остаточный фон, которым на данном этапе рассмотрения будем пренебрегать. Если же d < V, то есть исполнитель не может выполнить всех работ в рамках предъявляемых требований безопасности, то следует сократить число мегаполисов, подлежащих обслуживанию, выбирая из N мегполисов первоначальной задачи некоторое меньшее их число Ni, Ni < N .В этом случае решая по методу ДП задачу о демонтаже Ni источников мы сталкиваемся с той ситуацией, что после выполнения работ по демонтажу этих N1 источников при перемещении к пункту эвакуации ж* исполнитель будет находиться под ощутимым радиационным воздействием тех источников, которые соответствуют мегаполисам, исключенным из списка заданий в связи с неравенством d < V. Это повлечет изменение терминальной компоненты f нашего аддитивного критерия: «нулевая функция» для случая выполнения всех работ по демонтажу N мегаполисов должна быть заменена функцией, оценивающей дозу, получаемую исполнителем на этапе эвакуации. Новое значение глобального экстремума V1 с ненулевой терминальной функцией не должно превосходить d. Это условие должно, конечно, обеспечиваться рациональным выбором Ni мегаполисов.
Приведенные выше рассуждения показывают, насколько важна задача (2.20) на этапе планирования работ: к построению решения задачи (2.17) следует (по большому счету) приступать лишь тогда, когда удается обеспечить неравенство V ^ d, а это можно проверить, используя только процедуру отмеченную в замечании 2 и реализуемую с некоторой
экономией ресурсов памяти вычислителя. Итак, данная процедура может быть весьма полезной на этапе планирования работ и, в частности, для целей отбраковки недопустимых (по степени радиационного воздействия на работников) вариантов набора заданий.
§ 6. Некоторые особенности системы целевых функций в задаче о демонтаже
Отметим, что конструкция решения задачи маршрутизации на основе ДП, изложенная в §§ 2-4, является весьма общей и применимой для самых различных инженерных задач. В частности, это отражено в [1,2,25], где рассматривалось применение данной конструкции для решения задачи управления инструментом при листовой резке даталей на машинах с ЧПУ. Возможны и другие применения упомянутой схемы решения на основе ДП.
В предыдущем параграфе, однако, было обращено внимание на одну особенность, присущую именно задаче о последовательном демонтаже радиационно опасных элементов. Сейчас мы продолжаем детализацию общих построений, ориентируясь на последнюю задачу. В дальнейшем полагаем, что X есть конечномерное арифметическое пространство фиксированной размерности (два или три), выпуклые оболочки мегаполисов М1,..., М* попарно дизъюктны, и при каждом ; е 1, N источник с номером ; есть точка из выпуклой оболочки М,. В этой связи напомним прежде всего (2.12), учитывая, что в данной конкретизации общих положений значения функций с, с1,..., с*, f имеют смысл доз радиации, получаемых исполнителем при выполнении операций по демонтажу. Отметим, что в части, определяющей зависимость от списка заданий (то есть множества семейства имеет место следующее свойство аддитивности:
(с(х, у, к) = ^ с(х,у, {к}) Vx е X ^ е X VK е N (6.1)
йек
&(с,(*, к) = ^ с,(*, {к}) V; е ^ V* е М, VK е
йек
Иными словами, первые N+1 функции в (2.12) полностью определяются той «частью» этих функций, которая отвечает воздействию одиночных источников. Следовательно, в данной задаче определяющими являются функции
(х, у, к) ^ с(х, у, {к}): X х X х ^ Е+, (6.2)
(*,к) ^ с1(*, {к}): М1 х 1, N ^ Е+, ..., (*, к) ^ с* (*, {к}): М* х ^ Е+. (6.3)
Итак (см. (6.2), (6.3)), определяющими являются функции, характеризующие всякий раз радиационное воздействие того или иного одиночного источника (реализуется естественный вариант декомпозиции). Заметим, что подробный анализ таких функций, включая их построения, был проведен в [26], куда мы и отсылаем за подробностями. Сейчас же отметим на содержательном уровне некоторые особенности функций (6.3) в сравнении с (6.2). Функция (6.2) оценивает внешние перемещения без выполнения каких-либо работ; ее значения при фиксации к е получаются интегрированием нелинейной зависимости (имеется в виду обратная пропорциональность квадрату расстояния до источника с номером к) при перемещении с достаточно большой скоростью; имеются некоторые особенности типа прохождения через источник, которые учитываются введением соответствующих штрафов.
При ; е 1, N значения функции
(*, к) ^ с,(*, {к}): М, х ^ (6.4)
имеют более сложную структуру. При к € 1, N \ (3} и г = (х,у) € М, значение с,(г, {к}) реализуется в виде суммы трех слагаемых, отвечающих «далекому» воздействию источника с номером к на исполнителя, выполняющего работы, связанные с посещением М,. Речь идет о дозе радиации, получаемой при перемещении из х € М, к точке вблизи источника с номером 3, дозе, получаемой при выполнении собственно демонтажа без каких-либо существенных перемещений, и, наконец, о дозе, получаемой в результате воздействия источника с номером к при перемещении исполнителя от демонтированного 3-го источника к точке у. Если же к = 3, то с,(г, {к}) = с,(г, (3}) реализуется в виде суммы дозы, получаемой при перемещении от х к точке вблизи источника («близкое» воздействие) и дозы, получаемой при осуществлении демонтажа самого источника с номером 3. В рассматриваемой ниже модели предполагается, что скорость перемещений при выполнении внутренних работ существенно меньше скорости внешних перемещений.
Терминальная компонента критерия f в простейшем случае может полагаться тождественно равной нулю. Это соответствует случаю, когда V ^ в, и все задания по демонтажу могут быть выполнены (см. § 5). Если же это не так, то после выполнения части работ по демонтажу остаются излучающие элементы. Будем сейчас для простоты считать, что эта часть работ исчерпывается нашими мегаполисами М1,..., Мм (в противном случае следует сделать необходимые переобозначения). Вместе с тем остаются источники (их влияние, конечно, должно учитываться и при построении с, с1,..., см так же, как и влияние источников, связанных с М1,..., Мм), которые продолжают оказывать воздействие на исполнителя при его движении к пункту эвакуации. Это воздействие оценивается функцией f, построение которой подобно в идейном отношении построению с; имеется в виду определение конкретной дозы радиации при перемещении от последнего из обслуживаемых мегаполисов к пункту эвакуации.
Заметим, что в последнем случае мы фактически рассматриваем задачу с большим количеством радиационно опасных источников, часть из которых, допускающую обслуживание с непревышением дозы в, мы выделяем для демонтажа, но в построении функций (2.12) учитываем при этом влияние уже всех радиационно опасных объектов. С этой постановкой связана, конечно, следующая проблема: какие именно источники из имеющейся полной совокупности следует выбрать для осуществления демонтажа? Решение этого вопроса можно связать с построениями по схеме замечания 2, выполняемыми для различных вариантов заданий по демонтажу. В последующих построениях мы данную детализацию f не рассматриваем, ориентируясь в примерах на вариант, связанный с выполнением всех заданий, то есть по сути на ситуацию V ^ в.
Вернемся к (6.1)-(6.4), используя построения [5, раздел 3]. Совсем кратко обсудим первое выражение в (6.1). Точнее, при к € 1, N полагаем, что с(к) € ^+[Х х X] определяется следующими условиями: при х € X и у € X
Итак, в (6.5) мы рассматриваем дозу, создаваемую (одним) источником с номером к при перемещении исполнителя из х в у. Конкретное значение (6.5) вычисляется по формулам, приведенным в [26, раздел 6] (отметим, что в [26] исследовался вопрос о построении оптимизирующих вставок в эвристические решения, но представления [26, раздел 6] касаются нахождения значений (6.5) с выделением регулярного и аномального случаев). Итак (см. (6.1), (6.5)), при х € X, у € X и К € N
с(к)(х,у) = с(х,у, {к}).
(6.5)
Функции стоимости внутренних работ.
Рассмотрим построение функции с,, фиксируя для сокращения формулировок ; е 1,N. Мы исходим при этом из второго соотношения в (6.1). В этой связи полагаем, что при
к е и ; е
е гс+ДО,-])&(, е №+)&(, е ]) (6.7)
таковы, что реализуется представление
с,(*,{к}) = с^о*))++с(;з)(рг2(*)) V* е м,. (6.8)
Кроме того, полагаем, что при к е 1,N, ; е и х е М, значение с(д(х) есть доза радиации, определяемая действием источника с номером к на исполнителя при перемещении последнего из пункта х в точку (демонтажа), близкую к источнику с номером ;. Эту последнюю точку фиксируем на сегменте прямой, проходящей через х и ;-й источник и отстоящую от последнего на заданном малом расстоянии, зависящим от Конкретные значения с(А:1) определяются по аналогии со значениями функций с(1),..., с(*) посредством представлений [26, раздел 6] (см. регулярный случай внешнего перемещения в [26, с. 99]).
„„00„, „„„ -7М ™ о™ △
Далее, при к е 1, N и ; е 1, N константа с()2 имеет вид [5, (3.7)], то есть с(;2 = 7^ А . Здесь есть время, необходимое для демонтажа ;-го источника, а 7^ е \ {0} при ; = к есть коэффициент, зависящий от интенсивности источника с номером к и расстояния между источниками; при ; = к в виде 7^ = 7,, имеем константу, зависящую как от интенсивности ;-го источника, так и от расстояния между точкой вблизи этого источника (точка демонтажа), в которую прибывает исполнитель из пункта х, и точкой, где расположен сам источник. Итак, с, 2 е есть доза радиации, получаемая при демонтаже; исполнитель останавливается, не доходя до источника и, с применением соответствующих технических средств «выключает» источник, затрачивая на это время А^.
Значение ^(у), где к е , ; е и у е М,, определяется по разному для
случаев, когда ; = к и ; = к. Если ; = к, то данное значение находится подобно значениям функций с(1),..., с(*) и соответствует регулярному случаю внешнего перемещения в [26, раздел 6] (случай «дальнего» воздействия на исполнителя). Если же ; = к, то
с$5(у) = <й(у) = 0.
В последующем изложении мы ограничиваемся случаем, когда f (у) = 0, считая, что удается демонтировать все источники, и пренебрегая остаточным фоном при перемещении исполнителя к пункту эвакуации.
Таким образом, получаем задачу маршрутизации с достаточно сложными функциями стоимости. Ниже рассматриваются варианты решения модельных задач; в каждом из них найден глобальный экстремум V и оптимальный маршрутный процесс в ситуации, когда терминальная функция тождественно равна нулю.
§ 7. Вычислительный эксперимент
В этом параграфе сначала мы продолжим детализацию общей постановки. Прежде всего условимся, что в рассматриваемых ниже задачах X = М х М, то есть рассматривается плоский случай. Мы следуем при этом соглашениям предыдущего параграфа. В рассматриваемых вариантах были выбраны «большие» мегаполисы (первый пример) и «малые» мегаполисы (второй пример). Было обнаружено, что во втором примере время счета уменьшалось более чем в два раза, что косвенно говорит о том, что в вопросах, связанных с размерностью, допускающей эффективное решение, существенно не только значение N, но и мощности |М1|,..., |М*| (прочие основные численные параметры совпадают в обоих примерах).
Рассматриваемые в данной работе алгоритмические конструкции были реализованы в виде программы для ПЭВМ (написана на языке программирования C++), работающей под управлением 64-хразрядной операционной системы Windows (начиная с Windows 7). Вычислительная часть программы выполняется в отдельном от интерфейса пользователя потоке. Исходные данные и результаты счета хранятся в текстовом файле специальной структуры. Для случая решения задачи на плоскости имеется возможность графического отображения «мегаполисов», а также маршрута и трассы их посещения, при этом масштаб выбранного участка графика может быть увеличен; возможно сохранение графика в файл графического формата bmp в цветном или черно-белом виде.
Вычислительный эксперимент проводился на ПЭВМ с центральным процессором Intel Core i7, объемом ОЗУ 64 ГБ с установленной операционной системой Windows 7 Максимальная SP1. Рассматривалась задача оптимизации обхода «мегаполисов», представляющих собой равномерные «сетки», получаемые размещением точек на равных угловых расстояниях на окружностях. Источники излучения, ассоциированные с «мегаполисами», расположены в пределах контуров вышеупомянутых окружностей. Задано несколько начальных точек траектории для выбора оптимальной точки старта. Предполагается, что в нашей модели после выполнения всех работ исполнитель остается в последней точке траектории. это соответствует предположению о том, что перемещения в точку эвакуации происходят при нулевом фоне и потому не оцениваются. Количество мегаполисов и количество адресных пар, накладывающих ограничения на порядок посещения «мегаполисов», во всех примерах равно 33. Предполагается, что скорость перемещения вне контуров «мегаполисов» (внешние перемещения) в 4 раза больше скорости перемещения исполнителя внутри «мегаполисов» (внутренние работы). Интенсивности излучения, длительности выполнения внутренних работ непосредственно возле источников излучения, а также радиусы окрестностей вокруг источников для выбора точек производства работ по демонтажу задаются в качестве параметров задачи независимо для каждого «мегаполиса».
Пусть в рассматриваемых примерах интенсивности излучения варьируются в пределах от 1.3 до 5.5, радиусы окрестностей источников излучения от 1.1 до 1.4, длительности выполнения работ по демонтажу источников находятся в пределах от 1.1 до 1.7. На рисунках значком * отмечены источники излучения, а ▽ _ точки выполнения работ по демонтажу возле источников. Полагаем далее, что M = Mj х Mj при j £ 1, N.
Рассмотрим пример с «мегаполисами», каждый из которых состоит из 12 точек на окружности. Множество допустимых начальных точек траектории:
{(90, 40); (0, 0); (-60, -5); (-70, -100); (70, -50)}. Получены следующие результаты:
Величина совокупных затрат (суммарная доза) составила 226.797. Выбрана начальная точка траектории (-60, -5).
Первая упорядоченная пара траектории ((-85, -28.66); (-98.66, -25)) £ М2з х М2з. Финальная упорядоченная пара перемещений ((-24, -96.93); (-12, -90)) £ M19 х M19. Время счета составило 22 час. 58 мин. 3 сек. График маршрута и трасса приведены на рис. 1.
В следующем примере «мегаполисы» содержат 6 точек на окружности. Задано следующее множество допустимых начальных точек траектории:
{(90, 45); (0, 0); (-70, -100); (40, 0); (-80, 35)}.
Получены следующие результаты:
Величина совокупных затрат составила 227.755.
Выбрана начальная точка траектории (-80,35).
60 40 20 0 -20 -40 -60 -80 -100
-100 -80 -60 -40 -20 0 20 40 60 80
Рис. 1. Траектория движения по мегаполисам (пример 1)
60 40 20 0 -20 -40 -60 -80 -100
-100 -80 -60 -40 -20 0 20 40 60 80 Рис. 2. Траектория движения по мегаполисам (пример 2)
Первая упорядоченная пара траектории ((-100,10); (-100,10)) G M31 х M31. Финальная упорядоченная пара перемещений ((-85, 46.34); (-80, 55)) G Mi6 х Mi6. Время счета составило 5 час. 53 мин. 43 сек. График маршрута и трасса приведены на рис. 2.
§8. Заключение
В статье построены метод и, на его основе, алгоритм нахождения оптимального решения задачи маршрутизации, ориентированной на инженерные приложения. В основе подхода находится широко понимамое динамическое программирование (ДП); по постановке допускаются ограничения в виде условий предшествования и функции стоимости, зависящие от списка заданий (здесь имеется в виду список заданий, не выполненных на текущий момент). Объектами посещения являются непустые конечные множества — мегаполисы. При построении расширения основной задачи (для последующего использования аппарата ДП) допустимость маршрутов по предшествованию подменяется допустимостью в смысле вычеркивания заданий из списка. В задаче с условиями предшествования построение всего массива значений функции Беллмана подменяется построением системы ее слоев, что приводит к увеличению размерности задач, допускающих эффективное решение. Построенный оптимальный алгоритм реализован на ПЭВМ.
Финансирование. Работа выполнена в рамках исследований, проводимых в Уральском математическом центре при финансовой поддержке Министерства науки и высшего образования Российской Федерации (номер соглашения 075-02-2021-1383).
СПИСОК ЛИТЕРАТУРЫ
1. Петунин А. А. О некоторых стратегиях формирования маршрута инструмента при разработке управляющих программ для машин термической резки материала // Вестник Уфимского государственного авиационного технического университета. 2009. Т. 13. № 2 (35). С. 280-286. https://elibrary.ru/item.asp?id=15134316
2. Chentsov A. G., Chentsov P. A., Petunin A. A., Sesekin A. N. Model of megalopolises in the tool path optimisation for CNC plate cutting machines // International Journal of Production Research. 2018. Vol. 56. Issue 14. P. 4819-4830. https://doi.org/10.1080/00207543.2017.1421784
3. Ченцов А. Г., Ченцов А. А., Сесекин А. Н. Задачи маршрутизации перемещений с неаддитивным агрегированием затрат. М.: URSS, 2020.
4. Коробкин В. В., Сесекин А. Н., Ташлыков О. Л., Ченцов А. Г. Методы маршрутизации и их приложения в задачах повышения безопасности и эффективности эксплуатации атомных станций. М.: Новые технологии, 2002.
5. Ченцов А. Г., Ченцов А. А., Сесекин А. Н. О задаче последовательного обхода мегаполисов с условиями предшествования и функциями стоимости с зависимостью от списка заданий // Труды Института математики и механики УрО РАН. 2020. Т. 26. № 3. С. 219-234. https://doi.org/10.21538/0134-4889-2020-26-3-219-234
6. Ченцов А. Г. Экстремальные задачи маршрутизации и распределения заданий: вопросы теории. Ижевск: Институт компьютерных исследований, 2008.
7. Gutin G., Punnen A. P. The traveling salesman problem and its variations. Boston: Springer, 2007. https://doi.org/10.1007/b101971
8. Cook W. J. In pursuit of the traveling salesman: mathematics at the limits of computation. Princeton: Princeton University, 2012.
9. Гимади Э. Х., Хачай М. Ю. Экстремальные задачи на множествах перестановок. Екатеринбург: УМЦ УПИ, 2016.
10. Литл Дж., Мурти К., Суини Д., Кэрел К. Алгоритм для решения задачи о коммивояжере // Экономика и математические методы. 1965. Т. 1. Вып. 1. С. 94-107.
11. Беллман Р. Применение динамического программирования к задаче о коммивояжере // Кибернетический сборник. М.: Мир, 1964. Вып. 9. С. 219-228.
12. Хелд М., Карп Р. М. Применение динамического программирования к задачам упорядочения // Кибернетический сборник. М.: Мир, 1964. Вып. 9. С. 202-218.
13. Меламед И. И., Сергеев С. И., Сигал И. Х. Задача коммивояжера. Вопросы теории // Автоматика и телемеханика. 1989. Вып. 9. С. 3-33. http://mi.mathnet.ru/at6414
14. Меламед И. И., Сергеев С. И., Сигал И. Х. Задача коммивояжера. Точные алгоритмы // Автоматика и телемеханика. 1989. Вып. 10. С. 3-29. http://mi.mathnet.ru/at6433
15. Меламед И. И., Сергеев С. И., Сигал И. Х. Задача коммивояжера. Приближенные алгоритмы // Автоматика и телемеханика. 1989. Вып. 11. С. 3-26. http://mi.mathnet.ru/at6463
16. Куратовский К., Мостовский А. Теория множеств. М.: Мир, 1970.
17. Дьедонне Ж. Основы современного анализа. М.: Мир, 1964.
18. Кормен Т., Лейзерсон Ч., Ривест Р. Алгоритмы: построение и анализ. М.: МЦНМО, 2002.
19. Ченцов А. Г. К вопросу о маршрутизации комплексов работ // Вестник Удмуртского университета. Математика. Механика. Компьютерные науки. 2013. Вып. 1. C. 59-82. https://doi.org/10.20537/vm130107
20. Чеблоков И. Б., Ченцов А. Г. Об одной задаче маршрутизации с внутренними работами // Вестник Удмуртского университета. Математика. Механика. Компьютерные науки. 2012. Вып. 1. C. 96-119. https://doi.org/10.20537/vm120109
21. Ченцов А. Г. Одна параллельная процедура построения функции Беллмана в обобщенной задаче курьера с внутренними работами // Вестник Южно-Уральского государственного университета. Сер. «Математическое моделирование и программирование». 2012. Вып. 12. C. 53-76. http://mi.mathnet.ru/vyuru57
22. Chentsov A.A., Chentsov A.G., Sesekin A.N. An extremal routing problem with constraints and complicated cost functions // Advanced Control Systems: Theory and Applications. River Publishers, 2021. P. 21-52.
23. Lawler E. L. Efficient implementation of dynamic programming algorithms for sequencing problems. Stichting Mathematisch Centrum, 1979.
24. Ченцов А. Г., Ченцов А. А. К вопросу о нахождении значения маршрутной задачи с ограничениями // Проблемы управления и информатики. 2016. № 1. С. 41-54.
25. Петунин А. А., Ченцов А. Г., Ченцов П. А. Оптимальная маршрутизация инструмента машин фигурной листовой резки с числовым программным управлением. Математические модели и алгоритмы. Екатеринбург: Изд-во Урал. ун-та, 2020.
26. Ченцов А. Г., Ченцов А. А. Модельный вариант задачи о последовательной утилизации источников излучения (итерации на основе оптимизирующих вставок) // Известия Института математики и информатики Удмуртского государственного университета. 2017. Т. 50. C. 83-109. https://doi.org/10.20537/2226-3594-2017-50-08
Поступила в редакцию 05.07.2021
Ченцов Александр Георгиевич, д. ф.-м. н., член-корреспондент РАН, главный научный сотрудник, отдел управляемых систем, Институт математики и механики им. Н. Н. Красовского УрО РАН, 620108, Россия, г. Екатеринбург, ул. С. Ковалевской, 16;
профессор, кафедра информационных технологий и систем управления, Уральский федеральный университет, 620002, Россия, г. Екатеринбург, ул. Мира, 19. ORCID: https://orcid.org/0000-0001-6568-0703 E-mail: chentsov@imm.uran.ru
Ченцов Алексей Александрович, к. ф.-м. н., научный сотрудник, Институт математики и механики им. Н. Н. Красовского УрО РАН, 620108, Россия, г. Екатеринбург, ул. С. Ковалевской, 16. ORCID: https://orcid.org/0000-0002-0646-9147 E-mail: chentsov.a@binsys.ru
Сесекин Александр Николаевич, д. ф.-м. н., профессор, заведующий кафедрой прикладной математики и механики, Уральский федеральный университет, 620002, Россия, г. Екатеринбург, ул. Мира, 19; ведущий научный сотрудник, Институт математики и механики им. Н. Н. Красовского УрО РАН, 620108, Россия, г. Екатеринбург, ул. С. Ковалевской, 16. ORCID: https://orcid.org/0000-0002-1339-9044 E-mail: a.n.sesekin@urfu.ru
Цитирование: А. Г. Ченцов, А. A. Ченцов, А. Н. Сесекин. Одна задача маршрутизации работ в условиях повышенной радиации // Известия Института математики и информатики Удмуртского государственного университета. 2021. Т. 58. С. 94-126.
A. G. Chentsov, A. A. Chentsov, A. N. Sesekin
One task of routing jobs in high radiation conditions
Keywords: dynamic programming, dismantling, route.
MSC2020: 49L20, 90C39
DOI: 10.35634/2226-3594-2021-58-06
The problem of sequential bypass of megalopolises is investigated, focused on the problem of dismantling a system of radiation hazardous objects under constraints in the form of precedence conditions. The radiation impact on the performers is assessed by the doses received during movements and during the performance of dismantling works. The route problem of minimizing the dose load of workers carrying out dismantling in one or another sequence of operations is considered. The procedure for constructing an optimal solution using a variant of dynamic programming is investigated. On this basis, an algorithm is built, implemented on a PC. Examples of the numerical solution of a model problem for the minimum dose load are given.
Funding. The work was performed as a part of the research conducted in the Ural Mathematical Center with the financial support of the Ministry of Science and Higher Education of the Russian Federation (Agreement number 075-02-2021-1383).
REFERENCES
1. Petunin A. A. About some strategies of the programming of tool route by developing of control programs for thermal cutting machines, Vestnik Ufimskogo Gosudarstvennogo Aviatsionnogo Tekhnich-eskogo Universiteta, 2009, vol. 13, no. 2 (35), pp. 280-286 (in Russian). https://elibrary.ru/item.asp?id=15134316
2. Chentsov A. G., Chentsov P. A., Petunin A. A., Sesekin A. N. Model of megalopolises in the tool path optimisation for CNC plate cutting machines, International Journal of Production Research, 2018, vol. 56, issue 14, pp. 4819-4830. https://doi.org/10.1080/00207543.2017.1421784
3. Chentsov A.G., Chentsov A. A., Sesekin A.N. Zadachi marshrutizatsii peremeshchenii s neaddi-tivnym agregirovaniem zatrat (Move routing problems with non-additive cost aggregation), Moscow: URSS, 2020.
4. Korobkin V.V., Sesekin A.N., Tashlykov O.L., Chentsov A. G. Metody marshrutizatsii i ikh prilo-zheniya v zadachakh povysheniya bezopasnosti i effektivnosti ekspluatatsii atomnykh stantsii (Routing methods and their applications in improving the safety and efficiency of nuclear power plants operation), Moscow: Novye Tekhnologii, 2002.
5. Chentsov A.G., Chentsov A. A., Sesekin A.N. On the problem of sequential traversal of megalopolises with precedence conditions and cost functions depending on a list of tasks, Trudy Instituta Matematiki i Mekhaniki UrO RAN, 2020, vol. 26, no. 3, pp. 219-234. https://doi.org/10.21538/0134-4889-2020-26-3-219-234
6. Chentsov A. G. Ekstremal'nye zadachi marshrutizatsii: voprosy teorii (Extreme problems of routing and distribution of tasks: theoretical questions), Izhevsk: Institute of Computer Science, 2008.
7. Gutin G., Punnen A. P. The traveling salesman problem and its variations, Boston: Springer, 2007. https://doi.org/10.1007/b101971
8. Cook W. J. In pursuit of the traveling salesman: mathematics at the limits of computation, Princeton: Princeton University, 2012.
9. Gimadi E. Kh., Khachai M.Yu. Ekstremal'nye zadachi na mnozhestvakh perestanovok (Extreme problems on sets of permutations), Yekaterinburg: UMC UPI, 2016.
10. Little G. D. C., Murty K. G., Sweeney D. W., Karel C. An algorithm for the traveling salesman problem, Operations Research, 1963, vol. 11, issue 6, pp. 972-989. https://doi.org/10.1287/opre.1L6.972
11. Bellman R. Dynamic programming treatment of the travelling salesman problem, Journal of the ACM, 1962, vol. 9, issue 1, pp. 61-63. https://doi.org/10.1145/321105.321111
12. 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. https://doi.org/10.1137/0110015
13. Melamed 1.1., Sergeev S. I., Sigal I.Kh. The traveling salesman problem. I: Theoretical issues, Automation and Remote Control, 1989, vol. 50, no. 9, pp. 1147-1173. https://zbmath.org/?q=an:0705.90070
14. Melamed 1.1., Sergeev S. I., Sigal I. Kh. The traveling salesman problem. II: Exact methods, Automation and Remote Control, 1989, vol. 50, no. 10, pp. 1303-1324. https://zbmath.org/?q=an:0705.90071
15. Melamed I. I., Sergeev S. I., Sigal I. Kh. The traveling salesman problem. Approximate algorithms, Automation and Remote Control, 1989, vol. 50, no. 11, pp. 1459-1479. https://zbmath.org/?q=an:0704.90095
16. Kuratowski K., Mostowski A. Set theory, Amsterdam: North-Holland, 1967.
17. Dieudonne J. Foundations of modern analysis, New York: Academic Press, 1960.
18. Cormen T. H., Leiserson C. E., Rivest R. L. Introduction to algorithms, Cambridge, MA: MIT Press, 1990.
19. Chentsov A. G. To question of routing of works complexes, Vestnik Udmurtskogo Universiteta. Matematika. Mekhanika. Komp'yuternye Nauki, 2013, issue 1, pp. 59-82 (in Russian). https://doi.org/10.20537/vm130107
20. Cheblokov I. B., Chentsov A. G. About one route problem with interior works, Vestnik Udmurtskogo Universiteta. Matematika. Mekhanika. Komp'yuternye Nauki, 2012, issue 1, pp. 96-119 (in Russian). https://doi.org/10.20537/vm120109
21. Chentsov A. G. A parallel procedure of constructing Bellman function in the generalized courier problem with interior works, Vestnik Yuzhno-Ural'skogo Universiteta. Ser. Matematicheskoe Modeli-rovanie i Programmirovanie, 2012, issue 12, pp. 53-76. http://mi.mathnet.ru/eng/vyuru57
22. Chentsov A.A., Chentsov A.G., Sesekin A.N. An extremal routing problem with constraints and complicated cost functions, Advanced Control Systems: Theory and Applications, River Publishers, 2021, pp. 21-52.
23. Lawler E. L. Efficient implementation of dynamic programming algorithms for sequencing problems, Stichting Mathematisch Centrum, 1979.
24. Chentsov A. G., Chentsov A. A. On the question of finding the value of the route problem with restrictions, Problemy upravleniya i informatiki, 2016, issue 1, pp. 41-54 (in Russian).
25. Petunin A. A., Chentsov A. G., Chentsov P. A. Optimal'naya marshrutizatsiya instrumenta mashin figurnoi listovoi rezki s chislovym programmnym upravleniem. Matematicheskie modeli i algoritmy (Optimal tool routing of CNC shaped sheet cutting machines. Mathematical models and algorithms), Yekaterinburg: Ural Federal University, 2020.
26. Chentsov A. G., Chentsov A. A. A model variant of the problem about radiation sources utilization (iterations based on optimization insertions), Izvestiya Instituta Matematiki i Informatiki Udmurtskogo Gosudarstvennogo Universiteta, 2017, vol. 50, pp. 83-109 (in Russian). https://doi.org/10.20537/2226-3594-2017-50-08
Received 05.07.2021
Aleksandr Georgievich Chentsov, Doctor of Physics and Mathematics, Corresponding Member, Russian Academy of Science, Chief Researcher, N. N. Krasovskii Institute of Mathematics and Mechanics, Ural Branch of the Russian Academy of Sciences, ul. S. Kovalevskoi, 16, Yekaterinburg, 620108, Russia; Professor, Ural Federal University, ul. Mira, 19, Yekaterinburg, 620002, Russia. ORCID: https://orcid.org/0000-0001-6568-0703 E-mail: chentsov@imm.uran.ru
Aleksei Aleksandrovich Chentsov, Candidate of Physics and Mathematics, Researcher, N. N. Krasovskii Institute of Mathematics and Mechanics, Ural Branch of the Russian Academy of Sciences, ul. S. Ko-valevskoi, 16, Yekaterinburg, 620108, Russia. ORCID: https://orcid.org/0000-0002-0646-9147 E-mail: chentsov.a@binsys.ru
Aleksandr Nikolaevich Sesekin, Doctor of Physics and Mathematics, Professor, Head of Department of Applied Mathematics and Mechanics, Ural Federal University, ul. Mira, 19, Yekaterinburg, 620002, Russia;
Leading Researcher, N. N. Krasovskii Institute of Mathematics and Mechanics, Ural Branch of the Russian Academy of Sciences, ul. S. Kovalevskoi, 16, Yekaterinburg, 620108, Russia. ORCID: https://orcid.org/0000-0002-1339-9044 E-mail: a.n.sesekin@urfu.ru
Citation: A.G. Chentsov, A.A. Chentsov, A.N. Sesekin. One task of routing jobs in high radiation conditions, Izvestiya Instituta Matematiki i Informatiki Udmurtskogo Gosudarstvennogo Universiteta, 2021, vol. 58, pp. 94-126.
