УДК 519.87 Вестник СПбГУ. Прикладная математика. Информатика... 2018. Т. 14. Вып. 2
MSC 90B85, 15A80, 65K05, 90C48
Прямое решение минимаксной задачи размещения в прямоугольной области на плоскости с прямоугольной метрикой*
П. В. Плотников, Н. К. Кривулин
Санкт-Петербургский государственный университет, Российская Федерация, 199034, Санкт-Петербург, Университетская наб., 7—9
Для цитирования: Плотников П. В., Кривулин Н. К. Прямое решение минимаксной задачи размещения в прямоугольной области на плоскости с прямоугольной метрикой // Вестник Санкт-Петербургского университета. Прикладная математика. Информатика. Процессы управления. 2018. Т. 14. Вып. 2. С. 116-130. https://doi.org/10.21638/11702/spbu10.2018.204
Рассматривается минимаксная задача размещения точечного объекта на плоскости с прямоугольной (манхэттенской) метрикой с ограничениями на допустимую область размещения и предлагается ее прямое аналитическое решение при помощи методов тропической (идемпотентной) математики. При отсутствии ограничений эта задача, которая также известна как задача Ролса или задача посыльного, имеет известные геометрическое и алгебраическое решения. Исследуется решение данной задачи с учетом дополнительных ограничений на область размещения, которая задана в форме прямоугольника. Сначала задача записывается в терминах тропической математики как задача тропической оптимизации, вводится параметр для обозначения минимума целевой функции и задача сводится к решению параметризованной системы неравенств. Эта система решается относительно одной из переменных, а условия существования решений используются для нахождения оптимальных значений другой переменной с помощью вспомогательной задачи оптимизации. Затем полученное общее решение преобразуется в набор прямых решений, записанных в компактной замкнутой форме для различных случаев соотношений между исходными параметрами задачи. Приведены графические примеры решения задачи для разных вариантов расположения допустимой области размещения на плоскости.
Ключевые слова: задача размещения Ролса, размещение с ограничениями, прямоугольная метрика, идемпотентное полуполе, тропическая оптимизация, полное решение.
1. Введение. Одной из перспективных и быстро развивающихся областей теоретической и прикладной математики является тропическая (идемпотентная) математика, которая связана с изучением теории и приложений полуколец с идемпотент-ным сложением [1—9]. В технике, экономике и управлении достаточно часто можно встретить оптимизационные задачи, которые формулируются и эффективно решаются в терминах тропической математики. Краткий обзор таких задач представлен, например, в работе [10]. Если традиционно использование операций max и min при формулировке математических постановок прикладных задач часто усложняет их решение, то на языке идемпотентной алгебры работа с такими операциями упрощается. Вследствие этого методы тропической математики могут оказаться удобными при решении минимаксных задач. Следует отметить важную особенность, заключающуюся в возможности преобразовывать нелинейные задачи в линейные в терминах тропической математики.
* Работа выполнена при финансовой поддержке Российского фонда фундаментальных исследований (грант № 18-010-00723).
(¡5 Санкт-Петербургский государственный университет, 2018
Существует класс задач оптимизации, в которых требуется найти наилучший способ размещения объектов на плоскости с дополнительными ограничениями на допустимую область исследования. Значительную роль при решении этих задач играет выбор метрики, при помощи которой осуществляется подсчет расстояний. Решения для некоторого класса задач размещения с чебышёвской и прямоугольной метриками в терминах тропической математики представлены в [11—13]. Решение в явном виде минимаксной задачи размещения точечного объекта на плоскости с прямоугольной метрикой без ограничений на область исследования изучено в статьях [14, 15]. Такого рода минимаксная задача обычно называется задачей Ролса или задачей посыльного. Хорошо известно [16, 17] геометрическое решение данной задачи. Изучение задачи с ограничением в виде отрезка прямой линии проведено в работе [18].
Настоящее исследование направлено на дальнейшее развитие и применение методов тропической оптимизации для решения задач размещения. Рассматривается минимаксная задача размещения точечного объекта на плоскости с прямоугольной (манхэттенской) метрикой с ограничениями на допустимую область. Такая задача может встречаться, например, при оптимальном размещении объектов экстренной помощи населению (противопожарные службы, службы скорой помощи и др.) в городах с дорожной сетью в виде параллельных и перпендикулярных между собой улиц, при котором минимизируется расстояние от объекта до самого дальнего возможного очага возгорания или самого удаленного места проживания больного. Другие примеры включают оптимальное размещение сервера локальной сети (центрального пульта систем видеонаблюдения, охранной и пожарной сигнализаций и т. п.) в зданиях и сооружениях с прямоугольными и перпендикулярными друг другу маршрутами прокладки проводных линий связи по критерию минимизации потерь при передаче сигнала.
В рассматриваемой задаче ограничения на область размещения представлены в виде произвольного прямоугольника со сторонами, параллельными осям координат (линиям дорожной сети, линиям связи). При решении применяется подход, развитый в работе [10], который заключается в замене исходной задачи оптимизации, записанной в терминах тропической математики, на параметризованную систему неравенств. На основании этого подхода найдено прямое полное решение задачи, представленное в терминах тропической математики, а также в обычной форме. В отличие от известных алгоритмических решений полученный результат описывает все множество решений в явном виде с помощью простых расчетных формул, удобных как для дальнейшего анализа решений, так и для непосредственных вычислений. Это позволяет эффективно решать задачу размещения в случаях, когда алгоритмическое решение оказывается по тем или иным причинам невозможным или нецелесообразным.
2. Задача размещения на плоскости с ограничениями. Рассмотрим минимаксную задачу размещения точечного объекта на плоскости с прямоугольной метрикой, которая заключается в поиске на основе уже имеющегося набора объектов новой точки из Б С К2, расстояние от которой до самого дальнего объекта из набора с учетом некоторого дополнительного слагаемого было бы минимальным.
Расстояние в прямоугольной метрике между двумя точками х = (х\, Х2)т и у = (У1,У2)Т на плоскости К2 может быть вычислено по формуле
Р(х, у) = |Х1 - У11 + |Х2 - У2 | • (1)
Пусть задан набор точек гг = (тц, Г2г)т € К2 и чисел Ьг € К для всех % = 1,^,ш. Числа Ьг могут отражать дополнительные затраты, связанные с перемещением между объектами. Например, в задачах, где точки описывают положение объектов на
горизонтальной плоскости, эти числа могут иметь смысл высоты (вертикальной координаты) объекта. Введем функцию ф(х), которая определяет максимальное расстояние от точки х = (xi ,Х2 )T до набора точек vi следующим образом:
ф(х) = max (р(Vi, x) + hi)= max (|rii - xi| + |r2i - X2I + hi) = ф(х1 ,x2). (2)
l^i^m l^i^m ^ '
Предположим также, что заданы числа f,g,p,q G R такие, что f ^ g и p ^ q. Определим допустимую область размещения как прямоугольник
S = {(xi,X2)T | f < xi < g, p < X2 < q}. (3)
Рассматриваемая минимаксная задача размещения точечного объекта с ограничениями в форме прямоугольника S формулируется так:
mn ф(х). (4)
Для решения такого рода задачи могут быть применены различные подходы. Например, с помощью методов линейного программирования можно получить алгоритмическое итерационное решение задачи, но этот подход не гарантирует прямого решения в явном виде. Для построения такого аналитического решения в п. 6 будет использован аппарат тропической математики.
3. Элементы тропической математики. Введем числовое множество X, на котором определены ассоциативные и коммутативные операции сложения ф и умножения Обозначим через (X,0,1, ф, ®) заданное на X при помощи этих операций коммутативное полукольцо с нулем 0 и единицей 1. Сложение будем считать идем-потентным (т. е. для любого числа x G X выполняется x ф x = x), а умножение обратимым (т. е. для каждого x = 0 существует обратный элемент x-l такой, что x ® x-l = 1). Так как (X \ {©}, 1, образует коммутативную группу по умножению, то описанную структуру (X,0,1, ф, ®) принято называть идемпотентным полуполем.
Операция возведения в степень с целым показателем вводится стандартным образом. Для любого ненулевого числа x G X и натурального числа n определим x0 = 1, xn = x ® xn-i, x-n = (x-l)n и О" = ©. Будем считать, что операция возведения в целую степень ненулевого числа x может быть распространена на случай произвольного показателя степени.
Далее для упрощения математических выкладок знак умножения ® в алгебраических выражениях, как обычно, будет опускаться.
В силу того, что сложение идемпотентно, можно задать отношение частичного порядка ^ так: x ^ y тогда и только тогда, когда x ф y = y. Из этого определения следуют неравенства x ^ xфy и y ^ xфу, а также равносильность неравенства xфy ^ z и системы неравенств x ^ z и y ^ z для любых x,y,z G X. Кроме того, нетрудно проверить свойства монотонности операций сложения и умножения, по которым при условии x ^ y для любого z выполняются неравенства x^z ^ yфz и xz ^ yz. Наконец, для любых x,y = © из неравенства x ^ y следует x-l ^ y-1. В дальнейшем будем дополнительно предполагать, что введенный частичный порядок — линейный. Кроме того, несложно проверить неравенство xl/'2yi/2 ^ xфy, которое является тропическим аналогом неравенства между геометрическим и арифметическим средними.
В качестве примеров алгебраических структур описываемого типа можно привести вещественные идемпотентные полуполя
^max,+ = ^ и {-сю}, -то, 0, max, +}, Rmin,+ = ^ U {+ю}, +ю, 0, min, +}, RmaX,x = (R + U {0}, 0,1, max, x}, Rmin,x = (R + U {+ю}, +ю, 1, min, x},
где R — множество вещественных чисел и R+ = {x G R| x > 0}.
Нулевым элементом для полуполя RmaXj+, которое называют (max, +)-алгеброй, служит —ю, единичным — число 0. В этом полуполе любому числу x G R можно сопоставить обратный x-1, который совпадает с противоположным числом —x в обычной алгебре. Для любой пары чисел x,y G R можно определить степень xy, которая равна арифметическому произведению xy. Порядок, порожденный операцией идемпотент-ного сложения, совпадает с обычным линейным порядком на множестве .
4. Задачи тропической оптимизации с ограничениями. Рассмотрим задачи оптимизации, сформулированные в терминах тропической математики (задачи тропической оптимизации), с одной и двумя переменными.
Сначала исследуется задача оптимизации с одной переменной, решение которой было получено в работе [18]. Будет показано, как это решение может быть записано в другой форме, более удобной для последующего применения. Затем формулируется новая задача тропической оптимизации с двумя переменными.
Предположим, что необходимо найти ненулевые решения t G X задачи
min at-1 ф bt ф c,
f < t < g, (5)
где a, b, c, f, g — числа из X, удовлетворяющие условию a,b,c, g > (D. Решение этой задачи дает следующий результат [18].
Лемма. Минимум в .задаче (5) равен л = a1/2b1/2 ф ag-1 ф bf ф c и достигается тогда и только тогда, когда
t = (л-1 a ф f )1-a(^-1b ф g-1)-a, 0 < а < 1. (6)
Преобразуем утверждение леммы к другому, более удобному для дальнейшего использования виду.
Следствие 1. Минимум в задаче (5) равен л = a1/2b1/2 ф ag-1 ф bf ф c и справедливы следующие утверждения:
1) если л = a1/2b1/2, то t = a1/2b-1/2;
2) если л = ag-1, то t = g;
3) если л = bf, то t = f;
4) если л = c, то t = (ac-1 ф f )1-a(bc-1 ф g-1)-a для любого 0 ^ а ^ 1. Доказательство. Представим t из формулы (6) в виде двойного неравенства
л-1a ф f < t < (л-1Ь ф g-1)-1.
Рассмотрим возможные значения л. Пусть л = a1/2b1/2. После подстановки в предыдущую формулу и применения свойств монотонности идемпотентного сложения получим
a1/2b-1/2 < a1/2b-1/2 ф f < t < (a-1/2b1/2 ф g-1)-1 < a1/2b-1/2,
откуда можно заключить, что t = a1/2b-1/2.
Аналогично проверяем, что при л = ag-1 будем иметь t = g, а при л = bf значение t будет равно f.
Если л = c, то сохраним представление для t в параметрической форме. □
Полученные выражения для решения представим в более компактной форме.
Следствие 2. Минимум в задаче (5) равен л = a1/2b1/2 ф ag 1 ф bf ф c и справедливы такие утверждения:
1) если л = a1/2b1/2 ф ag-1 ф bf, то t = (a-1/2b1/2 ф g-1)-1 ф f;
2) если л = c, то t = (ac-1 ф f )1-a(bc-1 ф g-1)-a, 0 < а < 1.
Доказательство. Объединим результаты первых трех пунктов в формулировке следствия 1 в один. Заметим, что общее значение минимума л может быть записано как сумма л = a1/2b1/2 ф ag-1 ф bf. Покажем, что все величины, при которых достигается минимум, вычисляются по формуле t = (a~ -1/2b1/2 ф g-1)-1 ф f.
Предположим, что л = a1/2b1/2. Тогда выполняются неравенства a1/2b1/2 ^ ag-1 и a1/2b1/2 ^ bf, которые после преобразования могут быть записаны так: a-1/2b1/2 ^ g-1 и a1/2b-1/2 > f. Тогда t = (a-1/2b1/2 ф g-1)-1 ф f = a1/2b-1/2, что соответствует результату 1) следствия 1.
Аналогичным образом проверяем, что при л = ag-1 выполняется равенство t = (a-1/2b1/2 фg-1)-1 ф f = g, а при л = bf — равенство t = (a-1/2b1/2 фg-1)-1 ф f = f. □
Теперь рассмотрим новую задачу тропической оптимизации с двумя переменными.
Пусть заданы числа a, b, c, d, f,p, g, q > 0. Требуется найти решения t, s > 0 задачи
min at-1s-1 ф bts-1 ф ct-1s ф dts, f < t < g, p < s < q.
В п. 5 дадим прямое полное решение этой задачи.
5. Решение задачи с двумя переменными. Общая схема решения задачи (7) следующая. Сначала вводится параметр для обозначения минимума целевой функции, и задача сводится к решению параметризованной системы неравенств. Эта система решается относительно одной из переменных, а условия существования решений используются для нахождения оптимальных значений другой переменной с помощью вспомогательной задачи оптимизации. Затем полученное общее решение задачи преобразуется в набор прямых решений, записанных в компактной замкнутой форме для различных случаев соотношений между исходными параметрами задачи.
Все решения рассматриваемой задачи описывает следующий результат.
Теорема 1. Введем обозначения
u = a1/2c1/2 ф aq-1 ф cp, v = b1/2d1/2 ф bq-1 ф dp.
Минимум в .задаче (7) равен л = u1/2v1/2 ф ug-1 ф vf ф a1/2d1/2 ф b1/2c1/2 и справедливы следующие утверждения:
1) если л = u1/2v1/2 при u = a1/2c}/2 и v = bq-1 или при u = aq-1, то
t = (a1/4c1/4 ф a1/2q-1/2)(b1/2q-1/2 ф v1/2)-1, s = q;
2) если л = u1/2v1/2 при u = a1/2c}/2 и v = dp или при u = cp, то
t = (a1/4c1/4 ф c1/2p1/2)(d1/2p1/2 ф v1/2)-1, s = p;
3) если л = ug-1, то t = g и s = (a-1/2c1/2 ф q-1)-1 ф p;
4) если л = vf, то t = f и s = (b-1/2d1/2 ф q-1)-1 ф p;
5) если л = a1/2d1/2, то
t = (a-1/2d-1/2u ф f )1-a(a-1/2d-1/2v ф g-1)-a, s = a1/2d-1/2t-1;
6) если л = ь1/2с1/2, то
г = (ь-1/2с-1/2и е /)1-а(ь-1/2с-1/2у е д-1)-а, в = ь1/2с-1/2г,
где а — вещественное число, удовлетворяющее условию 0 ^ а ^ 1.
Доказательство. Обозначим минимум целевой функции в задаче (7) через Тогда все решения задачи определяются системой, состоящей из уравнения и двух неравенств:
аг-1в-1 е Ые-1 е сг-1в е ¿Ьз = л, / < г < д, р < в < д.
В силу того, что л — минимальное значение целевой функции, равенство в первом уравнении можно заменить на неравенство и записать систему так:
аг-1в-1 е Ыв-1 е сг-1в е ¿Ьз < л, / < г < д, р < в < д.
Полученная система неравенств эквивалентна системе
аг-1в-1 ^ л, Ыв-1 ^ л, сЬ-1 в ^ л, сгв ^ л, / ^ г, г ^ д, р ^ в, в ^ д.
Заметим, что перемножение соответствующих частей, например первого и четвертого неравенств, дает аС ^ /л2, откуда с учетом условия леммы следует, что л ^ а1'2С1'2 > ®.
Решим систему относительно неизвестного в, считая г параметром. Решения неравенств аг-1в-1 ^ л и Ыв-1 ^ л относительно в имеют вид в ^ л-1 аг-1 и в ^ л-1 Ьг. Объединив их с условием в ^ р, получим в ^ л-1 аг-1 е л-1Ьг ер.
Рассмотрим неравенства сг-1в ^ л и ¿гв ^ л, которые решим относительно в-1, чтобы получить в-1 ^ л-1сЪ-1 и в-1 ^ л-1СЬ. Найденные решения вместе с условием в-1 ^ д-1 приводят к неравенству в-1 ^ л-1сЪ-1 е л-1 ¿г е д-1 или к эквивалентному ему неравенству в ^ (л-1 сЬ-1 е л-1СЬ е д-1)-1.
Записывая вместе неравенства для в и граничные условия на г, получим систему
л-1аг-1 е л-1ьг е р < в < (л-1сг-1 е л-1сг е д-1)-1, / < г < д. (8)
Множество значений в, удовлетворяющих первому неравенству, непусто, если выполняется условие л-1аг-1 е л-1ЬЬ е р ^ (л-1сЬ-1 е л-1СЬ е д-1 )-1, которое равносильно условию
(л-1аг-1 е л-1ьг ер)(л-1сг-1 е л-1сг е д-1) < 1.
Раскроем скобки в левой части, а затем заменим это неравенство эквивалентной системой неравенств, включая рд-1 ^ 1, которое прямо следует из условия р ^ д. Решение остальных неравенств системы относительно л дает
л > а1/2с1/2г-1, л > а1/2С1/2, л > ь1/2с1/2, л > ь1/2С1/2г, л ^ ад-1г-1, л ^ ¿рг, л ^ срг-1, л ^ ьд-1г.
Построенная система с учетом граничных условий на г равносильна паре неравенств
л > а1/2с112 е а112с}/2г-1 е ь1/2с1/2г е ь1/2с1/2 е ад-1г-1 е ¿рг е срг-1 е ьд-1г,
/ < t < д.
(9)
Заметим, что переход от исходной задачи к полученной системе, в которой л обозначает минимум целевой функции, включал только эквивалентные преобразования. Из этого вытекает, что первое неравенство системы задает точную нижнюю границу для л, выраженную через t, а потому теперь необходимо решить задачу оптимизации следующего вида:
min (a1/2c1/2 ф aq-1 ф cp)t-^(b1/2d1/2 ф bq-1 ф dp)t ф a1/2d1/2 ф b1/2c1/2,
f < t < g.
Введем такие обозначения:
u = al/2c}/2 ф aq-1 ф cp, v = b1/2d1/2 ф bq-1 ф dp, w = a1/2d1/2 ф b1/2c1/2. (10)
Задача (9) может быть записана в форме (5), где a заменяется на u, b на v и c на w. Применение вышеизложенной леммы для решения задачи дает минимальное значение
л = u1/2v1/2 ф ug-1 ф vf ф w,
причем, согласно следствию 1, выполняются следующие условия:
1) если л = u1/2v1/2, то t = u1/2v-1/2;
2) если л = ug-1, то t = g;
3) если л = vf, то t = f;
4) если л = w , то t = (w-1 u ф f )1-a(w-1v ф g-1)-a, 0 < а < 1.
Теперь уточним полученное решение, рассматривая различные значения, которые может принимать величина л.
5.1. Случай л = u1/2v1/2. В этом случае t = u1/2v-1/2, а двойное неравенство (8) для s можно записать в виде
u-1a ф v-1b ф p < s < (u-1c ф v-1 d ф q-1)-1.
Рассмотрим все значения, которые могут принимать u и v в соответствии с формулами (10). Пусть u = ai/2ci/2, тогда л = a1/Ac1/Av1/2, а t = a1/Ac1/Av-1/2. Проверим, какие значения может принимать v, и найдем соответствующие величины t и s. Исследуем случай, когда выполняется равенство v =
b1/2d1/2
. В силу тропического аналога неравенства между средним арифметическим и средним геометрическим выполняется неравенство
л = u1/2v1/2 = a1/Ab1/Ac1/Ad1/A < a1/2d1/2 ф b1/2c1/2 = w.
С учетом условия л = u1/2v1/2 ^ w оно приводит к равенству л = w, а потому решение при v =
b1/2d1/2
можно рассматривать как частный случай для общего случая л = w, который будет изучен ниже.
Предположим, что v = bq-1. Тогда выполняются равенства л = a1/Ab1/2c1/Aq-1/2 и t = a1/Ab-1/2c1/Aq1/2. Из двойного неравенства для s получим
q < a1/2c-1/2 ф q ф p < s < (a-1/2c1/2 ф b-1qd ф q-1)-1 < q,
откуда следует, что s = q.
Если v = dp, то л = a1/Ac1/Ad1/2p1/2 и t = a1/Ac1/Ad-1/2p-1/2, а неравенство для s приводит к результату
p < a-1/2c-1/2a ф d-1p-1b ф p < s < (a-1/2c-1/2c ф p-1 ф q-1)-1 < p,
который показывает, что s = p.
Пусть и = ад-1. Тогда л = а1/2д-1/2у1/2 и £ = а1/2д-1/2у-1/2. Кроме того, д < а-1 да 0 V-1 Ь 0 р < в < (а-1 до 0 у-1с! 0 д-1)-1 < д,
а значит, в = д.
При и = ор выполняются л = о1'2p1'2v'1'2 и £ = оí'2pí'2v-1'2. Имеем неравенство
р < о-1р-1а 0 v-1Ь 0 р < в < (р-1 0 v-1d 0 д-1)-1 < р,
из которого вытекает, что в = р.
Записывая вместе предыдущие результаты для л = u1'2v1'2, находим, что
• при и = а1'2о1/2 и v = Ьд-1 выполняются £ = а1'4Ь-1'2о1'4д1'2 и в = д;
• при и = ад-1 выполняются £ = a1'2д-1'2v-1'2 и в = д;
• при и = а1'2о1'2 и v = dp выполняются £ = a1'4о1'4d-1'2р-1'2 и в = р;
• при и = ор выполняются £ = оí'2pí'2v-1'2 и в = р.
Покажем, что выражения для которые соответствуют значениям и = а1'2о1'2 и v = Ьд-1 или значению и = ад-1, можно представить с помощью одной формулы: г = (а1'4о1'4 0 а1'2д-1/2)(Ь1/2д-1/2 0 v1/2)-1.
Предположим, что выполнены условия и = а1'2о1'2 и v = Ьд-1. Из первого условия следует, что и = а1'2о1'2 ^ ад-1, а значит а1'4о1'4 ^ а1'2д-1'2. Учитывая, что v1'2 = Ь1'2д-1'2, получим I = (а1'4о1'4 0 а1'2д-1'2)(Ь1'2д-1'2 0 v1'2)-1 = а1'4Ь-1'2о1'4д1'2.
В случае, если и = ад-1, выполняется неравенство ад-1 ^ а1'2о1'2, из которого вытекает, что а1'2д-1'2 ^ а1'4о1'4. Кроме того, всегда выполняется неравенство v ^ Ьд-1. Следовательно, I = (а1'4о1'4 0 а1'2д-1'2)(Ь1'2д-1'2 0 v1'2)-1 = а1'2д-1'^-1'2.
Аналогичным образом проверяем, что при и = а1'2о1'2 и v = dp и при и = ор справедливо общее выражение для £ в виде £ = (а1'4о1'4 0 о}'2р1'2)(ё1'2р1'2 0 ví'2)-1. 5.2. Случай л = ид-1. В этом случае Ь = д, а неравенство (8) принимает форму
1
а 0 и 1Ьд2 0 р < в < (и 1о 0 и 1dg2 0 д
и а 0 и Ьд 0 р ^ в ^ (и
Найдем соответствующие значения для в и л при различных величинах и. Предположим, что и = а1'2о1'2. Тогда л = а1'2о1'2д-1. Учитывая, что £ = д, из двойного неравенства для в получаем
а1'2о-1'2 < а1'2о-1'2 0 а-1'2Ьо-1'2д2 0 р < в < < (а-1'2о1'2 0 а-1'2о-1'Чд2 0 д-1)-1 < а1'2о-1'2,
откуда вытекает равенство в = а1'2о-1'2.
Аналогичным образом находим, что при и = ад-1 выполняются равенства л = ад-1д-1 и в = д, а при и = ор — равенства л = орд-1 и в = р.
Нетрудно проверить, что все значения в, при которых достигается минимум в этом случае, вычисляются по формуле в = (а-1'2 о1'2 0 д ) 0 р.
Пусть, например, и = а1'2о1'2. В этом случае выполняются условия а1'2о1'2 ^ ад-1 и а1'2о1'2 ^ ор. Они эквивалентны неравенствам а-1'2о1'2 ^ д-1 и а1'2о-1'2 ^ р, откуда следует, что в = (а-1'2о1'2 0 д-1)-1 0 р = а1'2о-1'2.
Случаи, когда и = ад-1 и и = ор, проверяются аналогично.
5.3. Случай л = vf. Учитывая, что здесь £ = ], неравенство (8) имеет вид
v-1af-2 0 v-1 Ь 0 р < в < (v-1оf-2 0 v-1d 0 д-1)-1.
Если V = то л = /. Из двойного неравенства для в вытекает, что
ь1'2вг1'2 < v-1a/-2 е ь1'2вг1'2 ер < в < (V-1 с/-2 е ь-1'2^'2 е д-1)-1 < ь1/2а-1/2,
а тогда выполняется равенство в = ь1'2Л~'1'2.
Аналогичным образом можно проверить, что при V = ьд-1 выполняются равенства л = ьд-1/ и в = д, а при V = dр — равенства л = ¿р/ и в = р. Все значения в, при которых достигается минимум в этом случае, вычисляются по формуле в = (ь-1'2ё}'2 е д-1)-1 е р. Если, например, V = ь1/2¿1/2, то справедливы неравенства
ь1'2с11'2
> ьд 1 и ь1'2й1'2 > ¿р, которые равносильны условиям ь 1'2й1'2 > д 1 и ь1'2в,-1'2 > р. Это означает, что в = (ь-1'2^'2 е д-1)-1 е р = ь1'2й-1'2. Аналогично можно проверить случаи V = ьд-1 и V = ¿р.
5.4. Случай л = В этом случае выражение для г включает параметр а. Покажем, что при этом величина в может принимать только два значения. Пусть сначала л = а1'2^1'2. Тогда неравенство (8) для в имеет вид
а1'2а-1'2г-1 е а-1'2ьа-1'Ч ер < в < (а,-1'2ся-1'2г-1 е а-1'2а1'2г е д-1 )-1. Для выражений слева и справа выполняются неравенства
а1'2а-1'21-1 < а1'2а-1'2г-1 е а-1'2ьа-1'Ч ер, (а,-1'2сЯ-1'Ч-1 е а-1'2а1'21 е д-1)-1 < (а-1'2¿1'2г)-1 = а1'2¿-1'Ч-1, откуда можно сделать вывод, что в = а1'23-1'2Ь-1, где
г = (а-1'2с!-1'2и е /)1-a(a-1'2d-1'2v е д-1)-а, 0 < а < 1. При л = Ъ1'2с1'2 таким же образом получим, что в = ь1'2с-1'2г, где
г = (ь-1'2с-1'2и е /)1-а(ь-1'2с-1'\ е д-1)-а, 0 < а < 1.
Исследование этого случая завершает доказательство теоремы. □
Сформулируем полученный результат в более компактной форме. Следствие 3. Введем обозначения
и = а1'2с1'2 е ад-1 е ср, V = Ъ1'2с}1'2 е ьд-1 е ¿р, " = а1'2сI1'2 е ь1'2с1'2.
Минимум в задаче (7) равен л = u1'2v1'2 е ид-1 е V/ е " и достигается тогда и только тогда, когда
г = (л-1и е /)1-а(л-^ е д-1)-а, (11)
в = (л-1аг-1 е л-1ьг е р)1-а(л-1 сг-1 е л-^г е д-1)-а, (12)
где а — вещественное число, удовлетворяющее условию 0 ^ а ^ 1.
Доказательство. При доказательстве теоремы 1 было показано, что множество решений (г, в) задачи (7) описывается при помощи двойного неравенства
л-1аг-1 е л-1ьг е р < в < (л-1сг-1 е л-1dг е д-1)-1,
где величина г является решением задачи (9), которое может быть получено с помощью леммы в параметрической форме (11).
Ясно, что двойное неравенство для в также может быть записано как равенство
с использованием некоторого параметра. Заметим, что в каждом из случаев решения, исследованных в теореме 1, по крайней мере одна из величин t и s определяется однозначно без помощи параметров. В силу этого параметрическое представление для s можно также записать с использованием параметра а в виде (12). □
6. Приложение к задаче размещения. Рассмотрим задачу размещения (4) и представим ее в терминах тропической математики. Сначала, применив операции идемпотентного полуполя IRmax,+, определим расстояние (1) между двумя векторами:
p(x, y) = (x-xyi © y-lxi)(x-1y2 © у-гх2).
Для заданного набора точек v, = (r1i, r2i)T € IR2 и чисел hi € IR, где i = l,...,m, целевую функцию (2) запишем следующим образом:
m
Ф(х1 ,Х2) = ф hi(x-1 rii © r-i1x1 )(x-1T2i © r2i1X2). i=l
Чтобы упростить полученное выражение, введем дополнительные обозначения
m m m m
a = ф hirur2i, b = ф hir-/ r2i, c = ф Ыгцг-^, d = ф hr-г-,1. i=1 i=1 i=1 i=1
Раскрыв скобки и перегруппировав слагаемые, запишем целевую функцию так:
Ф(^1, x2) = ax-1x-1 © bx^x-1 © cx-1 x2 © dx1 x2.
С учетом ограничений (3) задача (4) принимает форму задачи (7), где t = x1 и s = x2. Применение теоремы 1 дает следующий результат. Теорема 2. Введем обозначения
u = a1/2c1/2 © aq-1 © cp, v = b1/2d1/2 © bq-1 © dp.
Минимум в задаче (4) равен л = u1/2v1/2©ug-1 ©vf ©a1/2d1/2©b1/2c1/2 и справедливы следующие утверждения:
1) если л = u1/2v1/2 при u = a1/2c}/2 и v = bq-1 или при u = aq-1, то
x1 = (a1/4c1/4 © a1/2q-1/2)(b1/2q-1/2 © v1/2)-1, x2 = q;
2) если л = u1/2v1/2 при u = a1/2c}/2 и v = dp или при u = cp, то
x1 = (a1/4c1/4 © c1/2p1/2)(d1/2p1/2 © v1/2)-1, x2 = p;
3) если л = ug-1, то x1 = g и x2 = (a-1/2c1/2 © q-1)-1 © p;
4) если л = vf, то x1 = f и x2 = (b-1/2d1/2 © q-1)-1 © p;
5) если л = a1/2d1/2, то
x1 = (a-1/2d-1/2u © f )1-a(a-1/2d-1/2v © g-1)-a, x2 = a1/2d-1/2x-1;
6) если л = b1/2c1/2, то
x1 = (b-1/2c-1/2u © f )1-a(b-1/2c-1/2v © g-1)-a, x2 = b1/2c-1/2x1,
где а — вещественное число, удовлетворяющее условию 0 ^ а ^ 1.
Опишем найденное решение с использованием обычных обозначений.
Следствие 4. Введем обозначения
a = max (hi + ru + r^i), b = max (hi - тц + r^i),
i^i^m i^i^m
c = max (hi + ru - r2i), d = max (hi - ru - r2i),
l^i^m l^i^m
u = max((a + c)/2, a - q,c + p), v = max((b + d)/2, b - q,d + p).
Минимум в задаче (4) равен ц = max((u + v)/2, u- g,v + f, (a + d)/2, (b + c)/2) и справедливы следующие утверждения:
1) если ц = (u + v)/2 при u = (a + c)/2 и v = b - q или при u = a - q, то
x1 = max((a + c)/4, (a - q)/2) - max((b - q)/2, v/2), x2 = q;
2) если ц = (u + v)/2 при u = (a + c)/2 и v = d + p или при u = c + p, то
x1 = max((a + c)/4, (c + p)/2) - max((d + p)/2,v/2), x2 = p;
3) если ¡л = u - g, то xi = g и x2 = max(- max((-a + c)/2, -q),p);
4) если ¡л = v + f, то xi = f и x2 = max(- max((-b + d)/2, -q),p);
5) если ц = (a + d)/2, то
xi = (1 - a) max((-a - d)/2 + u, f) - a max((-a - d)/2 + v, -g), x2 = (a - d)/2 - xi;
6) если ц = (b + c)/2, то
xi = (1 - a) max((-b - c)/2 + u, f) - a max((-b - c)/2 + v, -g), x2 = (b - c)/2 + xi,
где a — вещественное число, удовлетворяющее условию 0 ^ a ^ 1.
Заметим, что вычислительная сложность решения задачи размещения в соответствии со следствием 4 является линейной по числу точек т.
Приведем примеры применения полученных результатов. Пусть имеется n =11 исходных точек на плоскости с координатами
ri = (3, 7)T, Г2 = (4, 4)T, r3 = (5, 3)T, Г4 = (5, 9)T, r5 = (6, 6)T, r6 = (7, 3)T, r7 = (7, 8)T, r8 = (8, 4)T, r9 = (9, 2)T, rw = (10, 7)T, щ = (11, 4)T.
Для простоты будем считать, что hi =0 для всех i. Зададим четыре различных области размещения:
51 = {(xi,x2)T | 4 < xi < 10, 3 < x2 < 9},
52 = {(xi,x2)T | 2 < xi < 12, 1 < x2 < 5.5},
53 = {(xi,x2)T I 2 < xi < 5, 1 < x2 < 11},
54 = {(xi,x2)T | 2 < xi < 12, 8 < x2 < 11}.
В результате применения расчетных формул из следствия 4 для задач с областями размещения Si, S2, S3 и S4 получим решения, в которых оптимальное размещение описывается следующими парами координат:
Х1 = а + 6.5, Х1 = 0.5а + 6.5, Х1 = 5, Х1 = 7.5,
Х2 = а + 5, Х2 = 0.5а + 7, Х2 = 5, Х2 = 8,
где параметр а удовлетворяет условию 0 ^ а ^ 1.
На рисунке, а, б, в, г представлены результаты решения задач размещения в прямоугольных областях. Кружками обозначены элементы исходного множества объектов, а пунктирным прямоугольником — допустимая область размещения. Множество точек оптимального размещения изображено отрезком жирной линии или жирной точкой.
*2
/
О о , о
- - о----О........
о
4 6 8 10 12 в
, о
о ! о
о о
о 2 4 6 8 10 12
10 8 6 4
2
10 8 6 4
2
4 6 8 10 12
г
от
о о о
о о
о 2 4 6
10 12 X,
Рисунок. Оптимальное размещение в областях Я1 (а), Я2 (б), Я3 (в) и Я4 (г)
6. Заключение. В работе рассмотрена задача размещения в прямоугольной области на плоскости с прямоугольной метрикой. На основе применения методов тропической оптимизации получено прямое полное решение задачи, представленное в форме, которая является удобной для дальнейшего анализа решений и пригодной для непосредственных расчетов с невысокой вычислительной сложностью.
Дальнейшие исследования могут включать решение задач размещения на плоскости с допустимой областью, отличной от прямоугольной, а также задач размещения в трехмерном пространстве.
Литература
1. Baccelli F., Cohen G., Olsder G. J., Quadrat J.-P. Synchronization and linearity. Chichester: Wiley, 1993. 514 p. (Wiley Series in Probability and Statistics)
2. Маслов В. П., Колокольцев В. Н. Идемпотентный анализ и его применение в оптимальном управлении. М.: Физматлит, 1994. 144 с.
3. Cuninghame-Green R. A. Minimax algebra and applications // Advances in Imaging and Electron Physics. Vol. 90 / ed. by P. W. Hawkes. San Diego: Academic Press, 1994. P. 1-121. DOI: 10.1016/S1076-5670(08)70083-1
4. Golan J. S. Semirings and affine equations over them. New York: Springer, 2003. 256 p. (Vol. 556 of Mathematics and its Applications). DOI: 10.1007/978-94-017-0383-3
5. Heidergott B., Olsder G. J., van der Woude J. Max plus at work. Princeton: Princeton University Press, 2006. 226 p. (Princeton Series in Applied Mathematics)
6. McEneaney W. M. Max-plus methods for nonlinear control and estimation. Boston: Birkhauser, 2006. 241 p. (Systems and Control: Foundations and Applications). DOI: 10.1007/0-8176-4453-9
7. Itenberg I., Mikhalkin G., Shustin E. I. Tropical algebraic geometry. Basel: Birkhauser, 2007. 104 p. (Vol. 35 of Oberwolfach Seminars). DOI: 10.1007/978-3-7643-8310-7
8. Gondran M., Minoux M. Graphs, dioids and semirings. New York: Springer, 2008. 383 p. (Vol. 41 of Operations Research/ Computer Science Interfaces). DOI: 10.1007/978-0-387-75450-5
9. Butkovic P. Max-linear systems. London: Springer, 2010. 272 p. (Springer Monographs in Mathematics). DOI: 10.1007/978-1-84996-299-5
10. Krivulin N. K. A multidimensional tropical optimization problem with nonlinear objective function and linear constraints // Optimization. 2015. Vol. 64, N 5. P. 1107-1129. DOI: 10.1080/ 02331934.2013.84062
11. Krivulin N. K. Algebraic solutions to multidimensional minimax location problems with Chebyshev distance // Recent Researches in Applied and Computational Mathematics: Intern. conference on Applied and Computational Mathematics (ICACM'11). WSEAS, 2011. P. 157-162.
12. Krivulin N. K. Algebraic solution to a constrained rectilinear minimax location problem on the plane // Multimedia Technology (ICMT), 2011 Intern. conference on. IEEE, 2011. P. 6216-6220. DOI: 10.1109/ICMT.2011.6002526
13. Krivulin N. K. A new algebraic solution to multidimensional minimax location problems with Chebyshev distance // WSEAS Transactions on Mathematics. 2012. Vol. 11, N 7. P. 605-614.
14. Krivulin N. K. An extremal property of the eigenvalue of irreducible matrices in idempotent algebra and solution of the Rawls location problem // Vestnik of Saint Petersburg University. Mathematics. 2011. Vol. 44, iss. 4. P. 272-281. DOI: 10.3103/S1063454111040078
15. Krivulin N. K., Plotnikov P. V. On an algebraic solution of the Rawls location problem in the plane with rectilinear metric // Vestnik of Saint Petersburg University. Mathematics. 2015. Vol. 48, iss. 2. P. 75-81. DOI: 10.3103/S1063454115020065
16. Elzinga J., Hearn D. W. Geometrical solutions for some minimax location problems // Transport. Sci. 1972. Vol. 6, N 4. P. 379-394. DOI: 10.1287/trsc.6.4.379
17. Francis R. L. A geometrical solution procedure for a rectilinear distance minimax location problem // AIIE Trans. 1972. Vol. 4, N 4. P. 328-332. DOI: 10.1080/05695557208974870
18. Krivulin N. K., Plotnikov P. V. Using tropical optimization to solve minimax location problems with rectilinear metric on the line // Vestnik of Saint Petersburg University. Mathematics. 2016. Vol. 3(61), iss. 4. P. 602-614. DOI: 10.3103/S1063454116040087
Статья рекомендована к печати проф. Л. А. Петросяном.
Статья поступила в редакцию 27 декабря 2017 г.; принята к печати 15 марта 2018 г.
Контактная информация:
Плотников Павел Владимирович — аспирант; pavplot@gmail.com
Кривулин Николай Кимович — докт. физ.-мат. наук, профессор; nkk@math.spbu.ru
Direct solution of a minimax location problem on the plane with rectilinear metric in a rectangular area
P. V. Plotnikov, N. K. Krivulin
St. Petersburg State University, 7—9, Universitetskaya nab., St. Petersburg, 199034, Russian Federation
For citation: Plotnikov P. V., Krivulin N. K. Direct solution of a minimax location problem on the plane with rectilinear metric in a rectangular area. Vestnik of Saint Petersburg University. Applied Mathematics. Computer Science. Control Processes, 2018, vol. 14, iss. 2, pp. 116-130. https://doi.org/10.21638/11702/spbu10.2018.204
A minimax single-facility location problem with rectilinear (Manhattan) metric is examined under constraints on the feasible location region, and a direct, explicit solution of the problem is suggested using methods of tropical (idempotent) mathematics. When no constraints are imposed, this problem, which is also known as the Rawls problem or the messenger boy problem, has known geometric and algebraic solutions. In the present article, a solution to the problem is investigated subject to constraints on the feasible region, which is given by a rectangular area. At first, the problem is represented in terms of tropical mathematics as a tropical optimization problem, a parameter is introduced to represent the minimum value of the objective function, and the problem is reduced to a parametrized system of inequalities. This system is solved for one variable, and the existence conditions of solution are used to obtain optimal values of the second parameter by using an auxiliary optimization problem. Then, the obtained general solution is transformed into a set of direct solutions, written in a compact closed form for different cases of relations between the initial parameters of the problem. Graphical illustrations of the solution are given for several positions of the feasible location region on the plane.
Keywords : Rawls location problem, constrained location, rectilinear metric, idempotent semifield, tropical optimization, complete solution.
References
1. Baccelli F., Cohen G., Olsder G. J., Quadrat J.-P. Synchronization and linearity. Chichester, Wiley Press, 1993, 514 p. (Wiley Series in Probability and Statistics)
2. Maslov V. P., Kolokol'tsov V. N. Idempotentnyj analiz i ego primenenie v optimal'nom upravlenii [Idempotent analysis and its applications in optimal control]. Moscow, Fizmatlit Publ., 1994, 144 p. (In Russian)
3. Cuninghame-Green R. A. Minimax algebra and applications. Advances in Imaging and Electron Physics, vol. 90 / ed. by P. W. Hawkes. San Diego, Academic Press, 1994, pp. 1-121. DOI: 10.1016/S1076-5670(08)70083-1
4. Golan J. S. Semirings and affine equations over them. New York, Springer Press, 2003, 256 p. (Vol. 556 of Mathematics and its Applications). DOI: 10.1007/978-94-017-0383-3
5. Heidergott B., Olsder G. J., van der Woude J. Max plus at work. Princeton, Princeton University Press, 2006, 226 p. (Princeton Series in Applied Mathematics)
6. McEneaney W. M. Max-plus methods for nonlinear control and estimation. Boston, Birkhauser Press, 2006, 241 p. (Systems and Control: Foundations and Applications). DOI: 10.1007/0-8176-4453-9
7. Itenberg I., Mikhalkin G., Shustin E. I. Tropical algebraic geometry. Basel, Birkhauser Press, 2007, 104 p. (Vol. 35 of Oberwolfach Seminars). DOI: 10.1007/978-3-7643-8310-7
8. Gondran M., Minoux M. Graphs, dioids and semirings. New York, Springer Press, 2008, 383 p. (Vol. 41 of Operations Research/ Computer Science Interfaces). DOI: 10.1007/978-0-387-75450-5
9. Butkovic P. Max-linear systems. London, Springer Press, 2010, 272 p. (Springer Monographs in Mathematics). DOI: 10.1007/978-1-84996-299-5
10. Krivulin N. K. A multidimensional tropical optimization problem with nonlinear objective function and linear constraints. Optimization, 2015, vol. 64, no. 5, pp. 1107-1129. DOI: 10.1080/ 02331934.2013.84062
11. Krivulin N. K. Algebraic solutions to multidimensional minimax location problems with Chebyshev distance. Recent Researches in Applied and Computational Mathematics. Intern. conference on Applied and Computational Mathematics (ICACM'11), WSEAS, 2011, pp. 157-162.
12. Krivulin N. K. Algebraic solution to a constrained rectilinear minimax location problem on the plane. Multimedia Technology (ICMT), 2011 Intern. conference on. IEEE, 2011, pp. 6216-6220. DOI: 10.1109/ICMT.2011.6002526
13. Krivulin N. K. A new algebraic solution to multidimensional minimax location problems with Chebyshev distance. WSEAS Transactions on Mathematics, 2012, vol. 11, no. 7, pp. 605-614.
14. Krivulin N. K. An extremal property of the eigenvalue of irreducible matrices in idempotent
algebra and solution of the Rawls location problem. Vestnik of Saint Petersburg University. Mathematics, 2011, vol. 44, iss. 4, pp. 272-281. DOI: 10.3103/S1063454111040078
15. Krivulin N. K., Plotnikov P. V. On an algebraic solution of the Rawls location problem in the plane with rectilinear metric. Vestnik of Saint Petersburg University. Mathematics, 2015, vol. 48, iss. 2, pp. 75-81. DOI: 10.3103/S1063454115020065
16. Elzinga J., Hearn D. W. Geometrical solutions for some minimax location problems. Transport. Sci., 1972, vol. 6, no. 4, pp. 379-394. DOI: 10.1287/trsc.6.4.379
17. Francis R. L. A geometrical solution procedure for a rectilinear distance minimax location problem. AIIE Trans, 1972, vol. 4, no. 4, pp. 328-332. DOI: 10.1080/05695557208974870
18. Krivulin N. K., Plotnikov P. V. Using tropical optimization to solve minimax location problems with rectilinear metric on the line. Vestnik of Saint Petersburg University. Mathematics, 2016, vol. 3(61), iss. 4, pp. 602-614. DOI: 10.3103/S1063454116040087
Author's Information:
Plotnikov Pavel V. — postgraduate student; pavplot@gmail.com
Krivulin Nikolai K. — Dr. Sci. in physics and mathematics, professor; nkk@math.spbu.ru