CC BY
69
Математическое моделирование. Оптимальное управление Вестник Нижегородского университета им. Н.И. Лобачевского, 2012, № 3 (1), с. 163-169
УДК 519.852
ТРЕХ- И ЧЕТЫРЕХИНДЕКСНЫЕ ЗАДАЧИ С ВЛОЖЕННОЙ СТРУКТУРОЙ © 2012 г. Л.Г. Афраймович, А.С. Катеров
Нижегородский госуниверситет им. Н.И. Лобачевского levafraimovich@gmail.com
Поступила в редакцию 20.01.2012
Рассматриваются многоиндексные задачи линейного программирования транспортного типа. Один из подходов к решению таких задач - сведение их к задаче поиска потока в сети. Описана концепция такого сведения, в основу которой положено соответствие между переменными исходной задачи и дугами вспомогательной сети. Ранее было показано, что условие 2-вложенности многоиндексных задач является достаточным условием сводимости. В работе доказано, что (в рамках рассматриваемой концепции сводимости) условие 2-вложенности является необходимым и достаточным условием сводимости в трех- и четырехиндексном случае.
Ключевые слова: линейное программирование, многоиндексные задачи, задачи транспортного типа, сводимость, потоковые алгоритмы.
Введение
Существует широкий класс прикладных задач распределения ресурсов, формализуемых в виде многоиндексных задач линейного программирования транспортного типа. Примерами таких задач являются (см. [1-3]) задача объемно-календарного планирования для подразделений предприятия, задача формирования портфеля заказов, задача планирования и управления процессом переработки газового конденсата и др. Таким образом, разработка эффективных методов решения специального подкласса задач линейного программирования - многоиндексных задач линейного программирования транспортного типа делает возможным сокращение вычислительных затрат при решении ряда важных большеразмерных прикладных задач.
Для решения рассматриваемых многоиндексных транспортных задач могут быть применены общие методы решения задач линейного программирования - симплекс-метод, алгоритм Кармаркара [4]. Также существует ряд работ, посвященных непосредственно методам решения многоиндексных задач линейного программирования. Наиболее изученным является класс двухиндексных задач [5]. В общей постановке многоиндексные задачи рассматриваются в [6]. Условия, при которых удается понизить размерность и (или) сократить количество индексов многоиндексных транспортных задач, обсуждаются в [7]. Геометрические свойства множества допустимых решений многоиндексных систем линейных неравенств обсуждаются в [8].
Интерес представляет решение целочисленных многоиндексных задач линейного программирования транспортного типа. Как известно, двухиндексная задача целочисленного линейного программирования разрешима за полиномиальное время [5]. Однако данная задача является КР-трудной уже в трехиндексном случае [9]. Целочисленные многоиндексные транспортные задачи исследуются также в работах [10-12].
Одним из перспективных направлений при разработке эффективных алгоритмов исследования многоиндексных задач линейного программирования является нахождение подклассов задач, для решения которых применимы потоковые методы. Рассматриваемая в работе концепция сводимости класса задач линейного программирования к классу задач поиска потока минимальной стоимости была сформулирована в [13]. Важной особенностью данной концепции является существование соответствия между переменными исходной задачи и дугами вспомогательной сети. Предлагаемая схема сведения гарантирует, что произвольный оптимальный поток вспомогательной сети будет определять такое оптимальное решение исходной задачи, в котором переменным присваивается значение потока вдоль соответствующих дуг сети. В [1] было показано, что 2-вложенность многоиндексной задачи является достаточным условием сводимости. В случае выполнения условий 2-вложенности решение многоиндексной задачи линейного программирования транспортного типа сводится к поиску потока минимальной стоимости в сети с 0(т+п) вершинами и 0(т+п) дугами, где п - количество перемен-
ных, т - количество неравенств системы ограничений исходной задачи. В [13] показано, что 2-вложенность является необходимым и достаточным условием сводимости (в рамках рассматриваемой схемы сведения) в трехиндекс-ном случае, иначе Р=КР.
Данная работа является продолжением исследований, начатых в [1, 13]. В работе показано, что условие 2-вложенности является необходимым и достаточным условием сводимости (в рамках рассматриваемой схемы сведения) многоиндексных задач к классу задач поиска потока минимальной стоимости в трех- и четы-рехиндексном случае.
Многоиндексные транспортные задачи
При постановке многоиндексных задач линейного программирования транспортного типа воспользуемся формализацией, предложенной в [6]. Пусть 5 е N и N(5) = {1 ,2,..., s}. Каждому числу I поставим в соответствие параметр у называемый индексом, который принимает значения из множества Jl = {1, 2,., п/}, где п1 > 2, I е N(5). Пусть f = {k1,k2,...,k(}сN(5); k¡ <kм,
Г ^
1 ГІ1
22ргр7 2 2 ... 2^/¿7 , ¿/ Є Е/ •
Ff е Е- , / е М ; {с^( } - заданная 5-индексная
матрица коэффициентов целевой функции; {х^ } - 5-индексная матрица неизвестных.
Тогда многоиндексная задача линейного программирования транспортного типа формализуется следующим образом:
< Ь., е Е-, / еМ, (1)
г/
> 0, є Е
N ( і)
N (б ) N (ї )
(2)
(3)
Далее класс всех многоиндексных задач вида (1)-(3) при фиксированном множестве М будем обозначать Щ(М). Матрицу системы ограничений (1) обозначим через Маг(М, п1,.,пі), так как она однозначно определяется множеством М, а также количеством принимаемых значений для каждого индекса: п1,п2,.,Пц.
Для описания подматриц Ма^(М, п1,.,пі) будем использовать следующую схему. Ограничение вида (1), соответствующее фиксированному множеству / є М и фиксированному
і = 1, t -1. Набор значений индексов Ff = = (Л,, Л2,. .,Л,)/ будем называть t-индексом. Множество всех ^индексов, соответствующих / определим как Е/ = {¿/ | F є Jк^ х Jk^ х... х Jkt} .
Там где это не вызывает неоднозначности, будем опускать нижний индекс / в записи t-индексов. Если , є Ег, , є Е^„
набору F/ є Е/,
обозначим ё(/, ). Строку
где
/" , / "С N (5) и /" П / " = 0 , то через ,/;,/; обозначим такой набор, что е Егиг . Да-
лее определим / = N(5) \ / , тогда согласно введенному обозначению FN(5) = .
Каждому набору Ff поставим в соответствие действительное число zFf, е Е/ . Данное
отображение множества ¿-индексов Е/ в множество действительных чисел назовем (как и в [6]) ¿-индексной матрицей и обозначим через } .
Рассмотрим 5-индексную матрицу {zN(5)} и введем следующее обозначение
,/еЕ/ Л, Л2 еЛ2
Введенное обозначение подсумм 5-индекс-ной матрицы будем использовать при формализации многоиндексных транспортных задач.
Пусть М - заданное множество, М с 2N(5);
{а„ }, {Ь„ } - заданные | /1 -индексные матри-
р/ р/
цы свободных коэффициентов, 0 < а„ < Ь,_,
/ /
матрицы Ма?г(М, п1,...,п5), определяемую двусторонним неравенством d(/,), обозначим
row(/,), е Е-, / е М. Столбец матрицы Ма?г(М), соответствующий переменной х, ,
обозначим сol(FN(5)), FN(5) е EN(5). Пусть заданы последовательности /(1),/(2),...,/№) еМ, е Е-^,е Е-^,...,F-(kL) е Е
F(1)
ГN ( і) 5
F(2) F (к2) є Е
1 N(і)т--і 1 N(і) ^ J^N(і):
/ ^ /(kl)
тогда подматрицу, образованную элементами матрицы Ма?г(М), находящимися на пересечении строк гон^/ ), ,
/(і
1 = 1,k1, и столбцов со/^щ^)) , ] = 1,k2, будем обозначать Мо1г((/(,), F(¡^s), 1 = 1, FN(S), У = 1, ^).
Концепция сводимости
Формализация общей концепции сводимости дана в работе [13]. Нас будет интересовать схема сведения многоиндексных задач линейного программирования транспортного типа к классу задач поиска потока минимальной стоимости, при которой решение многоиндексной задачи может быть определено как подмножество компонент решения соответствующей потоковой задачи. Ниже приводится формализация соответствующей схемы сведения.
єЕ
X
F
N ( і )
,еЕ
N (і N ( і )
Пусть А є Япхт, Ъ,Ъ-,Ъ+є Я", с є Ят - заданные параметры, х є Ят - вектор неизвестных. Через м(А, Ъ, с) будем обозначать задачу линейного программирования тіп{(с, х)| Ах < Ъ, х > 0} ; через м(А, Ъ~, Ъ+, с) - задачу линейного программирования тіп{(с,х) | Ъ~ < Ах< Ъ+,х > 0}. Для удобства через гом(А) и со1(А) будем обозначать количество строк и столбцов матрицы А, соответственно. Отметим, что задача м(А, Ъ~, Ъ , с) может быть описана с использованием обозначения вида м(А, Ъ, с). Тем не менее будем использовать обозначение м(А, Ъ , Ъ , с) в случае, когда хотим подчеркнуть, что система ограничений задачи представляет собой систему двусторонних неравенств. Пусть Щ - произвольный класс задач линейного программирования.
Определение 1. Будем говорить, что класс Щ является t1 - і1| t2 - і2| t3 - і3 сводимым к классу Щ”, если для любой задачи м = м(А',Ъ',с) є Щ' можно за время О(^) построить матрицу А”, за время О(}2) построить векторы Ъ”,с”, что м" = м(А",Ъ",с")єЩ" и при этом:
- задача м' совместна (ограничена) тогда и только тогда, когда совместна (ограничена) задача м”;
- если известно оптимальное (допустимое) решение х” задачи м”, то оптимальное (допустимое) решение х' задачи м' может быть построено за время О(^).
Здесь (- і1), (- і2), (- і3) - опциональные строковые обозначения вычислительных процедур, связанных с построением матрицы системы ограничений, свободных коэффициентов и коэффициентов целевой функции и с построением решения задачи соответственно.
Задачу м” (см. определение 1) будем называть задачей, соответствующей задаче м'. Иногда для удобства функции оценки вычислительной сложности t1, ^, ^ будем заменять на обозначения L или Р, подразумевая линейные или полиномиальные от размера индивидуальной задачи м' функции, соответственно.
В данной работе исследуется возможность сведения класса многоиндексных задач линейного программирования транспортного типа к классу задач поиска потока минимальной стоимости, который определяется следующим образом. Рассмотрим ориентированный граф G =
=(УО, АО), АО с УО2, здесь Уо и А0 - множество вершин и дуг графа О соответственно. Пусть /у, щ - пропускные способности дуги (і, л); еу -стоимость дуги (і,у); ху - неизвестная величина потока вдоль дуги (і, у), (і,у) є АО . Тогда через v(О■; /Лу, иу, еу ,(і, Л) є АО) обозначим следующую задачу поиска потока минимальной стоимости:
Z xj - Z xJ‘ = °’ i gVg ,
(i,j)eAa (j,i')EAa
lj ^ xij ^ uij , (Л j) G AG ,
Xj ^ °, (j, j) G AG ,
Z evxv ^ min
O', j)<A
Обозначим через Graph множество всех ориентированных графов. Класс задач поиска потока минимальной стоимости определим как
WGraph = lj , Uij , ej , (І, j) G Ag )| lj ,
uy g Z +,ey g Z, (i, j) g Ag ,G g Graph}.
Определение 2. Пусть W - класс задач линейного программирования с двусторонней системой линейных неравенств. Будем говорить, что класс Wявляется ti|t2 - equal\t3 - edge сводимым к классу WGraph, если класс W является ti|t2|t3 сводимым к классу WGraph; и если
v = v(G; lj, uj, ej, (U j) G Ag ) G WGraph является задачей, соответствующей задаче w = w(A,b~,b+, c) g W, то выполняются следующие условия. Найдутся такие инъективные функции а:{1,2,...,row(A)}^Ag, в:{1,2,...,col(A)}^Ag, что:
- Li) = K,ua(i) = b+, i = 1,r0w(A) ; l(u,v) = °, uuv =
= b*, (u,v) g Ag \{a(i)| i = 1,row(A)}, где b - некоторая достаточно большая величина, для
row( A)
определенности b* = Zb+ ;
k=1
- eK0 = Ci, i = 1, col(A); eUv = 0, (u, v) g Ag \ {P(i) | i = 1, col(A)};
- если xtj, (i, j) g Ag, является оптимальным (допустимым) решением задачи v, то y = (xpm,xp(2),...,xp(col(A))) будет являться оптимальным (допустимым) решением задачи w.
Таким образом, согласно определению 2, в случае ti|t2 - equal\t3 - edge сводимости класса W к классу WGraph гарантируется, что если w gW ,
v = v(G; lij , uj , Cj ,(І, j) G AG ) G WGraph и v является
задачей, соответствующей задаче w, то при построении задачи поиска потока минимальной стоимости v пропускные способности и стоимости дуг в задаче определяются через коэффициенты задачи w, а решение задачи w находится через подмножество компонент решения задачи v. Тогда можно предложить алгоритм решения задачи w, основанный на решении соответствующей задачи v и имеющий вычислительную сложность O(tx +12 +13 + ц(| Vg |,| Ag |)), где ц(^, m) - вычислительная сложность алгоритма ре-
шения задачи поиска потока минимальной стоимости в сети с n вершинами и m дугами. Далее в работе рассматриваются условия, при которых класс W(M) является t^t2 - equal\t3 -edge сводимым к классу WGraph.
Исследование сводимости
Так как система ограничений (1) задач класса W(M) определяется заданием множества M,M с 2N(s), будем решать задачу поиска таких множеств M, при которых класс задач W(M) является t1\t2 - equal\t3 - edge сводимым к классу
WGraph.
Определение 3. Множество M называется k-вложенным, если существует разбиение множества M на k подмножеств
Mt = {/i(°, /2),...,4° }, i = U, что j е f%, j = 1, mt -1, i = 1, k .
Теорема 1 [1]. Пусть M с 2N(s). Для того чтобы класс W(M) являлся L\L - equal\L - edge сводимым к классу WGraph, достаточно, чтобы множество M было 2-вложенным.
Теорема 2 [14]. Пусть M с 2N(s). Для того чтобы класс W(M) являлся t1|t2 - equal\t3 - edge сводимым к классу WGraph, необходимо, чтобы Matr(M, nb...,ns) для любых n1,.,ns была абсолютно унимодулярна.
Утверждение 1 [15]. Пусть s = 3, тогда множество M является 2-вложенным тогда и только тогда, когда |{/1| /1= k,/ є M}|< 2,
k = 0, s.
Теорема 3 [15]. Множество M является 2-вложенным тогда и только тогда, когда для любых {/1, /2, /3} є M {/1, /2, /3} - 2-вложенное.
Проверка условий 2-вложенности в общем случае связана с перебором всех возможных разбиений множества M на два подмножества и имеет сложность 0(2^). Однако возможно исследовать 2-вложенность множества M проверив выполнение условий теоремы 3, т.е. перебрав все тройки элементов из M. Соответствующая проверка требует 0(|M|3) вычислительных операций.
Определение 4. Кодом множества M = {/1, /:,.../} M с 2N(s), будем называть вектор K =
= (lbl2,...,lk), li = /\, i = 1, k .
Лемма 1. Пусть s=3. Если множество M = = {/1,/:,•••/:} M с 2N(s), не является 2-вложен-ным, то найдутся n1, n2, n3, такие, что Matr(M, n1,n2,n3) не абсолютно унимодулярна.
Доказательство:
По утверждению 1 множество М может быть не 2-вложенным только для множеств, соответствующих следующим кодам:
1. К = (1,1,1), тогда М = {{1},{2},{3}}.
2. К = (2,2,2), тогда М = {{1,2},{2,3},{1,3}}. Было установлено существование п\,п2,п3,
таких, что Ма?г(М, пьп2,п3) не абсолютно уни-модулярна для следующих множеств М: М = {{1},{2},{3}}, М = {{1,2},{2,3},{ 1,3}} (см. приложение). Лемма доказана.
Лемма 2. Пусть 5 = 4. Если множество М = = {/1,/:,• ••/:}, М с 2N(5), не является 2-вложен-ным, то найдутся п\,...,п4, такие, что Ма?г(М, пь... ,п4) не абсолютно унимодулярна. Доказательство:
Так как М и все его элементы множества, выполняются следующие условия:
/ ф/у для любых 1,у е {1,2,3}, 1 ф у, (6)
/у = -[А^0, k2j),..., ¿/[К у =1,3, тогда
^у) ф kl(у) для любых
у е {1,2,3}, 1,I е{1,...,|/у |}, 1 ф I. (7)
Рассмотрим все возможные коды множеств, состоящих из трёх элементов при 5 = 4, для каждого кода построим все возможные не 2-вложенные множества М.
• К = (1,1,1): _
/1 = {11},/2 = {12},/3 = {13}, 1у е {1,...,4},у = 1,3. Положим М = {{1},{2},{3}}.
• К = (2,1,1): _
/1 = {11, 12},/2 = {13},/3 = {14}, 1у е {1,...,4},у = 1,4.
Пусть /1 = {1, 4}, тогда /3,/4 е{1,2,3,4}\ / = = {2,3}, иначе М - 2-вложенное. Получаем М = = {{1,4},{2},{3}}.
• К = (2,2,1):
/1 = {11, 12}, /2 = {13, 14}, /3 = {15}, е{1,...,4},
у = 1,5.
Возможны 2 случая:
а) /1 п /2 =0, тогда / и /2 = {1,2,3,4} и для любого /5 е {1,...,4} М- 2-вложенное.
б) | /1 п /2 |= 1, допустим/1 = {1,4},/2 = {3,4}, тогда 15 е {1,2,3,4} \ (/ и /2) = {2}, иначе М - 2-вложенное. Получаем М = {{1,4},{3,4},{2}}.
• К = (2,2,2):
/1 = {11, 12},/2 = {13, 14},/3 = {15, 1б}, 1у е{1,...,4},
у=1,6.
Возможны 2 случая:
а) Найдутся I, т е {1,2,3} , что /1 п /т =0, предположим/1 = {1,4},/2 = {2,3}, тогда /5 е /1, 16 е /2, иначе нарушится условие (6). Пусть М = = {{1,4},{2,3},{ 1,2}}.
б) Для любых I, т е {1,2,3}, I фт, / п/т ф0 . Но тогда | /1 п /т |= 1, I, т е {1,2,3}, I ф т . Пусть / = {1,4}, /2 = {2,4}. Допустимы 2 под-случая:
61) {1,2,3,4} \ (/ и /2) = {3} с /,, то есть /5 = 3, в этом случае / п/2 с /3, тогда М = {{1,4},
{2,4},{3,4}}. 1 2 3
62) {1,2,3,4}\С/1и/2)={3}£^ тогда / п./3 = = {Л }, /2 п /3 = {1''2 }, 1''1 ф 1''2, иначе нарушается условие (6) либо {3} с /3. Единственное возможное множество М = {{1,4},{2,4},{1,2}}. Данное множество идентично с точностью до перестановки множеству М = {{1,2},{2,3},{1,3}}, поэтому не нарушим общности рассуждений, если положим М = М.
• К = (3,1,1):
У! = {/\, ^ /2 = {1Л}, /3 = ОбК
^ е {1,...,4}, у = 15.
В таком случае/2 = {1,. ,4}\/1 = {1'}, иначеМ
- 2-вложенное. Тогда / и /2 = {1,2,3,4} и для любого 15 е{1,...,4} либо М - 2-вложенное, либо /2 = /3. Таким образом, для данного кода не существует не 2-вложенного множества.
• К = (3,2,1):
/1 = {1Ь ^ /2 = 04 1"5}, /3 = ОбК
^ е{1,...,4}, у = 1,6.
Аналогично случаю К = (3,1,1).
• К = (3,2,2):
/1 = {1Ь ^ /3}, /2 = 04 /5}, /3 = 06 1"7},
е {1,...,4}, у = 17.
Пусть {1,2,3,4}\/1 = {х^}, тогда 1 е /2,/3, иначе М - 2-вложенное. Пусть 14 = 1'1 и 16 = 1'1. При этом 15,17 е /, иначе нарушается условие (7), и 15 ф 17 , иначе нарушается условие (6).
Например, положим / = {1,2,4}, при этом 14 = = 16 = 3, а 15,17 е /. Возможным множеством будет М = { {1,2,4},{2,3},{ 1,3}}.
• К = (3,3,1):
/1 = {1Ь ^ /3}, У2 = 04 ^ /6}, /3 = ^7^
^ е {1,...,4}, у = 1,7.
Аналогично случаю К = (3,1,1).
• К = (3,3,2):
/1 = {1Ь ^ 1"3}, /2 = ^ ^ /3 = ОЧ 1"8},
е{1,...,4},у = 1,8.
Тогда 11 ф 12, где {1^} = {1,2,3,4}\ /, {1"2 } = = {1,2,3,4}\ /2, при этом для любого 1 е{1,...,4},
1 ф /,\ , 1 ф 12 выполняется 1 е / п /2. Тогда
/3 = {1 1 ,1 2 } , иначе /3 с / или /3 с /2 .
Например, пусть /1 = {1,2,4}, /2 = {2,3,4}, тогда /3 = {3,4}. M = { {1,2,4},{2,3,4},{ 1,3}}.
• K = (3,3,3):
/1 = {ib i2, 13}, /2 = {i4, i5, i6}, /3 = {i7, i8, i9}i
ij є {1,...,4}, j = 1,9.
Тогда i1 ф/2, где {i'1} = {1,2,3,4}\ /1, {i’2 } = = {1,2,3,4}\ /2, при этом для любого i є{1,...,4},
1 ф i,1, i ф i'2 выполняется i є/1 п/,. Тогда i7 = i'1,
і" 8 = 2 , иначе /3 = /1 или /3 = /2, ^9 є /1 ^ ./2 .
Например, пусть /1 = {1,2,4}, /2 = {2,3,4}, тогда /3 = {i,1,3}. Пусть i = 4, M = {{2,3,4}, {1,2,4},{ 1,3,4}}.
Таким образом, возможны только следующие не 2-вложенные множества M = {/1, /2, /3}, M с 2N(s) при s = 4: {{1},{2},{3}}; {{1,4},{2}, {3}}; {{1,4},{3,4},{2}}; {{1,4},{2,3},{1,2}};
{{1,4},{2,4},{3,4}}; {{1,2},{2,3},{1,3}}; {{1,2,4}, {2,3},{ 1,3}}; {{1,2,4},{2,3,4},{1,3}}; {{2,3,4}, {1,2,4}, {1,3,4}}. Все остальные совпадают с приведенными с точностью до перенумерации элементов N(s). Для каждого из приведённых множеств были найдены n1,.,n4, при которых Matr(M, n1,.,n4) не абсолютно унимодулярна (см. приложение). Для случаев {{1},{2},{3}} и {{1,2},{2,3},{ 1,3}} были найдены n1',n2',n3', при которых Matr (M, n1', n2', n3') не абсолютно унимодулярна, а следовательно, не абсолютно унимодулярна будет и Matr(M, n1', n2', n3', n4) при любых n4, так как Matr (M, n1', n2', n3') является её подматрицей. Лемма доказана.
Теорема 4. Пусть M с 2N(s) ,s < 4. Для того чтобы класс W(M) являлся t1|t2 - equal\t3 - edge сводимым к классу WGraph, где t1(n), t2(n), t3(n) > n, n є N, необходимо и достаточно, чтобы множество M было 2-вложенным. Доказательство:
Необходимость. При s = {1,2} множество
n.N(s) п.
2 является 2-вложенным, следовательно, любое его подмножество является 2-вложенным. Тогда, согласно теореме 1, класс W(M) является L \L - equal\L - edge сводимым к классу WGraph.
Далее пусть s є {3,4}. Если множество M не 2-вложенное, то по теореме 3 существует не 2-вложенное множество M' = {/1,/2,/3}, M'сM . Тогда найдутся n1,.,ns (согласно лемме 1 при s = = 3 и лемме 2 при s = 4), такие, что Matr(M',n1,...,ns) не является абсолютно уни-модулярной. Так как Matr(M',n1,...,ns) - подматрица Matr(M, n1,.,ns), то Matr(M, n1,...,ns) также не является абсолютно унимодулярной. Отсюда по теореме 2 класс задач W(M) не является t1|t2 - equal\t3 - edge сводимым к классу
WGraph. Если множество M 2-вложенное, то, согласно теореме 1, класс W(M) является L \L -equal\L - edge сводимым к классу WGraph.
Достаточность. Автоматически следует из теоремы 1. Теорема доказана.
Таким образом, доказательство теоремы 4 удалось описать путём перебора всех возможных множеств вида {/1, /2, /3}. Однако такой способ доказательства будет неэффективен при s > 5 в связи с ростом числа множеств вида {/1, /2, /3}.
Заключение
В представленной работе исследован класс трех- и четырехиндексных задач линейного программирования транспортного типа. Рассмотрены вопросы сводимости многоиндексных задач к классу задач поиска потока в сети. В рамках исследуемой схемы сведения показано, что для класса многоиндексных задач с числом индексов, не превосходящим 4, условие 2-вложенности множества, определяющего систему ограничений, является одновременно необходимым и достаточным условием сводимости к классу задач поиска потока минимальной стоимости.
Данные результаты обобщают исследования, проведённые в [1, 13]. В дальнейшем планируется исследовать сводимость многоиндексных транспортных задач линейного программирования к классу задач поиска потока в древовидных сетях.
Приложение
В рамках проведённых исследований была разработана параллельная программная система, которая по заданным множеству M, а также n1,.,ns строит Matr(M, n1,...,ns). Далее проверяется абсолютная унимодулярность построенной матрицы, таким образом определяется выполнение необходимого условия сводимости теоремы 2. Данная программная система запускалась на суперЭВМ МП-20 (РФЯЦ-ВНИИЭФ, г. Саров).
M = { {1,2},{2,3}, {1,3} }
H = Matr( ({1,2},(1){3}),({1,3},(1){2}),({2,3},(1){1});
(1,1,2)N(3), (1,2,1)N(3) ,(2,1,1)N(3) )
' 0 1 1 ^
H1 =1 0 1
V1 1
det( H1) = 2;
M = {{1},{2},{3}}
H 2 = Matr(({ 1}, (1,3){2,3} ),({1}, (2,2){2,3} ),({1}, (2,3){2,3) ),({3},(1,1){1,2} ),({3},(2,2){1Д})),
({2},(1,2){13}), ({2},(3,3){13});
(1,1,2) N (3),(1,1,3) N (3), (1,2,2) N (3) , (2,2,2) N (3), (2,2,3) N (3), (3,1,3) N (3), (3,2,3) N (3))
С 0 1 0 0 0 1 0'
0 0 1 1 0 0 0
0 0 0 0 1 0 1
1 1 0 0 0 0 0
0 0 0 1 1 0 0
1 0 1 0 0 0 0
V 0 0 0 0 0 1 1,
H 2 =
det( H 2) = 2;
M = {{1,2},{2,3},{1,4}}
H 3 = Matr(({1,2},(1,1){34}), ({2,3}, (1,1){1,4}), ({2,3},(1,2){14}),({1,4},(1,1 ){2,3} ),({1,4},(1,2){2,3});
(1,1,1,2) N (4),(1,1,2,1) N (4) , (1,1,2,2) N (4),
(1,2,1,1) N (4),(2,1,1,1) N (4))
С 00011 ^
0 10 10 10 10 0
1 0 0 0 1
H 3 =
V 0 110 0 , det H 3 =-2 M = { {1,4} ,{2}, {3} }
H = Matr(({1,4}, (1,3){2,3}), ({1,4},(2,2){23}), ({1,4} ,(2,3 ){23}),({3},(1,1,1);1,2,4}), ({3} ,(2,2,1);1,2,4})),
({2},(1,2,1){134} ),({2},(3,3,1){134});
(1.1.2.1)n (4),(1,1,3,1)n (4), (1,2,2,1)n w,(2,2,2,1)n (4),
(2.2.3.1)n (4),(3,1,3,1)n (4), (3,2,3,1)n (4))
H = H 2, det( H) = 2
M = {{1,4},{3,4},{2}}
H = Matr(({1,4}, (1,3){2,3}), ({1,4},(2,2){23}),
({1,4},(2,3){23} ),({3,4},(1,1){12> ),({3,4},(2,2){12>)),
({2},(1,2,1){134} ),({2},(3,3,1){134});
(1.1.2.1)n (4),(1,1,3,1)n (4), (1,2,2,1)n (4),(2,2,2,1)n (4),
(2.2.3.1)n (4),(3,1,3,1)n (4), (3,2,3,1)n (4))
H = H 2, det( H) = 2
M = {{1,4},{2,4},{3,4}}
H = Matr(({1,4}, (1,3){2,3}), ({1,4},(2,2){23}),
({1,4},(2,3){23} ),({3,4},(1,1){12> ),({3,4},(2,2){12>)),
({2,4},(1,2){13} ),({2},(3,3){13});
(1.1.2.1)n (4),(1,1,3,1)n (4), (1,2,2,1)n (4),(2,2,2,1)n (4),
(2.2.3.1)n (4),(3,1,3,1)n (4), (3,2,3,1)n (4))
Н = Н 2, det(H) = 2 М = {{1,2,4},{2,3},{1,3}}
Н = Ма?г(({1,2,4},(1){3}),({1,3},(1,1){24}), ({2,3},(1,1){м});
(1,1,2,1)N(4) , (1,2,1,1)N(4), (2,1,1,1)N(4) )
Н = Н1, det(Н) = 2 М = { {1,2,4},{2,3,4},{1,3} }
Н = Ма?г(({1,2,4},(1){3}),({1,3},(1,1){24}),({2,3,4}, (1){1} ); (1,1,2,1) N (4),(1,2,1,1) N (4),(2,1,1,1) N (4))
Н = Н1, det( Н) = 2 М = { {1,2,4},{2,3,4},{1,3,4} }
Н = Маг(({1,2,4},(1){3}), ({1,3,4},(1){2}), ({2,3,4}, (1){1} ); (1,1,2,1) N (4),(1,2,1,1) N (4),(2,1,1,1) N (4))
Н = Н1, det(Н) = 2
Список литературы
1. Афраймович Л.Г., Прилуцкий М.Х. Многоиндексные задачи распределения ресурсов в иерархических системах // Автоматика и телемеханика. 2006. №. 6. С. 194-205.
2. Афраймович Л.Г., Прилуцкий М.Х. Многоиндексные задачи оптимального планирования производства // Автоматика и телемеханика. 2010. № 10. С. 148-155.
3. Прилуцкий М.Х., Костюков В.Е. Оптимизационные задачи объёмно-календарного планирования для нефтеперерабатывающих предприятий // Системы управления и информационные технологии. 2007. №2.1 (28). С. 188-192.
4. Схрейвер А. Теория линейного и целочисленного программирования. В 2-х т. М.: Мир, 1991.
5. Гейл Д. Теория линейных экономических моделей. М.: Мир, 1969.
6. Раскин Л.Г., Кириченко И.О. Многоиндексные задачи линейного программирования. М.: Радио и связь, 1982.
7. Junginger W. On representatives of multi-index transportation problems // Europ J. of Operational Research. 1993. V. 66 (3). P. 353-371.
8. Емеличев В.А., Ковалев М.М., Кравцов М.К. Многогранники, графы, оптимизация. М.: Наука, 1981.
9. Гэри М., Джонсон Д. Вычислительные машины и труднорешаемые задачи. М.: Мир, 1982.
10. Spieksma F.C.R. Multi-index assignment problems: complexity, approximation, applications / in: P.M. Pardalos, L.S. Pitsoulis (Eds.). Nonlinear Assignment Problems: Algorithms and Applications. Dordrecht: Kluwer Academic Publishers, 2000. P. 1-11.
11. Сергеев С.И. Новые нижние границы для три-планарной задачи назначения. Использование классической модели // Автоматика и телемеханика. 2008. № 12. С. 53-75.
12. Гимади Э.Х., Глазков Ю.В. Об асимптотически точном алгоритме решения одной модификации трехиндексной планарной задачи о назначениях // Дискретный анализ и исследование операций. Сер. 2. 2006. Т. 13. № 1. С. 10-26.
13. Афраймович Л.Г. Трехиндексные задачи линейного программирования с вложенной структурой // Автоматика и телемеханика. 2011. № 8. С. 109-120.
14. Афраймович Л.Г., Прилуцкий М.Х. Многопродуктовые потоки в древовидных сетях // Известия РАН. Теория и системы управления. 2008. № 2. С. 57-63.
15. Катеров А.С. Исследование сводимости трёх-индексных транспортных задач к поиску потока в сети с применением параллельных вычислений // Прикладная информатика и математическое моделирование: Межвузовский сборник научных трудов. М.: МГУП им. Ивана Фёдорова, 2011. С. 47-57.
THREE- AND FOUR-INDEX PROBLEMS WITH NESTED STRUCTURE L.G. Afraimovich, A. S. Katerov
Multi-index linear programming transport type problems are considered. One of the approaches to solving such problems is to reduce them to flow algorithms. We describe a concept of this reduction which is based on correspondence between the variables of the original problem and the arcs of an auxiliary network. It was previously shown that two-nested structure of multi-index problems is a sufficient condition for reducibility. In this paper, we prove that the two-nested structure (in the framework of the reducibility concept considered) is a necessary and sufficient condition for reducibility in the cases of three- and four-index problems.
Keywords: linear programming, multi-index problems, transport type problems, reducibility, flow algorithms.
