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

УДК 519.8

Вестник СПбГУ. Сер. 1. Т. 2(60). 2015. Вып. 2

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

Н. К. Кривулин, П. В. Плотников

Санкт-Петербургский государственный университет,

Российская Федерация, 199034, Санкт-Петербург, Университетская наб., 7-9

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

Ключевые слова: идемпотентное полуполе, спектральный радиус матрицы, полное решение, прямоугольная метрика, задача Ролса о размещении.

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

В работе [8] изучается минимаксная задача размещения одиночного объекта на плоскости с прямоугольной метрикой (метрикой Li), которая возникает, например, при размещении пунктов экстренной помощи населению (больниц, пожарных депо) в городах с прямолинейными перпендикулярными друг другу улицами. Эта задача известна как задача Ролса, названная так по имени автора теории справедливости, американского философа Дж. Ролса (J. Rawls, 1921-2002), а также как задача посыльного.

Задача Ролса сводится в терминах тропической математики к нахождению таких векторов x, при которых достигается минимум в задаче

min x-Ax, (1)

x

где x- обозначает мультипликативно сопряженный вектор и A — некоторая матрица.

Известно [11], что минимум в задаче (1) равен спектральному радиусу А матрицы A и достигается на любом собственном векторе, который отвечает А. Для решения задачи Ролса в работе [8] потребовалось показать, что на самом деле множество решений задачи (1) шире, чем множество собственных векторов матрицы A. Было построено решение задачи Ролса в явном виде, однако проверить, что таким образом получены все решения и других решений нет, при этом не удалось.

* Работа выполнена при финансовой поддержке РГНФ в рамках научного проекта № 13-02-00338.

Новое общее решение задачи (1) при весьма общих условиях было затем найдено в работах [12-14], что позволяет теперь уточнить полученные в [8] результаты. Цель настоящей работы состоит в развитии на основе найденного общего решения задачи (1) нового метода исследования задачи Ролса, который гарантирует ее полное решение.

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

2. Вводные замечания и предварительные результаты. В этом разделе приведен краткий обзор основных понятий и результатов тропической математики, на которые опирается последующее исследование и решение задачи Ролса. Подробное изложение теории и методов тропической математики, а также примеры разнообразных приложений, представлены, например, в работах [1-6].

