Русская Википедия:Последовательность Сильвестра

Материал из Онлайн справочника
Перейти к навигацииПерейти к поиску

Файл:Sylvester-square.svg
Графическая демонстрация сходимости суммы 1/2 + 1/3 + 1/7 + 1/43 + … к 1. Каждая строка <math>k</math> квадратов со стороной <math>1/k</math> имеет общую площадь <math>1/k</math>, и все квадраты вместе в точности покрывают большой квадрат с площадью 1. Квадраты со стороной 1/1807 или меньше слишком малы, чтобы их увидеть на рисунке, и они не показаны.

Последовательность Сильвестра — Шаблон:Не переведено 5, в которой каждый очередной член равен произведению предыдущих членов плюс единица. Первые несколько членов последовательности:

Шаблон:Nums, … (Шаблон:OEIS long).

Названа по имени Джеймса Сильвестра, который первым исследовал её в 1880 году. Значения её членов растут как Шаблон:Не переведено 5, а сумма обратных членов образует ряд долей единицы, который сходится к 1 быстрее, чем любой другой ряд дробей единицы с тем же числом членов. Рекуррентное соотношение, которое определяет члены последовательности, позволяет числам в последовательности быть разложенными на множители проще, нежели другие числа того же порядка, но ввиду очень быстрого роста членов ряда полное разложение на простые множители известно только для некоторых членов этой последовательности. Значения, полученные с использованием этой последовательности, используются для образования конечного представления 1 в виде египетской дроби, Шаблон:Не переведено 5 многообразий Эйнштейна и как источник данных для Шаблон:Не переведено 5.

Формальные определения

Формально последовательность Сильвестра можно определить формулой:

<math>s_n = 1 + \prod_{i = 0}^{n - 1} s_i</math>.

Шаблон:Не переведено 5 равно 1, так что <math>s_0 = 2</math>.

Другим образом можно определить последовательность с помощью рекуррентного соотношения:

<math>\displaystyle s_i = s_{i-1}(s_{i-1}-1)+1</math>, где <math>s_0 = 2</math>.

Эквивалентность определений доказывается прямой индукцией.

Общая формула и асимптотики

Члены последовательности Сильвестра растут со скоростью Шаблон:Не переведено 5. В частности, можно показать, что:

<math>s_n = \left\lfloor E^{2^{n+1}}+\frac12 \right\rfloor,</math>

где число <math>E</math> приблизительно равно Шаблон:Num[1]. Эта формула приводит к следующему алгоритму:

s0 — ближайшее целое к E2; s1 — ближайшее целое к E4; s2 — ближайшее целое к E8; для sn, берём E2, возводим в квадрат n раз и берём ближайшее целое.

Этот алгоритм был бы приемлемым, если бы мы имели лучший путь вычисления E вместо вычисления чисел sn с последующим вычислением квадратных корней.

Рост элементов последовательности Сильвестра со скоростью двойной экспоненты совершенно не удивителен, если сравнивать с последовательностью чисел Ферма Fn. Числа Ферма часто задаются по формуле двойной экспоненты <math>2^{2^n}+1</math>, но их можно также задать по формулам умножения, подобным формулам последовательности Сильвестра:

<math>F_n = 2 + \prod_{i = 0}^{n - 1} F_i.</math>

Связь с египетскими дробями

Доли единицы, образованные числами, обратными значениям последовательности Сильвестра, образуют бесконечный ряд:

<math>\sum_{i=0}^{\infty} \frac1{s_i} = \frac12 + \frac13 + \frac17 + \frac1{43} + \frac1{1807} + \cdots.</math>

Частичные суммы этого ряда имеют простую форму

<math>\sum_{i=0}^{j-1} \frac1{s_i} = 1 - \frac{1}{s_j-1} = \frac{s_j-2}{s_j-1}.</math>

Это можно доказать по индукции или прямо, если заметить, что из рекурсии вытекает

<math>\frac{1}{s_i-1}-\frac{1}{s_{i+1}-1}=\frac{1}{s_i},</math>

Таким образом, сумма телескопического ряда будет равна

<math>\sum_{i=0}^{j-1} \frac{1}{s_i} = \sum_{i=0}^{j-1} \left( \frac{1}{s_i-1}-\frac{1}{s_{i+1}-1} \right) = \frac{1}{s_0-1} - \frac{1}{s_j-1} = 1 - \frac{1}{s_j-1}.</math>

