CC BY
33
Математика
УДК 519.6 DOI: 10.14529/mmph170101
О ДИНАМИЧЕСКОЙ ЗАДАЧЕ ПОСТРОЕНИЯ ОСТОВА ПОЛИЭДРАЛЬНОГО КОНУСА
С.И. Бастраков, Н.Ю. Золотых
Нижегородский государственный университет, Нижний Новгород, Российская Федерация E-mail: sergey.bastrakov@gmail.com
Рассматривается динамическая задача построения остова полиэдрального конуса. Задача состоит в последовательном выполнении операций добавления или удаления неравенств из фасетного описания полиэдрального конуса с соответствующим перестроением остова. Обсуждается возможность применения метода двойного описания для выполнения обеих операций, приводятся оценки трудоемкости. Для операции удаления неравенства анализируется зависимость размера выхода от размера входа.
Ключевые слова: система линейных неравенств; полиэдральный конус; построение двойственного описания, метод двойного описания.
Основные понятия и обозначения
Любой выпуклый полиэдр (многогранник) в Qd может быть представлен в двух видах: Вершинное описание - как сумма по-Минковскому выпуклой и конической оболочек конечных систем векторов:
P = conv{vb v2, ..., vk} + cone{ub u2, ..., up}, v, e Qd, i = 1, 2, ..., k, Uj e Qd, j = 1, 2, ...,p. Фасетное описание - как множество решений конечной системы линейных неравенств:
P = {х e Qd: Ax > b}, A e Qmxd, b e Qm. Ограниченный выпуклый полиэдр называется политопом.
Частным случаем выпуклого полиэдра является полиэдральный (многогранный) конус, который может быть представлен как коническая оболочка конечной системы векторов и как множество решений конечной однородной системы линейных неравенств. Неприводимая порождающая система векторов многогранного конуса называется его остовом. Обозначим U(C) остов конуса C.
Задача перехода между вершинным и фасетным описаниями выпуклого полиэдра (построение одного описания по заданному другому описанию) называется задачей построения двойственного описания полиэдра. Частным случаем данной задачи являются задачи построения выпуклой оболочки и нахождения общего решения системы линейных неравенств. Задача построения двойственного описания полиэдра является одной из центральных в теории систем линейных неравенств [1-4]. Известно множество алгоритмов решения данной задачи [5-10] и близких к ней задач, включая вопросы трудоемкости и организации параллельных вычислений [11-14]. Задача
построения двойственного описания выпуклого полиэдра в Qd сводится к аналогичной задаче
для полиэдрального конуса в Qd+1 (см., например, [4]). В дальнейшем изложении для краткости под полиэдрами и конусами понимаются, соответственно, выпуклые полиэдры и полиэдральные конусы.
Пусть A e Qm d. Обозначим A(I), где I с {1, 2, ..., m}, подматрицу матрицы A, составленную из строк с номерами из I.
Постановка задачи
В работе рассматривается следующая задача, тесно связанная с задачей построения двойственного описания. Заданы неприводимые вершинное и фасетное описания полиэдрального конуса Co в Qd:
С0 = (х е О*: Ах > 0} = сопе{иь и2, ..., ир}, А е Отх*, щ е О*, г = 1, 2, ...,р.
Требуется выполнить последовательность из N операций над конусом С0, обозначим Ск конус, полученный в результате выполнения операции с номером к (к = 1, 2, ..., Л). Рассматриваются операции двух типов:
1. Пересечение текущего конуса с заданным полупространством:
Ск = (х е Ск_ь акх > 0}. (1)
2. Удаление неравенства из фасетного описания текущего конуса [15]. Пусть Ск_1 = (х е О*: Ак_1х > 0}, где матрица Ак_! имеет тк строк. Задается номер удаляемого неравенства /к (1 < /к < тк). Результатом операции является конус
Ск = (х е О*: Ак х > 0}, Ак = А^Ш, 2, ..., тк} \ {/к}), (2)
где матрица Ак состоит из всех строк матрицы Ак_г, за исключением строки с номером /к.
В результате выполнения каждой из указанных операций требуется получить неприводимые вершинное и фасетное описания конуса Ск. Таким образом, операция (1) соответствует добавлению неравенства к фасетному описанию конуса, а операция (2) - удалению неравенства из фасетного описания, в обоих случаях производится перестройка остова. Отметим, что требование неприводимости описаний конуса может быть ослаблено для операции (1), однако является существенно важным для операции (2), как показано в [15]. Предполагается, что порядок выполнения операций фиксирован.
Необходимость выполнения подобных операций может возникать в некоторых задачах автоматического анализа и верификации программ для ЭВМ с использованием полиэдральных подходов [15-17]. При этом отметим, что в настоящей работе задача исследуется только с теоретической точки зрения.
Аналогичная задача для выпуклого полиэдра в О сводится к случаю конуса в О . Рассматриваемая задача построения остова конуса в дальнейшем называется динамической по аналогии с динамической задачей построения выпуклой оболочки точек, множество которых может меняться путем добавления или удаления точек в процессе работы. Дальнейшее изложение производится для случая конуса. Для простоты предполагается, что исходный конус и все промежуточные конусы являются острыми: гапк(Ак) = к = 1, 2, ..., Л, в этом случае векторы остова определяются единственным образом с точностью до умножения на положительные скаляры (см., например,
[4]).
В разделе 2 рассматривается операция добавления неравенства к фасетному описанию конуса с помощью метода двойного описания, приводятся оценки трудоемкости. Раздел 3 посвящен операции удаления неравенства из фасетного описания и анализу зависимости размера выхода от размера входа.
Операция добавления неравенства к фасетному описанию конуса
Рассмотрим выполнение операции (1), состоящей в добавлении нового неравенства к фасет-ному описанию многогранного конуса с перестройкой остова. Данная операция возникает в ин-крементных (по классификации [18]) алгоритмах построения двойственного описания полиэдра. В связи с тем, что для всех промежуточных конусов известны неприводимые вершинное и фа-сетное описания, целесообразным представляется использование метода двойного описания [5], иногда также называемого алгоритмом Моцкина-Бургера [19] или алгоритмом Черниковой [6]. Существует множество модификаций метода, например, [20-22]. В настоящем разделе приводится краткое описание метода в контексте выполнения операции добавления неравенства. В связи с этим не рассматривается начальный шаг метода, а также порядок обработки неравенств, который может существенно влиять на трудоемкость метода двойного описания в общем случае (см., например, [18, 20]).
На каждом шаге метода двойного описания выполняется пересечение текущего конуса с полупространством, определяемым очередным обрабатываемым неравенством. В начале шага известны вершинное и фасетное описания текущего конуса Ск_г, а также множество Е(Ск_г) пар смежных векторов его остова. В конце шага те же данные определяются для конуса Ск. Таким образом, шаг метода двойного описания соответствует операции (1).
Бастраков С.И., О динамической задаче построения
Золотых Н.Ю. остова полиэдрального конуса
Рассмотрим выполнение одного шага метода. Остов конуса Ск-1 разбивается на три подмножества в зависимости от того, выполняется ли для них добавляемое неравенство:
и+ = {и е С/(С-1): ак и > 0}, и = {и е ЩС^): ак и = 0}, и- = {и е ЩС^): ак и < 0}. Строятся комбинации пар смежных векторов из и+ х и-:
и± = {(аи)у + (-ау)и: и е и+, V е и-, {и, V} е Е(Ск-1)}. (3)
Искомый остов конуса С' определяется следующим образом:
и(Ск) = и+ и и, и и±. (4)
Информация о смежности векторов остова конуса Ск определяется следующим образом. Все пары векторов из множества и+, являющиеся смежными в конусе Ск-1, являются смежными и в конусе Ск. Таким образом, смежность необходимо определить для векторов из множества и0 и и±. Для каждого вектора и хранится множество 1пс(и) номеров неравенств, определяющих конус, которые удовлетворяются этим вектором как равенства.
Известны несколько критериев того, является ли пара векторов остова {и, V} смежной:
1. Алгебраический критерий [5]: гапк(Л(/пс(и) п 1пс(у))) = ё — 2.
2. Комбинаторный критерий [17]:
$ м> е и0 и и±: ^ Ф и, ^ Ф V, 1пс(и) п /nc(v) с 1пс(м>).
3. Графовый критерий (графовая модификация комбинаторного критерия) [22]. Определим неориентированный граф О = (V, Е) следующим образом: V = и0 и и±, Е = {{и, V}: и, V е V, |/пс(и) п /nc(v)| > ё — 2}. Тогда критерий смежности пары имеет вид:
$w е и0 и и±: w Ф и, w Ф V, {и, w} е Е, {V, w} е Е, 1пс(и) п /nc(v) с 1пс(м>).
4. Критерий на основе к-мерных деревьев [21].
Условие |/пс(и) п /nc(v)| > ё — 2 является необходимым (но в общем случае не достаточным) для смежности пары {и, V}.
Обозначим Тац(Ск-1) трудоемкость выполнения шага метода двойного описания для текущего конуса Ск-1. В качестве меры трудоемкости будем использовать количество арифметических операций над целыми числами, вопросы о битовой сложности и росте коэффициентов рассматриваются, например, в [3]. Для изложенного алгоритма справедливы следующие верхние оценки трудоемкости.
Утверждение 1 [23]. Трудоемкость шага метода двойного описания для текущего острого конуса Ск-1 составляет:
• Таёё(Ск-1) е 0(|и(Ск_1)|-(|и(Ск_1)| + ё) + ЩСк-^-т—ё2) при использовании алгебраического критерия смежности;
• Таёё(Ск-1) е 0(|и(Ск_1)|-(|и(Ск_1)| + ё) + |и(Ск_1)|-тк_1-(|и(Ск_1)|2 + ё)) при использовании комбинаторного и графового критериев смежности.
Как показывает утверждение 1, зависимость верхней оценки трудоемкости шага метода двойного описания от размера остова текущего конуса квадратичная при использовании алгебраического критерия и кубическая при использовании комбинаторного и графового критериев смежности. Отметим, что максимальный размер выхода шага метода двойного описания асимптотически ограничен сверху величиной | и(Ск-1)|2.
Для анализа трудоемкости в следующих разделах введем также обозначение для общей трудоемкости метода двойного описания. Обозначим Тоом(т, ё) трудоемкость метода двойного описания для построения остова конуса в Qё, заданного системой из т неравенств. Для описываемого варианта метода двойного описания справедливы следующие оценки:
Утверждение 2 [23]. Трудоемкость метода двойного описания для построения остова острого конуса при любом фиксированном ё составляет:
• Тоом(т, ё) е 0(тё+2) при использовании алгебраического критерия смежности;
• Тоом(т, ё) е 0(тзё2+1) при использовании комбинаторного и графового критериев смежности.
Операция удаления неравенства из фасетного описания конуса
Настоящий раздел посвящен операции (2) удаления неравенства из фасетного описания конуса. Отметим, что «наивным» способом решения задачи является прямое построение остова конуса Ск по системе неравенств, без использования информации о конусе Ск-1. В этом случае, при использовании метода двойного описания, трудоемкость составляет Тоом(тк-1—1, ё).
В [15] предложен алгоритм удаления неравенств из фасетного описания, названный авторами инкрементным. Обозначим — фасету, соответствующую удаляемому неравенству. Строится множество фасет смежных с С помощью метода двойного описания вычисляется остов конуса Саф определяемого неравенствами, соответствующими фасетами из В остове конуса Са^ находится подмножество векторов Щпе№, которые не принадлежат конусу Ск_1. Множество векторов Щ(Ск_1) и Щпе„ порождает конус Ск_1, однако не обязательно является неприводимым. Для завершения операции достаточно удалить избыточные элементы из множества Щ(Ск_1) и Щпе№. Таким образом, трудоемкость инкрементного алгоритма удаления неравенства составляет ?7гетоуе(Ск_1) е 0(7тм( l-Fad.il, *)). При этом в худшем случае полученная оценка совпадает с оценкой трудоемкости «наивного» алгоритма. Оценка в терминах количества фасет конуса Ск_1 может быть получена с помощью утверждения 2.
Отметим, что алгоритм [15] позволяет одновременное удаление нескольких неравенств, отличие от случая одного неравенства состоит лишь в том, что в добавляются фасеты, смежные хотя бы с одной из удаляемых.
В общем случае задача удаления неравенства из фасетного описания является существенно более трудной по сравнению с задачей добавления неравенства, как показывает следующий пример. Рассмотрим семейство политопов СС25(5), являющихся произведением 5 циклических политопов С25(5), конструкция и некоторые свойства данного семейства приведены в [18, 24]. Как показано в [24], при удалении любой вершины V политопа Р = СС25(я) количество фасет многогранника сопу(Р \ {V}), не являющихся фасетами Р, составляет 0(5 (5_1)5). Таким образом, удаление любой вершины такого политопа приводит к образованию очень большого количества новых фасет. При сведении к случаю удаления неравенства из фасетного описания конуса, из этого вытекает следующее утверждение:
Утверждение 3. Размер выхода задачи удаления одного неравенства из фасетного описания конуса не ограничен полиномом от длины входа.
Следующая теорема устанавливает один из подклассов, для которых размер выхода полиномиален от длины входа.
Теорема 1. Пусть острый конус Ск получен путем удаления неравенства ах > 0 из фасетного описания острого конуса Ск_1 и выполнено условие
"и е Щ(Ск) $ V е Щ(С*_0: (и, V} е Е(Ск), av > 0. (5)
Тогда справедливо неравенство
|Щ(Ск)| < |Щ(Ск_1)| + |Е(Ск_1)|. Доказательство. Рассмотрим обратную операцию: конус Ск_1 может быть получен из Ск путем добавления неравенства ах > 0 с помощью метода двойного описания, рассмотренного в разделе 2. При этом, ввиду соотношения (4), векторы из ЩСк), лежащие в полупространстве ах > 0, присутствуют и в и(Ск_\). Таким образом, справедливо неравенство
Щ(Ск)| < | ЩСк^) + | Щ(Ск) \ ЩСк^). (6)
Для оценки величины |Щ(Ск) \ Щ(Ск_1)| отметим, что векторы из этого множества создают комбинации со смежными им векторами из открытого полупространства ах > 0, согласно (3). Из условия (5) следует, что каждый вектор из ЩСк) \ Щ(Ск_1) создаст хотя бы одну комбинацию. Таким образом, каждый вектор из этого множества лежит на пересечении продолжения ребра конуса Ск_1, составленного парой векторов из Щ(Ск_\), один из которых являлся его «парой» при комбинировании, а другой - результатом комбинирования. Таким образом, количество ребер конуса Ск_1 дает верхнюю оценку величины |Щ(Ск) \ Щ(Ск_1)|. Подстановка этой оценки в соотношение (6) и дает доказываемое неравенство.
Так как величина |Е(Ск_г)| ограничена количеством (неупорядоченных) пар векторов остова и всегда меньше, чем | Щ(Ск_а)|2 / 2, справедливо следующее
Следствие 1. В условиях теоремы справедливо неравенство
|Щ(Ск)| < |Щ(Ск_1)| + |Щ(Ск_1)|2 / 2. Таким образом, в условиях теоремы 1 зависимость размера выхода от размера входа не более, чем квадратичная. Отметим, что доказательство теоремы дает полиномиальный от длины входа (в условиях теоремы) алгоритм построения ЩСк) с помощью продолжения ребер конуса Ск_1, составленных их пары векторов, один из которых лежит в удаляемой фасете, а другой - нет, до пересечения с фасетами. Смежность векторов может быть определена с помощью одного из
Бастраков С.И., О динамической задаче построения
Золотых Н.Ю. остова полиэдрального конуса
критериев, описанных в разделе 2. Принимая во внимание то, что оценки утверждения 1 также
являются полиномиальными от длины входа, справедлива следующая теорема:
Теорема 2. Пусть выполняется последовательность операций добавления и удаления неравенств вида (1), (2) и для каждой операции удаления выполнено условие (5). Тогда существует алгоритм выполнения последовательности операций с трудоемкостью, ограниченной полиномом от длины входа при любом фиксированном числе операций.
Заключение
В статье рассматривается задача построения остова полиэдрального конуса, множество неравенств которого динамически изменяется путем добавления новых и удаления существующих неравенств. Рассматривается применение метода двойного описания для операции добавления нового неравенства, соответствующей итерации метода, приводятся оценки трудоемкости. Для задачи удаления неравенства из фасетного описания конуса рассматривается «наивный» и ин-крементный алгоритмы, приводятся оценки трудоемкости. Приводится пример, показывающий, что размер выхода при удалении неравенства в общем случае не ограничен полиномом от длины входа. Выделяется подкласс задач, для него доказывается, что размер выхода ограничен квадратом от размера входа. Показано, что для данного класса существует алгоритм выполнения последовательного выполнения операций с трудоемкостью, ограниченной полиномом от длины входа при любом фиксированном числе операций.
Авторы благодарят В.Н. Шевченко и А.Н. Максименко за полезные обсуждения.
Литература
1. Черников, С.Н. Линейные неравенства / С.Н. Черников. - М.: Наука, 1968. - 488 с.
2. Емеличев, В.А. Многогранники, графы, оптимизация / В.А. Емеличев, М.М. Ковалев, М.К. Кравцов. - М.: Наука, 1981. - 344 с.
3. Схрейвер, А. Теория линейного и целочисленного программирования / А. Схрейвер. - М.: Мир, 1991. - Том 1. - 360 с.
4. Циглер, Г. Теория многогранников / Г. Циглер. - М.: МЦНМО, 2014. - 586 с.
5. Метод двойного описания / Т.С. Моцкин, Х. Райфа, Д. Л. Томпсон, Р.М. Тролл // Матричные игры: сб. науч. тр. - М.: Физматгиз, 1961. - С. 81-109.
6. Черникова, Н.Б. Алгоритм для нахождения общей формулы неотрицательных решений системы линейных неравенств / Н.Б. Черникова // Журн. вычисл. математики и мат. физики. - 1965. - Т. 5, № 2. - С. 334-337.
7. Avis, D. A Pivoting Algorithm for Convex Hull and Vertex Enumeration of Arrangements and Polyhedra / D. Avis, K. Fukuda // Discrete and Computational Geometry. - 1992. - Vol. 8. - Issue 3. -P.295-313.
8. Chazelle, B. An Optimal Convex Hull Algorithm in Any Fixed Dimension / B. Chazelle // Discrete and Computational Geom. - 1993. - Vol. 10. - Issue 4. - P. 377-409.
9. Barber, C.B. The Quickhull Algorithm for Convex Hulls / C.B. Barber, D.P. Dobkin, H. Huhdanpaa // ACM Transactions on Mathematical Software. - 1996. - Vol. 22, no. 4. - P. 469-483.
10. Bremner, D. Primal-Dual Methods for Vertex and Facet Enumeration / D. Bremner, K. Fukuda,
A. Marzetta // Discrete and Computational Geometry. - 1998. - Vol. 20. - Issue 3. - P. 333-357.
11. Шевченко, В.Н. Модификация алгоритма Фурье-Моцкина для построения триангуляции /
B.Н. Шевченко, Д.В. Груздев // Дискретный анализ и исследование операций. Сер. 2. - 2003. -Т. 10, № 1. - С. 53-64.
12. Панюков, А.В. Параллельные реализации симплекс-метода для безошибочного решения задач линейного программирования / А.В. Панюков, В.В. Горбик // Вестник ЮУрГУ. Серия «Математическое моделирование и программирование». - 2011. - Вып. 9. - № 25(242). - C. 107-118.
13. Панюков, А.В. Представление суммы Минковского для двух полиэдров системой линейных неравенств / А.В. Панюков // Вестник ЮУрГУ. Серия «Математическое моделирование и программирование». - 2012. - Вып. 14, № 40(299). - C. 108-119.
14. Панюков, А.В. Подход к решению систем линейных алгебраических уравнений с интервальной неопределенностью в исходных данных / А.В. Панюков, В.А. Голодов // Вестник ЮУрГУ. Серия «Математическое моделирование и программирование». - 2013. - T. 6, № 2 - C. 108119.
15. Amato, G. Efficient Constraint/Generator Removal from Double Description of Polyhedra / G. Amato, F. Scozzari, E. Zaffanella // Electronic Notes in Theoretical Computer Science. - 2014. -Vol. 307. - P. 3-15.
16. Cousot, P. Automatic discovery of linear restraints among variables of a program / P. Cousot, N. Halbwachs // Conference Record of the Fifth Annual ACM Symposium on Principles of Programming Languages. - 1978. - P. 84-96.
17. Bagnara, R. Applications of polyhedral computations to the analysis and verification of hardware and software systems / R. Bagnara, P.M. Hill, E. Zaffanella // Theoretical Computer Science. -2009. - Vol. 410. - P. 4672-4691.
18. Avis, D. How Good are Convex Hull Algorithms? / D. Avis, D. Bremner, R. Seidel // Computational Geometry. - 1997. - Vol. 7, Issues 5-6. - P. 265-301.
19. Burger, E. Uber homogene lineare Ungleihungssysteme / E. Burger // Zeitschrift Angewandte Math. Mehanik. - 1956. - Vol. 36. - Issue 3-4. - P. 135-139.
20. Fukuda, K. Double Desription Method Revisited / K. Fukuda, A. Prodon // Lecture Notes in Computer Science. - 1996. - Vol. 1120. - P. 91-111.
21. Terzer, M. Accelerating the Computation of Elementary Modes Using Pattern Trees / M. Terzer, J. Stelling // Lecture Notes in Computer Science. - 2006. - Vol. 4175. - P. 333-343.
22. Золотых, Н.Ю. Новая модификация метода двойного описания для построения остова многогранного конуса / Н.Ю. Золотых // Журн. вычисл. математики и мат. физики. - 2012. - Т. 52, № 1. - С. 153-163.
23. Бастраков, С.И. Алгоритмические вопросы построения двойственного описания выпуклого полиэдра: дис. ... канд. физ.-мат. наук / С.И. Бастраков. - Н. Новгород, 2016. - 100 с.
24. Bremner, D. Incremental Convex Hull Algorithms Are Not Output Sensitive / D. Bremner // Discrete and Computational Geometry. - 1999. - Vol. 21. - Issue 1. - P. 57-68.
Поступила в редакцию 11 октября 2016 г.
Bulletin of the South Ural State University Series "Mathematics. Mechanics. Physics" _2017, vol. 9, no. 1, pp. 5-12
MSC 52B55 DOI: 10.14529/mmph170101
ON THE DYNAMIC PROBLEM OF COMPUTING GENERATORS OF A POLYHEDRAL CONE
S.I. Bastrakov, N.Yu. Zolotykh
Lobachevsky State University of Nizhni Novgorod, Nizhniy Novgorod, Russian Federation
This paper considers a dynamic problem of computing generators of a polyhedral cone. The problem is to sequentially perform operations of adding and removing inequalities from a facet description of the polyhedral cone with a corresponding re-computation of generators. The application of a double description method for both operations is discussed and complexity estimation is given in the paper.
Adding a new inequality corresponds to a single step of the double description method. It can be performed with time complexity being quadratic or cubic of the input size for the current step, depending on the modification of the method and adjacency tests chosen. We give complexity bounds for adding a single inequality with widely used algebraic and combinatorial adjacency tests.
The problem of removing inequalities is intrinsically much harder, compared to adding inequalities. We briefly describe the naive and incremental algorithms and show an example with output size being superpolynomial of the input size in case of removing a single inequality. A subclass of problems with certain adjacency properties is investigated, for this subclass we prove that the output size is bounded by a quadratic function of the input size. Finally, we prove that for the distinguished subclass any finite sequence of adding and removing inequalities can be performed in polynomial time of the input size.
Keywords: system of linear inequalities; polyhedral cone; computing dual description; double description method.
Бастраков С.И., Золотых Н.Ю.
О динамической задаче построения остова полиэдрального конуса
References
1. Chernikov S.N. Lineynye neravenstva (Linear inequalities). Moscow, Nauka Publ., 1968, 488 p. (in Russ.).
2. Emelichev V.A., Kovalev M.M., Kravtsov M.K. Mnogogranniki, grafy, optimizatsiya (Polyhedrons, graphs, optimization). Moscow, Nauka Publ., 1981, 344 p. (in Russ.).
3. Schrijver A. Theory of Linear and Integer Programming. New York, Wiley, 1986. 471 p. DOI: 10.1002/net.3230200608
4. Ziegler G.M. Lectures on Polytopes. Graduate Texts in Mathematics, vol. 152. Berlin, New York: Springer-Verlag, 1995. 373 p. DOI: 10.1007/978-1-4613-8431-1
5. Motzkin T.S., Raiffa H., Thompson G.L., Trall R.M. The Double Description Method. Contributions to the Theory of Games. Princeton University Press, 1953, Vol. 2, pp. 51-74. DOI (Chapter): 10.1515/9781400881970-004
6. Chernikova N.V. Algorithm for finding a general formula for the non-negative solutions of a system of linear inequalities. USSR Computational Mathematics and Mathematical Physics, Vol. 5, Issue 2, 1965, pp. 228-233. D0I:10.1016/0041-5553(65)90045-5
7. Avis D., Fukuda K. A Pivoting Algorithm for Convex Hull and Vertex Enumeration of Arrangements and Polyhedra. Discrete and Computational Geometry, 1992, Vol. 8, Issue 3, pp. 295-313. DOI: 10.1007/BF02293050
8. Chazelle B. An Optimal Convex Hull Algorithm in Any Fixed Dimension. Discrete and Computational Geom., 1993, Vol. 10, Issue 4, pp. 377-409. DOI: 10.1007/BF02573985
9. Barber C.B., Dobkin D.P., Huhdanpaa H. The Quickhull Algorithm for Convex Hulls. ACM Transactions on Mathematical Software, 1996, Vol. 22, no. 4, pp. 469-483. DOI: 10.1145/235815.235821
10. Bremner D., Fukuda K., Marzetta A. Primal-Dual Methods for Vertex and Facet Enumeration. Discrete and Computational Geometry, 1998, Vol. 20, Issue 3, pp. 333-357.
11. Shevchenko V. N., Gruzdev D. V. Modifikatsiya algoritma Fur'e-Motskina dlya postroeniya triangulyatsii (A modification of the Fourier-Motskin algorithm for constructing a triangulation). Diskretnyy analiz i issledovanie operatsiy, Ser. 2, 2003, Vol. 10, no. 1, pp. 53-64. (in Russ.).
12. Panyukov A.V., Gorbik V.V. The parallel simplex-method achievements for errorless solving of linear programming problems. Bulletin of South Ural State University. Series: Mathematical Modelling, Programming & Computer Software, 2011, Issue 9, no. 25(242), pp. 107-118. (in Russ.).
13. Panyukov A.V. The Linear Inequalities Set Representation of Minkovski's Sum for Two Polyhedrons. Bulletin of South Ural State University. Series: Mathematical Modelling, Programming & Computer Software, 2012, Issue 14, no. 40(299), pp. 108-119. (in Russ.).
14. Panyukov A.V., Golodov V.A. Approach to Solve the Set of Linear Algrebraic Equations with Interval Uncertainty of Data Given. Bulletin of South Ural State University. Series: Mathematical Modelling, Programming & Computer Software, 2013, Vol. 6, no. 2, pp. 108-119. (in Russ.).
15. Amato G., Scozzari F., Zaffanella E. Efficient Constraint/Generator Removal from Double Description of Polyhedra. Electronic Notes in Theoretical Computer Science, 2014, Vol. 307, pp. 3-15. DOI: 10.1016/j.entcs.2014.08.002
16. Cousot P., Halbwachs N. Automatic discovery of linear restraints among variables of a program. Conference Record of the Fifth Annual ACM Symposium on Principles of Programming Languages, 1978, pp. 84-96. DOI: 10.1145/512760.512770
17. Bagnara R., Hill P.M., Zaffanella E. Applications of polyhedral computations to the analysis and verification of hardware and software systems. Theoretical Computer Science, 2009, Vol. 410, Issue 46, pp. 4672-4691. DOI: 10.1016/j .tcs.2009.07.033
18. Avis D., Bremner D., Seidel R. How Good are Convex Hull Algorithms? Computational Geometry, 1997, Vol. 7, Issues 5-6, pp. 265-301. DOI:10.1016/S0925-7721(96)00023-5
19. Burger E. Über homogene lineare Ungleichungssysteme. Zeitschrift Angewandte Math. Me-hanik, 1956, Vol. 36, Issue 3-4, pp. 135-139. DOI: 10.1002/zamm.19560360308
20. Fukuda K., Prodon A. Double Desription Method Revisited. Lecture Notes in Computer Science, 1996, Vol. 1120, pp. 91-111. DOI: 10.1007/3-540-61576-8_77
21. Terzer M., Stelling J. Accelerating the Computation of Elementary Modes Using Pattern Trees. Lecture Notes in Computer Science, 2006, Vol. 4175, pp. 333-343. DOI: 10.1007/11851561_31
22. Zolotykh N.Yu. New modification of the double description method for constructing the skeleton of a polyhedral cone. Journal of Computational mathematics and Mathematical Physics, 2012, Vol. 52, Issue 1, pp. 146-156. DOI: 10.1134/S0965542512010162
23. Bastrakov S.I. Algoritmicheskie voprosy postroeniya dvoistvennogo opisaniya vypuklogo polye-dra: dis. ... kand. fiz.-mat. nauk (Algorithmic issues of double description of a convex polyhedron. Cand. phys. and math. sci. diss.). Nizhniy Novgorod, 2016, 100 p. (in Russ.).
24. Bremner D. Incremental Convex Hull Algorithms are Not Output Sensitive. Discrete and Computational Geometry, 1999, Vol. 21, Issue 1, pp. 57-68. DOI: 10.1007/PL00009410
Received October 11, 2016
