Алгоритм поиска целочисленных соотношений
Целочисленное соотношение между набором вещественных чисел — это набор целых чисел , не все из которых равны нулю, таких, что
Алгоритм поиска целочисленных соотношений — это алгоритм поиска целочисленных соотношений между числами. В частности, если заданы вещественные числа с определённой точностью, алгоритм либо находит целочисленное соотношение между ними, либо определяет, что такой связи нет для коэффициентов, абсолютная величина которых меньше некоторой верхней границы[англ.] [1].
История
Для случая n = 2 расширение алгоритма Евклида может определить, существует ли целочисленное соотношение между двумя вещественными числами и . Алгоритм образует последовательность элементов разложения в непрерывную дробь. Если существует целочисленная взаимосвязь между числами, то их частное рационально и алгоритм, в конечном счёте, завершится.
- Алгоритм Фергюсона — Форкейда опубликовали в 1979 Фергюсон[англ.] и Форкейд[2]. Хотя в статье идёт речь о любом n, не совсем ясно, решает ли статья полностью задачу, поскольку в ней отсутствуют некоторые детали алгоритма, нет доказательств и точных границ, что существенно для достоверной имплементации.
- Первым алгоритмом с полным доказательством был ЛЛЛ-алгоритм, который разработали Арьен Ленстра, Хендрик Ленстра и Ласло Ловас в 1982[3].
- Юхан Хостад, Беттина Джаст, Джефри Лагариас и Клаус-Петер Шнорр разработали алгоритм HJLS в 1986[4][5].
- Фергюсон разработал в 1988 алгоритм PSOS[6]
- Фергюсон и Бейли разработали алгоритм PSLQ в 1992 и в 1999 в значительной степени упростили (вместе с Арно)[7][8]. В 2000 Джек Донгарра и Фрэнсис Салливан [9] включили алгоритм PSLQ в «Первую десятку алгоритмов столетия», хотя и признано, что большей частью он эквивалентен алгоритму HJLS[10][11].
- Алгоритм ЛЛЛ улучшали множество авторов. Современная имплементация алгоритма ЛЛЛ может решать задачи поиска целочисленных соотношений с n, большим 500.
Приложения
Алгоритм поиска целочисленных соотношений имеет многочисленные приложения. Первое приложение — определение, не является ли заданное вещественное число x алгебраическим, для чего производится поиск целочисленного соотношения между множеством степеней числа x {1, x, x2, ..., xn}. Второе приложение — поиск целочисленной связи между вещественным числом x и набором математических констант, таких как e, и ln(2), что приводит к выражению x в виде линейной комбинации этих констант.
Типичным подходом в экспериментальной математике является применение численных методов и арифметики произвольной точности для поиска приближённого значения бесконечного ряда, бесконечного произведения или интеграла с большой точностью (обычно берётся по меньшей мере 100 значащих цифр), а затем используется алгоритм поиска целочисленного соотношения для определения целочисленной связи между этим значением и набором математических констант. Если целочисленное соотношение найдено, оно говорит о возможном выражении в замкнутой форме[англ.] исходного ряда, произведения или интеграла. Полученная гипотеза затем может быть проверена с помощью формальных алгебраических методов. Чем выше используемая точность вычислений, тем выше уверенность, что найденное целочисленное соотношение не является просто числовым артефактом[англ.].
Заметным успехом такого подхода было использование алгоритма PSLQ для поиска целочисленного соотношения, которое привело к формуле Бэйли — Боруэйна — Плаффа для числа . Алгоритм PSLQ также помог найти новые тождества, в которые входит многомерная дзета-функция[англ.], и их появление в квантовой теория поля. Алгоритм PSLQ помог в распознании точек бифуркации логистического отображения. Например, если B4 является четвёртой точкой бифуркации логистического отображения, константа α=-B4(B4-2) является корнем многочлена 120-й степени с максимальным коэффициентом 25730[12][13]. Алгоритмы поиска целочисленных соотношений комбинируются с таблицами математических констант высокой точности и эвристическими методами поиска, такими как обратный символьный калькулятор[англ.] или инвертор Плуффера[англ.].
Поиск целочисленных соотношений может быть использован для разложения многочленов высокой степени[14].
Примечания
- ↑ Поскольку вещественные числа могут быть заданы только с конечной точностью, алгоритм без задания границы коэффициентов всегда найдёт целочисленную связь при достаточно больших коэффициентах. Результат интересен, только когда величина коэффициентов в целочисленном соотношении мала по сравнению с точностью задания вещественных чисел.
- ↑ Weisstein, Eric W. Integer Relation (англ.) на сайте Wolfram MathWorld.
- ↑ Weisstein, Eric W. LLL Algorithm (англ.) на сайте Wolfram MathWorld.
- ↑ Weisstein, Eric W. HJLS Algorithm (англ.) на сайте Wolfram MathWorld.
- ↑ Håstad, Just, Lagarias, Schnorr, 1989, с. 859–881.
- ↑ Weisstein, Eric W. PSOS Algorithm (англ.) на сайте Wolfram MathWorld.
- ↑ Weisstein, Eric W. PSLQ Algorithm (англ.) на сайте Wolfram MathWorld.
- ↑ Ferguson, Bailey, 1992.
- ↑ Cipra, 2000.
- ↑ Chen, Stehlé, Villard.
- ↑ Ferguson, Bailey, Arno.
- ↑ Bailey, Broadhurst, 2000, с. 1719-1736.
- ↑ Kotsireas, Karamanos, 2004, с. 2417-2423.
- ↑ van Hoeij, 2002, с. 167-189.
Литература
- Johan Håstad, Bettina Just, Jeffrey Lagarias, Claus-Peter Schnorr. Polynomial time algorithms for finding integer relations among real numbers. (англ.) // SIAM J. Computing. — 1989. — Vol. 18. — P. 859–881.
- Helaman R. P. Ferguson, David H. Bailey. A Polynomial Time, Numerically Stable Integer Relation Algorithm (англ.) // RNR Technical Report RNR-91-032. — 1992.
- Barry A. Cipra. The Best of the 20th Century: Editors Name Top 10 Algorithms (англ.) // SIAM News. — 2000. — Vol. 33, iss. 4.
- David H. Bailey, David J. Broadhurst. Parallel Integer Relation Detection: Techniques and Applications (англ.) // Mathematics of Computation. — 2000. — Vol. 70, iss. 236. — P. 1719-1736.
- I. S. Kotsireas, K. Karamanos. Exact Computation of the bifurcation Point B4 of the logistic Map and the Bailey-broadhurst Conjectures (англ.) // I. J. Bifurcation and Chaos. — 2004. — Iss. 14 (7).
- M. van Hoeij. Factoring polynomials and the knapsack problem (англ.) // J. of Number Theory. — 2002. — Iss. 95. — P. 167-189.
- Jingwei Chen, Damien Stehlé, Gilles Villard. A New View on HJLS and PSLQ: Sums and Projections of Lattices. In the proceedings of ISSAC'13 (англ.). ISSAC'13 (June 26-29, 2013). Дата обращения: 24 февраля 2017.
- Helaman R. P. Ferguson, David H. Bailey, Steve Arno. ANALYSIS OF PSLQ, AN INTEGER RELATION FINDING ALGORITHM (англ.) (PDF) (3 июля 1997). Дата обращения: 24 февраля 2017.
- M. van Hoeij. Factoring polynomials and the knapsack problem (англ.) // J. of Number Theory. — Vol. 2002, iss. 95.
Ссылки
- Recognizing Numerical Constants by David H. Bailey and Simon Plouffe
- Ten Problems in Experimental Mathematics Архивная копия от 10 июня 2011 на Wayback Machine by David H. Bailey, Jonathan M. Borwein, Vishaal Kapoor, and Eric W. Weisstein