CC BY
30
ЧЕБЫШЕВСКИЙ СБОРНИК
Том 17 Выпуск 3
УДК 517.36
ОТ ДИОФАНТОВЫХ ПРИБЛИЖЕНИЙ ДО ДИОФАНТОВЫХ УРАВНЕНИЙ
А. Д. Брюно (г, Москва)
Аннотация
Пусть в вещественном n-мерном пространстве Rn = {X} задан о m однородных вещественных форм /¿(X), i = l,...,m, 2 ^ m ^ n. Выпуклая оболочка множества значений G(X) = (|/i(X)|,..., |/m(X)|) G М^дляцелочисленпых X G Zn во многих случаях является выпуклым многогранным множеством, граница которого для ||X|| < const вычисляется с помощью стандартной программы. Точки X G Zn, для которых значения G(X) лежат на этой границе, названы граничными. Они являются наилучшими диофантовыми приближениями для корневых множеств указанных форм. IIx вычисление даёт глобальное обобщение цепной дроби. Для n = 3 обобщить цепную дробь безуспешно пытались Эйлер, Якоби, Дирихле, Эрмит, Пуанкаре, Гурвиц, Клейн, Минковский, Врун, Арнольд и многие другие.
Пусть p(£) — целый неприводимый в Q многочлен степени иА — его корень. Набор основных единиц кольца Z[A] можно вычислить по граничным точкам некоторой совокупности линейных и квадратичных форм, построенных по корням многочлена p(£). До сих
n=2
n=3
граничных точек в Мп и автоморфизм их обр азов в М^. В логарифмической проекции М!^ на Мт-1 можно найти фундаментальную область для группы вторых автоморфизмов, соответствующих единицам.
С помощью этих конструкций можно находить целочисленные решения диофантовых уравнений специального вида. Аналогично вычисляются все указанные объекты для других колец поля Q(A). Приведены примеры.
Наш подход обобщает цепную дробь, позволяет вычислить наилучшие совместные приближения, основные единицы алгебраических колец поля Q(A) и все решения некоторого
n
Ключевые слова: обобщение цепной дроби, диофантовы приближения, набор основных единиц, фундаментальная область, диофантово уравнение.
Библиография: 16 названий.
FROM DIOPHANTINE APPROXIMATIONS TO DIOPHANTINE EQUATIONS
A, D, Bruno (Moscow) Abstract
Let in the real n-dimensional space Rn = {X} be given m red homogeneous forms /¿(X), i = 1,..., m, 2 < m < n. The convex hull of the set of points G(X) = (|/i(X)|,..., |/m(X)|) for integer X G Zn in many cases is a convex polyhedral set. Its boundary for ||X || < const can be computed by means of the standard program. The points X G Zn are called boundary-points if G(X) lay on the boundary. They correspond to the best Diophantine approximations X for the given forms. That gives the global generalization of the continued fraction. For n = 3
Euler, Jacobi, Dirichlet, Hermite, Poincare, Hurwitz, Klein, Minkowski, Brun, Arnold and a lot of others tried to generalize the continued fraction, but without a succes.
Let p(£) be an integer real irreducible in Q polynomial of the order n and A be its root. The set of fundamental units of the ring Z[A] can be computed using boundary points of some set of linear and quadratic forms, constructed by means of the roots of the polynomial p(£). Similary one can compute a set of fundamental units of other rings of the field Q(A). Up today such sets of fundamental units were computed only for n = 2 (using usual continued fractions) and n = 3 (using the Voronoi algorithms).
Our approach generalizes the continued fraction, gives the best rational simultaneous approximations, fundamental units of algebraic rings of the field Q(A) and all solutions of a
n
Keywords: generalization of continued fraction, Diophantine approximations, set of fundamental units, fundamental domain, Diophantine equation.
Bibliography: 16 titles.
1. Цепная дробь
Пусть ао и а\ — натуральные числа. Для нахождения их наибольшего общего делителя используется алгоритм Евклида последовательного деления с остатком:
а0 = а0а1 + а2, а1 = а^2 + а3, а2 = а2а3 + а4,...
где натуральные числа ао, а1, а2,... суть неполные частные. Это алгоритм разложения числа а = ао/а1 в правильную цепную дробь [1], и он применим к любым вещественным числам а. При этом а0 = [а], где [а] — целая часть чпела а, а1 = [1/(а — ао)], • • •, т.е.
а = ао +----, (1)
а1 +--
а2 +
afc+Л = (0 1 W ak \
ak+2j \1 -ak) \ak+i)
ak = [ak/ак+i].
Если разложение (1) оборвать на а к и свернуть эту оборванную цепную дробь в рациональное число рк/Цк, то получается подходящая дробь, которая даёт наилучшее рациональное
а
Рк Рк-Л = /рк-1 Рк-2\ (ак 1 Цк <7к-1/ VЦк-1 <?к-2/ \1 0
ак Ц-1 = (0 1 1, а*(Рк Рк-11 = ±1,
1 0) —ак) \Як Цк-1)
т.е. векторы (ак, ак+1) и (рк,Цк) принадлежат сопряжённым плоскостям, и пара векторов (Рк,Цк) (Рк-1,Цк-1) может служить базисом в одной из них. Лагранж [1, § 10] доказал, что
а
то есть последовательность неполных частных ао, а1, а2, аз,..., начиная с какого-то номера состоит из повторяющегося отрезка ак, ак+1,..., ак+г-
Итак, разложение числа в цепную дробь: просто; дает наилучшие рациональные приближения к числу; конечно для рационального числа; периодично для квадратичных иррацио-нальностей [1, § 10]; устроено как для почти всех чисел [1, гл. III] для кубических иррацио-нальностей [2]. Кроме того, оно обладает ещё рядом замечательных свойств.
2. Глобальное обобщение цепной дроби и наилучшие диофантовы приближения
Обобщить цепную дробь для векторов безуспешно пытались Эйлер, Якоби, Дирихле, Эр-мит, Пуанкаре, Гурвиц, Клейн, Минковский, Брун, Арнольд и многие другие [3,4], [5, п. 1.2]. Только пошаговые алгоритмы Вороного [6] безотказны, но сложны.
В [5,7,8] предложено следующее обобщение цепной дроби.
Пусть в п-мерном вещественном пространстве Мп с координатами Х = (ж1,..., жп) заданы т однородных вещественных форм (т.е. многочленов от переменных) Д(Х),..., /т(Х), 2 ^ ^ т ^ п.
Модули дг(Х) = |/г(Х )| форм /¿(Х), г = 1,..., т, задают отображение С(Х) = = (д1(Х),...,дт(Х)) пространства Мп в положительный ортант Б == М^ в т-мерном пространстве Мт с координатами Б = (в1,..., вт): вг = дг(Х) = |/г(Х)|, г = 1,..., т. При этом целочисленная решётка Ъп С Мп отображается в некоторое множество Z С Б. Замыкание выпуклой оболочки Н множества Z\0 является выпуклым множеством. Все целочисленные точки Х £ Zn\0, отображающиеся на границу дН множества Н, назовём граничными.
Х
Решение задачи 1. В дальнейшем ограничимся случаями, когда выпуклое множество H является многогранным, т.е. его граница dH состоит из вершин, рёбер, граней различных размерностей и не содержит непрерывных «кривых» частей. В этих случаях граница dH вычисляется с помощью стандартных программ для вычисления выпуклых многограных оболочек [9,10]. Это и даёт алгоритмическое обобщение цепной дроби на любую размерность. Примеры см.
в [5].
В частности, это даёт возможность вычислить наилучшие совместные рациональные приближения Qi/Qo, • • • j Qm/Qo к вещественным числам Д,..., ßm, где Qo, Qi,..., Qm G Z и /¿(q0, Qi) = Q0ßi — Qi, i = 1,..., m. Здесь m = m и n = m + 1.
Пример 1. Пусть /1 = xia — ж2; /2 = ж1; где a G R а > 0. Здесь n = m = 2. Каждой вершине ломаной dH с ж1 = p ж2 = Q G Z+ соответствует подходящая дробь q/p цепной дроби числа, а. Эта точка (ж1,ж2) является граничной. Но, вообще говоря, не каждой подходящей дроби s/r, r, s G Z+ цепной дроби числа а соответствует вершина ж1 = ^ x2 = s ломаной dH.
Гипотеза. Если все /i,..., /m суть линеиные и квадратичные формы, то граница dH не имеет непрерывных кривых участков, т. е. является многогранной.
Более того, до сих пор неизвестно ни одного набора форм /i,..., /m, для которого граница dH не была бы многогранной.
3. Основные единицы кольца Z[A]
Пусть дан целый неприводимый в Q вещественный многочлен
ко = г+ыга-1 +...+bn-ie + bn (2)
с целыми коэффициентами bi, т.е. он не разлагается в произведение двух нетривиальных многочленов с коэффициентами из Q. Ему соответствует кольцо Z[A] чисел вида
{(X) = ж1 + ж2 A + ... + жпАп-1
(3)
с целыми коэффициентами Хг, где Л — корень многочлена (2) и X = (х1,..., хп) € Каждому числу (3) соответствует квадратная матрица ^(£) = (с!^):
п— 1
Л^С(Х ) = ^ Л, г = 0, 1,...,п — 1.
3=0
Определитель det ^(£) называется нормой числа (3) и об означается N (£). Норма произведения чисел равна произведению их норм: N (£1 ■ £2) = N(£1) ■ N(£2)- Те числа (3), у которых норма N(£) = ±1, называются единицами [11, гл. II]. В дальнейшем предполагаем, что среди корней многочлена р(£) нет единиц. Существует такой набор единиц £ = (е1,...,ег), что всякая единица е € ^[Л] однозначно представляется в виде
е = ±еа1 ••• еагт, (4)
где аг — целые числа. Эти единицы е1,... ,ег называются основными.
Задача 2. Для фиксированного многочлена (2) найти набор основных единиц кольца Z[Л].
Решение задачи 2. Пусть неприводимый в ^ (2) имеет I вещественных корней Л1,
..., Лг и к пар комплексно сопряжённых корней Лг+1, .. ., Лг+&, Аг+ь .. ., Аг+&, I + 2к = п. Здесь I ^ 0 к ^ 0. Рассмотрим т = к +1 форм
/¿(X) = (Ьг,Х), г = 1,...,1, /г+з (X) = (К1+3, X) (Кг+,-, X>, э = 1,..., к,
где
Ьг = (1, Лг, Л2, . . . , ЛП—^ , (Ьг, X) = Х1 + ЛгХ2 + ... + ЛП-1Хп,
= (1, Лг+3, л21+з ,..., Лп—/) , К1+3 = (1, Аг+з, Аг+з,..., Ап+/) .
По теореме Дирихле [11, гл. II, § 4, п. 3] для многочлена (2) число основных единиц г = к+1 — 1. Далее предполагаем, что т = к + I ^ 2. Ибо, если к + I ^ 1, то г ^ 0 и по теореме Дирихле основные единицы отсутствуют.
Теорема 1 ( [11, гл. II, § 1, п. 2]). Для чисел (3) с X = (х1,... ,хп) € Мп
N(£) = /(X) = /1 (X) .../„(X). (5)
Поэтому для всех единиц вида (3)
/(X) = ±1п ) = |/(X)| = 1. (6)
Пусть °Жп — множество точек X € Ъп со свойством (6). Рассмотрим для него (т. е. для X €
конструкции раздела 1: множество Ъ значений
),...,^)) с б=мт,
гДе 9г(X) = |/¿(X)|, г = 1,..., т, выпуклую оболочку Н множества 2 и её границу дЙ. Граница дЙ имеет размерность т — 1 = г, не имеет кривых участков и состоит из вершин, рёбер и граней.
Теорема 2. Все грани границы, dH являются симплексами, а значение G0 = (1,1,..., 1) является её вершиной.
Пусть А — некоторая (m — ^^ь границы дН, содержащая вершину Go =
(1,1,..., 1), a Ri,..., Rra-i — её рёбра, содержащие G0.
Теорема 3. Пусть Gi — вторая вершина ребра Ri, отличная от вершины G0; i = 1,..., m — 1. Числа, (3), у которых G(X) = Gi; i = 1,..., m — 1; образуют набор основных единиц кольца, Z[A|.
Следовательно, для вычисления основных единиц надо на некотором ограниченном множестве ||X|| < const X g Zn вычислить кусок границы dH, содержащий (m — 1)-мерную грань А
Каждому числу (3) соответствует матрица
T (£) = xiE + X2 B + ... + xnBn-i, Где e — единичная, a B — это матрица, сопровождающая многочлен (2):
B=
/ 0 1 0 ••• 0 0 \
0 0 1 ••• 0 0
0 0 0 ••• 0 1
\—Ьп —bn-1 —bn-2 ••• —b2 — bij
Если число (3) является единицей, то матрица Т(£) унимодулярна и линейное преобразование X* = Т(£)^ в Мп является автоморфизмом множества Й и индуцирует автоморфизм
в** = д^)вг, г = 1,..., т, (7)
множества Й в Б = М™. Следовательно, каждой единице е соответствует период Т(е) обобщён-
т—1
Лагранжа [1, § 10], доказанной для п = I = 2, к = 0, т. е. т = к + I = 2.
4. Фундаментальная область
В поле О(Л) всякая целая степень Ь ^ п числа £ из (3) однозначно записывается в виде многочлена от Л степени п — 1, ибо
Лп = — (Ьп + Ьп- 1Л +... + Ь1Лп—1) .
Поэтому отношение двух многочленов от А однозначно записывается в виде многочлена от А п — 1
Теорема 4. Пусть X = (хь ... ,Хп), У = (ш,... ,Уп), 2 = (л,... ,2п) € О1 и £ (X )-£(У) = £(2), тогда /.(X) ■ /¿(У) = /¿(2), г = 1,..., т.
Следствие 1. В условиях теоремы 4 д»^) ■ дг(У) = $¿(2), г = 1,..., т.
Логарифмическая замена
Н^) = 1п gг(X), г = 1,..., т, Я = (Ль...,Лт)
взаимно однозначно переводит S = R™ в Rm. При этом (m — 1)-мерная граница dH многогранного множества H переходит в (m — 1)-мерную поверхность, которая взаимно однозначно
проектируется на Rm-i = {Я7}, где H' =f (hi,..., hm-i). На Rm-i автоморфизм (7) принимает вид
h* = lngi(X)+ hi, i = 1,..., m — 1, (8)
т. е. является параллельным переносом. Единицы кольца Z[A] образуют абелеву группу по умножению. По теореме 4 их логарифмы H образуют абелеву группу по сложению. В Rm-i имеется фундаментальная область F относительно сдвигов (8) этой группы. Пусть m-мерные векторы Gi, i = 1,...,m — 1, соответствующие основным единицам теоремы 3, имеют вид Gi = (gii, . . . ,gmi). Положим
ri = (ln gii,..., ln gm-i,i), i = 1,..., m — 1. (9)
Теорема 5. В Rm-i фундаментальная область относительно сдвигов (8), (9) — это (m — 1)
F = {Я' = ^iri + ... + ^m-iГт-i, 0 < ^ < 1, i = 1,..., m — 1} . (10)
При вычислении границы выпуклой оболочки некоторого множества точек трудности возрастают вместе с ростом количества точек. Чтобы уменьшить эти трудности можно вычисления разбить на следующие 6 шагов.
Шаг 1. Сначала в кольце Z[A] находим все единицы с X £ Zn из области ||X|| < const, вычисляя значения g(X) в этих X.
Шаг 2. Затем на множестве единиц {X} надо вычислить границу dH их выпуклой оболочки
H.
Шаг 3. По теореме 3 из dH выделяем набор основных единиц, представленных в S вершинами Gi,..., Gm-i.
Шаг 4. По теореме 5 находим фундаментальную область (10).
Шаг 5. Теперь выпуклая оболочка значений G(X) с £(X) £ Z[A] вычисляется только по тем X, у которых H'(X) попадают в фундаментальную область (10) и её близкую окрестность.
Шаг 6. По этой части границы dH восстанавливается вся граница dH с помощью периодов Gi i = 1, . . . , m — 1
5. Диофантовы уравнения
Многочлену p(£) степени n из (2) соответствует форма f (X) из (5) степени n от n переменных X = (xi,..., xn). Её коэффициенты являются многочленами от коэффициентов bi многочлена (2). Так, при n = 2
f (X) = xi — biXiX2 + ,
при n = 3
f (X) = xi — biX2X2 + b2Xix2 — Ьзх2 + (bi — 2b^ xix3 + (3Ьз — bib2)xiX2X3 —
— bib2X2X3 + (b2 — 2bib^ xix2 — b2b3X2X2 + b^.
Для любого п при хз = ... = хп = 0 имеем / (X) = хп — Ь1хп-1х2 + Ь2хп-2х2 — ...
... + (—1)п-1Ьп-1х1хп-1 + (—1)пЬпхп. Задача 3. Для заданного многочлена (2) найти все решения (3) с £ € ^[Л] уравнения
/ (X )= в, (И)
где число в рационально, в = ±1-Решение задачи 3.
Теорема 6 ( [11, гл. II, § 5, теорема 1]). Все решения (3) с £ € ^[Л] уравнения (11) имеют
вид
£ = £°е^ ...еаг, з = 1,...,7, (12)
где £0,..., £0 — конечное множество выделенных решений и а1,..., аг — любые целые числа. Если выделенных решений £0 нет, то уравнение (11) не имеет решений.
Далее даётся набросок конструктивного доказательства этой теоремы, позволяющий вычислять соответствующие постоянные и выделенные решения £0,..., £0.
Согласно предыдущим разделам находим набор основных единиц £ = (е1,..., ет-1) кольца Z[Л] и по ним строим фундаментальную область Т в координатах Я' = (Л-1,..., Лт-1).
Лемма 1. Для всех точек X € Мп; у которых логарифмические проекции
Я' = (Л1,...,Лт-1) лежат в фундаментальной области Т, справедливы, оценки
Цг ^ д^) ^ V», г = 1,..., т — 1, (13)
где 0 < < V — вещественные числа,.
Лемма 2. Для всех точек X, у которых Я' € Т и выполнено равенство д^) = |в справедливы оценки
Цт ^ дт^) ^ V™, (14)
где 0 < цт < V™ — вещественные числа,.
Нижние оценки в (13) нужны для получения верхней оценки в (14).
Поскольку (К, X) = (КК, X) + г (ЭК, X), то
(К, X) (К, X) = (КК, X)2 + (ЭК, X)2 .
Лемма 3. Если
7 = ёе^Ль..., Л,, КК,+1, ЭКг+ь..., КК^, ЭК^) = 0, (15)
мп
Поскольку
7 = (_1г)к det (Л1,..., Лг,Кг+ьКг+1,... ) ,
и последний определитель отличается от определителя Вандермонда Ш ненулевым множите-
п
лем, а Ш = Л (Аг _ Аj), то условие (15) эквивалентно условию, что у многочлена (2) нет
1<%<3
кратных корней. Но это так по условию, что многочлен (2) неприводим в
Неравенства (13), (14) выделяют в Мп всего 2т ограниченных областей вида
^г ^ Кг/г(Х) ^ Кг = ±1, г = 1, . . . , т.
В каждой из них количество целочисленных точек Х £ Zn конечно. В каждой из этих точек
/ ( Х ) Хг
среди этих точек Хг оставляем толь ко те Х0,..., Х0, для которых отно шения £(X,0)/£(Xj,) по 2 = г не лежат в Z[А]. Тогда = £(Х0), 2 = 1,..., являются выделенными решениями уравнения (11) и все решения этого уравнения имеют вид (12).
6. Обобщения
6.1. Единицы с положительной нормой
Для единицы е норма N(е) = ±1. Иногда нужны только единицы, у которых норма
п
£ = (е1,..., ет-1), надо в описанной процедуре раздела 3 оставлять только те точки Х £ Zn, для которых /(Х) = +1, и по ним указанным выше способом выделить мультипликативный базис. При нечётном п всякой единице е соответствует единица е' с N(е') = 1: это либо е, либо _е Х _Х
6.2. Произвольный порядок
Согласно [11, гл. II,§ 2] полный модуль в поле 0>(А), содержащий число 1 и являющийся кольцом, называется порядком поля 0>(А). Очевидно, что кольцо Z[А] является порядком поля 0>(А). Но в этом поле могут быть и другие порядки. Например, если в записи (3) все х2 — чётные, то получим подкольцо кольца Z[А]. Все результаты разделов 3-5, доказанные для порядка Z[А], справедливы для любого порядка О тол я 0>(А). Пусть ,..., шп — базис порядка О, т.е. все чиела а £ О имеют вид
а = + ^2^2 + ... + Уп^п, Уг £ Z. (16)
При записи этих чисел в виде (3) коэффициенты Хг могут быть рациональными числами.
О
могут иметь рациональные коэффициенты Хг- Существует такой набор единиц е1,... , ег £ О, что все единицы поля имеют вид (4). Назовём эти единицы основными. Для них справедливы все конструкции и теоремы разделов 3-5. Только матрица периода Т(е) может иметь рациональные элементы, но det Т(е) = N(е) = ±1. Поэтому для отыскания основных единиц
/(Х)
ными Хг. Дальнейшие вычисления такие же, как для кольца Z[А]. Мультипликативный базис единиц с положительной нормой образует набор квазиосновных единиц £ = (£1,..., ет-1) с положительной нормой. Здесь также можно найти фундаментальные области Т и Т", соответствующие наборам £ и £.
6.3. Максимальный порядок
В поле 0>(А) имеется максимальный порядок О. Его базис ¿1,..., ¿V называется фундаментальным, о его вычислении см. [11, гл. II, § 2]. Всё сказанное для порядка ^[А] справедливо и для максимального порядка. В частности, он имеет набор основных единиц £ = (е1,... , ег), набор квазиосновных единиц £ = (¿?1,... , еГ) с положительными нормами и соответствующие им фундаментальные области Т и Т.
7. Пример 2
Пусть р(£) = — 5. Тогда п = 2, а корни многочлена р(£) суть А = ±л/5 и ±2.23605. Поэтому к = 0 I = 2 т = к +1 = 2 и г = т —1 = 1. Основная единица максимального порядка 0 есть е = (1 + А)/2 и 1.61803, т. е. в записи (3) ж1 = ж2 = 1/2. Поскольку N(X) = ж2 — 5x2, то N(е) = — 1 и квазиосновная единица максимального порядка П есть е = е2 = (3 + А)/2 и 2.61803. Основная единица кольца Z[А] есть е3 = 2 + А и 4.23605 с нормой N (е3) = —1. Наконец, квазиосновная единица этого кольца есть е = е6 = 9 + 4А и 17.94421.
Уравнение (11) здесь имеет вид ж2 — 5ж2 = в- Положим в = 4 и найдем все целочисленные решения (ж1,ж2) уравнения
ж? — 5ж2 = 4, (17)
т.е. £(ж!,ж2) = ж1 + ж^л/б € Z [ . Ограничимся решениями ж1,ж2 ^ 0, остальные решения получаются изменением знаков. Поэтому наша основная единица — это е = е6 = 9 + 4\/5 17.94421 < 18 Фундаментальная область Т в координатах ж1,ж2 — это
1 < /1 = ж1 + ж^л/5 < е < 18. На кривой (17) над Т выполнены неравенства 4/е ^ /2 = ж1 — ж^\/5 ^ 4, т. е.
(18)
2 г
ж1 — ж^5 ^ 4.
9
(19)
На плоскости (ж1,ж2) € М2 неравенства (18), (19) выделяют четырёхугольник, ограниченный /1 = 1 /1 = 18 /2 = 2/9 /2 = 4
Х2 к
4-
1 ■""04 4 Х\
Рис. 2: Область в М2, где содержатся все выделенные решения показана штриховкой.
Вершины этого четырёхугольника суть
$1 =
$2 =
$3 =
$4 =
111 1 ,2 + 9,275 _ 975
' 1 9
9 + -, 7 _
1
9' у/5 9л/5
V , \/5 'б ^/5' 2, 10~
<0.61111, 0.1/392), ¿(9.11111, 3.9/524),
(11, 3.13051),
<2.5, -0.6/082).
Теперь для каждого целого х2 : 0 ^ х2 ^ 3, переберём все целые х1 : 1 ^ х1 < 11 и выберем среди таких точек (х1, х2) решения уравнения (17). Получаем точки (2,0), (3,1),
/ -+ 3 д/5 3 -+ у/5
(/, 3). Поскольку отношения (3 + у/Щ /2, (/ + 3л/5) /2, -^ = —й— не лежат в ^\/5],
" 3 + у5 2
то получим три выделенных решения е0 = 2, е0 = 3 + \/5, е0 = / + 3\/5. По теореме 6 все
решения с х1,х2 ^ 0 имеют вид
е = ео (9+4\5)а
г = 1, 2, 3, 0 < а £ Z.
8. Пример 3
Пусть р(е) = е3 _ /е _ 2, все его корни вещественны:
А1 и _2.489288, А2 - _0.289168, А3 - 2.//845/.
Здесь п = т = I = 3, к = 0 г = т _ 1=2. Фундаментальный базис максимального порядка
О есть 1, А, (А + А2) /2. Вычисления проведём в 6 шагов раздела 4.
Шаг 1. Вычисляя значения $(У) на точках е = У1 + У2А + Уз (А + А2) /2 с целыми уг, находим единицы ег = (У1,У2,Уз) е1 = (0, 0,1) е2 = (1, 2, 2) ез = (_2, 0,1) е4 = (_10, _2, 3), е5 = (5,2, _2), ее = (0, 2, _1).
Шаг 2. Вычисляем выпуклую оболочку соответствующих точек Со, Ог = О(У) и получаем V неё 6 двумерных граней. Их логарифмические проекции на плоскость Л-1, Л-2 показаны на рис. 3. На нём логарифмические проекции рёбер показаны прямыми отрезками, хотя они являются криволинейными. Заметим, что ег+з = е-1, г = 1, 2, 3.
Шаг 3. Здесь любая пара ег и £j = ег±1 (г, 2 = 1,...,6) образует набор фундаментальных единиц.
е1 е3 .
0,£1,£2,£З.
Шаг 5. Логарифмическая проекция границы выпуклой оболочки значений О (У) по У £ Z3 с Н'(У) £ Т показана на рис. 4.
Тут имеются две новые вершины: ¿1 = (0,1,1) и ¿2 = (1,1,1). На них $(У) = 2. Имеется четырёхугольная грань с вершинами 0,51, £2, ¿2-
Рис. 3: Логарифмическая проекция вершин, рёбер и граней многогранника дН. Показаны проекции единиц, близкие к нулю. Проекции рёбер выпрямлены.
Рис. 4: Логарифмическая проекция многогранника дН на фундаментальную область.
Шаг 6. Сдвигая фундаментальную область рис. 4 на целочисленные линейные комбинации логарифмов известных единиц, получаем схематическую проекцию всего многогранника дН на плоскость ^1,^2, показанную на рис. 5. На рис. 6 показана точная логарифмическая проекция многогранника дН на плоскость Л-1, Л,2; проекции рёбер криволинейны. Отрезки в ромбах — это ошибки: их не должно быть. Этот рисунок взят из [5], куда он попал из [12].
Рис. 6: Аналог рис. 5 с точными проекциями рёбер. Рёбра, разделяющие ромбы на треугольники, ошибочны.
Этот пример взят у Вороного [6, § 59, пример]. Там найдены две пары основных единиц £2,£з и ^2, £4, но нет аналогов наших рисунков. На самом деле, в этом примере граница dH вычисляется сразу как выпуклая оболочка значений G(Y) по Y £ Z3, ибо размерность задачи n = 3 невелика. Но здесь показано разбиение на шаги, которое может быть полезным при больших размерностях n и m.
9. Предшественники
Для n = 2, k = 0 l = 2 когда m = 2 и r = 1, т. е. для вещественных квадратичных полей способ вычисления основной единицы максимального порядка Q, основанный на разложении в цепную дробь, описан в книге [11, гл. II, § 7]. В конце этой книги в табл. 1 приведены значения основных единиц £ > 1 максимальных порядков полей Q ^ д/d^ для 2 ^ d ^ 101 d £ Z.
Для n = 3 m = 2 (r = 1) и n = 3 m = 3 (r = 2) основные единицы максимальных порядков вычислял Вороной [6] с помощью своего пошагового обобщения цепной дроби. В [3, 5,13] вычислены многоугольники и многогранники dH.
Для n = 4, k = 2 l = 0, т. е. m = 2 (r = 1), Парусников [14] вычислил единицы максимальных порядков полей Q(A) для 41 многочлена (2) с помощью пошагового алгоритма, основанного на выпуклом многоугольнике. Но большинство найденных им единиц не являются основными, а являются лишь их целыми степенями.
Предварительные версии этой статьи — это препринт [15] и статья [16].
СПИСОК ЦИТИРОВАННОЙ ЛИТЕРАТУРЫ
1. Хинчин А. Я. Цепные дроби. 3-е изд. М.: Физматгиз, 1961.
2. Брюно А. Д. Разложения алгебраических чисел в цепные дроби // Журнал вычислительной матем. и матем. физики. 1964. Т. 4, № 2. С. 211-221.
3. Bruno A. D. New generalizations of continued fraction. I // Functiones et Approximatio. 2010. vol. 43, no. 1. Pp. 55-104.
4. Bruno A. D. On geometric methods in works by V. I. Arnold and V. V. Kozlov. Preprint of arXiv, No 1401.6320.
5. Брюно А. Д. Универсальное обобщение алгоритма цепной дроби // Чебышевский сборник (Тула), 2015, том 16, выпуск 2. С. 35-65.
6. Вороной Г. Ф. Об одном обобщении алгорифма непрерывных дробей. Варшава: Из-во Варш. Ун-та, 1896. Также: Собр. соч. в 3-х томах. Киев: Из-во АН УССР, 1952. Т. 1. С. 197-391.
7. Брюно А. Д. Структура многомерных диофантовых приближений// ДАН, 2010. Т. 433, № 5. С. 587-589. *
8. Bruno A. D. Structure of the best diophantine approximations and multidimensional generalizations of the continued fraction // Чебышевский сборник (Тула), 2010. том 11, вып. 1. С. 68-73.
9. Fukuda К. 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."
10. Barber С. В., Dobkin D. P., Huhdanpaa Н. Т. The Quickhull algorithm for convex hulls // ACM Trans, on Mathematical Software, 22(4):469-483, Dec. 1996, http://www.qhull.org.
11. Боревич 3. И., Шафаревич И. Р. Теория чисел. 2-е изд. М.: Наука, 1972.
12. Брюио А. Д., Парусников В. И. Многогранники модулей троек линейных форм // Препринты 1111 \! им. М.В.Келдыша. 2003. № 93. 20 с.
URL: http://library.keldysh.ru/preprint.asp?id=2003-93
13. Брюно А. Д. Обобщения цепной дроби // Чебышевский сборник (Тула), 2006, том 7, вып. 3. С. 4-71.
14. Парусников В. И. Четырёхмерное обобщение алгоритма цепных дробей // Препринты 1111 \! им. М.В.Келдыша. 2011. № 78. 16 с. URL: http://library.keldysh.ru/preprint. asp?id=2011-78.
15. Брюно А. Д. От диофантовых приближений к диофантовым уравнениям // Препринты 1111 \! им. М.В.Келдыша. 2016. № 1. 20 с. URL: http://library.keldysh.ru/preprint. asp?id=2016-l
16. Брюно А. Д. Вычисление наилучших диофантовых приближений и основных единиц алгебраических полей // ДАН, 2016. Т. 468, № 1. С. 7-11.
REFERENCES
1. Khinchin, A. Ya. 1963, Continued fractions, Noordhoff, Groningen.
2. Bruno, A. D. 1964, "The expansion of algebraic numbers into continued fractions", USSR Сотр. Math. Math. Phys., Vol. 4, no. 2, Pp. 1-15.
3. Bruno, A. D. 2010, "New generalizations of continued fraction. I", Functiones et Approximatio. vol. 43, no. 1. Pp. 55-104.
4. Bruno, A. D. 2014, On geometric m,et,hods in works by V. I. Arnold and V. V.Kozlov, Preprint of arXiv, No 1401.6320.
5. Bruno, A. D. 2015, "Universal generalization of the continued fraction algorithm", Chebyshevsky sbornik, vol. 16, no. 2, pp. 35-65.
6. Voronoi, G. F. 1896, On Generalization of the Algorithm of Continued Fraction, Warsawa University.
7. Bruno, A. D. 2010, "The structure of multidimensional Diophantine approximations", Doklady Mathematics, vol. 82, no. 1. Pp. 587-589.
8. Bruno, A. D. 2010, "Structure of the best diophantine approximations and multidimensional generalizations of the continued fraction", Chebyshevsky sbornik, vol. 11, no. 1, pp. 68-73.
9. Fukuda, K. 2008, "Exact algorithms and software in optimization and polyhedral computation", Proceed. ISSAC'08 of XXI International Sym,posiu,m, on Symbolic and Algebraic Commutations, ACM NY, USA, Pp. 333-334.
10. Barber, С. В. k, Dobkin, D. P. k, Huhdanpaa, H. T. 1996, "The Quickhull algorithm for convex hulls", A CM Trans, on Mathematical Software, 22(4):469-483, http://www.qhull.org.
11. Borevich, ZI к Shafarevich, IR 1966, Number Theory, Academic Press.
12. Bruno, A. D. k, Parusnikov, V. I. 2003, "Polvhedra of absolute values for triple of linear forms", Preprint no. 93 of the Keldysh Inst, of Applied Math., Moscow. URL: http://library.keldvsh.ru/ preprint.asp?id=2003-93.
13. Bruno, A. D. 2006, "Generalization of continued fraction", Chebyshevsky sbornik, vol. 7, no. 3, pp. 4-71.
14. Parusnikov, V. I. 2011, "4-dimensional generalization of the continued fractions", Preprint no. 78 of the Keldysh Inst, of Applied Math., Moscow. URL: http://library.keldvsh.ru/preprint. asp?id=2011-78.
15. Bruno, A. D. 2016, "From Diophantine approximations to Diophantine equations", Preprint no. 1 of the Keldysh Inst, of Applied Math., Moscow. URL: http://library.keldvsh.ru/preprint. asp?id=2016-l
16. Bruno, A. D. 2016, "Computation of the best Diophantine approximations and the fundamental units of the algebraic fields", Doklady Mathematics, vol. 93, no. 3. Pp. 243-247.
Институт прикладной математики им. М. В. Келдыша РАН.
Получено 5.05.2016 г.
Принято в печать 12.09.2016 г.