Пусть числовое множество X замкнуто относительно операций сложения ф и умножения (, а также содержит нуль 0 и единицу 1. Будем считать, что (X, 0,1, ф, () является коммутативным полукольцом, в котором сложение идемпотентно (т.е. х ф х = х), а умножение обратимо (т.е. у каждого элемента х = 0 существует обратный элемент х-1 такой, что х ( х-1 = 1), и называть это полукольцо идемпотентным полуполем.

Целая степень числа в полуполе определяется обычным путем: для любого х = 0 и целого п > 0 выполняется х0 = 1, хп = хп-1 ( х, х-п = (х-1)" и 0" = 0. Считается, что операция возведения в степень может быть распространена на случай рационального показателя. В дальнейшем знак умножения ( в выражениях будет опускаться. Обозначение степени будет использоваться в смысле идемпотентного полуполя.

Идемпотентное сложение определяет на X частичный порядок по правилу х < у тогда и только тогда, когда х ф у = у. Операции сложения и умножения в полуполе, очевидно, являются монотонными относительно этого порядка. Предполагается, что частичный порядок может быть продолжен до линейного порядка на X. Заметим, что тогда X образует делимую упорядоченную абелеву группу по умножению.

Примером идемпотентного полуполя рассматриваемого типа является полуполе вещественных чисел Ктах,+ = (К и {-то}, —то, 0, тах, +).

В Ктах,+ нулевым элементом является —то, а единичным — число 0. Каждому элементу х € К сопоставляется обратный элемент х-1, равный —х в обычной алгебре. Для любой пары х, у € К определена степень ху, значение которой соответствует арифметическому произведению ху. Порядок, порожденный операцией идемпотентного сложения, совпадает с обычным линейным порядком на множестве К.

Множество столбцов (векторов) порядка п с элементами из X обозначается X". Вектор, все элементы которого равны 0, называется нулевым, а вектор, который вообще не имеет нулевых элементов, называется регулярным.

Сложение двух векторов и умножение вектора на число выполняется обычным путем с использованием операций ф и ( .

Вектор Ь линейно зависит от векторов а1,..., ат, если Ь = и^а1 ф • • • ф штап, где и>1,..., и>т € X. Как обычно, множество всех линейных комбинаций векторов а1,..., ат образует линейную оболочку этих векторов.

Для ненулевого вектора а = (а) € X" определен мультипликативно сопряженный вектор-строка а- = (а-), где а- = а-1, если а^ = 0, и а- = 0 — в противном случае.

Обозначим через Xmxn множество матриц из m строк и n столбцов, состоящих из элементов множества X. Сложение и умножение матриц, а также умножение на число выполняются по стандартным правилам с заменой обычных операций на ф и

Рассмотрим квадратные матрицы из множества Xnxn. Матрица, диагональные элементы которой равны 1, а остальные элементы равны 0, называется единичной и обозначается I. Для любой матрицы A и целого m > 0 выполняется A0 = I,

Am = Am-1A.

След матрицы A = (aij) вычисляется по формуле trA = ац ф • • • ф ann. Пусть задана некоторая матрица A G Xnxn. Рассмотрим задачу нахождения регулярных векторов x G Xn, при которых достигается минимум в задаче (1).

Чтобы описать общее решение задачи, потребуются следующие определения. Спектральным радиусом матрицы A называется величина

n

Л = 0tr1/m(Am).

m= 1

Матрице A можно сопоставить матрицу (известную также как матрица Клини)

A* = I ф A ф---ф An-1.

В работах [12-14] получено общее решение задачи в следующей форме.

Теорема 1. Пусть A —матрица со спектральным радиусом Л > 0 и A\ = Л-1 A. Тогда минимум в задаче (1) равен Л, а общее регулярное решение имеет вид

x = A*xv, v G Xn.

На основе этого результата ниже будет найдено полное решение минимаксной задачи размещения одиночного объекта на плоскости с прямоугольной метрикой.

3. Задача размещения с прямоугольной метрикой. Пусть на плоскости R2 задана прямоугольная метрика, которая для любых двух векторов x = (ж1,ж2)т и У = (y 1, У2 )T вычисляется в обычной арифметике по формуле

p(x, у) = |Х1 - У11 + |Х2 - У2 |.

Рассмотрим следующую минимаксную задачу размещения. Предположим, что для всех i = 1,..., m заданы векторы vi = (r^, r'2i)T G R2 и числа wi G R. Требуется найти все векторы x = (x ,X2)T G R2, при которых достигается минимум

min max (p(vi, x) + wi). (2)

x£R2 1<i<m

Существует геометрическое решение этой задачи [15, 16], известной как задача Ролса или задача посыльного. Алгебраическое решение в терминах тропической математики было дано в работе [8], однако оставалось неясным, является ли оно полным.

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

Сначала так же, как в работе [8], представим задачу в терминах идемпотентного полуполя Ктах,+ - Прямоугольная метрика записывается следующим образом:

у) = (х- ф У- Ж1)(Ж2 У2 ф У2 Х2). Целевая функция в задаче (2) может быть представлена в форме

^J^ Wi(x- 1 rij ф r- 1 xi)(x2 1 r2i © r2i 1Ж2) = ax- ^2 ф bx 1 x- 1 ф cx- 1 x- 1 Ф dx 1 Ж2,

i=1

где после группировки членов использованы обозначения

m m m m

a = ф Wjr1ir2i1, b = ф Wjr1i1r2i, c = ф Wjr1ir2i, d =ф Wir1i1 r^1. i=1 i=1 i=1 i=1

Тогда задача размещения (2) принимает вид

min (ax-1 x2 Ф bx1x-1 ф cx-^x-1 ф dx1x2). (3)

xGR2

Чтобы перейти к векторной форме записи, в работе [8] вводятся вектор и матрица:

/ x1 \ /0 a 0

У = ( x2 I , A = I b 0 c

x12 1 0 d 0

Затем исходная двумерная задача сводится к расширенной задаче оптимизации на множестве трехмерных векторов

min y-Ay, (4)

yeR3

из числа решений которой следует выбрать те векторы, у которых первый и последний элементы — взаимно-обратные (противоположные — в терминах обычной арифметики) .

4. Анализ и новое решение задачи. Пусть А — спектральный радиус матрицы A и A\ = А-1 A. По теореме 1 минимум в задаче (4) равен А. Все решения имеют вид

y = A>, v G R3, (5)

и поэтому образуют линейную оболочку столбцов матрицы A^.

Для определения спектрального радиуса А найдем степени матрицы A

ab 0 ac \ / 0 a2b ф acd 0

A2 = | 0 ab ф cd 0 ! , A3 = I ab2 © bcd 0 abc © c2d bd 0 cd / \ 0 abd ф cd2 0

а затем вычислим

А = tr(A) ф tr1/2(A2) ф tr1/3(A3) = (ab ф cd)1/2. Заметим, что при этом выполняются неравенства ab < А2, cd < А2 и abcd < А4.

Чтобы записать решения расширенной задачи (4), необходимо вычислить мат-

рщу

( 1 А-1а А-2ае \

Л*х = I ® Ал ® Л\ = ( А-1Ъ 1 А-1с ) .

\ А-2Ъа А-1а 1 /

Нетрудно проверить, что второй столбец матрицы ЛЛ является линейной комбинацией двух других. Действительно,

( 1 \ ( А-2ае \ ( А-1а \

А-1а ( А-1Ъ ) ® А-1 а ( А-1с ) = ( 1 ) . \ А-2Ъа / \ 1 ) \ А-Ч )

Тогда для представления всех решений (5) достаточно построить линейную оболочку только первого и последнего столбцов Лл в виде

/ 1 А-2ае \ у = ( А-1Ъ А-1с ) V, V е К2.

\ А-2ъа 1 /

Переходя к скалярной форме записи, получим три равенства:

У1 = «1 ® А-2ае«2, у2 = А-1Ъ«1 ® А-1 е«2, уз = А-2Ъ^У1 ® «2.

Для решения исходной задачи сначала нужно найти все векторы у = (у1, у2, Уз)Т, которые удовлетворяют этим равенствам и условию для координат

У1 = у-1.

После этого останется определить все векторы х = (х , Х2)т, координаты которых определяются равенствами

Х1 = «1 ® А-2ае«2, Х2 = А-1Ъ«1 ® А-1с«2. (6)

Условие для первой и третьей координаты вектора у записывается в виде

® А-2ае«2 = (А-2Ъ^«1 ® г>2)-1.

После умножения на А-2® «2 имеем равенство А-2Ъ^«2 ® «1«2 ® А-2ае«2 = 1. Полученное равенство равносильно системе неравенств

< АЪ-1/2а-1/2, «2 < 1, «2 < Аа-1/2с-1/2, (7)

в которой хотя бы одно неравенство должно выполняться как равенство.

Рассмотрим три случая, в которых одно из неравенств заменяется равенством.

4-1- Случай 1. Первое неравенство выполняется как равенство. Заменим первое неравенство в (7) равенством, которое будет однозначно определять значение «1. После подстановки этого значения во второе неравенство имеем систему

«1 = АЪ-1/2а-1/2, «2 < А-1Ъ1/2а1/2, «2 < Аа-1/2е-1/2.

В силу условия абей < Л4 выполняется неравенство Л 161/2й1/2 < Ла 1/2с 1/2, а тогда полученная система сводится к системе

«1 = Л6-1/2й-1/2, «2 < Л-161/2й1/2.

Подставим найденные значения в решение (6). Применяя неравенства аб < Л2 и ей < Л2, находим, что Л-2ае«2 < Л-3а61/2ей1/2 < и Л-1е«2 < Л-261/2ей1/2 < Л-1

Тогда решение (6) принимает вид

X = = Л6-1/2й-1/2, х2 = Л-16«1 = 61/2й-1/2.

Нетрудно проверить, что при условии а = 1 решение можно записать в форме

_ ( Л2а-1 (а1-а6-ае1-ай-а)1/2 \ Х = ^ (аа — 1ЬаС1 —а)1/2 ) . (8)

4-2. Случай 2. Второе неравенство выполняется как равенство. Предположим, что второе неравенство в системе (7) является равенством. В этом случае система принимает вид

< Л6-1/2й-1/2, «2 = 1, «2 < Ла-1/2е-1/2.

Воспользуемся равенством «2 = «-1. Первое и третье условия системы можно объединить в виде двойного неравенства

Л-1 а1/2е1/2 < «1 < Л6-1/2й-1/2.

Чтобы записать решение исходной задачи в форме (6), сначала заметим, что

Л-2ае«2 < Л-1 а1/2е1/2 < «1.

При этом, если Л2 = аб, то имеем неравенство Л-1 е«2 < Л-161/2е1/2 < Л-16«1, откуда следует, что решение можно записать в виде

X = «1, х2 = Л-16«1, Л-1а1/2е1/2 < < Л6-1/2й-1/2.

Введем параметр 0 < а < 1 и запишем решение в параметрической форме. Заменим двойное неравенство для «1 выражением

«1 = аа/26(а-1)/2е(1-а)/2й-а/2 = Л2а-1(а1-а6-ае1-ай-а)1/2.

Заметим, что при а = 0 это выражение совпадает с левой границей двойного неравенства для «1, а при а =1 оно равно правой границе.

После подстановки приходим к результату, который обобщает решение (8) в форме

Л2а-1 (а1-а6-ае1-^-а)1/2 \ 1 (9)

(аа-1б«с1-ай-а)1/2 ^ 0 < а < 1. (9)

В случае, когда Л2 = ей, справедливо неравенство Л-16«1 < Л-161/2е1/2 < Л-1 е«2. С учетом этого неравенства имеем

Х1 = «1, Х2 = Л-1е«2 = е1/2<Г1/2«-1, Л-1 а1/2е1/2 < «1 < Л6-1/2й-1/2.

Записывая двойное неравенство для «1 в параметрическом виде

VI = а(1-а)/2Ъ-а/2са/2 й(а-1)/2 = Л2а-1(а1-аЪ-ас1-ай-а)1/2, снова получаем решение задачи, представленное как (9).

4-3. Случай 3. Третье неравенство выполняется как равенство. Заменив третье неравенство в системе (7) на равенство, получаем

V1 < ЛЪ-1/2й-1/2, «1 < Л-1а1/2с1/2,

V2

Ла-1/2с-1/2.

В силу неравенства ЛЪ 1/2й 1/2 > Л 1а1/2с1/2 система принимает вид

«1 < Л-1а1/2с1/2,

«2

= Ла-1/2с-1/2.

Учитывая, что < Л 1а1/2с1/2 = Л 2ас«2 и Л 1Ъ«1 < Л 2а1/2Ъс1/2 < Л 1

С«2,

имеем

Х1

= Л-2 ас«2 = Л-1а1/2 с1/2,

Х2

= Л-1 с«2 = а-1/2с1/2.

Осталось проверить, что это решение совпадает с (8) при а = 0.

4-4- Формулировка основного результата. Результаты анализа задачи (3) для всех рассмотренных случаев можно объединить в виде следующего утверждения, которое уточняет решение, полученное в [8].

Теорема 2. Минимум Л = (аЪ ф ей)1/2 в задаче (3) достигается тогда и только тогда, когда

Л2а-1 (а1-аЪ-ас1-ай-а )1/2 (аа-1Ъас1-ай-а)1/2

0 < а < 1,

где

е

¿=1

^ ГН Г2-\

Ъ = 0 Гн1^

¿=1

С = ф ШгГн^, ¿=1

й =0

¿=1

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

Следствие 1. Минимум Л = тах(а + Ъ, с + й)/2 в задаче (2) достигается тогда и только тогда, когда

(2а - 1)Л + (1 - а)(а + с)/2 - а(Ъ + й)/2 (1 - а)(с - а)/2 - а(й - Ъ)/2

0 < а < 1,

где

а = тах (шг + ги - ^), с = тах (^ + гЦ + ^),

1< ¿<т

Ъ = тах (и^ - Гн + ^),

л, = тах (^ - Гц - ^).

1< ¿<т

а

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

Литература

1. Cuninghame-Green R. A. Minimax algebra and applications // Advances in Imaging and Electron Physics / ed. by Peter W. Hawkes. San Diego, CA: Academic Press, 1994. Vol. 90 of Advances in Imaging and Electron Physics. P. 1-121.

2. Маслов В. П., Колокольцов В.Н. Идемпотентный анализ и его применение в оптимальном управлении. М.: Физматлит, 1994. 144 с.

3. Golan J. S. Semirings and Affine Equations Over Them: Theory and Applications. New York: Springer, 2003. Vol. 556 of Mathematics and Its Applications. 256 p.

4. Heidergott B., Olsder G. J., van der Woude J. Max-plus at Work: Modeling and Analysis of Synchronized Systems. Princeton Series in Applied Mathematics. Princeton, NJ: Princeton University Press, 2006. 226 p.

5. Кривулин Н. К. Методы идемпотентной алгебры в задачах моделирования и анализа сложных систем. СПб.: Изд-во С.-Петерб. ун-та, 2009. 256 с.

6. Butkovic P. Max-linear Systems: Theory and Algorithms. Springer Monographs in Mathematics. London: Springer, 2010. 272 p.

7. Zimmermann K. Disjunctive optimization, max-separable problems and extremal algebras // Theoret. Comput. Sci. 2003. Vol. 293, N 1. P. 45-54.

8. Krivulin N. K. An extremal property of the eigenvalue of irreducible matrices in idempotent algebra and solution of the Rawls location problem // Vestnik St. Petersburg Univ. Math. 2011. Vol. 44, N4. P. 272-281.

9. Krivulin N. A new algebraic solution to multidimensional minimax location problems with Chebyshev distance // WSEAS Trans. Math. 2012. Vol. 11, N7. P. 605-614.

10. Krivulin N. Complete solution of a constrained tropical optimization problem with application to location analysis // Relational and Algebraic Methods in Computer Science / ed. by P. Hofner, P. Jipsen, W. Kahl, M. E. Müller. Springer, 2014. Vol. 8428 of Lecture Notes in Computer Science. P. 362-378.

11. Cuninghame-Green R. A. Describing industrial processes with interference and approximating their steady-state behaviour // Oper. Res. Quart. 1962. Vol. 13, N 1. P. 95-100.

12. Krivulin N. A multidimensional tropical optimization problem with nonlinear objective function and linear constraints // Optimization. 2015. Vol.64, N5. P. 1107-1129.

13. Krivulin N. A constrained tropical optimization problem: Complete solution and application example // Tropical and Idempotent Mathematics and Applications / ed. by G. L. Litvinov, S. N. Sergeev. Providence, RI: AMS, 2014. Vol. 616 of Contemp. Math. P. 163-177.

14. Krivulin N. Extremal properties of tropical eigenvalues and solutions to tropical optimization problems // Linear Algebra Appl. 2015. Vol.468. P. 211-232.

15. Elzinga J., Hearn D. W. Geometrical solutions for some minimax location problems // Transport. Sci. 1972. Vol. 6, N4. P. 379-394.

16. Francis R. L. A geometrical solution procedure for a rectilinear distance minimax location problem // AIIE Trans. 1972. Vol. 4, N4. P. 328-332.

Статья поступила в редакцию 25 декабря 2014 г. Сведения об авторах

Кривулин Николай Кимович — доктор физико-математических наук, доцент, профессор; nkk@math.spbu.ru

Плотников Павел Владимирович — аспирант; pavplot@gmail.com

ON ALGEBRAIC SOLUTION OF THE RAWLS LOCATION PROBLEM ON THE PLANE WITH RECTILINEAR METRIC

Nikolay K. Krivulin, Pavel V. Plotnikov

St.Petersburg State University, Universitetskaya nab., 7-9, St.Petersburg, 199034, Russian Federation; nkk@math.spbu.ru, pavplot@gmail.com

A minimax single-facility location problem with rectilinear metric is examined by using methods of tropical (idempotent) mathematics. This problem, known as the Rawls problem or the messenger problem, appears

in the location of emergency services (hospitals, fire stations) in the cities with straight rectangular streets. In terms of idempotent algebra, the problem is reduced to minimizing a functional, given on the set of three-dimensional vectors by an appropriately constructed matrix and calculated by using the multiplicative conjugate transposition. The minimum of the objective function is found subject to constraints in the form of a relation that holds between components of the vectors. A new result of the spectral theory of matrices in idempotent algebra is applied, which offers a general solution to the problem of minimizing such functionals without additional constraints. Based on the result, a general solution to the problem with constraints on the feasible solution is given in terms of the tropical algebra. The solution obtained is then used to derive a complete solution to the Rawls location problem, which extends a known particular solution of the problem under consideration. Refs 16.

Keywords: idempotent semifield, spectral radious of matrix, complete solution, rectilinear metric, the Rawls location problem.

References

1. Cuninghame-Green R. A. "Minimax algebra and applications", Advances in Imaging and Electron Physics / ed. by Peter W. Hawkes. San Diego, CA: Academic Press, 1994. Vol. 90 of Advances in Imaging and Electron Physics. P. 1-121.

2. Maslov V. P., Kolokoltsov V. N., Idempotent Analysis and Its Applications to Optimal Control Theory (Moscow, Nauka, 1994) [in Russian].

3. Golan J. S. Semirings and Affine Equations Over Them: Theory and Applications (New York, Springer, 2003. Vol. 556 of Mathematics and Its Applications) 256 p.

4. Heidergott B., Olsder G.J., van der Woude J. "Max-plus at Work: Modeling and Analysis of Synchronized Systems", Princeton Series in Applied Mathematics (Princeton, NJ, Princeton University Press, 2006) 226 p.

5. Krivulin N. K., Methods of Idempotent Algebra for Problems in Modeling and Analysis of Complex Systems (St. Petersburg, Saint Petersburg Univ. Press, 2009) [in Russian].

6. Butkovic P. Max-linear Systems: Theory and Algorithms, Springer Monographs in Mathematics (London, Springer, 2010) 272 p.

7. Zimmermann K. "Disjunctive optimization, max-separable problems and extremal algebras", Theoret. Comput. Sci. 293(1), 45-54 (2003).

8. Krivulin N. K. "An extremal property of the eigenvalue of irreducible matrices in idempotent algebra and solution of the Rawls location problem", Vestnik St. Petersburg Univ. Math. 44(4), 272-281 (2011).

9. Krivulin N. "A new algebraic solution to multidimensional minimax location problems with Chebyshev distance", WSEAS Trans. Math. 11(7), 605-614 (2012).

10. Krivulin N. "Complete solution of a constrained tropical optimization problem with application to location analysis", Relational and Algebraic Methods in Computer Science / ed. by P. Höfner, P. Jipsen, W. Kahl, M. E. Möller (Switzerland, Springer, 2014). Vol. 8428 of Lecture Notes in Computer Science. P. 362-378.

11. Cuninghame-Green R. A. "Describing industrial processes with interference and approximating their steady-state behaviour", Oper. Res. Quart. 13(1), 95-100 (1962).

12. Krivulin N. "A multidimensional tropical optimization problem with nonlinear objective function and linear constraints", Optimization 1-23 (2013).

13. Krivulin N. "A constrained tropical optimization problem: Complete solution and application example", Tropical and Idempotent Mathematics and Applications / ed. by G. L. Litvinov, S. N. Sergeev (Providence, RI, AMS, 2014). Vol. 616 of Contemp. Math. P. 163-177.

14. Krivulin N. "Extremal properties of tropical eigenvalues and solutions to tropical optimization problems", Linear Algebra Appl. 1-22 (2014).

15. Elzinga J., Hearn D. W. "Geometrical solutions for some minimax location problems", Transport. Sci. 6(4), 379-394 (1972).

16. Francis R. L. "A geometrical solution procedure for a rectilinear distance minimax location problem", AIIE Trans. 4(4), 328-332 (1972).