2019, Т. 161, кн. 1 С. 127-144
УЧЕНЫЕ ЗАПИСКИ КАЗАНСКОГО УНИВЕРСИТЕТА. СЕРИЯ ФИЗИКО-МАТЕМАТИЧЕСКИЕ НАУКИ
ISSN 2541-7746 (Print) ISSN 2500-2198 (Online)
УДК 514.174.5+514.172.45+519.667+519.853 <М: 10.26907/2541-7746.2019.1.127-144
СИМПЛИЦИАЛЬНЫЙ АЛГОРИТМ ПОИСКА НЕПОДВИЖНЫХ ТОЧЕК С 2а ЦЕЛОЧИСЛЕННЫМИ
МЕТКАМИ
М.Н. Матвеев
Московский физико-технический институт, г. Долгопрудный, 141701, Россия
Аннотация
Симплициальные алгоритмы, применяемые для поиска неподвижных точек непрерывных функций, могут использовать или векторные, или целочисленные метки. Алгоритмы с целочисленными метками являются более простыми и в силу дискретности целочисленных меток устойчивыми к ошибкам округления. Однако вычислительная эффективность существующих алгоритмов с целочисленными метками ограничивается недостаточной гибкостью их конструкции, в частности, в пространстве КМ эти алгоритмы могут использовать только от <¿+1 до 2й меток. Не более чем сопоставимое с размерностью пространства количество меток делает алгоритмы с целочисленными метками недостаточно быстрыми, особенно в задачах большой размерности. Настоящая работа преодолевает эту трудность и представляет новый симплициальный алгоритм поиска неподвижных точек с 2 целочисленными метками.
Ключевые слова: триангуляции, многогранники, вееры, симплициальные алгоритмы, алгоритмы поиска неподвижных точек
Введение
Симплициальные алгоритмы названы так, потому что в процессе работы, в частности в процессе поиска неподвижных точек х* = /(х*) непрерывных функций /(х) в пространстве М^, они перебирают симплексы некоторой триангуляции. Для того чтобы перебор симплексов был организован должным образом, каждой вершине триангуляции присваивается некоторая метка. Метки могут быть или векторами (как правило, с действительными координатами), или (скалярными) целыми числами и называются соответственно векторными или целочисленными. Векторные метки чаще всего являются векторами вида /(х) — х, где х - вершина триангуляции. В общем же случае и те и другие метки присваиваются в соответствии с некоторым правилом, называемым правилом помечивания.
С самого своего возникновения симплициальные алгоритмы используют как векторные, так и целочисленные метки. Первый алгоритм поиска неподвижных точек, разработанный в 1967 г. Скарфом в серии работ [1-3], может применяться также и в значительно более широком контексте [4] (имеется русский перевод [5]), вместо симплексов перебирает так называемые примитивные множества и использует векторные метки. Однако в том же году Коэн [6] заметил, что уже доказательство леммы Шпернера [7] может быть интерпретировано как симплициальный алгоритм с целочисленными метками, а годом позднее Кун [8] описал и модификацию алгоритма Скарфа, в котором перебираются симплексы триангуляции и используются целочисленные метки.
В целом векторные метки применяются для решения более широкого класса задач, а целочисленные значительно упрощают вычисления. Использование и тех
и других меток прошло достаточно длительный путь развития и сформировало свою собственную философию. Хороший обзор этих вопросов, а также достаточно полное описание симплициальных алгоритмов можно найти в [9] (имеется русский перевод [10]) и [11]. Симплициальные алгоритмы как с векторными, так и с целочисленными метками различают главным образом по количеству лучей, вдоль которых эти алгоритмы могут покидать точку начала поиска. Для алгоритмов с целочисленными метками, за исключением 2 -лучевого алгоритма Ямамото [12], количество таких лучей совпадает с количеством меток.
Большее количество реализаций симплициальных алгоритмов поиска неподвижных точек использует векторные метки. Более того, описаны общие способы построения алгоритмов с векторными метками для произвольного количества лучей: сначала в контексте поиска решений нелинейных уравнений такой способ предложили Коджима и Ямамото [13], а затем в несколько более широком контексте поиска стационарных точек - Тальман и Ямамото [14]. Замечательными отдельными реализациями симплициальных алгоритмов с векторными метками являются 2й-лучевой, или октаэдральный, алгоритм Райта [15] и (3й — 1)-лучевой алгоритм Коджимы и Ямамото [13].
Симплициальные алгоритмы поиска неподвижных точек с целочисленными метками менее разнообразны. Реализации этих алгоритмов в пространстве мй по существу представлены единственной конструкцией ван дер Лаана и Тальмана [16], основанной на декартовом произведении симплексов и использующей от с! +1 до 2с! меток. Аналогичная конструкция с 2! метками для декартова произведения ! одномерных симплексов или куба описана Рэйзером [17]. Помимо этого существует только единственный в своем роде алгоритм с целочисленными метками, число лучей которого не совпадает с числом меток, - это 2-лучевой алгоритм Ямамото [12], который использует с! +1 целочисленную метку.
Дискретность целочисленных меток значительно упрощает симплициальные алгоритмы и, что еще более важно, делает эти алгоритмы устойчивыми к ошибкам округления. Поэтому с вычислительной точки зрения меньшая распространенность целочисленных меток выглядит достаточно ограничительной. Это послужило причиной для разработки общего способа [18] построения симплициальных алгоримов с произвольным количеством целочисленных меток. Используя общий способ [18], настоящая работа представляет новый симплициальный алгоритм поиска неподвижных точек с 2й целочисленными метками, являющийся аналогом 2й-лучевого или октаэдральногго алгоритма Райта [15] с векторными метками.
При построении нового алгоритма работа опирается на два основных принципа. Первый из них состоит в том, что чем больше меток или лучей, вдоль которых алгоритм может покидать точку начала поиска, тем выше его вычислительная эффективность. Теоретически этот принцип достаточно хорошо обоснован в [18], где приведены также подтверждающие его результаты численных экспериментов в пространстве м4. Косвенно подтверждающими этот принцип являются также результаты численных экспериментов в пространствах мй, ! ^ 50, более высокой размерности, приведенные в [13].
Большое количество меток имеет, однако, и обратную сторону: чем больше меток, тем, вообще говоря, больше времени нужно для их вычисления. Для того чтобы симплициальный алгоритм был достаточно быстрым, вычисления используемых им меток должны быть эффективными. Второй принцип настоящей работы состоит в том, что эффективность вычисления меток может быть обеспечена регулярностью расположения лучей, вдоль которых алгоритм может покидать точку начала поиска. На основании этого принципа в качестве исходной геометрической конструкции для построения нового алгоритма берется !-мерный куб, имеющий
2й регулярно расположенных вершин. Количество вершин куба и определяет количество меток алгоритма.
Содержание работы самодостаточно для ее понимания. В то же время она не имеет целью и не может заменить ни статью [18], где приведена, обоснована и проиллюстрирована общая схема построения алгоритмов с целочисленными метками, ни статью [16] и тем более книгу [11], которые содержат разные, но одинаково хорошие доказательства сходимости алгоритмов с целочисленными метками. Оба доказательства являются верными и для алгоритма, который описывается в настоящей работе, и здесь не дублируются.
Вместо этого работа акцентирует внимание на чертах, присущих исключительно новому симплициальному алгоритму. Сначала доказываются свойства триангуляции К\ пространства мй, необходимые для применения общей схемы [18]. Затем рассматривается, каким образом триангуляция К\ может быть использована для построения симплициального многогранника из куба. После этого описывается разделение приближений неподвижных точек, найденных новым алгоритмом, на слабые и сильные. В заключении приводится ряд замечаний, касающихся триангуляции К, которая существенно используется в работе.
1. Настройка триангуляции
В основе любого симлициального алгоритма лежит триангуляция. Множество, будь то целое пространство или только его часть, будучи разделенным на симплексы, предоставляет дополнительные возможности для отыскания неподвижных точек. Алгоритму удается организовать перебор симплексов триангуляции некоторым рациональным образом, таким что после достаточно малого по сравнению с возможным числа шагов оказывается найден симплекс, являющийся приближением неподвижной точки.
Понятие триангуляции определяется следующим образом. Пусть Б - произвольное множество в мй. Набор Т выпуклых многогранных множеств называется разбиением Б, если
(1) Б является объединением всех множеств из Т;
(2) пересечение любых двух множеств из Т или пусто, или является их общей собственной гранью;
(3)1 любая точка Б имеет окрестность, пересекающуюся только с конечным числом множеств из Т.
Если разбиение Т состоит исключительно из симплексов, то есть выпуклых оболочек аффинно-независимых точек, то такое разбиение называется симплици-альным разбиением или триангуляцией.
Существует много различных триангуляций, описание и сравнение некоторых из них можно найти, например, в [9] (или [10]). Для алгоритма настоящей работы является подходящей и будет использоваться триангуляция К1 (пространства мй), открытая Фройденталем [19]. Для построения К1 все пространство мй сначала разбивается на кубы, а затем каждый куб - на симплексы. Заметим, что для описания К1 удобно использовать единичные кубы, на практике же, как правило, применяется триангуляция К1 мелкости О - с уменьшенным в О раз относительно единичного размером кубов и, соответственно, симплексов.
В пространстве м3 триангуляция К1 иллюстрируется рис. 1, в общем случае ее описание выглядит следующим образом. Пусть щ, г = 1,...,!,- единичные векторы некоторого ортогонального базиса пространства мй и пусть zd - множество
1Свойство (3) не всегда используется при определении разбиений. Разбиения, удовлетворяющие (3), называются локально-конечными.
Рис. 1. Триангуляция К1 пространства к3: (а) единичный куб иг с разбиением двух его граней на симплексы; (б) разбиение куба Пг на 3! = 6 симплексов из множества Т(иг) с указанием перестановок; (в) разбиение конуса х3 ^ х2 ^ х1 ^ 0 на симплексы триангуляцией К1
!-мерных точек из мй с целыми координатами. Для произвольной точки г £ 1й рассмотрим с!-мерный единичный куб
= { х £ мй | г, < х, < г, + 1,г = 1,...,! }
в пространстве мй (см. рис. 1, а). Очевидно, что пространство мй совпадает с объединением кубов иг по всем г £ 1й.
Далее для любой перестановки п = (пх,...,п^) чисел 1,...,! определим симплекс тх(п) как выпуклую оболочку (! + 1)-й точки г, ух(п),..., уа(п) с целыми координатами в мй, где у1(п) = г+иП1 и у,(п) = уг_х(п) + для всех г = 2,...,!. Обозначим через Т(Пх) множество симплексов тг(п) для всех перестановок п чисел 1,..., с! (см. рис. 1, б) и через Кх объедиение множеств Т(Пг) по всем г £ 1й. Полное доказательство того, что Кх является триангуляцией пространства мй, можно найти в [9-11], далее потребуется только часть этого доказательства, представленная следующей леммой о симплексе тх (п) триангуляции Кх, содержащем некоторую точку х пространства мй.
Лемма 1. Пусть х £ мй, тогда х £ тг (п) £ Кх, где координаты г являются целыми частями координат х, то есть г, = [х,\ для всех г = 1, ..., с!, а перестановка п определяется порядком убывания координат хп — вектора х — г, то есть хп — ^ х^+1 — для всех г = 1, .. ., с! — 1.
Доказательство. Введем обозначения а, = хП. — гП. для всех г = 1,... ,с!, а также вй = ай и в% = а, — х для всех г = 1,... ,с! — 1. По построению в% ^ 0 для всех г = 1 , . . . , ! и
й
У^вг = (а^ — а^-х) + (а^-х — 2) +-----+ (а^— — а.а) + а.а = а^ < 1.
г=3
Легко убедиться, что выпуклая комбинация вершин 2, ... , Уй(п) симплекса
й
т2(п) с коэффициентами 1 — ^^А,вь... совпадает с х, действительно,
¿=1
(1 — ^ + в^¿(п) = г + ^ a¿Mп¿ = х. ¿=1 ¿=1 ¿=1
Таким образом, х € т2 (п).
□
Триангуляция К пространства сама по себе хорошо подходит для целого ряда существующих симплициальных алгоритмов [9, 10]. Однако для ее использования в алгоритме, который описывается в данной работе, требуется доказать, что одновременно со всем пространством мй она триангулирует также конусы СП = = { х € кй | хП1 ^ хП2 ^ • • • ^ хПа ^ 0 } для перестановок п чисел 1,..., й. Конусы СП естественным образом возникают при разбиении кубов и на симплексы множества Т(иг) (см. рис. 1, б, в), а триангуляции ТП конусов СП, образованные из симплексов К, описываются следующей леммой.
Лемма 2. Множество симплексов ТП гуляцией конуса СП .
{ т € К1 | т С СП } является триан-
Доказательство. Рассмотрим произвольную точку х € СП. Точка х однозначно определяет точку г € такую, что 2 = для всех i = 1,..., й. Так
как из хП1 ^ хП2 ^ • • • ^ хПа ^ 0 следует, что и 2П1 ^ 2П2 ^
> > 0, то
2 € СП. Точки х и 2, в свою очередь, однозначно определяют непересекающиеся
п
множества индексов /1,..., /п, У /¿- = {1,..., й}, таких что x¿ — 2 = х^- — ^ для
¿=1
всех ^ ^ € , к = 1,.. ., п, и x¿ — 2 > х^- — 2^- для всех i € /^, € /^ , k¿ > к^-. Выберем перестановку £ индексов 1,..., й таким образом, чтобы индексы каждого из множеств /^, к = 1, ...,п, следовали в ней в том же порядке, в котором они следуют в перестановке п, и индексы множеств /¿- предшествовали индексам множеств /^ для всех к = 1,. .., п, ] > к.
По лемме 1 точка х принадлежит симплексу т2 (£) триангуляции К. Докажем, что т2 (£) С СП. Идея доказательства состоит в том, что при последовательном переходе сначала от вершины 2 к вершине у1, а затем от вершины у к вершине ^¿+1 симплекса т2 (£) происходит последовательное увеличение координат точки 2 в порядке, определяемом перестановкой £, но это не приводит к выходу вершин симплекса из конуса СП. Например, если £ = п, то справедлива следующая цепочка вытекающих друг из друга неравенств:
2П1 > 2П2 > . . > > 0,
2П1 +1 > 2П2 > . . > > 0,
2П1 +1 > 2П2 +1 > . . > > 0,
2П1 +1 > 2П2 +1 > . . > + 1 > 0.
Само доказательство проведем методом от противного. Предположим, что вершина ут симплекса т2 (£) в первый раз выходит из конуса СП при т-м увеличении £т-й координаты точки 2, то есть предположим, что найдется индекс 1 ^ т ^ й, такой что ут / СП, но у € СП для всех 1 ^ i < т. Рассмотрим индексы п = £т и = . Случай к < т здесь невозможен, потому что при этом увеличение
-й координаты вектора 2 предшествовало бы увеличению его п-й координаты.
й
й
й
В случае т > к, в свою очередь, невозможно, чтобы индексы £т и принадлежали одному и тому же множеству /j, поскольку по построению это означало бы выполнение обратного неравенства т < к .
Предположим, наконец, что £т и £к принадлежат разным множествам /j. Тогда х^^ ^ хП€, но хП—1 — [х^^ ] < xПí — [хп ]. Такая комбинация условий возможна только при одновременном выполнении неравенств х^^ > хп. и [хП4-1 ] > [хп ] . Второе из этих неравенств означает, что г^^ > гп.. Но так как г £ 1й, то верно также и неравенство г^^ ^ г^ + 1 .В свою очередь, так как ут_х £ Сп (или г £ Сп, если т = 1), то из этого неравенства следует, что и ут £ Сп. Таким образом, т-е увеличение £т-й координаты точки г не может приводить к выходу вершины ут симплекса т2 (£) из конуса Сп и доказывает, что т2 (£) С Сп. Отсюда следует, что конус Сп содержится в объединении симплексов из Т-. Обратное включение очевидно. □
Иллюстрация доказательства леммы 2 с указанием перестановок п и £, расположения конуса Сп, точек х и г, а также симплекса т2 (£) приведена на рис. 1, в. На этом же рисунке штриховкой выделены аффинные плоскости хз = 1 и хз = 2, участвующие в образовании из конуса Сп симплексов, соответственно, 1 ^ хз ^ х2 ^ хх ^ 0 и 2 ^ хз ^ х2 ^ хх ^ 0. Выделенные гиперплоскости иллюстрируют, что триангуляция Кх выполняет симплициальное разбиение не только всего пространства и конусов Сп, но еще и симплексов = {х £ Сп | хпл ^ О}, где О £ 1+ . Верно следующее утверждение.
Следствие 1. Для любого О £ 1+ множество симплексов Т= {т £ Кх | т С } является триангуляцией симплекса .
Доказательство. Предположим, что в найдется точка у, не принадлежащая ни одному симплексу из ТП° . Поскольку симплексы замкнуты, а их количество в ТП° конечно, без ограничения общности можно считать, что у не принадлежит гиперплоскости хПл = О, то есть, уПл < О. В то же время, так как у £ Сп, по лемме 2 найдется симплекс т £7! такой, что у £ т. По построению пересечение любого симплекса триангуляции Кх с гиперплоскостью хПЛ = О или пусто или является его собственной гранью. Так как у £ т и уПЛ < О, это означает, что т С и, следовательно, т £ . Полученное противоречие доказывает, что БП° содержится в объединении симплексов из . Обратное включение очевидно. СИ
2. Симплициальный многогранник, основанный на триангуляции К-±
Согласно общей схеме [18], любой алгоритм с целочисленными метками строится на основе симплициального многогранника, содержащего во внутренности начало координат. Настоящая работа предполагает, что для того чтобы построенный алгоритм получился вычислительно эффективным, такой многогранник должен иметь также большое количество регулярно расположенных вершин. Большое количество регулярно расположенных вершин имеет куб, а триангуляция Кх , разбивая поверхность единичных кубов на симплексы таким образом, что множество вершин симплексов этого разбиения совпадает с множеством вершин куба, превращает кубы в потенциальные симплициальные многогранники.
Поэтому в качестве основы для построения симплициального многогранника с большим количеством регулярно расположенных вершин естественно взять именно куб. Пусть = { х £ мй | —1 ^ х, ^ 1, г = 1,...,!}. Вершинами куба являются точки в, координатами которых могут являться только ±1, поэтому далее будем называть такие точки в векторами (ненулевых) знаков. Подходящий
симплициальный многогранник будем искать в виде выпуклой оболочки векторов
ь
где в - вектор знаков, к - количество отрицательных координат вектора в, а Ь и 2 - некоторые параметры. Справедлива следующая лемма о коэффициентах 1 + , участвующих в образовании векторов в'.
Лемма 3. Пусть 2>Ь> 0, 2> 1, х, у ^ 0, к ^ 1 и х + у + к ^ 2, тогда
1
11
<
1
1 + ь/^у 1 + 1 + ь/2ж 1 + ь/
х+й
при условии, что у < х.
Доказательство. Рассмотрим
А =
В =
1
1
ь/2у+Й — ь/2у
ь (1/2У+Й — 1/2У)
1 + ь/2У 1 + (1 + ь/2У )(1 + ) 1 + ь/2У + + ь2/22у+й:
ь/2Ж+Й — ь/2
ь (1/2Ж+Й — 1/2Ж)
1 + ь/2ж 1+ ь/2
■х+й
(1 + ь/2х) (1 + ь/2ж+Й) 1 + ь/2х + ь/2ж+Й + ь2/22
2х+й
А — В
где С = Ас — Вс и
ьс
(1 + ь/2У) (1 + ь/2У+й) (1 + ь/2х) (1 + ь/2ж+й ) '
.1 1 \ / ь ь ь2 Ас = — — — 1 + — + +
1
+
ь
+
ь
+
ь2
1
ь
ь
ь2
У+й 2х+у+й 2х+у+2й 22х+у+2й 2х+у 2х+у+й 22х+у+й '
Вс =
1 п/ ь ь ь2 \
2х+й 2х ) ( + + + 22у+й )
1
2х
+
ь
+
ь
+
ь2
1
ь
ь
ь2
С = Ас — Вс =
1
1
— — +
ь2
ь2
1
1
+ — —
ь2
+
ь2
ь2
1
= -(- — и + ,
11
----+
2х
ь V 1
1—
1
1 — 1
ь2
ь2
1—1
1 — !|Ц —
2х
ь2
х+у+й
2
2й
Так как 2> 1, х, у ^ 0, у<х, то 1/2У — 1/2х > 0. Условия 2>ь, х + у + к ^ 2 означают, что 1 — ь2/2х+у+й > 0. Наконец, из неравенств 2 > 1, к ^ 1 вытекает, что 1/2й — 1 < 0. Таким образом, С < 0. Так как 2>ь> 0, то из С< 0 следует, что и А — В < 0. Доказательство леммы завершено. □
1
1
X
и
2
2
1
к
х
Г \2
Рис. 2. Построение многогранника и± 1 в пространстве к3 : (а) кубы и о и и±1, симплексы то (е) и т± , гиперграни т+ , , т_ , ; (б) и (в) расположение вершин многогранника и± 1 относительно вершин куба и±1; (г) поляра 0± 1 к многограннику и± 1
Построение симплициального многогранника из куба Ц±1 в пространстве К3 иллюстрируется на рис. 2, а, б, в, а его существование и свойства в общем случае устанавливаются следующей теоремой.
Теорема 1. Пусть в 1, .. ., - векторы знаков, являющиеся вершинами куба и±1, пусть, далее, вершины в0,.. ., куба Цо , соответствующие векторам в1,. .., в2'1, получаются путем замены отрицательных компонент в 1, .. ., на нулевые и пусть векторы . .., , соответствующие векторам 31,...,32^ и в0,. .., , определяются как
4 = (1+^, 2 = 1>...,
где Ь и г - скалярные величины такие, что г > Ь > 0, а к^ - число отрицательных компонент вектора знаков з^, 0 ^ kj ^ 1. Тогда
1) выпуклая оболочка и± векторов з^, . .., з2^ является симплициальным многогранником, то есть многогранником, все собственные грани которого являются симплексами;
2) выпуклая оболочка 1 -элементного множества векторов sj, 2 £ J, является гипергранью или (1 — 1) -мерной гранью многогранника и± тогда и только тогда, когда выпуклая оболочка векторов з0 , ] £ J, соответствующих векторам sj , 2 £ J, образует на поверхности куба и о гипергрань симплекса триангуляции К •
Доказательство. Без ограничения общности достаточно рассмотреть симплекс то(е) € Т(Цо) С К для тождественной перестановки е чисел 1,...,!. В кубе и±1 симплексу то (е) соответствует симплекс т±, являющийся выпуклой оболочкой векторов знаков
«о = (-1,-1,..., -1,-1),
51 = ( + 1, -1,..., -1, -1),
«2 = ( + 1, +1,..., -1,-1),
в^-1 = (+1, +1,..., +1,-1), = (+1, +1,..., +1, +1).
Симплекс т± имеет в точности две гиперграни, лежащие на поверхности куба и±1: это выпуклая оболочка т+ векторов знаков в1,..., в^ и выпуклая оболочка т- векторов знаков во,. , содержащиеся соответственно в гипергранях
= { х € и±1 | х1 = 1} и В- = { ж € и±1 | х^ = -1 } куба и±1. Согласно необходимым и достаточным условиям, полученным в [20, 21], для доказательства того, что выпуклая оболочка т^ векторов в1,..., в^, соответствующих векторам знаков в1,..., в^, и выпуклая оболочка т- векторов во,..., «¿-1, соответствующих векторам знаков во,..., в^-, являются гранями многогранника и±1, достаточно доказать, что сумма коэффициентов разложения любой вершины V' многогранника Ц±1 по векторам в^,..., в^ и во,..., меньше 1, если V' / т' , то есть V' = в^,...,в^, или V / т-, то есть V' = во,..., £¿-1.
Линейные комбинации, образующие разложения возможных вершин многогранника и±1, будем искать как линейные комбинации векторов знаков в1,..., в^ и во,..., , образующие разложение некоторой вершины V куба Ц±1. Причем при рассмотрении симплекса т+ С , то есть при анализе линейных комбинаций векторов в1,..., в^ будем считать, что V € , а при рассмотрении симплекса т- С В-, то есть при анализе линейных комбинаций векторов во,..., в^-, будем считать, что V € В-. В общем виде описанная вершина куба Ц±1 может быть представлена !-мерным вектором
Х1 X У
V = (+1,..., +1,-1,..., -1, +1,..., +1,...,'-!,.",-1, +1,..., +1,-1,..., -1).
Вектор V начинается, возможно, с цепочки (идущих подряд) +1, в частности, он обязательно начинается по крайней мере одной +1, если V € . Затем в векторе V располагается цепочка -1, затем цепочка +1 и т. д. Завершается вектор V, возможно, цепочкой -1, в частности, он обязательно завершается по крайней мере одной -1, если V € В-. Предположим, что всего вектор V имеет I цепочек +1, которым предшествуют цепочки -1. Число цепочек -1, за которыми следуют цепочки +1, тоже получается равным I, будем считать, что эти цепочки -1 имеют длины Ж1,..., ж;. Будем считать также, что длина заключительной I-й цепочки +1 равна к, а длина заключительной цепочки -1, которая может присутствовать в векторе V, равна у.
Доказательство проведем индукцией по количеству цепочек I вектора V. Условие I = 0 в случае V € равносильно тому, что V € т+ , то есть вектор V совпадает с одним из векторов в1,..., в^, а в случае V € В- - тому, что V € т- , то есть V совпадает с одним из векторов во,..., в^-1. В любом случае в разложении V участвует единственный вектор, а сумма коэффициентов разложения равна 1 , причем то же самое, очевидно, верно и для разложения вектора V', соответствующего вектору V, по векторам в^,..., в^ или во,..., 1.
Предположим, что для любого вектора знаков и £ У й-, количество цепочек I которого меньше или равно п, п ^ 0, сумма коэффициентов разложения вектора и', соответствующего и, по векторам ..., , если ад £ , и по векторам ..., в£г-1, если и £ , меньше или равна 1. Покажем, что в таком случае сумма коэффициентов разложения вектора V', соответствующего вектору знаков V £ У с количеством цепочек I = п +1, по тем же самым векторам в1,..., или ..., в£г-1 строго меньше 1.
Воспользуемся очевидным представлением вектора V в виде V = и + — - в^-у-й, где вектор
и = (+1,..., +1,—1,7. ."7—1 +1,..., +1,...,—1,.. ,7—1, —1,..., —1, —177. .7—1)
имеет ту же структуру, что и вектор V, только с меньшим количеством цепочек +1 , которым предшествуют цепочки —1 .
Коэффициенты разложения вектора V' по векторам и', и найдем
следующим образом.
V — и = в^-у — в^-у-й,
(1 + Ь/г*) (1 + Ь/г*+А;) =(1 + Ь/гУ ) (1 + )
(1 + Ь/г*)V — (1 + Ь/г*+й) и = (1 + Ь/гУ) — (1 + Ь/гУ+Й) ,
V'___и' = ___^-у-й
(1 + Ь/г*) (1 + Ь/,г*+Й) = (1 + Ь/гУ) (1 + ),
1 + Ь/г* N ' /1 + Ь/г* \ ' / 1 + Ь/г* \ '
и + К , ,/.,, ^-у — К , ^-у-й, (1)
1 + и \1 + Ь/г^ V1 + ,
где х = Х1 + • • •+х; + у. Поскольку I = п +1 > 0, то вектор V имеет по крайней мере одну цепочку +1, которой предшествует цепочка — 1. Это означает, что х, к ^ 1, х > у, учитывая, что у ^ 0, находим также, что х + у + к ^ 2. Таким образом, для чисел г, Ь, х, у и к выполнены все условия леммы 3, из которой следует, что сумма коэффициентов разложения вектора V ' по векторам и ', и
меньше 1 , то есть
1 + Ь/г* 1 + Ь/г* 1 + Ь/г*
1 + ь/,г*+й + 1 + ь/гу 1 + < 1 (2)
Пусть и обозначают суммы коэффициентов разложения векторов V' и и ' по векторам , если V, и £ , и по векторам ..., 1, если
V, и £ й-. Тогда из (1) имеем
1 + Ь/г* \ 1 + Ь/г* 1 + Ь/г* ( )
1 + + 1 + Ь/гУ — 1 + Ь/гУ+й . (3)
В свою очередь, поскольку по предположению индукции стш' ^ 1, то из (2) и (3)
г+(С ), V £
немедленно вытекает, что ст^ < 1. Для случаев т+(с ), V £ и т-(с й-),
V £ й- теорема доказана.
(с ^ £ +
ся к уже рассмотренным. Действительно, если V £ , то V = и — 2^1, где и £ и 2^1 = в1 + , а если V £ й-, то V = и + , где и £ й- и = —во — . В любом случае сумма коэффициентов разложения будет меньше суммы коэффициентов разложения ст^'. Но, если и £ т+ в случае и £ или и £ т- в случае и £ й- , то = 1, иначе, как только что было доказано, стш' < 1. Таким образом, а^' < < 1. □
Оставшиеся случаи т+ (с ), V £ и т-(с й-), V £ й- легко сводят-
3. Сильное и слабое приближение
Согласно общей схеме [18], симплициальный -алгоритм перемещается в пространстве, используя многогранник Ц±1, а управляет движением и проверяет, достигнуто ли приближение неподвижной точки, используя многогранник О±, являющийся полярой к и±1. По свойствам поляры [22] гиперграни О± находятся во взаимно-однозначном соответствии с вершинами Ц±1, а значит, и с векторами знаков, определяющими вершины и±1. Так как многогранник Ц±1 является сим-плициальным, из этого соответствия следует также, что многогранник О± является простым: в каждой его вершине пересекаются в точности ! гиперграней. В пространстве м3 многогранник О± изображен на рис. 2, г.
В процессе работы симплициальный 2^-алгоритм использует 2^ целочисленных меток, которые присваиваются вершинам симплексов триангуляции в соответствии с правилом помечивания. Эти метки задаются фактически векторами знаков, поэтому далее различия между векторами знаков и метками не делается. Используя многогранник О±, правило помечивания симплициального 2^ -алгоритма можно сформулировать следующим образом: вершине триангуляции х присваивается метка в, если вектор /(ж) - ж, где f (ж) - значение целевой функции в точке х , принадлежит конусу метки в , определяемому как коническая оболочка гиперграни многогранника О± 1, соответствующей метке в.
На рис. 2, г пунктирными линиями построены конус метки (+1, +1, +1), соответствующей гиперграни , и конус метки (-1, +1, -1), соответствующей гиперграни В, этот конус дополнительно выделен штриховкой, а также частично построены конусы меток (-1, +1, +1) и (+1, +1, -1) . Эти четыре построенных конуса демонстрируют важное отличие в степени точности найденного приближения, которое появляется при переходе от (относительно) малого количества меток к (относительно) большому. В частности, нахождение симплициальным 2^-алгоритмом всех четырех меток (+1, +1, +1), (-1, +1, -1), (-1, +1, +1) и (+1, +1,-1) означает, что достигнуто приближение неподвижной точки.
Приближение достигается всякий раз, когда пересечение гиперграней многогранника О±, соответствующих найденным меткам пусто, или, что то же самое, пересечение конических оболочек этих гиперграней совпадает с началом координат О, в том числе это так и для меток (+1, +1, +1), (-1, +1, -1), (-1, +1, +1) и (+1, +1, -1). Однако точность достигаемого при нахождении этих меток приближения определяется минимальным углом между векторами, лежащими в частично построенных конусах меток (-1, +1, +1) и (+1, +1, -1), — пересечение этих конусов само по себе совпадает с началом координат О . На рис. 2, г это угол между векторами V! и V2, выделенный штриховкой и являющийся достаточно малым. Точность приближения также будет здесь невысокой.
Для сравнения заменим одну из меток (-1, +1, +1) или (+1, +1, -1) на метку (-1, -1, -1). Конусом метки (-1, -1, -1) является коническая оболочка гиперграни многогранника О± , расположенная «под» гипергранью и, так же как и гипергрань , имеющая форму неправильного шестиугольника с параллельными противолежащими ребрами. Полученный набор меток по-прежнему дает приближение неподвижной точки, но теперь уже сами по себе пересекаются в начале координат О конусы меток (+1, +1, +1) и (-1, -1, -1). Понятие о минимальном угле между векторами, лежащими в этих конусах, дает «толщина» выделенного на рис. 2, г штриховкой конуса метки (-1, +1, -1), очевидно, что этот угол значительно превосходит угол между векторами Vl и V2.
Таким образом, в связи с появлением в симплициальном 2^-алгоритме большого (по сравнению с размерностью пространства) количества меток точность приближения неподвижной точки начинает значительно меняться в зависимости
(в)
Я
■-1 0 . .0 0
+1 -1 . .0 0
0 0 . . +1 -1
0 0 . .0 +1
Рис. 3. Различные объекты триангуляции в пространстве
(б) симплекс Яе при С = 4; (в) при С ^ триангуляция Яе переходит в триангуляцию конуса Се; (г) триангуляция стандартного симплекса Д2 мелкости С = 4 получается аффинным преобразованием с матрицей Я из симплекса в пространстве К2
от набора меток, на котором достигается это приближение. Будем называть приближение сильным, если среди найденного набора меток присутствуют метки и такие, что = — , и будем называть приближение слабым в противном случае. Понятие сильного и слабого приближения неподвижной точки влияет на организацию перезапусков или рестартов алгоритма из найденного приближения с целью нахождения приближения более высокой степени точности.
Если найдено сильное приближение (как в примере с векторами знаков (+1, +1, +1) и ( — 1, —1, —1)), то неподвижная точка, скорее всего, находится где-то недалеко, поэтому имеет смысл увеличить мелкость триангуляции при перезапуске алгоритма. Если найдено слабое приближение (как в примере с векторами знаков ( —1, +1, +1) и (+1, +1, —1)), то, скорее всего, вблизи этого приближения целевая функция всего лишь достаточно резко меняет свое поведение, поэтому имеет смысл перезапустить алгоритм из найденного приближения с прежней мелкостью разбиения и искать неподвижную точку в новом направлении.
В связи с появлением понятия сильного и слабого приближений следует обратить внимание также на пространственную асимметрию симплициального -алгоритма. По построению многогранник Ц±1 имеет выделенное направление от вектора знаков со всеми —1 до вектора знаков со всеми +1 .В пространстве К3 это направление задается диагональю ( — 1, —1, —1) - (+1, +1, +1) куба Ц±1. То же выделенное направление присутствует и в многограннике - на рис. 2, г гиперграни и В2 многогранника , соответствующие векторам знаков (+1, +1, +1) и ( — 1, —1, —1), имеют форму шестиугольников, не пересекаются между собой и пересекаются со всеми оставшимися гипергранями, имеющими форму четырехугольников (такими как, например, гипергрань В).
В результате такого расположения гиперграней многогранника симпли-
циальный 2^-алгоритм предположительно оказывается несколько более «предрасположенным» к нахождению слабого приближения, если в наборе его меток
присутствует вектор знаков, состоящий только из -1 или только из +1 . Однако геометрически все диагонали куба одинаковы. Это дает дополнительную возможность управлять поведением алгоритма в процессе работы: связывая перед стартом или рестартом алгоритма выделенное направление в пространстве с одной из диагоналей куба, можно в определенной степени повлиять на скорость нахождения и характер найденного приближения неподвижной точки.
В частности, изменение выделенного направления может быть использовано для уменьшения влияния пространственной асимметрии симплициального 2^-ал-горитма на результаты вычислений. В общем случае все расположения вектора /(х) — х в пространстве являются равновероятными, поэтому, скорее всего, он будет принадлежать конической оболочке некоторой гиперграни й многогранника 0± 1, отличной от и , например такой, как на рис. 2, г. Но в момент старта или рестарта алгоритма указанное равновероятное расположение вектора f (х) — х может быть практически в равной степени обеспечено просто выбором выделенного направления, что и было сделано в [18] при оценке скорости предварительной версии 2^-алгоритма в пространстве м4 .
4. Заключительные замечания
Триангуляция К была впервые описана Фройденталем [19] в 1942 г. Однако по замечанию самого Фройденталя это описание уже использовалось им [23] в 1939 г. и основано на его же статье [24] о триангуляции декартова произведения симплексов, вышедшей еще двумя годами ранее. Таким образом, триангуляция К пространства м^ известна уже около 80 лет - почти столетие. Далее К встречается в [25, с.140], [8, 9, 11, 26] и других работах, поэтому к настоящему времени она хорошо изучена. Обозначение К дал ей Тодд в книге [9], посвященной сим-плициальным алгоритмам. Согласно системе обозначений, принятой в [9], индекс 1 говорит о том, что К1 триангулирует все пространство м^.
Триангуляция К замечательна тем, что она триангулирует одновременно (1) единичные кубы V; (2) кубы Vа = {х £ м^ | 0 ^ х4 ^ С}, что может быть также интерпретировано как триангуляция единичных кубов V мелкости С; (3) все пространство м^, что также следует из (2) при С ^ , таким образом, (1) ^ (2) ^ (3); (4) симплексы , что может быть также интерпретировано как триангуляция симплекса ¿Л мелкости С; (5) стандартный симплекс
I 1 1
х £ 11 = 1, ^ 0, г = 1,..., й +1 > , получаемый из симплекса
¿=1 ) аффинным преобразованием и характерный для теоретико-игровых, например [27], и экономических, например [28], задач; (6) симплициальные конусы , что также следует из (4) при С ^ , таким образом, (4) ^ (5) ^ (6).
Все эти свойства являются, с одной стороны, «геометрически очевидным» (рис. 3), с другой стороны, достаточно трудоемкими для доказательства. Например, доказательство только (5) в [28] занимает 8 страниц текста, что сопоставимо по объему с отдельной статьей средних размеров. Поэтому каждая из работ (в том числе и оригинальная работа Фройденталя [19]), имеющая дело с триангуляцией К1, описывает ее лишь в той степени, в какой это необходимо для целей этой работы. Хорошая схема (3) ^ (6) ^ (4) доказательства всех описанных свойств триангуляции К была предложена в [9], однако реализации ее представлено не было. В результате до настоящего времени доказательства свойств триангуляции К являются разбросанными по различным источникам.
Так, в [9] и [11] доказано (3). В [26] доказано (1) и (2) .В [28] доказано (5). Лемма 2 данной работы доказывает переход (3) ^ (6) от группы свойств (1), (2),
(3), имеющих дело преимущественно с кубами, к группе свойств (4), (5), (6), имеющих дело преимущественно с симплексами. Что касается оригинальной работы Фройденталя [19], то она дает очень компактное описание перехода (1) ^ (4) для $П. По своей сути переход (1) ^ (4) в [19] является аналогом перехода (3) ^ (6) леммы 2, более того, он также использует перестановки, аналогичные ее перестановкам п и £. Однако имеются и существенные отличия, в частности, в лемме 2 перестановка £ определяется перестановкой п, а в работе [19], наоборот, в результате доказанными оказываются аналогичные по смыслу, но разные по способу его выражения утверждения.
Итак, лемма 2 является реализацией звена (3) ^ (6) (а вместе со следствием 1 и всей) схемы (3) ^ (6) ^ (4) доказательства, предложенной в [9], и интерпретацией части написанной в достаточно свободной манере работы Фройденталя [19] в более привычном для современной математики виде. Центральный результат работы -теорема 1 и лемма 2 дают полное обоснование нового симплициального алгоритма с 2d целочисленными метками, а описание понятий сильного и слабого приближения делает возможным создание программной реализации этого алгоритма в версии, хорошо зарекомендовавшей себя в [18] при вычислениях в пространстве R4 .
Литература
1. Scarf H. The approximation of fixed points of a continuous mapping // SIAM J. Appl. Math. - 1967. - V. 15, No 5. - P. 1328-1343. - doi: 10.1137/0115116.
2. Scarf H. The core of an N person game // Econometrica. - 1967. - V. 35, No 1. -P. 50-69. - doi: 10.2307/1909383.
3. Scarf H. On the computation of equilibrium prices // Ten Economic Studies in the Tradition of Irving Fisher. - N. Y., Sydney: John Wiley and Sons, 1967. - P. 207-230.
4. Hansen T., Scarf H. On the applications of a recent combinatorial algorithm: Cowles Foundation Discussion Papers No. 272. - New Haven: Cowles Foundation for Research in Economics, Yale University, 1969. - 43 p.
5. Хансен Т., Скарф Г. О приложениях нового комбинаторного алгоритма // Математическая экономика: Равновесные модели, оптимизационное планирование и управление / Ред. Б.С. Митягин. - М.: Мир, 1974. - С. 143-169.
6. Cohen D.I.A. On the Sperner lemma // J. Comb. Theory. - 1967. - V. 2. - P. 585-587.
7. Sperner E. Neuer beweis fur die invarianz der dimensionszahl und des gebietes // Abh. math. Seminar Univ. Hamburg. - 1928. - Bd. 6. - S. 265-272.
8. Kuhn H.W. Simplicial approximation of fixed points // Proc. Natl. Acad. Sci. U. S. A. -1968. - V. 61, No 4. - P. 1238-1242. - doi: 10.1073/pnas.61.4.1238.
9. Todd M.J. The computation of fixed points and applications // Lecture Notes in Economics and Mathematical Systems. - 1976. - V. 124. - vii, 129 p.
10. Тодд М.Дж. Вычисление неподвижных точек и приложения к экономике. - М.: Наука, 1983. - 112 c.
11. Yang Z. Computing equilibria and fixed points. - Boston: Kluwer Acad. Publ., 1999. - x, 344 p.
12. Yamamoto Y. A variable dimension fixed point algorithm and the orientation of simpli-ces // Math. Program. - 1984. - V. 30, No 3. - P. 301-312. - doi: 10.1007/BF02591935.
13. Kojima M., Yamamoto Y. A unified approach to the implementation of several restart fixed point algorithms and a new variable dimension algorithm // Math. Program. -1984. - V. 28. - P. 288-328. - doi: 10.1007/BF02612336.
14. Talman A.J.J., Yamamoto Y. A simplicial algorithm for stationary point problems on polytopes // Math. Oper. Res. - 1989. - V. 14, No 3. - P. 383-399. - doi: 10.1287/moor.14.3.383.
15. Wright A.H. The octahedral algorithm, a new simplicial fixed point algorithm // Mathematical Programming. - 1981. -V. 21, No 1. - P. 47-69. - doi: 10.1007/BF01584229.
16. van der Laan G., Talman A.J.J. A class of simplicial restart fixed point algorithms without an extra dimension // Math. Program. - 1981. - V. 20, No 1. - P. 33-48. -doi: 10.1007/BF01589331.
17. Reiser P.M. A modified integer labeling for complementarity algorithms // Math. Oper. Res. - 1981. - V. 6, No 1. - P. 129-139. - doi: 10.1287/moor.6.1.129.
18. Matveev M.N. An approximation of fixed points with integer labels // Proc. 18th IFAC World Congress (Milano). - 2011. - P. 10174-10179.
19. Freudenthal H. Simplizialzerlegungen von Beschrankter Flachheit // Ann. Math. - 1942. -V. 43. - P. 580-582.
20. Матвеев М.Н. Невидимые гиперграни и порождающие многогранники // Труды МФТИ. - 2011. - T. 3, № 1. - С. 102-106.
21. Матвеев М.Н. Минимальный непорождаемый веер // Учен. зап. Казан. ун-та. Сер. Физ.-матем. науки. - 2012. - T. 154, кн 1. - С. 202-207.
22. Бренстед А. Введение в теорию выпуклых многогранников. - М.: Мир, 1988. - 240 c.
23. Freudenthal H. Die Triangulation der differenzierbaren Mannigfaltigkeiten // Proc. Konink. Acad. Wetensch. Amsterdam. - 1939. - V. 42. - P. 880-901.
24. Freudenthal H. Eine Simplizialzerlegung des Cartesischen Produktes zweier Simplexe // Fund. Math. - 1937. - V. 29. - P. 138-144.
25. Lefschetz S. Introduction to topology. - Princeton: Princeton Univ. Press, 1949. - 228 p.
26. Kuhn H.W. Some combinatorial lemmas in topology // IBM J. Research and Development. - 1960. - V. 4, No 5. - P. 518-524. - doi: 10.1147/rd.45.0518.
27. Shapley L.S. On balanced games without side payments // Hu T.C., Robinson S.M. (Eds.) Mathematical Programming. - N. Y.: Acad. Press, 1972. - P. 261-290.
28. Scarf H., Hansen T. The computation of economic equilibria. Monograph 24, Cowles Foundation for Research in Economics. - New Haven and London: Yale Univ. Press, 1973. - x, 249 p.
Поступила в редакцию 17.09.18
Матвеев Михаил Николаевич, кандидат физико-математических наук, доцент кафедры высшей математики
Московский физико-технический институт
Институтский пер., д. 9, г. Долгопрудный, Московская обл., 141701, Россия E-mail: miklem@mail.mipt.ru
ISSN 2541-7746 (Print) ISSN 2500-2198 (Online) UCHENYE ZAPISKI KAZANSKOGO UNIVERSITETA. SERIYA FIZIKO-MATEMATICHESKIE NAUKI (Proceedings of Kazan University. Physics and Mathematics Series)
2019, vol. 161, no. 1, pp. 127-144
doi: 10.26907/2541-7746.2019.1.127-144
Simplicial Fixed Point Algorithm with 2d Integer Labels
M.N. Matveev
Moscow Institute of Physics and Technology, Dolgoprudny, 141701 Russia E-mail: miklem@mail.mipt.ru
Received September 17, 2018 Abstract
Simplicial fixed point algorithms can be based on either integer or vector labels. An apparent advantage of the algorithms with integer labels is their simplicity and stability under round-off errors due to the discrete nature of integer labels. At the same time, the application of all existing algorithms with integer labels is limited by some rigidness of their construction, in particular, in the d-dimensional space they can bear only from d + 1 to 2d labels. The numbers of labels comparable with the dimension of space make these algorithms not very fast, especially in high-dimensional tasks. This paper overcomes this rigidness and builds a new simplicial fixed point algorithm with 2d integer labels. To achieve this goal, the paper proves new properties of triangulation K1 and separates all approximations of fixed points into weak and strong ones. This separation, which has never been used until this moment, is caused by the high number of labels of the new algorithm.
Key words: triangulations, polytopes, fans, simplicial algorithms, fixed point algorithms
Figure Captions
Fig. 1. Some shapes of the triangulation K1 in r3: (a) a unit cube Uz with two facets triangulated; (b) all 3! = 6 simplices and their permutations in the triangulation T(Uz) of the cube Uz; (c) the subdivision of the cone x3 ^ x2 ^ x1 ^ 0 induced by the triangulation K1.
Fig. 2. The polytope U± 1 in r3: (a) the cubes UO and U±1, simplices to(e) and t± , facets t+ , F+1, t- , F- ; (b) and (c) vertices of the polytope U^ 1 with respect to vertices of the cube U±1; (d) the polar O± 1 to the polytope U± 1.
Fig. 3. Objects to be triangulated in r3: (a) the cube UG with G = 4; (b) the simplex
with G = 4; (c) the triangulation of the cone Ce obtained from the triangulations of SG when G ^ ; (d) the triangulation of the standard simplex A2 with the grid fineness G = 4 obtained by an affine map with the matrix Q from the simplex S1 in r2.
References
1. Scarf H. The approximation of fixed points of a continuous mapping. SIAM J. Appl. Math., 1967, vol. 15, no. 5, pp. 1328-1343. doi: 10.1137/0115116.
2. Scarf H. The core of an N person game. Econometrica, 1967, vol. 35, no. 1, pp. 50-69. doi: 10.2307/1909383.
3. Scarf H. On the computation of equilibrium prices. In: Ten Economic Studies in the Tradition of Irving Fisher. New York, Sydney, John Wiley and Sons, 1967, pp. 207-230.
4. Hansen T., Scarf H. On the applications of a recent combinatorial algorithm: Cowles Foundation Discussion Papers No. 272. New Haven, Cowles Found. Res. Econ., Yale Univ., 1969. 43 p.
5. Hansen T., Scarf H. On the applications of a recent combinatorial algorithm. In: Mitya-gin B.S. (Ed.) Matematicheskaya ekonomika: Ravnovesnye modeli, optimizatsionnoe planirovanie i upravlenie [Mathematical Economics: Equilibrium Models, Optimization Planning, and Control]. Moscow, Mir, 1974, pp. 143-169. (In Russian)
6. Cohen D.I.A. On the Sperner lemma. J. Comb. Theory, 1967, vol. 2, pp. 585-587.
7. Sperner E. Neuer beweis fur die invarianz der dimensionszahl und des gebietes. Abh. math. Seminar Univ. Hamburg, 1928, Bd. 6, S. 265-272. (In German)
8. Kuhn H.W. Simplicial approximation of fixed points. Proc. Natl. Acad. Sci. U. S. A., 1968, vol. 61, no. 4, pp. 1238-1242. doi: 10.1073/pnas.61.4.1238.
9. Todd M.J. The computation of fixed points and applications. Lect. Notes Econ. Math. Syst., 1976, vol. 124. vii, 129 p.
10. Todd M.J. Vychislenie nepodvizhnykh tochek i prilozheniya k ekonomike [The Computation of Fixed Points and Applications]. Moscow, Nauka, 1983. 112 p. (In Russian)
11. Yang Z. Computing Equilibria and Fixed Points. Boston, Kluwer Acad. Publ., 1999. x, 344 p.
12. Yamamoto Y. A variable dimension fixed point algorithm and the orientation of simplices. Math. Program., 1984, vol. 30, no. 3, pp. 301-312. doi: 10.1007/BF02591935.
13. Kojima M., Yamamoto Y. A unified approach to the implementation of several restart fixed point algorithms and a new variable dimension algorithm. Math. Program., 1984, vol. 28, pp. 288-328. doi: 10.1007/BF02612336.
14. Talman A.J.J., Yamamoto Y. A simplicial algorithm for stationary point problems on polytopes. Math. Oper. Res., 1989, vol. 14, no. 3, pp. 383-399. doi: 10.1287/moor.14.3.383.
15. Wright A.H. The octahedral algorithm, a new simplicial fixed point algorithm. Math. Program., 1981, vol. 21, no. 1, pp. 47-69. doi: 10.1007/BF01584229.
16. van der Laan G., Talman A.J.J. A class of simplicial restart fixed point algorithms without an extra dimension. Math. Program1981, vol. 20, no. 1, pp. 33-48. doi: 10.1007/BF01589331.
17. Reiser P.M. A modified integer labeling for complementarity algorithms. Math. Oper. Res., 1981, vol. 6, no. 1, pp. 129-139. doi: 10.1287/moor.6.1.129.
18. Matveev M.N. An approximation of fixed points with integer labels. Proc. 18th IFAC World Congress (Milano), 2011, pp. 10174-10179.
19. Freudenthal H. Simplizialzerlegungen von Beschränkter flachheit. Ann. Math., 1942, vol. 43, pp. 580-582. (In German)
20. Matveev M.N. Invisible faces and face polytopes. Tr. MFTI, 2011, vol. 3, no. 1, pp. 102106. (In Russian)
21. Matveev M.N. A minimal nonpolytopal fan. Uchenye Zapiski Kazanskogo Universiteta. Seriya Fiziko-Matematicheskie Nauki, 2012, vol. 154, no. 1. pp. 202-207. (In Russian)
22. Br0ndsted A. An introduction to convex polytopes. In: Graduate Texts in Mathematics, 1983, vol. 90. viii, 162 p.
23. Freudenthal H. Die Triangulation der differenzierbaren Mannigfaltigkeiten. Proc. Konink. Acad. Wetensch. Amsterdam, 1939, vol. 42, pp. 880-901. (In German)
24. Freudenthal H. Eine Simplizialzerlegung des Cartesischen Produktes zweier Simplexe. Fundam. Math., 1937, vol. 29, pp. 138-144. (In German)
25. Lefschetz S. Introduction to Topology. Princeton, Princeton Univ. Press, 1949. 228 p.
26. Kuhn H.W. Some combinatorial lemmas in topology. IBM J. Res. Dev., 1960, vol. 4, no. 5, pp. 518-524. doi: 10.1147/rd.45.0518.
27. Shapley L.S. On balanced games without side payments. In: Hu T.C., Robinson S.M. (Eds.) Mathematical Programming, New York, Acad. Press, 1972. pp. 261-290.
28. Scarf H., Hansen T. The Computation of Economic Equilibria. Monograph 24, Cowles Foundation for Research in Economics. New Haven and London, Yale Univ. Press, 1973. x, 249 p.
/ Для цитирования: Матвеев М.Н. Симплициальный алгоритм поиска неподвиж-( ных точек с 2d целочисленными метками // Учен. зап. Казан. ун-та. Сер. Физ.-матем. \ науки. - 2019. - Т. 161, кн. 1. - С. 127-144. - doi: 10.26907/2541-7746.2019.1.127-144.
I For citation: Matveev M.N. Simplicial fixed point algorithm with 2d integer labels. ( Uchenye Zapiski Kazanskogo Universiteta. Seriya Fiziko-Matematicheskie Nauki, 2019, \ vol. 161, no. 1, pp. 127-144. doi: 10.26907/2541-7746.2019.1.127-144. (In Russian)