Поскольку последовательность частичных сумм (sj−2)/(sj−1) сходятся к единице, весь ряд образует бесконечное представление единицы в виде египетской дроби:

<math>1 = \frac12 + \frac13 + \frac17 + \frac1{43} + \frac1{1807} + \cdots.</math>

Можно найти конечные представления единицы в виде египетской дроби любой длины путём обрывания этого ряда и вычитанием единицы из последнего знаменателя:

<math>1 = \tfrac12 + \tfrac13 + \tfrac16, \quad 1 = \tfrac12 + \tfrac13 + \tfrac17 + \tfrac1{42}, \quad 1 = \tfrac12 + \tfrac13 + \tfrac17 + \tfrac1{43} + \tfrac1{1806},\quad \dots.</math>

Сумма первых k членов бесконечного ряда даёт ближайшую нижнюю оценку единицы k-членными египетскими дробями.[2] Например, первые четыре члена добавляются к 1805/1806, а потому любая египетская дробь в открытом интервале (1805/1806,1) требует по меньшей мере пять членов.

Можно рассматривать последовательность Сильвестра как результат жадного алгоритма для египетских дробей, который на каждом шаге выбирает наименьший возможный делитель, оставляющий частичную сумму меньше единицы. Также члены ряда после первого можно рассматривать как делители Шаблон:Не переведено 5 числа 1/2.

Единственность быстрорастущих рядов с рациональными суммами

Как заметил сам Сильвестр, последовательность Сильвестра, похоже, единственная, которая имеет такую скорость роста, при этом сходясь к рациональному числу. Эта последовательность показывает пример того, что скорость роста двойной экспоненты недостаточна для того, чтобы последовательность целых чисел была рациональной последовательностьюШаблон:Sfn.

Из результата Бадеа[3] следует, что если последовательность целых чисел <math>a_n</math> растёт достаточно быстро, так, что

<math>a_n\ge a_{n-1}^2-a_{n-1}+1,</math>

и если ряд

<math>A=\sum\frac1{a_i}</math>

сходится к рациональному числу A, то для всех n, начиная с некоторого места, эта последовательность должна удовлетворять рекуррентному соотношению

<math>a_n= a_{n-1}^2-a_{n-1}+1</math>,

что можно использовать для определения последовательности Сильвестра.

Эрдёш[4] высказал гипотезу, что в результатах такого типа границу неравенства роста последовательности можно заменить на более слабое условие

<math>\lim_{n\rightarrow\infty} \frac{a_n}{a_{n-1}^2}=1.</math>

Делимость и разложение

Если i < j, из определения следует, что sj ≡ 1 (mod si). Таким образом, любые два члена последовательности Сильвестра взаимно просты. Последовательность можно использовать для доказательства бесконечности числа простых чисел, поскольку любое простое число может делить максимум одно число в последовательности. Никакой простой множитель числа в последовательности не может быть сравним с 5 (mod 6), и последовательность можно использовать для доказательства, что существует бесконечно много простых чисел, сравнимых с 7 (mod 12).Шаблон:Sfn

Шаблон:Unsolved Много остаётся неизвестного о разложении на множители членов последовательности Сильвестра. Например, неизвестно, являются ли все члены последовательности свободными от квадратов, хотя все члены, для которых известно разложение на простые, таковыми являются.

Как пишет Варди Шаблон:Harv, легко установить, какой из членов последовательности Сильвестра (если такой есть) делится на простое число p — просто вычисляем вычеты членов последовательности по модулю p согласно рекуррентной формуле, пока не встретится ноль (по модулю p) или встретится такой же остаток. Используя эту технику, Варди нашёл, что 1166 из первого миллиона простых чисел являются делителями чисел Сильвестра,[5] и ни один квадрат этих простых чисел не делит числа Сильвестра. Множество простых чисел, которые могут оказаться делителями членов ряда Сильвестра, имеет плотность ноль во множестве всех простых чисел. Более того, согласно Джонсу Шаблон:Sfn число таких простых меньших x равно <math>O(\pi(x) / \log\log\log x)</math>.Шаблон:Sfn

В следующей таблице приведены известные разложения этих чисел, (за исключением первых четырёх, являющихся простыми):[6]

