CC BY
29
УДК 519.87
Вестник СПбГУ. Прикладная математика... 2017. Т. 13. Вып. 4
Н. К. Кривулин1, С. А. Губанов2
ИСПОЛЬЗОВАНИЕ МЕТОДОВ ТРОПИЧЕСКОЙ ОПТИМИЗАЦИИ В ЗАДАЧАХ СЕТЕВОГО ПЛАНИРОВАНИЯ*
1 Санкт-Петербургский государственный университет, Российская Федерация, 199034, Санкт-Петербург, Университетская наб., 7—9
2 Санкт-Петербургский филиал АО «КБ "Луч"», Российская Федерация, 197376, Санкт-Петербург, ул. Академика Павлова, 14А
Статья посвящена решению задач сетевого планирования с использованием методов тропической оптимизации. Исследуются задачи составления оптимального плана сроков реализации проекта, который состоит в выполнении некоторого набора взаимосвязанных работ при заданных ограничениях на время их начала и завершения. Рассматриваются критерии оптимальности плана, при которых требуется максимизировать разброс времени начала или завершения работ. Такие задачи появляются при необходимости по тем или иным причинам (недостаток ресурсов, технологические ограничения, требования безопасности и т. п.) по возможности избежать одновременного начала или завершения всех работ проекта. Сначала формулируются задачи планирования в виде обычных задач оптимизации с ограничениями. Далее приводятся определения и результаты тропической математики, которые применяются в последующем анализе и решении задач тропической оптимизации. Описываются новые задачи тропической оптимизации с ограничениями и находятся их решения. Задачи планирования решаются путем их сведения к задачам тропической оптимизации. Построены новые решения задач планирования, которые в отличие от существующих алгоритмических решений дают результат в виде прямых расчетных формул в компактном векторном представлении, удобных для дальнейшего анализа и практических расчетов. Для иллюстрации полученных результатов представлены численные примеры. Библиогр. 15 назв.
Ключевые слова: тропическая математика, идемпотентное полуполе, тропическая оптимизация, управление проектами, сетевое планирование.
N. K. Krivulin1, S. A. Gubanov2
THE USE OF TROPICAL OPTIMIZATION METHODS IN PROBLEMS OF PROJECT SCHEDULING
1 St. Petersburg State University, 7—9, Universitetskaya nab., St. Petersburg, 199034, Russian Federation
2 St. Petersburg Office "KB Lutch", 14A, ul. Academician Pavlov, St. Petersburg, 197376, Russian Federation
The paper is devoted to the solution of problems in project scheduling by using methods of tropical optimization. Problems are examined that are to develop an optimal schedule for a project consisting in the execution of a set of interrelated tasks under given constraints on their initiation and completion time. The optimal schedule criteria are considered, which require the maximization of the deviation of the initiation or the deviation of the completion time of tasks. Such problems arise when, for some reason (such as a lack of resources, technical constraints,
Кривулин Николай Кимович — доктор физико-математических наук, профессор; nkk@math.spbu.ru
Губанов Сергей Александрович — программист; segubanov@mail.ru
Krivulin Nikolai Kimovich — doctor of physical and mathematical sciences, professor; nkk@math.spbu.ru
Gubanov Sergey Aleksandrovich — programmer; segubanov@mail.ru
* Работа выполнена при финансовой поддержке Российского гуманитарного научного фонда (проект № 16-02-00059).
© Санкт-Петербургский государственный университет, 2017
security requirements, and the like), there is a need to avoid a simultaneous start or finish for all tasks in the project. The paper begins with the formulation of scheduling problems in the form of usual optimization problems. Then, definitions and results of tropical mathematics are given, which are used in the subsequent analysis and solution of tropical optimization problems. New constrained problems of tropical optimization are considered, and their solutions are obtained. The scheduling problems are solved by reducing to tropical optimization problems. To illustrate the results obtained, numerical examples are presented. Refs 15.
Keywords : tropical mathematics, idempotent semifield, tropical optimization, project management, project scheduling.
1. Введение. Под управлением проектом обычно понимается согласование действий, выполняемых для достижения целей проекта при разумном распределении имеющихся ресурсов. Одним из подходов к управлению временем выполнения (составлению графика работ) проекта служит применение методов сетевого планирования, базирующихся на теории графов (например, метод критического пути или метод оценки и пересмотра планов [1]). Перспективным способом решения задач сетевого планирования является их сведение к задачам тропической оптимизации [2-4].
Тропическая (идемпотентная) математика изучает полукольца с идемпотентным сложением [5-9]. Задачи тропической оптимизации обычно задаются в терминах минимизации (максимизации) функций, определенных на векторах над идемпотентны-ми полуполями (идемпотентными полукольцами, в которых умножение обратимо) при ограничениях в виде тропических линейных уравнений и неравенств.
Цель настоящей работы — представить новые способы решения задач сетевого планирования на основе методов тропической математики, опирающиеся на результаты работ [10-14]. В частности, в статье [14] задача сетевого планирования заключалась в минимизации разброса времени завершения всех работ. Ниже рассматриваются задачи составления оптимального плана сроков выполнения работ в соответствии с критериями максимума разброса времени начала или завершения работ. Такие задачи появляются при необходимости по тем или иным причинам (недостаток ресурсов, технологические ограничения, требования безопасности и т. п.) по возможности избежать одновременного начала или завершения всех работ проекта.
Статья организована следующим образом. В п. 2 формулируются исследуемые задачи сетевого планирования, для которых определяются критерии оптимальности и задаются ограничения. В п. 3 приводятся основные определения, обозначения и предварительные результаты идемпотентной математики, необходимые для анализа и решения задач тропической оптимизации. В п. 4 описываются новые задачи тропической оптимизации с ограничениями, для которых предлагаются прямые решения в явном виде. В п. 5 полученные результаты применяются к задачам сетевого планирования, которые формулируются и решаются как задачи тропической оптимизации. Приводятся примеры численного решения рассматриваемых задач.
2. Задачи сетевого планирования. Рассмотрим задачи оптимизации, которые возникают при планировании сроков начала и завершения работ различных проектов и сводятся к максимизации некоторой целевой функции при ограничениях в форме равенств и неравенств. Целевая функция определяется критерием оптимальности плана, а ограничения задаются условиями на время начала и завершения работ проекта. В этих задачах целевая функция и ограничения выражаются при помощи операций максимума, сложения и вычитания.
Пусть имеется проект, состоящий в выполнении n работ при ограничениях типа «старт-финиш», «старт-старт», «поздний финиш» и «ранний старт». Ограничение «старт-финиш» не позволяет завершить работу, пока не прошло определенное время
после начала других работ. Минимальный допустимый интервал между временем начала одной работы и временем начала другой задается ограничением «старт-старт». Ограничение «ранний старт» обусловливает самое раннее возможное время начала, а ограничение «поздний финиш» указывает наиболее позднее возможное время завершения работы. Будем считать, что каждая работа заканчивается немедленно, как только будут выполнены все заданные ограничения на время ее завершения.
2.1. Ограничения на время выполнения работ. Запишем определенные выше ограничения в формальном виде. Для каждой работы i = 1,...,n введем обозначения: Xi — время начала работы; yi — время завершения работы; aj — минимальный промежуток времени между началом работы j = 1,...,n и завершением работы i, если величина интервала aij не задана, полагаем aij = —ж; bij — минимальный промежуток времени между началом работы j и началом i, если интервал bij не задан, то считаем, что bj = —ж; gi — наиболее раннее возможное время начала работы; hi — наиболее позднее возможное время завершения работы.
Ограничения типа «старт-финиш» выполняются для работы i, когда справедливы неравенства
aij + Xj < yi, j = 1,...,n.
Учитывая, что работа завершается немедленно при выполнении этих ограничений, хотя бы одно из неравенств выполняется как равенство. Тогда эти неравенства эквивалентны следующему условию:
max (aij + Xj) = yi.
l^j^n
Ограничения типа «старт-старт» описываются с помощью неравенств
bij + Xj ^ xi, j 1 ,...,n, которые можно объединить в одно неравенство
max (bij + Xj) ^ xi.
l^j^n
Ограничения «ранний старт» и «поздний финиш» для работы i определяются при помощи неравенств
Xi > gi, yi < hi.
2.2. Критерии оптимальности и задачи планирования. В задачах составления оптимального плана проекта при необходимости избежать одновременного начала или завершения работ можно рассматривать следующие критерии оптимальности, которые требуется максимизировать. Одним из критериев является максимальный разброс времени начала всех работ:
max Xi — min Xi = max Xi + max (—Xi).
l^i^n l^i^n l^i^n l^i^n
Другой критерий задает максимальный разброс времени завершения:
max yi — min yi = max yi + max (—yi).
l^i^n l^i^n l^i^n l^i^n
Рассмотрим примеры задач сетевого планирования с описанными выше условиями предшествования и критериями оптимальности. Подобные задачи изучались, например, в работе [15].
Предположим, что заданы ограничения «старт-старт» и «ранний старт». В ситуации нехватки ресурсов (например, людей меньше, чем задач) часто нужно распределить время начала выполнения работ на максимально большом интервале. Для составления графика работ при этих условиях необходимо решить задачу
max ( max xi + max (-xi)),
l^i^n l^i^n
max (bij + Xj) ^ xi, gi ^ xi, i = l,...,n.
l^j^n
Пусть заданы ограничения типа «старт-старт», «старт-финиш» и «поздний финиш». В случае, когда требуется по возможности избежать одновременного завершения всех работ, может оказаться полезным потребовать распределить время их завершения на максимально большом интервале. Такую задачу можно записать таким образом:
max ( max yi + max (-yi)),
l^i^n l^i^n
max (bij + Xj) ^ xmax (aij + Xj)
l^j^n l<,j<n
Заметим, что задачи (1) и (2) можно представить как задачи линейного программирования, а затем найти их численные решения путем применения одного из известных итерационных алгоритмов, например алгоритма симплекс-метода. На основе применения методов тропической математики будут получены прямые решения этих задач в явном виде в компактной векторной форме.
3. Элементы тропической математики. Приведем краткий обзор основных алгебраических определений, обозначений и предварительных результатов [8, 11-13], которые потребуются для последующей формулировки и решения задач тропической оптимизации. Дополнительную информацию можно найти, например, в работах [57, 9].
3.1. Идемпотентное полуполе. Рассмотрим множество X, которое замкнуто относительно ассоциативных и коммутативных операций сложения ф и умножения <8> и содержит их нейтральные элементы: нуль 0 и единицу 1. Сложение обладает свойством идемпотентности (x ф x = x для всех x G X), а умножение дистрибутивно относительно сложения и обратимо (каждый x G Х+, где Х+ = X \ {©}, имеет обратный x-1, удовлетворяющий равенству x-l ® x = 1). Так как X + образует группу по умножению, алгебраическую структуру (X, ®, 1, ф, обычно называют идемпо-тентным полуполем. Для упрощения записи знак умножения ® далее опускается.
Для любого x G Х+ и целого положительного p положим x0 = 1, xp = xp-lx, x-p = (x-l)p, (Dp = (D. Идемпотентное сложение задает следующий частичный порядок: x ^ y тогда и только тогда, когда x ф y = y. Символы отношения и задачи максимизации понимаются в контексте этого порядка. Операции полуполя — монотонные, т. е. из неравенства x ^ y следуют неравенства x ф z ^ y ф z и xz ^ yz для любого z G X. Считаем, что указанный частичный порядок дополнен до линейного порядка на .
Примерами рассматриваемого идемпотентного полуполя являются вещественные полуполя Rmax,+ = (RU{ —сю}, —ж, 0, max, +), Rmin,+ = (R и{+ю}, +ю, 0, min, +), Rmax,x = (R + U {0}, 0,1, max, х) и RminjX = (R + U {+ю}, +ю, 1, min, х), где R — множество вещественных чисел, R + = {x G R|x > 0}.
В полуполе RmaXj+ нуль 0 определен как —ж, а единица 1 — как 0. Любой элемент x G R имеет обратный x-l, которому соответствует —x в стандартных обозначе-
Уi, yi
^ hi
(2)
ниях. Для всех х,у € К степень ху совпадает с арифметическим произведением ху. Отношение порядка, определяемое идемпотентным сложением, отвечает обычному отношению линейного порядка на К.
3.2. Матрицы и векторы.. Пусть Xтхп — множество матриц, состоящих из т строк и п столбцов с элементами из X. Нулевой матрицей в Xтхп является матрица со всеми элементами, равными 0. Регулярной по строкам (по столбцам) называется матрица без нулевых строк (столбцов). Матрица регулярна, если она регулярна по строкам и по столбцам.
Сложение и умножение согласованных по размеру матриц, а также умножение матрицы на скаляр выполняются по обычным правилам с заменой арифметических сложения и умножения на операции ф и <. Эти операции обладают свойством монотонности, которое заключается в том, что из поэлементного неравенства А ^ В следуют поэлементные неравенства А ф С ^ В ф С, АО ^ ВО и хА ^ хВ для всех матриц А, В, С и О соответствующего размера и скаляра х.
Для любой матрицы А ее транспонированная матрица обозначается через АТ.
Для каждой матрицы А = (а^) € Xтхп введем мультипликативно сопряженную матрицу А- = (а-) £Жпхт, где а- = а-1, если а^ = 0, и а- = 0 — иначе.
Рассмотрим квадратные матрицы в Xпхп. Единичной называется матрица с элементами, равными 1 на диагонали и 0 — вне ее, которая обозначается через I.
След матрицы А = (а^) вычисляется по формуле ^ А = ац ф ■ ■■ ф апп.
Введем функцию
Тг( А) = ^ А ф - --ф^ Ап.
Если Тг(А) ^ 1, то определим матрицу, которую обычно называют матрицей Клини:
А* = I ф А ф---фАп-1.
Если все элементы матрицы выше или ниже диагонали равны , то матрица треугольная. Матрица является разложимой, если ее можно привести к блочно-тре-угольной форме путем перестановки строк вместе с такой же перестановкой столбцов. Иначе матрица называется неразложимой. Матрица с регулярными столбцами (строками) включает только ненулевые элементы и потому неразложимая. Нетрудно понять, что для любой неразложимой матрицы ее матрица Клини не имеет нулевых элементов.
Обозначим через п множество вектор-столбцов размерности п.
Вектор 0 = (©,..., 0)Т, все компоненты которого равны 0, нулевой. Вектор без нулевых компонент называется регулярным.
Вектор, состоящий из единиц, обозначается 1 = (1,..., 1 )Т.
Для любого ненулевого вектора х = (х^) € Хп определим мультипликативно сопряженную вектор-строку х- = (х-), где х- = х-1, если х^ = 0, и х- = 0 — иначе.
3.3. Решение векторных неравенств. Предположим, что заданы некоторые матрица А € Жтхп и вектор Ь € Хт. Рассмотрим задачу решения относительно неизвестного вектора х € п неравенства
Ах < Ь. (3)
Решение задачи дает следующее утверждение.
Лемма 1 [8]. Для регулярной по столбцам матрицы А и регулярного вектора Ь вектор х является решением неравенства (3) тогда и только тогда, когда
х < (Ь-А)-.
Пусть для заданных матрицы A G Xnxn и вектора b G Xn требуется найти регулярные векторы x G X, которые удовлетворяют неравенству
Ax ф b < x. (4)
Для решения этого неравенства можно применить следующий результат.
Лемма 2 [11]. Пусть x — общее регулярное решение неравенства (4). Тогда справедливы утверждения:
1) если Tr(A) ^ 1, то x = A*u для любого регулярного вектора u ^ b;
2) если Tr(A) > 1, то регулярных решений не существует.
При b = 0 неравенство (4) принимает вид
Ax ^ x,
решение которого находится как прямое следствие леммы 2.
4. Задачи тропической оптимизации. Методы тропической математики применяются при решении различных задач размещения объектов, сетевого планирования, принятия решений и др. В работах [2-4] приводится решение задач сетевого планирования посредством сведения их к задачам тропической оптимизации. Рассмотрим пример такой задачи оптимизации.
Даны матрица A G Xnxn и векторы p, q G Жп. Необходимо найти регулярные векторы x G n, которые решают задачу
max q-x(Ax)-p. (5)
В работе [13] предлагается такое решение.
Теорема 1. Предположим, что матрица A = (aj) имеет регулярные столбцы aj = (aij), а векторы p = (pj) и q = (qj) являются регулярными.
Тогда максимум в задаче (5) равен Д = q- A-p и достигается тогда и только тогда, когда вектор x = (xj) имеет компоненты
xk = aa-klps, xj < aa-ps, j = к, для всех a > 0 и индексов к и s, определяемых условиями
к = arg max q-la-p, s = arg max a-}pj.
lj,n j j ljn jk j
4.1. Задача с дополнительным ограничением-неравенством. Пусть имеются матрицы A, B G Xnxn и векторы p, q, g G Xn. Рассмотрим задачу, которая получена добавлением к задаче (5) ограничения Bx ф g ^ x и имеет вид
max q-x (Ax) p, Bx ф g < x.
Докажем следующую теорему.
Теорема 2. Пусть матрицы A и B такие, что Tr(B) и матрица D = AB* имеет регулярные столбцы dj = (dij), где B* = (b*) — матрица Клини. Предположим, что векторы p = (pj) и q = (qj) регулярны.
Тогда максимум в задаче (6) равен Д = q- B* D-p и достигается тогда и только тогда, когда x = B*u, где u = (uj) — любой вектор с компонентами
Uk = ad-klps, gj < uj < ad-ps, j = к,
при условии, что
а ^ max gj dsj q-1, к = arg max q-b*d-p, s = arg max d-}pj.
Доказательство. Сначала заметим, что для неравенства Bx фд ^ x выполняется условие леммы 2. Заменим это неравенство на его решение x = B*u, где u — любой регулярный вектор такой, что u ^ д.
В результате подстановки этого решения с заменой AB* на D получим задачу относительно неизвестного вектора u в следующем виде:
max q-B*u (Du) p, , s
ч \ ) и, (7)
u > g.
Рассмотрим задачу (7) без учета ограничения u ^ д. Матрица D = AB* имеет регулярные столбцы. Из регулярности вектора q следует, что вектор q-B* регулярен. Применяя теорему 1, получим, что максимум в такой задаче равен Д = q-B* D-p и достигается тогда и только тогда, когда вектор u имеет компоненты
uk = ad-^ps, Uj < ad-jps, j = к, (8)
для всех а > 0 и индексов к и s, определяемых условиями
к = arg max q-b*d-p, s = arg max d-}pj. (9)
Выясним, при каких ограничениях на параметр а вектор u будет удовлетворять неравенству u ^ д. Учитывая формулы (8), запишем условия существования такого вектора u в форме
ad-J1ps > gj, j = 1,...,n.
Поскольку матрица D имеет только регулярные столбцы, а вектор p является регулярным, скаляр dsjp-1 не равен нулю. Умножив на dsjp-1 обе части последнего неравенства с учетом монотонности умножения, получим неравенство а ^ gjdsjp-1.
Тогда при условии, что
а ^ max gj dsjp-1,
l^j^n
вектор u, компоненты которого находятся по формулам (8), будет удовлетворять неравенству u ^ д, а вектор x = B*u будет решением задачи (6). □
Рассмотрим матрицу D и заметим, что для неразложимой матрицы B матрица Клини B* не имеет нулевых элементов. Отсюда следует, что если матрица A регулярна по строкам, то матрица D = AB* будет иметь регулярные столбцы. В таком случае для выполнения условий теоремы 2 достаточно неразложимости матрицы B и регулярности по строкам матрицы A. Наконец, если матрица A имеет только регулярные столбцы, то столбцы матрицы D будут регулярными при любой матрице B.
4.2. Задача с двумя ограничениями-неравенствами. Предположим, что заданы матрицы A, B, C G Xnxn и векторы p, q, g, h G Xn. Добавим к задаче (5) два ограничения-неравенства Bx ^ x и Cx ^ h, и запишем полученную задачу так:
max q-x(Ax)-p,
v 7 (10) Bx ^ x, Cx ^ h.
Эту задачу можно решить при помощи следующей теоремы. 390 Вестник СПбГУ. Прикладная математика. Информатика... 2017. Т. 13. Вып. 4
Теорема 3. Предположим, что матрицы A и B такие, что Tr(B) ^ 1, матрица D = AB* имеет регулярные столбцы dj = (dij), где B* = (b*) — матрица Клини, а матрица C является регулярной по столбцам. Пусть векторы p = (pj), q и h регулярные, а вектор r = (rj) определен как r = (h-CB*)-.
Тогда максимум в задаче (10) равен Д = q-B*D-p и достигается тогда, когда x = B*u, где u = (uj) — любой вектор с компонентами
Uk = ad-klps, uj < ad-j1ps, j = к,
при условии, что
a < min dsjp-lrj, к = arg max q-b*d-p, s = arg max d-}pj.
Kjj<n s Kjj<n j j Kjj<n jk
Доказательство. Так же, как при доказательстве предыдущего утверждения, заменим неравенство Bx ^ x на его решение x = B*u, где u — любой регулярный вектор подходящего размера. Подставляя решение, получим задачу относительно u в виде
max q-B*u (Du) p, CB*u < h.
Затем применим теорему 1 для нахождения максимума в задаче без ограничений. Получаем минимум, который равен Д = q-B*D-q и достигается тогда и только тогда, когда вектор u имеет компоненты, которые определяются условиями (8), (9).
Определим, какие из полученных решений удовлетворяют неравенству CB*u ^ h. Поскольку матрица C является регулярной по столбцам и B* ^ I, матрица CB* также регулярна по столбцам. Тогда указанное неравенство удовлетворяет условию леммы 1, которая дает решение u ^ (h-CB*)- = r.
Подберем такие значения a, при которых любой вектор u, который задается условиями (8), будет удовлетворять неравенству u ^ r. Для этого достаточно, чтобы
для всех ] = 1,...,п выполнялись неравенства аЛ- рэ ^ гз. Решая их, получим для
3
всех ] неравенства а ^ ря Г, которые эквивалентны одному неравенству:
а < тт ¿эзр—гз.
Параметр а, удовлетворяющий этому условию, задает такие векторы и, при которых векторы х = В*и являются решениями задачи (10). □
5. Решение задач планирования. Приведем решение задач оптимального планирования (1) и (2), полученное путем их представления как задач тропической оптимизации.
Заметим, что запись задач (1) и (2) в обычных обозначениях включает только операции определения максимума, сложения и вычитания (вычисления, обратного по сложению). Поэтому представим эти задачи в терминах полуполя Ктах,+ .
5.1. Максимизация разброса времени начала работ. Рассмотрим задачу (1) максимизации разброса времени начала работ при ограничениях «старт-старт» и «ранний старт». С использованием операций полуполя КтаХ1+ задача принимает вид
max xi x
n
^^ bij xj ф gi ^ xi, i = 1
Определим соответствующие матрично-векторные обозначения B = (bij), g = (gi), x = (xi) и запишем задачу в векторной форме
max 1T xx-1,
(11)
Bx © g < x.
Сформулируем и докажем следующую лемму.
Лемма 3. Пусть матрица B = (bj) со столбцами bj = (bij) является неразложимой и Tr(B) ^ 1. Тогда максимум в задаче (11) равен Д = 1TB*(B*)-1, где B* — матрица Клини со столбцами b* = (bj), и достигается тогда и только тогда, когда x = B*u, где u = (uj) — любой вектор с компонентами
uk = a(bSkГ\ gj < uj < a(b*sj)-1, j = k, (12)
при условии, что
a ^ max gjb*sj, k = arg max 1Tb*(b*)-1, s = arg max (bjk)-1. (13)
l^j^n j l^j^n j j l^j^n j
Доказательство. Задача (11) имеет вид задачи (6), если положить p = 1, q = 1 и A = I. Проверим выполнение условий теоремы 2. Ясно, что векторы p и q — регулярные. Поскольку матрица B является неразложимой и Tr(B) ^ 1, матрица D = AB* = B* содержит только регулярные столбцы. Таким образом, условия теоремы выполнены.
Применяя теорему 2, получим, что максимум в задаче равен Д = 1TB*(B*)-1 и достигается тогда и только тогда, когда x = B*u, где u — вектор с компонентами (12), которые удовлетворяют условию (13). □
Рассмотрим проект, состоящий из n = 3 работ с ограничениями «старт-старт» и «ранний старт», заданными с использованием обозначения (D = -ж следующими матрицей и вектором:
-2 1 \ /2
B
-1
Убедимся, что Tr(B) и построим матрицу Клини B*. Сначала найдем матрицы
О 0 \ / -1 -2 1
B2 = | 1 -2 1 I , B3 = I 0 -1 2 -3 0 J \ -1 0-1
Затем вычислим Tr(B) = tr B © tr B2 © tr B3 = 0 = 1. Следовательно, условия леммы 3 выполнены. Матрица Клини имеет вид
0 -2 1
B* = I © B © B2 = I 1 0 2
-1 -3 0
Используя матрицу B*, построим матрицы
0 -1 1
(Б*)- = В*(В*)- = ( 2 0 3
-1 -2 0
а также вычислим максимум Д = 1ТВ*(В*)-1 = 3 целевой функции задачи.
Определим индексы к и в, а также нижнюю границу для параметра а, которые определяются условиями (13). Сначала найдем
1Т Ь*(Ь*)-1 = 2, 1Т Ь*(Ь*)-1 = 3, 1Т Ь*3(Ь*3 )-1 = 2,
откуда вытекает, что к = 2. Затем вычислим
(6*2)-1 =2, (6У-1 =0, (6*2)-1 = 3
и выберем в = 3.
Установим нижнюю границу для параметра а. После расчета
£1Ь31 = 1, 92Ъ*2 = -3, дзЪ*з = 0
находим максимум, равный 1, из чего следует, что а ^ 1.
Определим вектор х самых ранних сроков начала работ, соответствующих значению а = 1. Сначала вычислим границы (12) для элементов вектора и = (и1,и2,из)Т при условии, что а = 1. После соответствующих расчетов получаем
и1 = 2, и2 = 4, 0 < и3 < 1.
Положим из = 0 и по вектору и = (2, 4,0)Т найдем вектор
х = В и = 24
Учитывая, что при из = 1 решение остается прежним, заключаем, что наиболее раннее время начала работ задается равенствами
Х1 = 2, Х2 = 4, хз = 1.
5.2. Максимизация разброса времени завершения работ. Представим задачу (2) максимизации разброса времени завершения работ с ограничениями «старт-старт», «старт-финиш» и «поздний финиш» в терминах полуполя КтаХ1+ в форме
тах (&У-
1
^^ а^зхз — ^^ Ь3х3 ^ xi, Уг ^ ^г, ^ — 1,...,п. з=1 з=1
Введем обозначения
А = (аз), В = (Ьз ), х = (хг), у = (Уг), Ъ = (Ы).
Представим задачу в виде
max 1T yy-1,
УУ ' (14)
Ax = y, Bx ^ x, y ^ h.
Решение задачи (14) можно записать следующим образом.
Лемма 4. Предположим, что матрицы A и B такие, что матрица A является регулярной, Tr(B) ^ 1, а матрица D = AB* имеет регулярные столбцы dj = (dj), где B* — матрица Клини.
Тогда максимум в задаче (14) равен Д = 1TDD-1 и достигается тогда, когда x = B*u и y = Du, где u = (uj) — любой вектор с компонентами
Uk = ad-,1, Uj < ad-j, j = к, (15)
при условии, что
a < min dsj(h-dn)-, к = arg max 1Tdnd-1, s = arg max d-. (1«)
i<j<n j i<j<n j j i^j^n jk y J
Доказательство. Подставим равенство y = Ax в целевую функцию и первое ограничение. Получим задачу относительно x в виде
max 1T Ax(Ax)-1,
Bx ^ x, Ax ^ h.
Она сводится к задаче (10) при условии, что C = A, p = 1 и q = (1TA)-. Проверим, что выполняются условия теоремы 3. Для матрицы B справедливо неравенство Tr(B) ^ 1, все столбцы матрицы D = AB* регулярны и матрица C является регулярной. Определим вектор r = (h-AB*)- = (h-D)-. Учитывая, что векторы p и q регулярные, заключаем, что указанные условия выполняются.
Применяя теорему 3, получим, что максимум в задаче равен Д = 1TDD-1 и достигается тогда, когда x = B*u, где u — любой вектор с компонентами (15), удовлетворяющими условиям (16). □
Пусть имеется проект, который состоит из n = 3 работ. Предположим, что для этих работ выполняются ограничения «старт-старт», «старт-финиш» и «поздний финиш», заданные при помощи следующих матриц и векторов:
4 1 1 \ /®-2 1\ /5
A = | 2 2 0 I, B = I 0 0 2 I, h = I 4 0 13/ \ -1 00/ у 4
Сначала заметим, что матрица B определена так же, как в предыдущем примере. Из результатов вычислений для этого примера следует, что Tr(B) = 1. Построим матрицы
/0 -2 1 \ /4
B* = | 1 0 2 I, D = AB* = | 3 \ -1 -3 0 J \ 2
Затем вычислим матрицы
-4 -3 -2 \ /012
D = | -2 -2 -1 I, DD- = | 0 0 1 -5 -4 -3 / \ -1 -1 0
и определим максимум в задаче, который равен Д = 1TDD-1 = 2.
Найдем индексы к и в, а также нижнюю границу для параметра а, которые определяются условиями (16). Сначала вычислим
1Т = 2, 1Т = 1, 1Т ё3ё-1 = 2,
из чего вытекает, что к может принимать два значения: 1 и 3. Сначала положим к = 1. Учитывая, что
¿-11 = -4, ¿-11 = -3, Л-11 = -2,
выбираем в = 3. Аналогично проверяется, что при к = 3 следует снова взять в = 3.
Определим верхнюю границу для параметра а. Вычислим вспомогательный вектор
Н-В = ( -1 -2 0 )
и рассмотрим величины
¿31(Н-^)- = 3, ¿32 (Л.-^)- = 3, ¿зз(Н-йз)- = 3.
Учитывая, что минимум этих величин равен 3, заключаем, что а ^ 3.
Чтобы получить решение, которое соответствует наиболее позднему времени начала работ, положим а = 3. Определим компоненты вектора и = (и1,и2,из)Т при условии, что к =1. По формулам (15) находим
и1 = 1, и2 ^ 2, из ^ 0.
В случае, когда к = 3, компоненты вектора и удовлетворяют условиям
и1 ^ 1, и2 ^ 2, из = 0.
Выбирая вектор и = (1, 2,0)Т, получаем решение в виде
1 \ ( 5' х = В*и = ( 2 I , у = Аи = ( 4
03
6. Заключение. В статье изучены задачи оптимального планирования сроков выполнения работ проекта при наличии ограничений на время начала и завершения, а также порядок работ. В качестве критериев оптимальности плана берется минимум максимального разброса времени начала или завершения работ. Такие задачи появляются, например, при необходимости составления плана, который позволяет избежать одновременного начала или завершения всех работ проекта.
Для решения использованы методы тропической оптимизации (методы оптимизации функций, заданных на векторах над идемпотентными полуполями). Получены новые решения задач планирования, которые, в отличие от существующих алгоритмических решений, дают результат в виде прямых расчетных формул в компактном векторном представлении, удобных для дальнейшего анализа и практических расчетов.
Полученные результаты демонстрируют возможность применения методов тропической оптимизации в качестве удобного инструмента решения задач планирования. Разработка приложений тропической оптимизации для решения новых задач сетевого планирования, включая задачи с другими типами ограничений и критериями
оптимальности, представляется перспективным направлением дальнейших исследований.
Литература
1. Kerzner H. Project management. Hoboken: Wiley, 2013. 1296 p.
2. Cuninghame-Green R. A. Projections in minimax algebra // Math. Program. 1976. Vol. 10. P. 111-123. DOI: 10.1007/BF01580656
3. Zimmermann U. Linear and combinatorial optimization in ordered algebraic structures. Vol. 10 of Annals of Discrete Mathematics. Amsterdam: Elsevier, 1981. 390 p. DOI: 10.1002/net.3230130312
4. Butkovic P., Aminu A. Introduction to max-linear programming // IMA J. Manag. Math. 2009. Vol. 20, N 3. P. 233-249. DOI: 10.1093/imaman/dpn029
5. Маслов В. П., Колокольцев В. Н. Идемпотентный анализ и его применение в оптимальном управлении. М.: Физматлит, 1994. 144 с.
6. Golem J. S. Semirings and affine equations over them. Vol. 556 of Mathematics and Its Applications. New York: Springer, 2003. 256 p. DOI: 10.1007/978-94-017-0383-3
7. Heidergott B., Olsder G. J., van der Woude J. Max Plus at work. Princeton Series in Applied Mathematics. Princeton: Princeton University Press, 2006. 226 p.
8. Кривулин Н. К. Методы идемпотентной алгебры в задачах моделирования и анализа сложных систем. СПб.: Изд-во С.-Петерб. ун-та, 2009. 256 с.
9. Butkovic P. Max-linear systems. Springer Monographs in Mathematics. London: Springer, 2010. 272 p. DOI: 10.1007/978-1-84996-299-5
10. Krivulin N. Explicit solution of a tropical optimization problem with application to project scheduling // Mathematical Methods and Optimization Techniques in Engineering / eds: D. Biolek, H. Walter, I. Utu, C. von Lucken. Athens: WSEAS Press, 2013. P. 39-45.
11. Krivulin N. 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. 840624
12. Krivulin N. A constrained tropical optimization problem: complete solution and application example // Tropical and Idempotent Mathematics and Applications. Vol. 616 of Contemporary Mathematics / eds: G. L. Litvinov, S. N. Sergeev. Providence: AMS, 2014. P. 163-177. DOI: 10.1090/conm/616/12308
13. Krivulin N. A maximization problem in tropical mathematics: a complete solution and application examples // Informatica. 2016. Vol. 27, N 3. P. 587-606. DOI: 10.15388/Informatica.2016.101
14. Кривулин Н. К., Губанов С. А. Решение задачи сетевого планирования на основе методов тропической оптимизации // Вестн. С.-Петерб. ун-та. Сер. 10. Прикладная математика. Информатика. Процессы управления. 2016. Вып. 3. С. 62-72. DOI: 10.21638/11701/spbu10.2016.306
15. Butkovic P., Tam K. P. On some properties of the image set of a max-linear mapping // Tropical and Idempotent Mathematics. Vol. 495 of Contemporary Mathematics / eds: G. L. Litvinov, S. N. Sergeev. Providence: AMS, 2009. P. 115-126. DOI: 10.1090/conm/495/09694
Для цитирования: Кривулин Н. К., Губанов С. А. Использование методов тропической оптимизации в задачах сетевого планирования // Вестник Санкт-Петербургского университета. Прикладная математика. Информатика. Процессы управления. 2017. Т. 13. Вып. 4. С. 384-397. https://doi.org/10.21638/11701/ spbu10.2017.405
References
1. Kerzner H. Project management. Hoboken, Wiley, 2013, 1296 p.
2. Cuninghame-Green R. A. Projections in minimax algebra. Math. Program., 1976, vol. 10, pp. 111123. DOI: 10.1007/BF01580656
3. Zimmermann U. Linear and combinatorial optimization in ordered algebraic structures. Vol. 10 of Annals of Discrete Mathematics. Amsterdam, Elsevier, 1981, 390 p. DOI: 10.1002/net.3230130312
4. Butkovic P., Aminu A. Introduction to max-linear programming. IMA J. Manag. Math., 2009, vol. 20, no. 3, pp. 233-249. DOI: 10.1093/imaman/dpn029
5. 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 Press, 1994, 144 p. (In Russian)
6. Golan J. S. Semirings and affine equations over them. Vol. 556 of Mathematics and Its Applications. New York, Springer, 2003, 256 p. DOI: 10.1007/978-94-017-0383-3
7. Heidergott B, Olsder G. J., van der Woude J. Max Plus at work. Princeton Series in Applied Mathematics. Princeton, Princeton University Press, 2006, 226 p.
8. Krivulin N. K. Metody idempotentnoj algebry v zadachakh modelirovaniya i analiza slozhnykh sistem [Methods of idempotent algebra, for problems in modeling and analysis of complex systems]. Saint Petersburg, Saint Petersburg State University Publ., 2009, 256 p. (In Russian)
9. ButkoviC P. Max-linear systems. Springer Monographs in Mathematics. London, Springer, 2010, 272 p. DOI: 10.1007/978-1-84996-299-5
10. Krivulin N. Explicit solution of a tropical optimization problem with application to project scheduling. Mathematical Methods and Optimization Techniques in Engineering. Eds. by D. Biolek, H. Walter, I. Utu, C. von Lucken. Athens, WSEAS Press, 2013, pp. 39-45.
11. Krivulin N. 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. 840624
12. Krivulin N. A constrained tropical optimization problem: complete solution and application example. Tropical and Idempotent Mathematics and Applications. Vol. 616 of Contemporary Mathematics. Eds by G. L. Litvinov, S. N. Sergeev. Providence, AMS Publ., 2014, pp. 163-177. DOI: 10.1090/conm/616/12308
13. Krivulin N. A maximization problem in tropical mathematics: a complete solution and application examples. Informatica, 2016, vol. 27, no. 3, pp. 587-606. DOI: 10.15388/Informatica.2016.101
14. Krivulin N. K., Gubanov S. A. Reshenie zadachi setevogo planirovanija na osnove metodov tropicheskoj optimizacii [Solution of a project scheduling problem by using methods of tropical mathematics]. Vestnik of Saint Petersburg University. Series 10. Applied Mathematics. Computer Science. Control Processes, 2016, iss. 3, pp. 62-72. (In Russian) DOI: 10.21638/11701/spbu10.2016.306
15. Butkovic P., Tam K. P. On some properties of the image set of a max-linear mapping. Tropical and Idempotent Mathematics. Vol. 495 of Contemporary Mathematics. Eds. by G. L. Litvinov, S. N. Sergeev. Providence, AMS Publ., 2009, pp. 115-126. DOI: 10.1090/conm/495/09694
For citation: Krivulin N. K., Gubanov S. A. The use of tropical optimization methods in problems of project scheduling. Vestnik of Saint Petersburg University. Applied Mathematics. Computer Science. Control Processes, 2017, vol. 13, iss. 4, pp. 384-397. https://doi.org/10.21638/11701/spbu10.2017.405
Статья рекомендована к печати проф. Л. А. Петросяном. Статья поступила в редакцию 29 июня 2017 г. Статья принята к печати 12 октября 2017 г.
