УДК 519.8+330.4 Вестник СПбГУ. Математика. Механика. Астрономия. 2021. Т. 8 (66). Вып. 1 МБС 15А80, 90С24, 90С47, 90В35
Алгебраическое решение задачи оптимального планирования сроков проекта в управлении проектами*
Н. К. Кривулин, С. А. Губанов
Санкт-Петербургский государственный университет,
Российская Федерация, 199034, Санкт-Петербург, Университетская наб., 7—9
Для цитирования: Кривулин Н. К., Губанов С. А. Алгебраическое решение задачи оптимального планирования сроков проекта в управлении проектами // Вестник Санкт-Петербургского университета. Математика. Механика. Астрономия. 2021. Т. 8(66). Вып. 1. С. 73-87. https://doi.org/10.21638/spbu01.2021.107
Рассматривается задача временного планирования проекта, который состоит в выполнении некоторого набора работ при заданных ограничениях на время начала и завершения работ. В качестве критерия оптимальности плана берется разброс времени начала выполнения работ, который требуется минимизировать. Такие задачи возникают в управлении проектами при необходимости по технологическим, организационным, экономическим или иным причинам обеспечить по возможности одновременное начало выполнения всех работ проекта. Рассматриваемая задача планирования формулируется как минимаксная задача оптимизации с ограничениями, а затем решается при помощи методов тропической (идемпотентной) математики, которая занимается вопросами теории и приложений полуколец с идемпотентным сложением. Сначала изучается задача тропической оптимизации, заданная в терминах общего идемпотент-ного полуполя (идемпотентного полукольца с обратимым умножением), и находится полное аналитическое решение этой задачи. Полученный результат затем используется для построения прямого решения задачи планирования в компактной векторной форме, удобной для дальнейшего анализа решений и непосредственных вычислений. Приводится иллюстративный численный пример решения задачи оптимального планирования проекта, состоящего из четырех работ.
Ключевые слова: идемпотентное полуполе, тропическая оптимизация, минимаксная задача оптимизации, временное планирование проекта, управление проектом.
1. Введение. Управление проектом заключается в согласовании действий, выполняемых для достижения целей проекта при рациональном распределении имеющихся ресурсов [1, 2]. Для решения задач оптимального планирования сроков выполнения проекта (составления оптимального графика работ) в управлении проектами обычно применяют методы сетевого планирования, которые опираются на анализ проекта с помощью графовых и сетевых моделей (метод критического пути, метод оценки и пересмотра планов), методы линейного и смешанного целочисленного линейного программирования и другие [3-6]. Как правило, указанные методы предлагают алгоритмические решения, которые позволяют получить численные ре-
* Работа выполнена при финансовой поддержке РФФИ в рамках научного проекта (грант №20010-00145).
(¡5 Санкт-Петербургский государственный университет, 2021
зультаты при помощи итерационных вычислительных процедур или установить, что решений нет.
Еще один подход к решению задач временного планирования проектов состоит в применении методов и результатов тропической оптимизации. Задачи тропической оптимизации формулируются и решаются в терминах тропической (идемпотент-ной) математики, которая занимается вопросами теории и приложений полуколец с идемпотентным сложением [7-12]. В таких задачах обычно требуется минимизировать (максимизировать) функции, определенные на векторах над идемпотентными полуполями (идемпотентными полукольцами, в которых операция умножения обратима), при ограничениях в виде тропических векторных уравнений и неравенств. Примеры некоторых задач тропической оптимизации можно найти, например, в [13].
Методы тропической оптимизации применяются при решении различных прикладных задач в информатике, исследовании операций, математической экономике и в других областях. Для ряда задач, включая задачи оптимального планирования сроков выполнения проекта, использование указанных методов позволяет находить новые аналитические решения в форме, удобной для последующего анализа и непосредственных вычислений [14-18].
В задачах временного планирования обычно рассматривают проект, который состоит в параллельном выполнении определенного набора работ при заданных ограничениях в форме отношений предшествования и допустимых границ для времени начала и завершения работ. Используются различные критерии оптимальности, включая разброс времени начала или завершения работ, который требуется минимизировать (или, наоборот, максимизировать). Подобные критерии возникают в задачах планирования при необходимости по технологическим, организационным, экономическим или иным причинам, по возможности, обеспечить (или, наоборот, не допустить) одновременное начало или завершение всех работ проекта. Результаты решения таких задач с помощью методов тропической оптимизации представлены в [18-20].
В настоящей работе изучается задача оптимального планирования сроков выполнения работ проекта, для которых установлены отношения предшествования «старт-старт» и «старт-финиш», а также заданы нижние границы для времени начала работ (сроки выпуска) и верхние границы для времени завершения работ (директивные сроки). В качестве критерия оптимальности, который следует минимизировать, берется максимальный разброс времени начала работ.
Сначала рассматриваемая задача планирования формулируется как минимаксная задача оптимизации с ограничениями. Затем изучается некоторая задача тропической оптимизации, заданная в терминах общего идемпотентного полуполя, для которой находится полное аналитическое решение в компактной векторной форме. Полученный результат применяется для решения исходной задачи оптимального планирования, которая для этого представляется в виде задачи тропической оптимизации. Приводится иллюстративный численный пример решения задачи оптимального планирования проекта, состоящего из четырех работ.
2. Задача временного планирования сроков проекта. Рассмотрим следующую задачу временного планирования проекта (обзор задач планирования и методов их решения можно найти, например, в [3-5]). Пусть проект состоит из п работ, которые выполняются параллельно с учетом ограничений на время начала и завершения каждой работы в форме отношений предшествования «старт-старт»
и «старт-финиш», и границ для самого раннего допустимого времени начала работ (сроки выпуска) и самого позднего допустимого времени завершения работ (директивные сроки). Требуется определить время начала каждой работы так, чтобы с учетом ограничений минимизировать разброс времени начала работ.
Также, как в работах [14-18], введем обозначения, которые позволяют записать задачу планирования в виде, удобном для дальнейшего представления в форме задачи тропической оптимизации. Для каждой работы i = 1,..., n обозначим время начала и время завершения через xj и yj. Пусть gj указывает самое раннее допустимое время начала работы, а fj — наиболее позднее допустимое время завершения работы, которые задают для величин xj и yj нижнюю и верхнюю границы:
gi < Xj, yj < fj. (1)
Определим ограничения «старт-старт» для работы i в виде неравенств bj + x, < xj для всех j = 1,..., n, где bjj обозначает минимально допустимый временной интервал между началом работы i и началом работы j. Если величина bjj не задана, то положим bjj = —то. Объединение неравенств по всем j дает эквивалентное неравенство
max (bjj + ж,) < xj. (2)
1<j<n
Чтобы задать ограничения «старт-финиш», обозначим минимальный допустимый интервал между временем начала работы i и завершением работы j через cjj (положим cjj = —то, если интервал не задан). Тогда все отношения предшествования «старт-финиш» для работы i можно представить с помощью неравенства
j (cjj + ж,) < yj.
Учитывая, что должно выполняться неравенство yj < fj, получаем ограничение
max (cjj + ж,) < fj. (3)
1<j<n
Критерием оптимальности плана является максимальный разброс времени начала всех работ, который должен быть минимальным. Представим этот критерий в виде
max xj — min xj = max xj + max (—xj). (4)
1<j<n 1<j<n 1<j<n 1<j<n
Сформулируем задачу составления оптимального плана работ в соответствии с указанными ограничениями и введенным критерием: при заданных значениях параметров bjj, cjj, gj и fj, где i, j = 1,..., n, найти величины xj, которые обеспечивают
min max xj + max (—xj);
1<j<n 1<j<n
max (bj + x,) < xj, max (cjj + x,) < fj, (5)
1<j<n 1<j<n
gj < xj, i = 1,. .., n.
Нетрудно понять, что задача (5) может быть записана в форме задачи линейного программирования и решена в численном виде с помощью известных вычислительных процедур, таких как симплекс-алгоритм [21-23] и алгоритм Кармаркара
[24]. В следующих разделах для построения прямого аналитического решения рассматриваемой задачи будут использованы методы и результаты тропической оптимизации.
3. Элементы тропической алгебры. В этом разделе приводится краткий обзор основных алгебраических определений, обозначений и предварительных результатов, которые используются ниже для построения решения задачи тропической оптимизации и ее приложения к задаче планирования. Дальнейшие сведения по теории и приложениям тропической (идемпотентной) математики можно найти, например, в работах [7-12].
Пусть множество X замкнуто относительно ассоциативных и коммутативных операций сложения ф и умножения ( и включает их нейтральные элементы нуль 0 и единицу 1. Сложение обладает свойством идемпотентности, что означает выполнение равенства x фx = x для любого x G X, а умножение является дистрибутивным относительно сложения и обратимым, т. е. для любого ненулевого x существует элемент x-1 такой, что x ( x-1 = 1. Алгебраическую систему (X, 0,1, ф, () обычно называют идемпотентным полуполем. Знак умножения ( далее для удобства будет опускаться.
Операция сложения ф задает на X частичный порядок по правилу: x < y тогда и только тогда, когда xфy = y. В терминах этого порядка для сложения и умножения выполняются свойства монотонности, по которым для любых x, y, z G X при условии x < y справедливы неравенства xфz < yфz и xz < yz. Операция обращения является антитонной: для всех x, y = 0 из неравенства x < y следует x-1 > y-1.
Сложение имеет экстремальное свойство, которое состоит в том, что для любых x, y G X выполняются неравенства x < x ф y и y < x ф y. Кроме того, неравенство x ф y < z эквивалентно системе двух неравенств x < z и y < z. Далее будем считать, что указанный частичный порядок дополнен на X до линейного порядка.
Для каждого x = 0 и натурального p целая степень определяется обычным путем: x0 = 1, xp = xp-1x, x-p = (x-1 )p и 0p = 0. Предполагается, что операция возведения в степень естественным образом распространяется на случай рационального показателя.
Для любых x, y G X и рационального числа а > 0 справедливо биномиальное тождество (x ф y)a = xa ф ya, из которого вытекает тропический аналог неравенства между арифметическим и геометрическим средними в форме x ф y > (xy)1/2. Эти тождество и неравенство прямо обобщаются на случай произвольного числа слагаемых.
Примером вещественного полуполя, которое обычно называют (max, +)-алгеброй, является алгебраическая система Rmax,+ = (R U {-то}, —то, 0, max, +), где сложение ф определено как max, умножение ф как +, в роли нуля 0 выступает —то, в роли единицы 1 — число 0. Для любого числа x G R существует обратное x-1, которому соответствует противоположное число —x в обычной алгебре. Для любых x, y G R степени xy соответствует арифметическое произведение xy. Определяемое идемпотентным сложением отношение порядка соответствует обычному линейному порядку на R.
Обозначим через Xmxn множество матриц, которые состоят из m строк и n столбцов с элементами из X. Матрица, все элементы которой равны 0, является нулевой. Матрица, которая не имеет нулевых столбцов, называется регулярной по столбцам.
Сложение и умножение матриц, а также умножение матрицы на скаляр выполняются по обычным правилам с заменой арифметических операций + и х на операции ф и (. Отношение порядка, введенное выше, и связанные с ним свойства операций обобщаются на матрицы и понимаются покомпонентно.
Рассмотрим квадратные матрицы из множества Хпхп. Единичной является матрица I с элементами, равными 1 на главной диагонали и 0 вне ее.
Для квадратной матрицы А и натурального р неотрицательная целая степень определяется как обычно: А0 = I, Ар = Ар—1 А. Для степени суммы матриц справедливо биномиальное тождество
р
(А ф В)р = Вр ф0 0 В®0 (АВ®1 ••• АВ®к).
к=1 ®0+-----+ ®к =р—к
След матрицы А = (о,) € Хпхп вычисляется по формуле ^ А = а11 ф • • • ф а„„. Для любых подходящих по размеру матриц А и В справедливы равенства:
^(А ф В) = ^ А ф ^ В, ^(АВ) = ^(ВА), ^(жА) = ж tr А.
Для любой матрицы А € Хпхп можно ввести тропические аналоги определителя и спектрального радиуса матрицы, которые вычисляются соответственно по формулам
Тг(А) А ф---ф tr А", Л = ^ А ф — ф ^1/п (А").
Если выполняется условие Тг(А) < 1, то определена матрица Клини в виде
А* = I ф А ф---ф А"-1.
Применяя биномиальное тождество для матриц, можно представить матрицу Клини для суммы матриц А, В € Хпхп следующим образом:
п—1п—к—1
(А ф В)* = В* ф0 0 0 В®0 (АВ®1 ••• АВ®к). (6)
к=1 т=0 ®о +®1 +-----+®к =т
Обозначим множество векторов-столбцов из п элементов через X". Вектор, не имеющий нулевых компонент, называется регулярным.
Для любого вектора х транспонированный вектор обозначается через хт. Для ненулевого вектора-столбца х = (ж®) € X" определим его мультипликативно сопряженный вектор-строку х- = (ж—), где ж— = ж—1, если ж® = 0, и ж— = 0 в противном случае. Вектор, состоящий из единиц, обозначим через 1 = (1,..., 1)т.
Для вектора ж = (ж®) € X" и матрицы А = (о,) € Хтхп введем тропические аналоги нормы вектора и матрицы с помощью формул
п т п
УхУ = 0 ж® = 1Тх = хТ 1, ||А|| = 0 0^- = 1ТА1.
®=1 ®=1 ,7=1
4. Задачи тропической оптимизации. Задачи оптимизации, которые формулируются и решаются в терминах тропической алгебры (задачи тропической оптимизации), находят применение в различных областях техники, экономики и управления. С помощью методов тропической оптимизации нередко можно построить
новое прямое аналитическое решение задач, для которых известны только алгоритмические численные решения, включая задачи оптимального планирования сроков выполнения проектов (см., например, [13-18]).
Рассмотрим пример задачи тропической оптимизации и ее решения. Пусть заданы матрицы A, B, C G Жпхп и векторы g, f G Xn. Требуется найти регулярные векторы x G Xn, которые решают задачу
min x—Ax;
(7)
Bx < x, Cx < f, g < x.
Для описания решения задачи (7) введем обозначения:
Sk = 0 AB®1 ••• ABik, k = 1,..., n; (8)
0<ii +-----+ik <n —k
Tk = 0 Bi0 (AB®1 ••• ABik), k = 1,..., n - 1. (9)
0<®o +ii +-----+ifc <n—k—1
Решение задачи можно получить как прямое следствие результатов, полученных в [14, 16], и представить в следующей форме.
Теорема. Пусть A — матрица со спектральным радиусом Л > 0, а B — матрица, для которой Tr(B) < 1. Пусть C — регулярная по столбцам матрица, g — вектор, а f — регулярный вектор такие, что f—CB*g < 1.
Тогда минимальное значение целевой функции в задаче (7) равно
п п— 1
0 = 0 tr1/k(Sk) 0 0(f—CTkg)1/k, (10)
k=1 k=1
а все регулярные 'решения записываются в параметрической форме
x = (0—1A ф B)*u, (11)
где u — регулярный вектор параметров, который удовлетворяет условиям
g < u < (f—C(0—1A ф B)*)—. (12)
Заметим, что условие Tr(B) < 1 имеет смысл требования наличия регулярных решений x у неравенства Bx < x. Условие f— CB*g < 1 обеспечивает существование общих регулярных решений для всех неравенств системы ограничений задачи.
Чтобы применить результат теоремы для решения задачи планирования, потребуется изучить частный случай задачи (7). Положим A = 11T и рассмотрим задачу:
min x— 11T x;
(13)
Bx < x, Cx < f, g < x.
В следующем разделе для задачи (13) будет получено прямое решение, представленное в компактной векторной форме.
5. Решение задачи оптимизации. Для решения задачи тропической оптимизации (13) будем использовать результат вышеприведенной теоремы. Полученное решение затем представим в более компактном виде, опираясь на особую форму целевой функции задачи. Справедливо следующее утверждение.
Лемма. Пусть В — матрица, для которой Тг(В) < 1. Пусть С — регулярная по столбцам матрица, д — вектор, а f — регулярный вектор такие, что f— СВ*д < 1.
Тогда минимальное значение целевой функции в задаче (13) равно
в = ||В*||ф 0 — СВ®|| ||В7 д||, (14)
0<®+,<п —2
а все регулярные 'решения имеют вид
х = Ои, О = В* ф 0 в—1В®11тВ7, (15)
0<®+,<п—2
где и — любой регулярный вектор, который удовлетворяет условиям
д < и < (f—со)—. (16)
Доказательство. В силу того, что задача (13) имеет вид (7), применим теорему, приведенную выше. Для этого сначала найдем минимум в по формуле (10). Подставляя А = 11т в выражение АВ®1 • • • АВ®к и записывая скалярные множители в виде нормы матриц, получим
11тВ®1 ••• 11тВ®к = (ЦВ®11| • • • |В®к-1 ||)11тВ®к. (17)
Тогда равенство (8) принимает вид
^ = 0 (ЦВ®11| • • • |В®к-1 ||)11т в®к.
0<®1 +-----<п— к
Найдем след матрицы под знаком суммы. С учетом свойств следа имеем ^(ЦВ®11| • • • |В®к-1 ||)11тВ®к) = ЦВ®11| • • • |В®к||. С помощью полученного результата запишем сумму
пп
0tr1/k (5к ) = 0 0 (ЦВ®11| • • • ||В®к ||)1/к. к=1 к=1 0<®1+-----+®к <"—к
Рассмотрим первое слагаемое, соответствующее к = 1 в суммах слева и справа,
ОД)= 0 ЦВ®1||.
0<®1<п —1
Покажем, что другие слагаемые по величине не больше первого. Действительно, применяя неравенство между геометрическим и арифметическим средними, получим
(ЦВ®11| •••|В®к||)1/к < |в®1 || ф — ф ||В®к||,
откуда следует, что при всех к = 2,..., п выполняется
п—к п—1
0 (ЦВ®11| • • • ||В®к||)1/к <0 ||В®| < 0 ||В®| = ОД).
0<®1 +-----+®к<п —к ®=0 ®=0
В силу того, что первое слагаемое доминирует над остальными, в результате имеем
п п— 1
01г1/к(£)=ОД) = 0 |^|| = ||В*||.
к=1 ¿=0
Теперь, используя равенство (17), запишем (9) в виде
Тк = 0 (|В || — |^к-1 ||^011т .
0^о+¿1+-----^¿к<п—к—1
После умножения матрицы под знаком суммы слева на f—С и справа на д и замены скалярных множителей на норму соответствующих векторов получим равенство
f—С(|В || • • • |^к-1 ||^011тд = —КНВ41|| • • • |^к-1 ||)|^кд||. С учетом этого равенства составим сумму
п—1 п—1
0(f—СТ*д)1/к = 0 0 —CB¿0КНВ41|| • • • ||В^-1 ||)НВ4кд||)1/к
к=1 к=1 0^0 +¿1+-----<п—к—1
и рассмотрим первое слагаемое, которое соответствует к = 1 и имеет вид f—СТ1д = 0 —|||Я д||.
<п—2
Проверим, что величина других слагаемых не больше суммы ^(£1) ф f—СТ1д. Применяя неравенство между геометрическим и арифметическим средними, получим
—CB¿0 НСНВ41|| • • • НЯ^-1 ||)|В4кд||)1/к < (|Я ||ф- • •ф^-1 ^ф^—CB¿0|| д||. Следовательно, для всех к = 2,..., п будет выполняться неравенство
0 (Нf—CB¿0||(|Я|| • • • ||В^-11|)|д||)1/к <
0^0 +¿1 +-----Ик <п—к—1
п—1
<0Мф 0 — CB¿||||Bjд|| = ОД) ф f—СТ1д.
¿=0 0^+;?<п—2
С учетом последнего соотношения имеем двойное неравенство
п—1
f—СТ1д < 0(f— СТ*д)1/к < ОД) ф f—СТ1д.
к=1
Чтобы получить (14), осталось объединить полученные результаты:
п п—1
О = 01г1/к(£) ф0а—СТкд)1/к = ||В*||ф 0 ||f— CB¿||||Bjд||.
к=1 к=1 0^+з<п—2
Перейдем к параметрическому представлению всех решений задачи на основе формул (11) и (12). Составим матрицу Клини (в-Ы ф В)* = (в-111т ф В)* и воспользуемся тождеством (6). Применяя, как и раньше, равенство (17) для преобразования выражения под знаком суммы, окончательно получим
п-1п-к-1
(0-111тфВ)* = в*ф0 0 0 в-тв4° 11Тв^(||в41 н^нв^-1 ||).
к=1 т=0 ¿0 + ---+4к=т
Исследуем сумму в правой части и рассмотрим слагаемое при к =1 в виде
п—2
0 0 в- 1В4011т В*1 =0 в-^¿11^.
т=0 ¿о +41=т 0<г+^'<п —2
Заметим, что из найденного выражения (14) для минимума целевой функции задачи следует, что для всех г = 0,..., п — 1 выполняются неравенства в > ||В*|| > ||ВгН и, в частности, что в > ||71| = 1 > 0. Учитывая, что тогда для всех г будут справедливы неравенства в-11 Вг | < 1, получим неравенство
в-1в40 11тв4к(в-1 |ВН ••• в-1 ||в4к-1 Н) < в-1в40 11тв4к.
С учетом последнего неравенства для всех к = 2,..., п — 1 выполняется
п-к-1
0 0 в-тв4°11Тв4к(||в41 || • • • нв4^-1 Н) < 0 в-1 в®иТв^,
т=0 ¿0 + ---+4к=т 0<г+^'<п-2
откуда следует, что слагаемое суммы для к = 1 доминирует над остальными, которые поэтому можно отбросить. Тогда матрица Клини приобретает вид
(в-111т ф В)* = в* ф 0 в-1в411т в^ = а.
0<г+^<п-2
Подставляя полученное выражение в формулы (11) и (12), приходим к параметрическому описанию решений задачи в виде (15) и (16). □
Оценим вычислительную сложность построенного решения. Сначала заметим, что матрица В* находится путем прямых матричных вычислений за 0(п4) скалярных операций или с помощью алгоритма Флойда — Уоршелла [22] за 0(п3) операций.
При условии, что имеются все необходимые степени матрицы В, определение значений ||В*|| и Тг(В) дополнительно потребует 0(п2) скалярных операций. Вычисление векторов f- СВ4 и В^д и их норм имеет вычислительную сложность порядка 0(п3).
Вычисление всех матриц В411тВ^ можно осуществить следующим образом: сначала найти все векторы В41 и 1тВ^, выполнив 0(п3) операций, а затем перемножить соответствующие векторы за не более чем 0(п4) операций. При этом общий порядок сложности таких вычислений (и всего решения в целом) не будет превосходить 0(п4).
6. Приложение к задаче оптимального планирования. Вернемся к задаче оптимального планирования (5) и запишем ее в терминах (тах, +)-алгебры (полуполя Ктах,+) в виде задачи тропической оптимизации. Чтобы сформулировать задачу
в векторной форме, введем матрицы и векторы:
B = (bj), C = (cj), x = (xi), g = (gi), f = (fi).
Используя операции (max, +)-алгебры для того, чтобы представить целевую функцию (4) в скалярной, а затем в векторной формах, получим
Xi Ш
„-1
x-11T x.
1<i<n 1<j<n
Ограничения (2), (3) и (1) в скалярной форме принимают вид
(Q^) bij< Xi, СО cijxj < gi < Xi, i = 1, . . . , n.
1<j<n 1<j<n
После перехода к векторным обозначениям имеем неравенства Bx < x, Cx < f, g < x.
Записывая вместе целевую функцию и ограничения, приходим к задаче
min x—11T x;
Bx < x, Cx < f, g < x.
Для решения полученной задачи с помощью вышеприведенной леммы, необходимо проверить условия леммы. Заметим, что в контексте задач временного планирования в управлении проектами элементы неизвестного вектора x определяют время начала работ, а элементы заданного вектора f показывают самое позднее время завершения работ. Поскольку указанное время естественно считать определенным (т.е. большим, чем 0 = -то), требование регулярности векторов x и f оказывается выполненным. Условие Tr(B) < 1 обеспечивает непротиворечивость отношений предшествования «старт-старт», а условие f— CB*g < 1 — непротиворечивость отношений «старт-старт», «старт-финиш» вместе с заданными сроками выпуска и директивными сроками работ.
Пусть условия леммы выполнены. Тогда ее применение позволяет вычислить минимальный разброс 0 времени начала работ по формуле (14). Все векторы x, которые задают время начала работ с минимальным разбросом, записываются в параметрической форме с помощью соотношений (15) и (16).
Для иллюстрации полученного результата и демонстрации вычислительной техники приведем пример решения задачи планирования проекта, состоящего из n = 4 работ.
Пример. Рассмотрим проект, в котором ограничения заданы с использованием обозначения 0 = —то при помощи следующих матриц и векторов:
B
/ 0 -1 0 0 ^ ( 4 0 0 0 0 ( 6 \
1 0 0 0 , C = 1 3 0 2 0 6
-2 -1 0 0 4 0 5 0 , g = 0 , f = 6
\ 0 -3 2 0 V 5 1 4 3 0 1 6 )
Чтобы применить результат леммы, сначала проверим выполнение ее условий. Вычисление степеней матрицы В дает результат:
0 -1 -1 0 0 -1 -1 0
В2 = 1 0 0 0 , В3 = В4 = 1 0 0 0
0 -1 0 0 0 -1 0 0
0 1 20 2 1 2 0
Теперь найдем следы ^ В = ^ В2 = ^ В3 = ^ В4 = 0. После в
величины
Тг(В) = ^ В ф ^ В2 ф tr В3 ф tr В4 = 0 заключаем, что условие Тг(В) < 1 = 0 выполняется. Далее находим
В* = 7 ф В ф В2 ф В3 =
В*д
/0 -1 -1 0 \
10 0 0
0 -1 0 0
\ 2 1 2 0 У
f-С = ( -1 -3 -1 -3 ) , f-СВ*д = -1,
0 1
0
2
откуда следует, что условие f-СВ*д = -1 < 1 также выполнено.
Определим минимальный разброс в времени начала работ. Сначала построим вспомогательные векторы:
f-СВ = f-СВ2 = ( -1 -2 -1 -3 ), Вд = В2д
0 1
0
2
Вычислим нормы ||В*|| = 2, -С|| = -СВ|| = ||f-СВ2|| = -1, а также нормы ||д|| = 0, ||Вд|| = ||В2д|| = 2. Окончательно находим в по формуле
в = ||В*|| ф -С||||д|| ф -СВ||||д|| ф -С||||Вд||ф
ф -СВ2||||д|| ф -СВ||||Вд|| ф -С|||В2д|
2.
Опишем все решения х задачи в параметрической форме. Вычисление вспомогательных матриц дает результат:
В11Т
2ч 1 Т
В211
В11ТВ
0 0 0 0 \ / 1 0 2 0
Т = 1 1 1 1 , 11ТВ 1 0 2 0
0 0 0 0 1 0 2 0
2 2 2 2 ) \ 1 0 2 0
( 1 0 2 0 \ 1 1 2 0 \
2 1 3 0 11ТВ2 = 1 1 2 0
1 0 2 0 1 1 2 0
3 2 4 2 1 1 2 0 )
Найдем следующие матрицу и вектор:
g = в* e (—2)(llT e в11т e llTв e в211т e в11тв e llTв2)
0 -l о -2
1 о 1 -1 о -l о -2
y 2 12 о y
(f—CG)
l 2 l
V3/
Теперь можно записать все решения х = (хх,х2,хз,х4)т с помощью вектора параметров и = (их, и2, из, и4)т в форме
о -l о
1 о l
о l о
о о о
u
l 2 l
V2 1 2 0/ \0/ V3/
Учитывая, что все столбцы матрицы С коллинеарны, ее можно представить в виде
G
о -l о -2 \ о
1 о 1 -1 1
о -1 о -2 о
2l 2о 2
о -l о -2 .
Введем скалярный параметр и = м©( — 1)и2©из®( —2)м4. После замены вектора параметров и на новый параметр и и вычисления границ для и получим решение
/ 0 \
1
0
2
о < u< 1.
В терминах обычных арифметических операций координаты вектора х записываются в виде
Xl = u,
X2 = u +1,
хз = u,
X4 = u + 2,
о < u< 1.
Авторы благодарят рецензентов за ряд важных замечаний и предложений, которые были учтены при подготовке окончательного варианта статьи.
x
u
Литература
1. Троцкий М., Груча Б., Огонек К. Управление проектами, пер. с польск. Москва, Финансы и статистика (2006).
2. Kerzner H. Project Management. 10th ed. Hoboken, Wiley (2010).
3. Pinedo M. L. Scheduling. 3rd ed. New York, Springer (2008).
4. Vanhoucke M. Project Management with Dynamic Scheduling. 2nd ed. Berlin, Springer (2012). https://doi.org/10.1007/978-3-642-40438-2
5. Blazewicz J., Ecker K. H., Pesch E., Schmidt G., Sterna M., Weglarz J. Handbook on Scheduling. In: International Handbooks on Information Systems. Cham, Springer (2019). https://doi.org/10.1007/978-3-319-99849-7
6. Бурков В. Н., Буркова И. В., Засканов В. Г. Метод сетевого программирования в задачах календарного планирования. Автомат. ителемех., (6), 17—28 (2020). https://doi.org/10.31857 /S0005231020060025
7. Baccelli F. L., Cohen G., Olsder G. J., Quadrat J.-P. ¡Synchronization and Linearity. In: Wiley Series in Probability and Statistics. Chichester, Wiley (1993).
8. Маслов В. П., Колокольцов В. Н. Идемпотентный анализ и его применение в оптимальном управлении. Москва, Физматлит (1994).
9. Golan J. S. Semirings and Affine Equations Over Them. In: Mathematics and Its Applications, vol. 556. New York, Springer (2003). https://doi.org/10.1007/978-94-017-0383-3
10. Heidergott B., Olsder G. J., van der Woude J. Max Plus at Work. In: Princeton Series in Applied Mathematics. Princeton, Princeton Univ. Press (2006).
11. Gondran M., Minoux M., Graphs, Dioids and Semirings. In: Operations Research/Computer Science Interfaces, vol.41. New York, Springer (2008). https://doi.org/10.1007/978-0-387-75450-5
12. Butkovic P. Max-linear Systems. In: Springer Monographs in Mathematics. London, Springer (2010). https://doi.org/10.1007/978-1-84996-299-5
13. Krivulin N. A multidimensional tropical optimization problem with nonlinear objective function and linear constraints. Optimization 64 (5), 1107-1129 (2015). https://doi.org/10.1080/02331934.2013. 840624
14. Krivulin N. A constrained tropical optimization problem: complete solution and application example. In: G. L. Litvinov, S. N. Sergeev (eds). Tropical and Idempotent Mathematics and Applications. Contemporary Mathematics, vol.616, 163-177. Providence, AMS (2014). https://doi.org/10.1090/conm/616/12308
15. Krivulin N. Extremal properties of tropical eigenvalues and solutions to tropical optimization problems. Linear Algebra Appl. 468, 211-232 (2015). https://doi.org/10.1016/j.laa.2014.06.044
16. Krivulin N. Direct solution to constrained tropical optimization problems with application to project scheduling. Comput. Manag. Sci. 14 (1), 91-113 (2017). https://doi.org/10.1007/s10287-016-0259-0
17. Krivulin N. Tropical optimization problems with application to project scheduling with minimum makespan. Ann. Oper. Res. 256 (1), 75-92 (2017). https://doi.org/10.1007/s10479-015-1939-9
18. Krivulin N. Tropical optimization problems in time-constrained project scheduling. Optimization 66 (2), 205-224 (2017). https://doi.org/10.1080/02331934.2016.1264946
19. Кривулин Н.К., Губанов С. А. Решение задачи сетевого планирования на основе методов тропической оптимизации. Вестник Санкт-Петербургского университета. Сер. 10. Прикладная математика. Информатика. Процессы управления 12, вып. 3, 62-72 (2016). https://doi.org/10.21638/11701/spbu10.2016.306
20. Кривулин Н. К., Губанов С. А. Использование методов тропической оптимизации в задачах сетевого планирования. Вестник Санкт-Петербургского университета. Прикладная математика. Информатика. Процессы управления 13, вып. 4, 384-397 (2017). https://doi.org/10.21638/11701/ spbu10.2017.405
21. Kantorovich L. V. Mathematical methods of organizing and planning production. Management Sci. 6 (4), 366-422 (1960). https://doi.org/10.1287/mnsc.6A366
22. Романовский И. В. Алгоритмы решения экстремальных задач. Москва, Наука (1977).
23. Васильев Ф.П., Иваницкий А. Ю. Линейное программирование. 3-е изд. Москва, Факториал Пресс (2008).
24. Karmarkar N. A new polynomial-time algorithm for linear programming. Combinatorica 4 (4), 373-395 (1984). https://doi.org/10.1007/BF02579150
Статья поступила в редакцию 30 мая 2020 г.;
после доработки 7 сентября 2020 г.; рекомендована в печать 17 сентября 2020 г.
Контактная информация:
Кривулин Николай Кимович — д-р физ.-мат. наук, проф.; nkk@math.spbu.ru
Губанов Сергей Александрович — аспирант; segubanov@mail.ru
Algebraic solution of a problem of optimal project scheduling in project management*
N. K. Krivulin, S. A. Gubanov
St. Petersburg State University, 7—9, Universitetskaya nab., St. Petersburg, 199034, Russian Federation
For citation: Krivulin N. K., Gubanov S. A. Algebraic solution of a problem of optimal project scheduling in project management. Vestnik of Saint Petersburg University. Mathematics. Mechanics. Astronomy, 2021, vol. 8(66), issue 1, pp. 73-87. https://doi.org/10.21638/spbu01.2021.107 (In Russian)
A problem of optimal scheduling is considered for a project that consists of a certain set of works to be performed under given constraints on the times of start and finish of the works. As the optimality criterion for scheduling, the maximum deviation of the start time of works is taken to be minimized. Such problems arise in project management when it is required, according to technological, organizational, economic or other reasons, to provide, wherever possible, simultaneous start of all works. The scheduling problem under consideration is formulated as a constrained minimax optimization problem and then solved using methods of tropical (idempotent) mathematics which deals with the theory and applications of semirings with idempotent addition. First, a tropical optimization problem is investigated defined in terms of a general idempotent semifield (an idempotent semiring with invertible multiplication), and a complete analytical solution of the problem is derived. The result obtained is then applied to find a direct solution of the scheduling problem in a compact vector form ready for further analysis of solutions and straightforward computations. As an illustration, a numerical example of solving optimal scheduling problem is given for a project that consists of four works.
Keywords: idempotent semifield, tropical optimization, minimax optimization problem, project scheduling, project management.
References
1. Trocki M., Grucza B., Ogonek K. Zarazqdzanie Projectami. Warszawa, PWE (2003). (In Polish) [Russ. ed.: Upravlenie proektami. Moscow, Finansy i statistika Publ. (2006)].
2. Kerzner H. Project Management. 10th ed. Hoboken, Wiley (2010).
3. Pinedo M.L. Scheduling. 3rd ed. New York, Springer (2008).
4. Vanhoucke M. Project Management with Dynamic Scheduling. 2nd ed. Berlin, Springer (2012). https: //doi.org/10.1007/978-3-642-40438-2
5. Blazewicz J., Ecker K.H., Pesch E., Schmidt G., Sterna M., Weglarz J. Handbook on Scheduling. In: International Handbooks on Information Systems. Cham, Springer (2019). https://doi.org/10.1007/978-3-319-99849-7
6. Burkov V. N., Burkova I.V., Zaskanov V. G. The network programming method in project scheduling problems. Avtomat. i Telemekh., (6), 17—28 (2020). https://doi.org/10.31857 /S0005231020060025 (In Russian) [Engl. transl.: Autom. Remote Control 81 (6), 978-987 (2020). https: //doi.org/10.1134/S000511792006003X].
7. Baccelli F. L., Cohen G., Olsder G. J., Quadrat J.-P. Synchronization and Linearity. In: Wiley Series in Probability and Statistics. Chichester, Wiley (1993).
8. Maslov V. P., Kolokoltsov V.N. Idempotent analysis and its applications to optimal control theory. Moscow, Nauka Publ. (1994). (In Russian)
9. Golan J. S. Semirings and Affine Equations Over Them. In: Mathematics and Its Applications, vol.556. New York, Springer (2003). https://doi.org/10.1007/978-94-017-0383-3
10. Heidergott B., Olsder G. J., van der Woude J. Max Plus at Work. In: Princeton Series in Applied Mathematics. Princeton, Princeton Univ. Press (2006).
11. Gondran M., Minoux M., Graphs, Dioids and Semirings. In: Operations Research/Computer Science Interfaces, vol.41. New York, Springer (2008). https://doi.org/10.1007/978-0-387-75450-5
*This work is supported by Russian Foundation for Basic Research (grant No. 20-010-00145).
12. Butkovic P. Max-linear Systems. In: Springer Monographs in Mathematics. London, Springer (2010). https://doi.org/10.1007/978-1-84996-299-5
13. Krivulin N. A multidimensional tropical optimization problem with nonlinear objective function and linear constraints. Optimization 64 (5), 1107-1129 (2015). https://doi.org/10.1080/02331934.2013. 840624
14. Krivulin N. A constrained tropical optimization problem: complete solution and application example. In: G. L. Litvinov, S. N. Sergeev (eds). Tropical and Idempotent Mathematics and Applications. Contemporary Mathematics, vol.616, 163-177. Providence, AMS (2014). https://doi.org/10.1090/conm/616/12308
15. Krivulin N. Extremal properties of tropical eigenvalues and solutions to tropical optimization problems. Linear Algebra Appl. 468, 211-232 (2015). https://doi.org/10.1016/j.laa.2014.06.044
16. Krivulin N. Direct solution to constrained tropical optimization problems with application to project scheduling. Comput. Manag. Sci. 14 (1), 91-113 (2017). https://doi.org/10.1007/s10287-016-0259-0
17. Krivulin N. Tropical optimization problems with application to project scheduling with minimum makespan. Ann. Oper. Res. 256 (1), 75-92 (2017). https://doi.org/10.1007/s10479-015-1939-9
18. Krivulin N. Tropical optimization problems in time-constrained project scheduling. Optimization 66 (2), 205-224 (2017). https://doi.org/10.1080/02331934.2016.1264946
19. Krivulin N. K., Gubanov S. A. 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 12, iss. 3, 62-72 (2016). https://doi.org/10.21638/11701/spbu10.2016.306 (In Russian)
20. 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 13, iss. 4, 384-397 (2017). https://doi.org/10.21638/11701/spbu10.2017.405 (In Russian)
21. Kantorovich L. V. Mathematical methods of organizing and planning production. Management Sci. 6 (4), 366-422 (1960). https://doi.org/10.1287/mnsc.6A366
22. Romanovskiy I. V. Algorithms for Solving Extremal Problems. Moscow, Nauka Publ. (1977). (In Russian)
23. Vasilyev F. P., Ivanitskiy A. Y. Linear Programming. 3rd ed. Moscow, Factorial Press (2008). (In Russian)
24. Karmarkar N. A new polynomial-time algorithm for linear programming. Combinatorica 4 (4), 373-395 (1984). https://doi.org/10.1007/BF02579150
Received: May 30, 2020 Revised: September 7, 2020 Accepted: September 17, 2020
Authors' information:
Nikolay K. Krivulin — nkk@math.spbu.ru Sergey A. Gubanov —segubanov@mail.ru