n Множители sn
4 13 × 139
5 3263443, простое
6 547 × 607 × 1033 × 31051
7 29881 × 67003 × 9119521 × 6212157481
8 5295435634831 × 31401519357481261 × 77366930214021991992277
9 181 × 1987 × 112374829138729 × 114152531605972711 × 35874380272246624152764569191134894955972560447869169859142453622851
10 2287 × 2271427 × 21430986826194127130578627950810640891005487 × P156
11 73 × C416
12 2589377038614498251653 × 2872413602289671035947763837 × C785
13 52387 × 5020387 × 5783021473 × 401472621488821859737 × 287001545675964617409598279 × C1600
14 13999 × 74203 × 9638659 × 57218683 × 10861631274478494529 × C3293
15 17881 × 97822786011310111 × 54062008753544850522999875710411 × C6618
16 128551 × C13335
17 635263 × 1286773 × 21269959 × C26661
18 50201023123 × 139263586549 × 60466397701555612333765567 × C53313
19 775608719589345260583891023073879169 × C106685
20 352867 × 6210298470888313 × C213419
21 387347773 × 1620516511 × C426863
22 91798039513 × C853750

Здесь Pn и Cn обозначают простые и составные числа с n десятичными цифрами.

Приложения

Бойер, Галики и Коллар Шаблон:Harv использовали свойства последовательности Сильвестра для определения большого числа Шаблон:Не переведено 5 многообразий Эйнштейна, имеющих дифференциальную топологию сфер нечётных размерностей или экзотических сфер. Они показали, что число различных сасакианских эйнштейновских метрик на топологической сфере размерности 2n − 1 по меньшей мере пропорционально sn, а потому растёт со скоростью двойной экспоненты (от n).

Как пишут Галамбос и Вогингер Шаблон:Harv, Браун Шаблон:Harv и Лианг Шаблон:Harv использовали значения, полученные из последовательности Сильвестра, для построения примеров нижней границы для Шаблон:Не переведено 5 алгоритмов упаковки в контейнеры. Зайден и Вогингер Шаблон:Harv подобным же образом использовали последовательность для нижней границы производительности двумерного алгоритма раскроя[7].

Задача Знама касается множеств чисел, таких, что каждое число в множестве делит, но не равно произведению всех остальных множеств плюс единица. Без условия эквивалентности значения последовательности Сильвестра решают эту задачу. Если это условие ставится, имеется другое решение, получаемое из рекуррентного соотношения, подобного определению последовательности Сильвестра. Решения задачи Знама имеют приложения к классификации особых точек поверхностей (Brenton, Hill 1988) и теории недетерминированных конечных автоматов.Шаблон:Sfn

Куртис Шаблон:Harv описывает приложение ближайшего приближения к единице k-членными суммами долей единицы к нижней границе числа делителей любого совершенного числа, а Мюллер Шаблон:Harv использует то же самое свойство для Шаблон:Не переведено 5 числа некоторых групп.

См. также

Примечания

Шаблон:Примечания

Литература

Ссылки

Шаблон:Rq

  1. В книге Грэма, Кнута и Паташника Шаблон:Harv данное утверждение приведено в качестве упражнения. См. также Голомба Шаблон:Harv.
  2. Это утверждение обычно приписывается Кёртису Шаблон:Harv, но Миллер Шаблон:Harv сделал то же самое утверждение в более ранней работе. См. также Шаблон:Harvnb, Шаблон:Harvnb и Шаблон:Harv.
  3. Шаблон:Harv
  4. Шаблон:Harv, обзор работ, посвящённых этой гипотезе — Шаблон:Harv, см. также Шаблон:Harvnb.
  5. Здесь, похоже, возникла опечатка, поскольку Андерсен (Andersen) нашёл 1167 простых делителей в этом диапазоне.
  6. Все простые множители p чисел Сильвестра sn с p < 5Шаблон:E и n ≤ 200 перечислены Варди. Кен Такусагава (Ken Takusagawa) перечислил разложения вплоть до s9 Шаблон:Wayback и разложение s10 Шаблон:Wayback. Остальные разложения взяты из списка разложений последовательности Сильвестра Шаблон:Wayback, которую ведёт Йенс Круз Андерсен (Jens Kruse Andersen). По состоянию на 13/06/2014.
  7. В своей работе Зайден и Вогингер ссылаются на последовательность Сильвестра как на «последовательность Зальцера», опираясь на работу Зальцера Шаблон:Harv о ближайшей аппроксимации.