CC BY
133
ЧЕБЫШЕВСКИЙ СБОРНИК
Том 16 Выпуск 2 (2015)
УДК 511.36
УНИВЕРСАЛЬНОЕ ОБОБЩЕНИЕ АЛГОРИТМА ЦЕПНОЙ ДРОБИ
А. Д. Брюно (г. Москва)
Аннотация
1. Простое обобщение. Пусть в трехмерном вещественном пространстве заданы три вещественные однородные линейные формы. Их модули дают отображение этого пространства в другое. В нем рассматривается выпуклая оболочка образов всех целочисленных точек первого пространства, кроме его начала координат. Замыкание этой выпуклой оболочки названо модульным многогранником. Наилучшие целочисленные приближения к корневым подпространствам заданных форм дают точки, образы которых лежат на границе модульного многогранника. Граница модульного многогранника вычисляется любой стандартной программой вычисления выпуклых оболочек. Алгоритм дает также периодичность для кубических иррациональностей с положительным дискриминантом. Обобщить цепную дробь пытались Эйлер, Якоби, Дирихле, Эрмит, Пуанкаре, Гур-виц, Клейн, Минковский, Вороной и многие другие.
2. Универсальное обобщение. Пусть в п-мерном вещественном пространстве Мп заданы I линейных и к квадратичных форм (п = I + 2к). Модули этих форм задают отображение пространства Мп в положительный ортант Б = Мт т-мерного вещественного пространства Мт, т = I + к. При этом целочисленная решётка Ъп в Мп отображается в некоторое множество Z в Б. Замыкание выпуклой оболочки Н множества 2\0 является многогранным множеством. Целочисленные точки из Мп, отображающиеся на границу дН многогранника Н, дают наилучшие диофантовы приближения к совокупности корневых подпространств т заданных форм. В алгебраическом случае, когда заданные формы определённым образом связаны с корнями многочлена степени п, доказывается, что многогранник Н имеет т—1 независимый период. Это обобщение теоремы Лагранжа о периодичности цепной дроби квадратичной иррациональности. По теореме Дирихле соответствующее поле алгебраических чисел имеет ровно т — 1 фундаментальных единиц. Граница дН многогранника Н вычисляется стандартной программой вычисления выпуклых оболочек.
Ключевые слова: цепная дробь, модульный многогранник, программа вычисления выпуклого многогранника.
Библиография: 75 названий.
UNIVERSAL GENERALIZATION OF THE CONTINUED FRACTION ALGORITHM
A. D. Bruno (Moscow)
Abstract
1. Simple generalization. Let three homogeneous real linear forms be given in a three-dimensional real space. Their moduli give a mapping of the space into another space. In the second space, we consider the convex hull of images of all integer points of the first space except its origin. This convex hull is called the modular polyhedron. The best integer approximations to the root subspaces of these forms are given by the integer points whose images lie on the boundary of the modular polyhedron. For the concret three linear forms, any part of the boundary of the modular polyhedron can be computed by means of any standard program for computation of a convex hull. The algorithm gives the best approximations, and it is periodic for cubic irrationalities with positive discriminant. It also allows to understand why matrix algorithms proposed by Euler, Jacobi, Dirichlet, Hermite, Poincare, Hurwitz, Brun, Guting and others are not universal: proper algorithm is composed from several different matrix algorithms.
2. Universal generalization. Let l linear forms and k quadratic forms (n = l + 2k) be given in the n-dimensional real space Rn. Absolute values of the forms define a map of the space Rn into the positive orthant S of the m-dimensional real space Rm, where m = l + k. Here the integer lattice Zn in Rn is mapped into a set Z in S. The closure of the convex hull H of the set Z\0 is a polyhedral set. Integer points from Rn, which are mapped in the boundary dH of the polyhedron H, give the best Diophantine approximations to root subspaces of all given forms. In the algebraic case, when the given forms are connected with roots of a polynomial of degree n, we prove that the polyhedron H has m — 1 independent periods. It is a generalization of the Lagrange Theorem, that continued fractions of a square irrationality is periodic. For the certain set of the m forms, any part of the boundary dH of the polyhedron H can be computed by a program for computing convex hulls.
3. Main achievement. Best Diophantine approximations can be computed by a global algorithm using a standard program for computing convex hulls, instead of step-by-step computations as in the continued fraction algorithm. It gives a solution of the problem, that majority of main mathematicians of the XIX century tried to solve.
Keywords: continued fraction, modular polyhedron, program for computing convex hull.
Bibliography: 75 titles.
1. Введение 1.1. Цепная дробь
Пусть ао и а - натуральные числа. Для нахождения их наибольшего общего делителя используется алгоритм Евклида [1] последовательного деления с остатком:
а0 = а0а1 + а2, а1 = а^2 + а3, а2 = а2а3 + а4,...
где натуральные числа ао,Я1,Я2, • • • суть неполные частные. Это алгоритм разложения числа а = ао/а1 в правильную цепную дробь [2], и он применим к любым вещественным числам а. При этом а0 = [а], где [а] - целая часть числа а, а1 = [1/(а — ао)], . . . , т. е.
1
а = ао +--1-, (1)
а1 +--
1
а2 H--
аз + .
•■(::;;)=(? Л) (а:,) ■ -=а^
Если разложение (1) оборвать на а: и свернуть эту оборванную цепную дробь в рациональное число p:/q:, то получается подходящая дробь, которая дает наилучшее рациональное приближение к числу а. При этом
P: p:-i\ = (P:-i Р:-А /а: Л (3)
q: q:-i ) \q:-i q:-2/ V1 0/
? Il-1 =(? 1 ) ,detfP: M =±1,
1 0 J \1 -a:j Vq: q:-i J
т.е. векторы (а:, a:+i) и (p:,q:) принадлежат сопряженным плоскостям, и пара векторов (p:, q:), (p:-i, q:-i) может служить базисом в одной из них. Цепные дроби (1) и соотношения (2), (3) рассматривал Валлис [3] в 1655 г. В 1737 г. Эйлер [4] дал название "непрерывная дробь" (fractio continua). Лагранж [5] доказал, что для квадратичных иррациональностей а разложение в цепную дробь периодично (и обратно), т.е. последовательность неполных частных ао, ai, а2, аз,..., начиная с какого-то номера, состоит из повторяющегося отрезка а:, а:+!,..., а:+4. Итак, алгоритм разложения числа в цепную дробь:
1. прост;
2. дает наилучшие рациональные приближения к числу;
3. конечен для рационального числа;
4. периодичен для квадратичных иррациональностей [2];
5. устроен как для почти всех чисел [2] для кубических иррациональностей [37]. Кроме того, он обладает еще рядом замечательных свойств.
1.2. История обобщений
В 1775 г. Эйлер [6] сделал первую попытку обобщить алгоритм цепной дроби на векторы. Впоследствии его подход развивали Якоби [7, 8], Пуанкаре [9], Брун [10], Гютинг [67], Брюно [16, 17, 19, 21] (два алгоритма) и Парусников [16, 17], Пустыль-ников [13] и др. Они по аналогии с (2) строили матричные алгоритмы вида
Лк+1 = Ок Лк, к = 0,1,2,..., (4)
где Лк - п-мерный вектор и С к - квадратная п-матрица с целыми элементами, которая образована по вектору Лк, и ёе! Ск = ±1. Эти алгоритмы просты, но, вообще говоря, не дают наилучших рациональных приближений к вектору и не всегда при п = 3 обладают аналогом свойства 4.:
4'. периодичность для кубических иррациональностей.
Еще Эрмит [15] критиковал алгоритм Якоби. В работах [16, 17, 18, 19, 20, 21, 22, 23] было проведено сравнение качества матричных алгоритмов и было установлено, что ни один их них не обладает свойствами 2 и 4' для всех векторов Ло. При этом оказалось, что наихудшим является алгоритм Пуанкаре [9]. В 1842 г. Дирихле [24] при п = 3 предложил рассматривать не трехмерный вектор Л = Ло, а две линейные однородные формы 11(Х) и 12(Х) такие, что ^(Л) = 12(Л) = 0. В 1850 г. Эрмит [25], развивая этот подход, предложил свое обобщение цепной дроби. Наконец, в 189596 гг. Клейн [26], Минковский [27] и Вороной [28] независимо пришли к тому, что надо в М3 рассматривать тройку однородных линейных форм 11(Х),12(X),13(X) и предложили свои концепции обобщения цепной дроби. При этом Клейн ограничился общими геометрическими соображениями, а Минковский и Вороной предложили конкретные алгоритмы. Впоследствии подход Клейна переоткрывали Скубенко [29] и Арнольд [30] и называли многогранники Клейна парусами и многогранниками Арнольда [31, 32, 74] соответственно. И хотя в [16] был предложен алгоритм вычисления многогранников Клейна, в [17, 18, 19, 20, 21, 22, 23] было выяснено, что многогранники Клейна-Скубенко-Арнольда не дают основы для хорошего алгоритма, обобщающего цепную дробь. Только алгоритмы Минковского и Вороного обладают свойствами 2 и 4' но они весьма громоздки. Их применению и развитию было посвящено много работ (см. [33, 34, 35]). Свой подход к обобщению цепной дроби предложил Гурвиц [36] в 1894 г., но без алгоритма. Интерес автора к обобщению цепной дроби возник в связи с его работой [37] 1964 г., которая была повторена Ленгом [38] в 1972 г. (см. также Старк [39]). В 2003 г. автор [40] предложил новое обобщение цепной дроби. Оно состоит в следующем.
Пусть для X = (х1,х2,х3) € М3 заданы три однородные линейные формы /¿(X), г = 1, 2, 3. Тогда в М3 имеются три плоскости С% = {X : /¿(X) = 0} и три прямых
Сг П Сj. Модули форм /¿(X)| = ) задают отображение целочисленной решетки
й3 С К3 в некоторое множество Z неотрицательного октанта Б = М+ = {У : У ^ 0}. Пусть Н - выпуклая оболочка множества Z и дН - ее граница. Множество дН является двумерной многогранной поверхностью, состоящей из вершин, ребер и граней. При этом плоскостям Сг С К3 соответствуют координатные квадранты октанта Б, а прямым Сг П С^ С М3 соответствуют координатные лучи октанта Б. Наилучшим целочисленным приближением плоскостей С% и прямых С% П С в М3 соответствуют вершины многогранной поверхности дН в Б. Поэтому нахождение наилуч-
ших диофантовых приближений сводится к вычислению поверхности дН [73], которая была названа «модульным многогранником». Модульный многогранник является выпуклой оболочкой четырех многогранников Клейна. Для одиннадцати наборов форм /¿(X),г = 1, 2, 3 модульные многогранники вычислены в [44] и показаны их логарифмические проекции. Свойства модульных многогранников изучены в [42, 43, 44, 45]. Оказалось, что все грани модульного многогранника являются треугольниками с вершинами С(В),г = 1,2,3, где вектор С(Х) = ($1(Х), д2(Х), $з(Х)). При этом ёе^В^гВз) = 0 или ±1 или ±2. Если коэффициенты форм /¿(X) принадлежат одному кубическому полю, то поверхность дН дважды периодична. В [43, 69] был предложен алгоритм движения по поверхности дН. Однако он был довольно громоздким.
1.3. Постановка задачи
Рассмотрим более общую ситуацию.
Пусть в п-мерном вещественном пространстве Мп с координатами X = (ж1,..., жп) заданы т однородных вещественных форм /1(Х),... ,/т(Х), которые не имеют общего корня X = 0, 2 ^ т ^ п. Целочисленная ненулевая точка X € Zn\0 С Мп называется минимальной, если нет такой другой ненулевой целочисленной точки У € Zn\0 € Мп, У = —X, выполнены неравенства
|/г(У)| ^ /¿^)|, г = 1,..., т.
Задача 1. Найти все минимальные точки.
Частичное 'решение задачи 1. Модули форм /¿(X) задают отображение С^) пространства Мп в положительный ортант Б = М^" т-мерного пространства Мт с координатами Б = (в1,...,вт): вг = ¿'¿(X) = /¿^)|, г = 1,...,т. При этом целочисленная решетка Zn С Мп отображается в некоторое множество Z С Б. Замыкание выпуклой оболочки Н множества Z\0 является выпуклым множеством. Все целочисленные точки X € ^п\0, отображающиеся на границу дН множества Н, являются минимальными и дают наилучшие диофантовые приближения к совокупности корневых пространств т форм /¿(X). Возможны минимальные точки X, образы которых С^) = ),...,¿т^)) не лежат на границе дН. Однако будем искать погра-
ничные минимальные точки X, для которых С^) € дН.
Пример 1. Если т = п = 2, а /1^) и /2^) — линейные формы, то дН -выпуклая ломаная на плоскости М2, лежащая в Б = М+. Координаты х1 и х2 прообраза X = (ж1,ж2) ее вершины С = (¿1 , ¿2) - это числитель и знаменатель подходящей цепной дроби числа, равного котангенсу наклона одной из корневых прямых /1^) = 0 или /2^) = 0. Точки из Z7 лежащие на ребрах ломанной дН, отсутствуют. Промежуточным дробям этой цепной дроби соответствуют не минимальные пограничные точки.
В дальнейшем ограничимся случаями, когда выпуклое множество Н является многогранным, т.е. его граница дН состоит из вершин, ребер и граней различных размерностей и не содержит непрерывных «кривых» частей. В этих случаях граница дН вычисляется с помощью стандартных программ для вычисления выпуклых многогранников [68, 73].
Пусть иррациональные Л1 и Л2 € М - корни квадратного уравнения Л2 + аЛ + Ь = 0 с целыми коэффициентами. Тогда цепные дроби чисел Л1 и Л2 периодичны и ломаная дН для пары форм /1(Х) = х1 — Л1х2 и /2(Х) = х1 — Л2х2 также периодична.
Задача 2. Найти все периоды поверхности дН, т.е. ее линейные автоморфизмы X = ТУ.
1.4. О содержании статьи
Статья состоит из трех параграфов.
В § 2 рассматривается случай трех линейных форм от трех переменных. Для этого случая определяется модульный многогранник, изучаются его свойства. Подробно рассматривается подслучай, когда коэффициенты всех форм принадлежат одному кубическому полю. В нем модульный многогранник имеет два независимых периода, которым соответствуют две фундаментальные единицы поля. Наконец, приводится 11 примеров проекций модульных многогранников на плоскость.
В § 3 обсуждается случай нескольких линейных и нескольких квадратичных форм. Определяется модульный многогранник. В алгебраическом случае он имеет столько периодов, сколько фундаментальных единиц имеет соответствующее поле. Через периоды можно вычислить эти фундаментальные единицы.
2. Модульный многогранник и его глобальные свойства
2.1. Общие свойства [69]
Пусть в М3 заданы три вещественные независимые однородные линейные формы
/г(Х) = {Ьг,Х), г = 1, 2, 3, ёе1(Ь1,Ь2,^з) = 0, (5)
где X € М3, Ьг = (1г1,1г2,1г3) € М3, {■, ■) означает скалярное произведение, и пространство М3 двойственно пространству М3. Три линейные однородные формы (5) определяют отображение С(Х) = (#1(Х),д2(X),д3(X)), где дг(Х) = |/г(Х)|, г = 1,2,3.
В [43] предложена следующая конструкция. Вектор-функция С(Х) отображает пространство М3 с координатами X в пространство М3 с координатами О = (д1,д2,д33), точнее - в его неотрицательный октант М+ = Б. При этом множество Z3\0 всех целочисленных точек X кроме X = 0 отображается в некоторое множество Z в Б. Пусть Н - замыкание выпуклой оболочки точек множества Z и дН - его граница. Очевидно, что Н - это выпуклый многогранник (точнее - многогранное множество), который назван модульным, и дН - это выпуклая многогранная поверхность, состоящая из вершин Уг, ребер Кг и граней Гг. Пусть У1, У2, У3 € Z и
У- = О(Б,), В € 1 = 1, 2, 3. (6)
Положим ш(У1,У2,У3) = | ¿е1(В1В2В3)|. Очевидно ш принимает целые неотрицательные значения. Для грани Гг поверхности дН определим ш(Гг) как минимум
У2, У3) по всем тройкам точек У1, У2, У3 € Z, лежащих на грани Г и не лежащих на одной прямой, т. е. У1, У2, Уз € Z П Г и ёе^У1У2У3) = 0. Грань Г поверхности дН назовем простой, если она является треугольником с вершинами (6) и не содержит других точек из множества Z, и полупростой, если она является треугольником, содержащим внутри ровно одну точку множества Z, и имеет ш(Г) = 1. Для простой грани Г с вершинами (6) имеем ш(Г) = | ёе^В^2В3)|.
Теорема 1. Для граней Г поверхности дН всегда ш(Г) ^ 2.
Теорема 2. Грань Г с ) = 2 является простым треугольником.
Теорема 3. Если точки (6) лежат на одной грани Г и ш(У1,У2,У3) = 0, то одна из точек В1,В2,В3 является суммой двух других.
2.2. Периодичность [75, §7]
Определение 0. Унимодулярная 3x3 матрица Т называется периодом поверхности дН, если поверхность дН, соответствующая формам /¿(У) = (Ь,УТ*), г = 1, 2, 3, переводится в дН линейным преобразованием
9 = №5, г = 1,2,3, (7)
где № € М и №1№2№3 = 1.
Здесь звездочка означает транспонирование. Таким образом, периодичность дН означает, что для дН два линейных преобразования X * = ТУ * в М3 и (7) в Б дают один и тот же результат.
Пусть у неприводимого в Q кубического многочлена
р(А) = А3 + аА2 + ЬА + с (8)
дискриминант положительный. Тогда у него имеются три вещественных корня А1 < < А2 < А3. Пусть вектор Ь в форме (5) такой, что
(Ь*,Ь*,Ь3 ) = ЯШ, (9)
где Я - невырожденная матрица, а Ш - матрица Вандермонда для полинома (8).
Теорема 4. Если все коэффициенты многочлена (8) целые числа, и все элементы матрицы Я рациональные, то поверхность дН, соответствующая формам (7), двояко периодическая, т. е. ее логарифмическая проекция
П1 = ^ 1/1 (X)|, П2 = 1о§ |/2(X)| (10)
инвариантна относительно двух независимых параллельных переносов.
Положим, что модуль О состоит из чисел (Ь,Х), где X € Ъ3, и числа Г1,Г2,Г3 являются генераторами модуля О. В соответствии с теоремой Дирихле [66, гл. II] в этом случае имеется ровно две фундаментальные единицы и и V поля 0>(А) в модуле О, т.е. каждая единица поля имеет форму (и^п, где £к = 1 и к, 1,т € Ъ. Следовательно, поверхность дН имеет два независимых периода £>** и £>*. В логарифмической проекции каждому периоду соответствует свой параллельный перенос.
2.3. Примеры [44]
Здесь для т = п = 3 приводятся вычисления многогранных поверхностей дН для 11 тернарных кубических форм вида
) = (¿1, X)(£з,Х) , (11)
где fi(X) = ), г = 1, 2, 3. Десять из них, являющиеся экстремальными фор-
мами Девенпорта и Свиннертона-Дайера [62], рассматривались в [16] - [23], где для них были вычислены многогранники Клейна. Одиннадцатая форма (11) была взята из примера работы Вороного [28]. Результаты вычислений представлены в виде рисунков, где показаны логарифмические проекции (10) многогранников дН на плоскость п1 = 1п |Д(Х)|, п2 = 1п )|. Проекции вершин V(X) обозначены черными кружочками и соединены проекциями ребер. В каждой из этих точек V(X) = (|Д(Х)|, (X)|, (X)|) указано соответствующее значение вектора X и значение формы (11) для этого X.
Ниже каждая их этих 11 форм обсуждается в отдельном примере, и для нее приводятся начальные значения в виде многочлена (8) и таблицы
В1 Л(В1) f2(Bl) fз(Bl) Й1 £(В1) В2 f2(B2) fз(B2) Й2
Вз Л(Вз) f2(Bз) fз(Bз) аз £(Вз)
где В1, В2, Вз - исходный базис, А = (а1, а2, аз), Д(А) = ^(А) = 0, fi(Bj) точные значения. Во всех случаях поверхность дН просчитывается алгоритмом для вычисления выпуклого многогранника и имеет 2 независимых периода, которым в логарифмической проекции соответствуют параллельные переносы. Указываются фундаментальные области Ф относительно группы этих переносов.
Пример 2. Фо'рмаН1.
Согласно [16, 17], начальные данные суть р(А) = Лз + Л2 — 2А — 1,
(1,0,0) 11111 (0,1,0) -1.24697960 1.80193773 0.44504187 0.44504187 — А^2 (0,0,1) -0.55495813 -2.24697960 0.80193774 0.80193774 —1 — Ат
Логарифмическая проекция поверхности дН показана на рис. 1. Фундаментальная область Ф состоит из двух треугольников
Г1 = |С(1,0,0), С(0,1,0), С(0,0,1)}, Г2 = |С(1,0,0), С(0,0,1), С(1,0,1)}. (12) Пример 3. Форма Н2.
Согласно [16, 17], начальные данные суть р(А) = Аз — ЗА — 1,
Рис. 1: Логарифмическая проекция модульного многогранника для формы к1.
(1,0,0) 11111 (0,1,0) -1.87938524 1.53208889 0.34729636 0.34729636 —Лг+2 (0,0,1) -0.65270364 -2.87938524 0.53208889 0.53208889 —1 — Лг+1
Логарифмическая проекция поверхности дН показана на рис. 2. Фундаментальная область Ф состоит из двух треугольников (12).
Рис. 2: Логарифмическая проекция модульного многогранника для формы к2.
Пример 4. Форма Л,3.
Согласно [18], начальные данные суть р(Л) = Л3 — 3Л2 — Л + 1,
(1,0,0) 1 (0,1,0) 2.48119430
(0,0,1) -5.15632517
1
-1.17008647 -0.369102386
1
0.68889218 0.52542756
1 1 0.68889218
Лг
Л,2-3Л; + 1
0.52542756
Лг
Логарифмическая проекция поверхности дН показана на рис. 3. Фундаментальная область Ф состоит из двух треугольников
Г1 = {С(1,0,0), С(0,1,0), С(0,0,1)}, Г2 = {С(1,0,0), С(1,1,0), С(2,1,1)}. (13)
Рис. 3: Логарифмическая проекция модульного многогранника для формы к3.
Пример 5. Форма Ьд.
Согласно [20], ее начальные данные суть р(А) = Л3 + 9А2 + 6А — 1,
(1,0,0) 11111 (0,1,0) 0.87891770 -0.13776296 8.25884526 8.25884526 7 + 9А, + А2 (0,0,1) -7.25884526 0.121082302 1.13776296 1.13776296 А, + 1
Логарифмическая проекция поверхности дН показана на рис. 4. Фундаментальная область Ф состоит из двух треугольников
Г1 = |С(1, 0, 0), С(0,1, 0), С(0, 0,1)}, Г2 = |С(0,1, 0), С(0, 0,1), С(0,1,1)}. (14)
Рис. 4: Логарифмическая проекция модульного многогранника для формы Ьд.
Пример 6. Форма Л,4.
Согласно [19], ее начальные данные суть р(А) = А3 + 9А2 + 6А — 1,
(1,0,0) 1111 (0,1,0) 1.90687424 -2.92775456 -0.17911967 1.51261571 (0,0,1) -0.67598551 0.59666387 2.47932164 5.74622955
1
-3+15А;+2А? 5
11+10 А; +А?
Логарифмическая проекция поверхности дН показана на рис. 5. Фундаментальная область ограничена вершинами С(1,1, 4), С(0, 0,1), С(1, 0, 0), С(0, —1, 0), С(-5, —1, 2), С(3,1, —1), С(—1, 0,1), С(1,1, 3) и состоит из 9 треугольников.
Рис. 5: Логарифмическая проекция модульного многогранника для формы Н4.
Пример 7. Форма Н4, сопряженная к Н4.
Согласно [19], ее начальные данные суть р(А) = А3 + 9А2 + 6А — 1,
5
(1,0,0) 11111 (0,1,0) 0.26323621 -0.68494283 1.51261571 -0.17911967 13+26^+3Л'? (0,0,1) -0.38431852 -0.452820124 5.74622955 2.47932164 53+73^+8Л1
Логарифмическая проекция поверхности дН показана на рис. 6. Фундаментальная область поверхности дН состоит из 4 треугольников
Г1 = {С(2, —5,1), С(—1,1,0), С(—1, —3,1)}, Гг = {С(2, —5,1), С(—1,1,0), С(0,1,0)}, Г3 = {С(0,1,0), С(—1,1,0), С(—1, —3,1)}, Г 4 = {С(0,1,0), С(—1, —3,1), С(—1, —2,1)}.
Рис. 6: Логарифмическая проекция модульного многогранника для формы Н4.
Пример 8. Форма Ь5.
Согласно [21, 61], ее начальные данные суть р(А) = А3 + А2 — 4А + 1,
(1,0,0) 1 1 1 1 1 (0,1,0) 0.72610945 -0.37720285 3.65109341 3.65109341 —1 + 2Лг + Лг2 (0,0,1) -2.65109341 0.273890555 1.37720285 1.37720285 Лг
Логарифмическая проекция поверхности дН показана на рис. 7. Фундаментальная область Ф состоит из двух треугольников
Г1 = {О(1,0,0), 0(0,1,0), 0(0,0,1)}, Г2 = {0(0,1,0), 0(0,0,1), 0(0,1,1)}. (15)
Рис. 7: Логарифмическая проекция модульного многогранника для формы к5.
Пример 9. Форма Лб-
Согласно [22, 61], ее начальные данные суть р(Л) = Л3 + 4Л2 — 25Л — 1,
(1,0,0) 1 1 1 0.23382690 1
(0,1,0) 3.88530978 -0.24362329 1.35831350 0.12322591 -5+32Л;+2Л2 (0,0,1) -2.24362329 -0.64168650 1.88530978 0.31761024 -13+1221Л г+Л2
Логарифмическая проекция поверхности дН показана на рис. 8. Фундаментальная область состоит из 9 треугольников и ограничена вершинами С(5, —2, — 1), С(-3,1,1), С(2,0, —1), С(—1, —2,2), С(1, —2,1), С(0, —1,1), С(0,1,0), С(1,0,1), с(1,1,1), С(3,1,3), С(2,0,1), С(—2,1,1).
Рис. 8: Логарифмическая проекция модульного многогранника для формы к6.
Пример 10. Форма Л,7.
Согласно [23, 61], ее начальные данные суть р(А) = А3 — 2А2 — А + 1,
(1,0,0) 1 1 1 -0.08626792 1
(0,1,0) 0.643104132 0.307978528 5.04891734 0.21777981 А2 (0,0,1) 2.24697960 -0.801937736 0.55495813 -0.02393754 А2 — 2А,
Логарифмическая проекция поверхности дН показана на рис. 9. Фундаментальная область ограничена вершинами С(1, 0, 0), С(0, 0,1), С(—3,1, —3), С(1, 0,1), С(—1,1, —1), С(0,1, 0), С(—1, 2, 0), С(—1,1, 0) и состоит из 6 треугольников.
Рис. 9: Логарифмическая проекция модульного многогранника для формы к7.
Пример 11. Форма к7 сопряженная к к7.
Согласно [23], ее начальные данные суть р(Л) = Л3 — 2Л2 — Л + 1,
(1,0,0) 0.30141661 0.78485132 -0.08626792 1 (0,1,0) -0.09692113 -0.12085868 0.21777981 5.04891734 (0,0,1) 0.33863849 -0.31470094 -0.02393754 0.55495813
Здесь матрица {fi(Bj)} обратна матрице )} примера 10.
Логарифмическая проекция поверхности дН показана на рис. 10. Фундаментальная область Ф состоит из двух треугольников Г1 = {С(—1, 0,1), С(0, 0,1), С(0,1, 0)} и Г2 = {С(-1, 0,1), С(0,1, 0), С(1, 3, 0)}.
Рис. 10: Логарифмическая проекция модульного многогранника для формы Н7.
Пример 12. Форма Вороного.
В [28, § 53] рассмотрена форма (11), где Ьг = (1,Лг, (Лг + Л2)/2), г = 1,2,3, а Лг -корни уравнения Л3 — 7Л — 2 = 0. Начальные данные суть р(Л) = Л3 — 7Л — 2,
(1,0,0) 1 1 1 0.79184956 1
(0,1,0) -2.48928857 -0.28916855 2.77845712 1.95640956 Лг (0,0,1) 1.85363451 -0.10277505 5.24914054 2.20012002 (Лг + Лг2)/2
Логарифмическая проекция поверхности дН показана на рис. 11. Фундаментальная область Ф ограничена вершинами 0(1, 0, 0), 0(0,1, 0), 0(0, 0,1), 0(2, 5, 6), 0(1, 1, 2), 0(1, 1, 1) и состоит из 6 треугольников.
Рис. 11: Логарифмическая проекция модульного многогранника для формы Вороного.
Замечание 2. Формы Л1, Л2, Л3, Ьд и Л5 являются самосопряженными. Поэтому для них и их сопряженных форм Л многогранники дН и дН совпадают и устроены
просто. Формы Л,4,Л,6,Л,7 и форма Вороного не являются самосопряженными. Они отличаются от своих сопряженных форм; отличаются и их многогранники дН, дН, которые устроены сложнее. Здесь приведены многогранники для сопряженных форм Л,4 и Л,7, но отсутствуют многогранники для форм, сопряженных к Н6 и Вороного.
Замечание 2. Неуспех попыток найти единый матричный алгоритм вызван тем, что в сложных (несамосопряженных) случаях переходы от граней с ш = 1 к другим таким же граням дается последовательностью разных матричных алгоритмов. Только в самосопряженных случаях можно обойтись каким-то одним матричным алгоритмом.
3. Многомерный случай [70] 3.1. Модульный многогранник
Пусть в Мп = {X}, X = (х1, ..., жп) заданы I линейных форм
/¿(X) = X), г = 1, ...,1 (16)
и к квадратичных форм
/г+,-(X) = (К,, X)<К,, X), ; = 1, ..., к, (17)
I, к ^0, I + 2к = п, I + к = т ^2. (18)
Здесь - п-мерные вещественные векторы-строки, К, - п-мерные комплексные векторы-строки, X) означает скалярное произведение и черта сверху означает комплексное сопряжение. Очевидно, т ^ п. Предположим, что набор форм (16), (17) невырожден, т. е.
...,Ь*1,К*1,...,К*к,К*1,...,К*к) = 0. (19)
Здесь звездочка * означает транспонирование, т. е. Ь* и К* векторы-столбцы. Положим
дг(X) = /¿^)|, г = 1, ...,т, (20)
и
С^) = (gl(X), ..^д^)). (21)
Целочисленная точка X € называется минимальной или наилучшим (диофан-товым) приближением к корневым подпространствам
и = ^ /¿(X ) = 0} , г = 1, + к = т, X € Мп, (22)
совокупности форм (16), (17), если нет У € У = 0, У = —X:
С(У) < С(X). (23)
В частности, если п = 2,1 = 0 и к = 1, то этот случай не рассматривается, ибо т = к + 1 = 1 < 2, что противоречит (18).
Предлагается изучить структуру наилучших целочисленных приближений X € Ъп к корневым подпространствам (22) форм (16), (17) с помощью следующего подхода [69, 70].
Формулы (20), (21) задают отображение п-мерного пространства Мп = {X} в неотрицательный ортант Б = {5 ^ 0} т-мерного пространства Мт = {5}. При этом корневые подпространства (22) из Мп отображаются в координатные подпространства в Б, а целочисленные точки X = 0, т. е. X € Ъп \ {0}, отображаются в некоторое множество Z С Б. Пусть Н - замыкание выпуклой оболочки множества Z. Согласно определению 0 наилучшим целочисленным приближениям X € Ъп корневых подпространств (22) соответствуют точки ) множества Z, лежащие на границе дН многогранного множества Н.
Вершины, ребра и грани любых размерностей для дН вычисляются стандартной программой вычисления выпуклых оболочек [73] для т ^14 или [68] для любого т. Точки из Z П дН, отличные от вершин, также находятся. Они являются пограничными минимальными точками. Не пограничные минимальные точки этим подходом не отлавливаются.
3.2. Периодичность модульного многогранника
Рассмотрим линейное преобразование
X = ^ (24)
исходных координат, где В - неособая квадратная п-матрица. Напомним, что матрица В называется унимодулярной, если все её элементы - целые и ёе! В = ±1. Унимодулярная матрица В является периодом модульного многогранника Н, если линейному преобразованию (24) соответствует диагональное линейное преобразование
дг = №дг, г = 1, ...,т, = 1 (25)
между многогранниками Н и Н, где дг = дг ^XJ, или между их частями. Пусть неприводимый в Q многочлен
р(А) = Ап + ап-1Ап-1+ ... + Й1А + ао (26)
с целыми коэффициентами аг €Ъ имеет 1 вещественных корней А1, ..., Аг и к пар комплексно сопряжённых корней Аг+1, ..., Аг+к, Аг+ъ ..., Аг+к, 1 + 2к = п. Здесь 1 ^ 0, к ^ 0, но нам нужно чтобы к + 1 ^ 2. Пусть Ш — соответствующая матрица Вандермонда, К — неособая матрица с рациональными элементами и
(¿1, ¿2, КО = КШ, (27)
где звёздочка означает транспонирование, т. е. Ь\ — вектор-столбец. Рассмотрим т = к + 1 форм
¿(X) = (£г, X>, г = 1, ...,1, (28)
/г+,-(X) = (К, X>(К, X>, ; = 1, ...,к. ^
Теорема 5. Соответствующая формам (28) гиперповерхность дН размерности т — 1 имеет ровно т — 1 независимых периодов.
Согласно теореме Дирихле [66] поле 0>(Л) имеет ровно m — l независимых фундаментальных единиц. Следовательно, многогранная гиперповерхность дH размерности m — l имеет ровно m — l независимых периодов, произведение степеней которых не равно тождественному преобразованию. Алгоритмы, обобщающие цепную дробь, позволяют находить многогранники H. По многогранникам H находятся их периоды и по ним — фундаментальные единицы поля Q^) (см. [66]-[68], также [70, 71, 72, 7З]).
Теперь один из векторов Li, i = l, ..., l или Kj, j = l,..., k в (27) заменим произвольным вектором Li или Kj с сохранением вещественности или комплексности вектора Li или Kj и свойства невырожденности (19), и вместо формы fi(X) рассмотрим форму
f(X) = {Li, X), если i < l,
f(X) = {KU X){KU X), если i > l.
Пусть пересечение U корневых подпространств Lj всех нештрихованных форм (j = i) отлично от нуля, а пересечение подпространства U с корневым подпространством Li формы f'(X) состоит только из нуля X = О. Тогда поверхность д^, соответствующая m — l старым нештрихованным формам и одной штрихованной форме, будет m — l периодической вблизи подпространства U. Это обобщение теоремы Лагранжа о периодичности цепной дроби для квадратичной иррациональности.
4. Заключение
Наилучшие диофантовые приближения можно вычислять посредством глобального алгоритма, использующего стандартную программу вычисления выпуклых оболочек, вместо пошаговых вычислений как в алгоритме цепной дроби. Это дает решение проблемы, которую пыталось решить большинство ведущих математиков XIX века.
СПИСОК ЦИТИРОВАННОЙ ЛИТЕРАТУРЫ
1. Венков Б. А. Элементарная теория чисел. М.-Л.: ОНТИ, 19З7.
2. Хинчин А. Я. Цепные дроби. З-е изд. М.: Физматгиз, 1961.
3. Wallis J. A. Arithmetica infinitorum, 16ББ.
4. Euler L. De fractinibus continuis ^ Comm. Acad. Sci. Imper. Petropol. 17З7. V. 9.
Б. Lagrange J. L. Complement chez Elements d'algebre etc. par M. L. Euler. T. III. 1774.
6. Euler L. De relatione inter ternas pluresve quantitates instituenda УУ Petersburger Akademie Notiz. Exhib. August 14, 177Б У У Commentationes arithmeticae collectae. V. II. St. Petersburg. 1849. P. 99-104.
7. Jacobi C. G. J. Allgemeine Theorie der Kettenbruchänlichen Algorithmen, in welchen jede Zahl aus drei vorhergehenden gebildet wird ^ J. Reine Angew. Math. 1868. V. 69. P. 29-64. ^ Gesammelte Werke, Bd. IV. Berlin: Reimer, 1891. P. З8Б-426.
8. Jacobi C. G. J. Ueber die Auflösung der Gleichung ai^i + a2x2 + ... + anxn = // J. Reine Angew. Math. 1868. V. 69. P. 21-28.
9. Poincare H. Sur une generalization des fractiones continues // C.R. Acad. Sci. Paris. Ser. 1. 1884. V. 99. P. 1014-1016.
10. Brun V. En generalisation av Kjedebroken // Skrifter utgit av Videnskapsselskapeti Kristiania. I. Matematisk-Naturvidenskabelig Klasse 1919. 1920. № 6.
11. Perron O. Grundlagen för eine Theorie des Jacobischen Ketten-bruchalgorithmus // Math. Ann. 1907. V. 64. P. 1-76.
12. Bernstein L. The Jacobi-Perron algorithm — its theory and application. LNM 207. Berlin/Heidelberg/New York: Springer Verlag, 1971.
13. Пустыльников Л. Д. Обобщенные цепные дроби и эргодическая теория // УМН. 2003. Т. 58. № 1. С. 113-164.
14. Schweiger F. Multidimensional Continued Fractions. Oxford Univ. Press: New York, 2000.
15. Hermite Ch. Correspondance d'Hermite et de Stieltjes. T. II, lettres 232, 238, 408. Gauthier-Villars, Paris, 1905.
16. Брюно А. Д., Парусников В. И. Многогранники Клейна для двух кубических форм Давенпорта // Матем. заметки. 1994. Т. 56. № 4. С. 9-27.
17. Брюно А. Д., Парусников В. И. Сравнение разных обобщений цепных дробей // Матем. заметки. 1997. Т. 61. № 3. С. 339-348.
18. Parusnikov V. I. Klein polyhedra for complete decomposable forms // Number theory. Dvophantine, Computational and Algebraic Aspects. Editors: K. Gyo'ry, A. Petho' and V. T. Sos. De Gruyter. Berlin, New York. 1998. P. 453-463.
19. Парусников В. И. Многогранники Клейна для четвертой экстремальной кубической формы // Матем. заметки. 2000. Т. 67. № 1. С. 110-128.
20. Парусников В. И. Многогранники Клейна с большими гранями // Препринт № 93. М.: ИПМ им. М.В.Келдыша, 1997.
21. Парусников В. И. Многогранники Клейна для пятой экстремальной кубической формы // Препринт № 69. М.: ИПМ им. М.В.Келдыша, 1998.
22. Парусников В. И. Многогранники Клейна для шестой экстремальной кубической формы // Препринт № 69. М.: ИПМ им. М.В.Келдыша, 1999. 31 с.
23. Парусников В. И. Многогранники Клейна для седьмой экстремальной кубической формы // Препринт № 79. М.: ИПМ им. М.В.Келдыша, 1999. 32 с.
24. Lejeune Dirichlet G. P. Verallgemeinerung eines Satzes aus der Lehre von den Kettenbröchen nebst einigen Anwendungen auf die Theorie der Zahlen // S.- В. Press. Akad. Wiss. 1842. S. 93-95 // Werke. Bd. I. Berlin: Reimer, 1889. P. 635-638.
25. Hermite Ch. Lettres de M. Ch. Hermite â M. Jacobi sur différents objets de la theorie des nombres //J. Reine Angew. Math. 1850. Bd. 40. P. 261- 315 // Oeuvres, T. I, Paris: Gauther-Villares. 1905. P. 100-163 // Opuscule Mathematica de Jacobi. V. II.
26. Klein F. "Über eine geometrische Auffassung der gewöhnlichen Kettenbruchentwicklung // Nadir. Ges. Wiss. Göttingen Math.-Phys. Kl. 1895. № 3. P. 357-359.
27. Minkowski H. Generalisation de le theorie des fractions continues / Ann. Sci. Ec. Norm. Super, ser III, 1896, t. 13, p. 41-60. Also in: Gesamm. Abh. I. P. 278- 292.
28. Вороной Г. Ф. Об одном обобщении алгорифма непрерывных дробей. Варшава: Из-во Варш. Ун-та, 1896. Также: Собр. соч. в 3-х томах. Киев: Из-во АН УССР, 1952. Т. 1. С. 197-391.
29. Скубенко Б. Ф. Минимумы разложимых кубических форм от трех переменных // Записки научных семинаров Ленинградского отделения математич. ин-та им. Стеклова (ЛОМИ). 1988. Т. 168. С. 125-139.
30. Арнольд В. И. Цепные дроби. М.: МЦНМО, 2001.
31. Lachaud G. Polyedre d'Arnol'd et voile d'un câne simplicial: analogues du theoreme de Lagrange // C.R. Acad. Sei. Ser. 1. 1993. V. 317. P. 711-716.
32. Lachaud G. Polyedre d'Arnol'd et voile d'un câne simplicial, analogues du theoreme de Lagrange pour les irrationnels de degre quelconque. Pretirage N 93-17. Marseille: Laboratoire de Mathematiques Discretes du C.N.R.S., 1993.
33. Brentjes A. J. Multi-dimensional Continued Fraction Algorithms. Mathematical Centre Tracts 145, Amsterdam: Mathematisch Centrum, 1981.
34. Buchmann J. On the period length of the generalized Lagrange algorithm //J. Number Theory. 1987. V. 26. P. 8-37.
35. Быковский B. A. Теорема Валена для двумерных подходящих дробей // Матем. заметки. 1999. Т. 66. № 1. С. 30-37.
36. Hurwitz A. Üeber die angenaherte Darstellung der Zahlen durch rationale Bröche // Math. Ann. 1894. Bd. 44. P. 417-436.
37. Брюно А. Д. Разложения алгебраических чисел в цепные дроби // Журнал вычислительной матем. и матем. физики. 1964. Т. 4. № 2. С. 211-221.
38. Lang S., Trotter Н. Continued fractions of some algebraic numbers //J. Reine Angew. Math. 1972. Bd. 252. P. 112-134.
39. Старк X. M. Объяснение некоторых экзотических непрерывных дробей, найденных Бриллхартом // Вычисления в алгебре и теории чисел (ред. Б.Б. Венков и Д.К. Фаддеев). М.: Мир, 1976. С. 155-171.
40. Minkowski Н. Über die Annelierung an eine reele Grosse durch rationale Zahlen // Math. Annalen. 1901. B. 54. P. 91-124. Also in: Gesamm. Abh. I. P. 320-352.
41. Pipping N. Zur Theorie der Diagonalkettenbrüche // Acta Acad. Aboens., 1924. B. 3. 22 S.
42. Нечаев В. И. Диагональная цепная дробь // Математ. Энциклоп. М.: Советская Энциклоп. 1979. Т. 2. С. 123.
43. Брюно А. Д. Правильное обобщение цепной дроби // Препринт № 86. М.: ИПМ им. М.В.Келдыша, 2003. 17 с.
44. Брюно А. Д., Парусников В. И. Многогранники модулей троек линейных форм // Препринт N 93. М.: ИПМ им. М.В.Келдыша, 2003. 20 с.
45. Брюно А. Д. К обобщениям цепной дроби // Препринт N 10. М.: ИПМ им. М.В.Келдыша, 2004. 31 с.
46. Брюно А. Д. Алгоритм обобщенной цепной дроби // Препринт N 45. М.: ИПМ им. М.В.Келдыша, 2004. 32 с.
47. Брюно А. Д., Парусников В. И. Дальнейшее обобщение цепной дроби // Препринт N 40. М.: ИПМ им. М.В.Келдыша, 2005. 19 с.
48. Bruno A.D. and Parusnikov V. I. New generalizations of the continued fraction // Preprint no. 52 of the Keldysh Inst, of Applied Math.: Moscow, 2005. 20 p.
49. Брюно А. Д. Структура наилучших диофантовых приближений // ДАН. 2005. Т. 402. № 4. С. 439-444.
50. Брюно А. Д. Алгоритм обобщенной цепной дроби // ДАН. 2000. Т. 402. № 6. С. 732-736.
51. Parusnikov V. I. Comparison of several generalizations of the continued fraction // Чебышевский сборник (Тула). 2005. T. 5. № 4. P. 180-188.
52. Брюно А. Д. Свойства модульного многогранника // Препринт № 72. М.: ИПМ им. М.В.Келдыша, 2005. 31 с.
53. Klein F. Sur une representation geometrique du developpement, en fraction continue ordinare // Nouv. Ann. Math. (3), 1896, Bd. 15. P. 327-331.
54. Klein F. Ausgewahlte Kapitel der Zahlentheorie. Bd. I, Einleitung. Güttingen. 1896. P. 16-50.
55. Koksma J. F. Diophantische Approximationen. Berlin: Julius Springer, 1936.
56. Perron O. Die Lehre von den Kettenbrüchen. Teubner, Leipzig, 1913; Stuttgart, 1954, 1977.
57. Hurwitz A. Ueber eine besondere Art der Kettenbruch-Entwiklung relier Grossen // Acta math. 1889. В. 12. P. 367-405.
58. Коркина E. Двумерные цепные дроби. Самые простые примеры. // Труды Мат. ин-та им. Стеклова. 1995. Т. 203. С. 143-166.
59. Kontsevich M. L., Suhov Yu. M. Statistics of Klein polyhedra and multidimensional continued fractions // Amer. Math. Soc. Transi. (2). 1999. V. 197. P. 9-27.
60. Briggs K. Klein polyhedra (2013), Доступно по адресу: http://keithbriggs.info/klein-polyhedra.html.
61. Парусников В. И. Многогранники Клейна для трех экстремальных форм // Ма-тем. заметки. 2005. Т. 77. Na 4. С. 566-583.
62. Swinnerton-Dyer H. P. F. On the product of three homogeneous linear forms // Acta Arithmetica. 1971. V. 18. P. 371-385.
63. Делоне Б. Н., Фаддеев Д. К. Теория иррациональностей третьей степени. Труды МИ АН. Т. 11. М.-Л.: АН СССР, 1947.
64. Делоне Б. Н. Петербургская школа теории чисел. М.-Л.: АН СССР, 1947.
65. Касселс Дж. В. С. Введение в теорию диофантовых приближений. М.: Изд- во иностр. лит., пер. с англ., 1961.
66. Боревич З. И., Шафаревич И. Р. Теория чисел. Второе изд. М.: Наука, 1972.
67. Guting R. Zur Verallgemeinerung des Kettenbruchalgorithmus. I // J. Reine Angew. Math., 278/279 (1975).
68. Fukuda K. Exact algorithms and software in optimization and polyhedral computation // Proceed. ISSAC'08 of XXI International Symposium on Symbolic and algebraic computations, ACM NY, USA, 2008, 333-334.
69. Брюно А. Д. Обобщения цепной дроби // Чебышевский сборник, 7:3 (2006) 4-71.
70. Брюно А. Д. Структура многомерных диофантовых приближений // ДАН 433:5 (2010) 587-589.
71. Брюно А. Д., Парусников В. И. Двустороннее обобщение непрерывных дробей // ДАН 429:6 (2009) 727-730.
72. Брюно А. Д. Структура наилучших диофантовых приближений и многогмерное обобщение цепной дроби // Чебышевский сборник (Тула) 11:1 (2010) 68-73.
73. Maple 2015.0 // http://www.maplesoft.com/products/Maple.
74. Bruno A. D. On geometric methods in works by V. I. Arnold and V. V. Kozlov. Preprint of arXiv, No 1401.6320.
75. Bruno A. D. New generalization of continued fraction, I // Functiones et Approximate. 2010. vol. 43, no. 1. P. 55-104.
REFERENCES
1. Venkov, BA 1937, Elementary theory of numbers, ONTI, Moscow-Leningrad.
2. Khinchin, AYa 1963, Continued fractions, Noordhoff, Groningen.
3. Wallis, JA 1655, Arithmetica infinitorum.
4. Euler, L. 1737, "De fractinibus continuis", Comm. Acad. Sci. Imper. Petropol., vol. 9.
5. Lagrange, JL 1774, Complement chez Elements d'algebre etc. par M. L. Euler, T. III.
6. Euler, L. 1775, "De relatione inter ternas pluresve quantitates instituenda", Petersburger Akademie Notiz. Exhib., Commentationes arithmeticae collectae. V. II. St. Petersburg. 1849. pp. 99-104.
7. Jacobi, C. G. J. 1868, "Allgemeine Theorie der Kettenbruchänlichen Algorithmen, in welchen jede Zahl aus drei vorhergehenden gebildet wird", J. Reine Angew. Math., vol. 69, pp. 29-64. // Gesammelte Werke, Bd. IV. Berlin: Reimer, 1891. pp. 385-426.
8. Jacobi, C. G. J. 1868, "Ueber die Auflösung der Gleichung aixi+a2x2+.. .+anxn = fn", J. Reine Angew. Math., vol. 69, pp. 21-28.
9. Poincare, H. 1889, "Sur une generalization des fractiones continues", C.R. Acad. Sci. Paris, Ser. 1, vol. 99, pp. 1014-1016.
10. Brun, V. 1920, "En generalisation av Kjedebroken", Skrifter utgit av Videnskap-sselskapeti Kristiania. I, Matematisk-Naturvidenskabelig Klasse 1919, no. 6.
11. Perron, O. 1907, "Grundlagen för eine Theorie des Jacobischen Ketten- bruchalgorithmus", Math. Ann., vol. 64, pp. 1-76.
12. Bernstein, L 1971, The Jacobi-Perron algorithm — its theory and application, LNM 207, Springer Verlag, Berlin/Heidelberg/New York.
13. Pustylnikov, L.D. 2003, "Generalized continued fractions and the ergodic theory", Russian Math.-Surveys, vol. 58, no. 1.
14. Schweiger, F 2000, Multidimensional Continued Fractions, Oxford Univ. Press, New York.
15. Hermite, Ch. 1905, Correspondance d'Hermite et de Stieltjes, T. II, lettres 232, 238, 408, Gauthier-Villars, Paris.
16. Bruno, A.D. & Parusnikov, V.I. 1994, "Klein polyhedrals for two cubic Davenport forms", Math. Notes, vol. 56, no. 3-4, pp. 994-1007.
17. Bruno, A.D. & Parusnikov, V.I. 1997, "Comparison of various generalization of continued fractions", Math. Notes, vol. 61, no. 3, pp. 278-286.
18. Parusnikov, V.I. 1998, "Klein polyhedra for complete decomposable forms", Number theory. Dvophantine, Computational and Algebraic Aspects, Editors: K. Gyo'ry, A. Petho' and V. T. Sos, De Gruyter, Berlin, New York, pp. 453-463.
19. Parusnikov, V.I. 2000, "Klein polyhedra for the fourth extremal forms", Math. Notes, vol. 67, no. 1, pp. 87-102.
20. Parusnikov, V.I. 1997, "Klein's polyhedra with big faces", Preprint no. 93 of the Keldysh Inst. of Applied Math., Moscow.
21. Parusnikov, V.I. 1998, "Klein's polyhedra for the fifth extremal cubic form", Preprint no. 69 of the Keldysh Inst. of Applied Math., Moscow.
22. Parusnikov, V.I. 1999, "Klein's polyhedra for the sixth extremal cubic form", Preprint no. 69 of the Keldysh Inst. of Applied Math., Moscow.
23. Parusnikov, V.I. 1999, "Klein's polyhedra for the seventh extremal cubic form", Preprint no. 79 of the Keldysh Inst. of Applied Math., Moscow.
24. Lejeune Dirichlet, G.P. 1842, "Verallgemeinerung eines Satzes aus der Lehre von den Kettenbrüchen nebst einigen Anwendungen auf die Theorie der Zahlen", S.- B. Press. Akad. Wiss., S. 93-95 // Werke. Bd. I, Reimer, Berlin, 1889, pp. 635-638.
25. Hermite, Ch. 1850, "Lettres de M. Ch. Hermite â M. Jacobi sur differents objets de la theorie des nombres", J. Reine Angew. Math., Bd. 40, pp. 261-315 // Oeuvres, T. I, Gauther-Villares, Paris, 1905, pp. 100-163 // Opuscule Mathematica de Jacobi. V. II.
26. Klein, F. 1895, "Über eine geometrische Auffassung der gewöhnlichen Kettenbru-chentwicklung", Nadir. Ges. Wiss. Göttingen Math.-Phys. Kl., no. 3, pp. 357-359.
27. Minkowski, H. 1896, "Generalisation de le theorie des fractions continues", Ann. Sci. Ec. Norm. Super, ser III, t. 13, pp. 41-60. Also in: Gesamm. Abh. I. P. 278-292.
28. Voronoi, GF 1896, On Generalization of the Algorithm of Continued Fraction, Warsawa Üniversity.
29. Skubenko, B.F. 1991, "Minimum of decomposable cubic form of three variables", J. Sov. Math., vol. 53, no. 3, pp. 302-321.
30. Arnold, V.I. 1998, "Higher dimensional continued fraction", Regular and Chaotic Dynamics, vol. 3, no. 3, pp. 10-17.
31. Lachaud, G. 1993, "Polyedre d'Arnol'd et voile d'un cone simplicial: analogues du théorème de Lagrange", C.R. Acad. Sei., Ser. 1, vol. 317, pp. 711-716.
32. Lachaud, G. 1993, "Polyedre d'Arnol'd et voile d'un cone simplicial, analogues du theoreme de Lagrange pour les irrationnels de degre quelconque", Pretirage N 93-17, Laboratoire de Mathematiques Discretes du C.N.R.S., Marseille.
33. Brentjes, A.J. 1981, "Multi-dimensional Continued Fraction Algorithms", Mathematical Centre Tracts, 145, Amsterdam, Mathematisch Centrum.
34. Buchmann, J. 1987, "On the period length of the generalized Lagrange algorithm", J. Number Theory, vol. 26, pp. 8-37.
35. Bykovsky, V.A. 1999, "The Valen's theorem for two-dimensional convergent fraction", Math. Notes, vol. 66, no. 1.
36. Hurwitz, A. 1984, "Ueber die angenäherte Darstellung der Zahlen durch rationale Brüche", Math. Ann., Bd. 44, pp. 417-436.
37. Bruno, A.D. 1964, "Continued fraction expansion of algebraic numbers", USSR Comput. Math. and Math. Phys., vol. 4, no. 2, pp. 1-15.
38. Lang, S. & Trotter, H. 1972, "Continued fractions of some algebraic numbers", J. Reine Angew. Math., Bd. 252, pp. 112-134.
39. Stark, H.M. 1971, "An explanation of some exotic continued fractions found by Brillhart", Computers in Number Theory, Academic Press, London and New York, pp. 21-35.
40. Minkowski, H. 1901, "Uber die Annelierung an eine reele Grosse durch rationale Zahlen", Math. Anmalen. B. 54, pp. 91-124. Also in: Gesamm. Abh. I. pp. 320-352.
41. Pipping, N. 1924, "Zur Theorie der Diagonalkettenbrüche", Acta Acad. Aboens., B. 3, 22 S.
42. Nechaev, V.I. 1979, "Diagonal continued fraction", Math. Encyclop., Kluwer Acad. Publ.
43. Bruno, A.D. 2003, "The correct generalization of the continued fraction", Preprint no. 86 of the Keldysh Inst. of Applied Math., Moscow.
44. Bruno, A.D. & Parusnikov, V.I. 2003, "Polyhedra of absolute values for triple of linear forms", Preprint no. 93 of the Keldysh Inst. of Applied Math., Moscow.
45. Bruno, A.D. 2004, "On generalization of the continued fraction", Preprint no. 10 of the Keldysh Inst. of Applied Math., Moscow.
46. Bruno, A.D. 2004, "Algorithm of generalized continued fractions", Preprint no. 45 of the Keldysh Inst. of Applied Math., Moscow.
47. Bruno, A.D. & Parusnikov, V.I. 2005, "Further generalization of the continued fraction", Preprint no. 40 of the Keldysh Inst. of Applied Math., Moscow.
48. Bruno, A.D. & Parusnikov, V. I. "New generalizations of the continued fraction" Preprint no. 52 of the Keldysh Inst. of Applied Math., Moscow, 20 p.
49. Bruno, A.D. 2005, "Structure of the best Diophantine approximations", Doklady Mathematics, vol. 71, no. 3, pp. 396-400.
50. Bruno, A.D. 2000, "Generalized continued fraction algorithm", Doklady Mathematics, vol. 71, no. 3, pp. 446-450.
51. Parusnikov, V. I. 2005, "Comparison of several generalizations of the continued fraction", Chebyshevsky Sbornik (Tula), vol. 5. no. 4, pp. 180-188.
52. Bruno, A.D. 2005, "Properties of the modular polyhedron", Preprint no. 72 of the Keldysh Inst. of Applied Math., Moscow.
53. Klein, F. 1896, "Sur une representation geometrique du developpement, en fraction continue ordinäre", Nouv. Ann. Math. (3), Bd. 15, pp. 327-331.
54. Klein, F. 1896, "Ausgewählte Kapitel der Zahlentheorie", Bd. I, Einleitung. Göttingen. pp. 16-50.
55. Koksma, JF 1936, Diophantische Approximationen, Julius Springer, Berlin.
56. Perron, O 1913, Die Lehre von den KettenbrUchen, Teubner, Leipzig; Stuttgart, 1954, 1977.
57. Hurwitz, A. 1889, "Ueber eine besondere Art der Kettenbruch-Entwiklung relier Grossen", Acta math, B. 12, pp. 367-405.
58. Korkina, E.I. 1995, "Two-dimensional convergent fractions. The simplest examples", Proceedings of the Steklov Inst. of Math., vol. 209, pp. 124-144.
59. Kontsevich, M.L. & Suhov, Yu.M. 1999, "Statistics of Klein polyhedra and multidimensional continued fractions", Amer. Math. Soc. Transi. (2), vol. 197, pp. 9-27.
60. Briggs, K. Klein polyhedra (2013), Available at: http: //keithbriggs.info/klein-polyhedra.html.
61. Parusnikov, V.I. 2005, "Klein's polyhedra for three extremal forms", Math. Notes, vol. 77, no. 4, pp. 523-538.
62. Swinnerton-Dyer H.P.F. 1971, "On the product of three homogeneous linear forms", Acta Arithmetica, vol. 18, pp. 371-385.
63. Delone, B.N. & Faddeev, D.K. 1964, "The theory of irrationalities of the third degree", Am. Math. Soc. Transl. of Math. Monographs, vol. 10.
64. Delone, BN 1947, The Petersburg's School of Mathematics, AN SSSR, Moscow-Leningrad.
65. Cassels, JWS 1959, An Introduction to the Geometry of Numbers, Springer-Verlag, Berlin.
66. Borevich, ZI & Shafarevich, IR 1966, Number Theory, Academic Press.
67. Göting, R. 1975, "Zur Verallgemeinerung des Kettenbruchalgorithmus. I", J. Reine Angew. Math., 278/279.
68. Fukuda, K. "Exact algorithms and software in optimization and polyhedral computation", Proceed. ISSAC'08 of XXI International Symposium on Symbolic and algebraic computations, ACM NY, USA, 2008, pp. 333-334.
69. Bruno, A.D. 2006, "Generalization of continued fraction", Chebyshevsky sbornik, vol. 7, no. 3, pp. 4-71.
70. Bruno, A.D. 2010, "The structure of multidimensional Diophantine approximations", Doklady Mathematics, vol. 82, no. 1.
71. Bruno, A.D. & Parusnikov, V.I. 2009, "Two-way generalization of the continued fraction", Doklady Mathematics, vol. 80, no. 3, pp. 887-890.
72. Bruno, A.D. 2010, "Structure of the best Diophantine approximations and multidimensional generalizations of the continued fraction", Chebyshevskii Sbornik (Tula), vol. 11, no. 1. pp. 68-73.
73. Maple 2015.0 (2015), Available at: http://www.maplesoft.com/products/Maple.
74. Bruno, A.D. 2014, "On geometric methods in works by V. I. Arnold and V. V. Kozlov", Preprint of arXiv, No 1401.6320.
75. Bruno A. D. 2010 "New generalization of continued fraction, I", Functiones et Approximatio, vol. 43, no. 1. pp. 55-104.
Поступило 30.04.2015
