The Golomb--Keller formula expresses the next prime $p_{n+1}$ as a recurrence relation in terms of the first $n$ primes $p_1, \ldots, p_n$ using the Riemann zeta function and an Euler product, but requires taking a limit as $s \to \infty$, rendering it non-constructive. We transform this asymptotic formula into an effective recurrence by proving that a finite parameter $s \leq p_n$ suffices when combined with the ceiling function, establishing a constructive method valid for all $n \geq 1$.
The minimal integer parameter $s_n$ (OEIS A389650) reveals deep connections to prime constellations. We prove $\liminf_{n\to\infty} Ï_n = 0$ unconditionally, where $Ï_n = s_n/p_n$. The limit superior $C = \limsup Ï_n$ satisfies $\log Ï\lesssim C \leq 0.4332$, where $Ï\approx 1.46557$ is the supergolden ratio. The lower bound is conditional on the twin prime conjecture; the upper bound is unconditional. The constant $C$ relates to the densest admissible prime constellation, connecting to the Hardy--Littlewood conjectures.
The method extends to Dirichlet L-functions, yielding other effective formulas for calculating $p_{n+1}$ but also for predicting residues of $p_{n+1}$ modulo any integer with reduced precision requirements.
- ID статьи: 2508.02690
- Название: An effective analytic recurrence for prime numbers
- Автор: Benoit Cloitre
- Классификация: math.NT (Теория чисел), math.HO (История и обзор)
- Дата публикации: 12 октября 2025 г. (arXiv v2)
- Ссылка на статью: https://arxiv.org/abs/2508.02690v2
Формула Голомба-Келлера выражает следующее простое число pn+1 через первые n простых чисел p1,…,pn посредством рекуррентного соотношения, использующего дзета-функцию Римана и произведение Эйлера, но требует предельного перехода s→∞, что делает её неконструктивной. В данной работе доказано, что конечный параметр s≤pn в сочетании с функцией потолка достаточен для преобразования этой асимптотической формулы в эффективную рекуррентность, устанавливающую конструктивный метод, действительный для всех n≥1.
Минимальный целочисленный параметр sn (OEIS A389650) раскрывает глубокую связь с простыми созвездиями. Автор безусловно доказывает, что liminfn→∞σn=0, где σn=sn/pn. Верхний предел C=limsupσn удовлетворяет условию logψ≲C≤0.4332, где ψ≈1.46557 — супер-золотое сечение. Нижняя граница зависит от гипотезы о близнецах; верхняя граница безусловна. Константа C связана с наиболее плотными допустимыми простыми созвездиями, что соединяет с гипотезой Харди-Литтлвуда.
Метод распространяется на L-функции Дирихле, порождая другие эффективные формулы для вычисления pn+1 и позволяя предсказывать остаток pn+1 по модулю произвольного целого числа с пониженными требованиями к точности.
Поиск явных формул для простых чисел — классическая задача теории чисел. Хотя существуют прямые нерекурсивные формулы (такие как формула Виллана, формула Миллса), они вычислительно неосуществимы. Данная работа сосредоточена на рекуррентных соотношениях, то есть выражении pn+1 через p1,…,pn.
- Ганди 4 первым использовал примориал и функцию Мёбиуса для таких рекуррентностей
- Ванден Эйнден 19 упростил доказательство
- Якимчук 9 обобщил метод
- Голомб 5 независимо открыл формулу, используя аналитическую теорию чисел; позже переоткрыта Келлером 10
Классическая формула Голомба-Келлера имеет вид:
pn+1=lims→∞[(∏k=1n(1−pks1))ζ(s)−1]−1/s
Основная проблема этой формулы заключается в необходимости предельного перехода s→∞, что делает её непригодной для практических вычислений.
Данная работа применяет противоположный подход: сохраняет полное разложение дзета-функции, вычисляемое до рабочей точности, но использует конечное значение s. Таким образом, усекается показатель степени, а не ряд, что делает формулу конструктивной.
- Конструктивная рекуррентная формула: Доказано, что для всех n≥1 существует минимальное целое число sn такое, что:
pn+1=⌈(−1+ζ(sn)∏j=1n(1−pjsn1))−1/sn⌉
- Эффективные границы:
- Доказано sn≤2pn с использованием постулата Бертрана (теорема 10)
- Доказано sn≤pn с использованием теоремы Нагуры (теорема 12)
- Анализ асимптотического поведения:
- Безусловно доказано liminfn→∞σn=0 (предложение 13)
- Установлены границы для C:=limsupn→∞σn: 0.3823≲C≤0.4332
- Связь с простыми созвездиями: Обнаружена нижняя граница logψ≈0.3823 (зависит от гипотезы о близнецах), где ψ — супер-золотое сечение
- Расширение на L-функции Дирихле: Позволяет предсказывать свойства остатков pn+1(mod4) и т.д.
- Численные данные: Предоставлены значения sn для n=1 до 200 (OEIS A389650)
Дано первые n простых чисел p1,p2,…,pn; требуется конструктивно вычислить следующее простое число pn+1.
Определяется ключевая функция:
Dn(s)=∑k≥1gcd(k,Pn)=1k−s
где Pn=∏j=1npj — примориал первых n простых чисел.
Лемма 3: Для ℜ(s)>1,
Dn(s)=ζ(s)∏j=1n(1−pj−s)
Определяется h(s)=(Dn(s)−1)−1/s и доказывается:
- h(s)<pn+1 для всех s>1
- lims→∞h(s)=pn+1
- h(s) строго возрастает на (1,∞)
Предложение 6: Для каждого n≥1 существует единственное sn∗>1 такое, что h(sn∗)=pn+1−1.
Определяется sn=⌊sn∗⌋+1 как минимальный целочисленный параметр.
- Компромисс между точностью и усечением: В отличие от метода Келлера, усекающего сумму ряда, данная работа сохраняет полную дзета-функцию ζ(s), но использует конечное значение s
- Трюк с функцией потолка: Умело используется свойство ⌈h(s)⌉=pn+1 тогда и только тогда, когда s>sn∗
- Техника интегральных границ: Используются тонкие интегральные сравнения для контроля ошибок хвостовых членов
- Использованы PARI/GP и библиотека mpmath для Python с высокой точностью
- Для n≤200 требуется 100-битная точность
- Для n≈500 требуется примерно 2500-битная точность (из-за усиления эффектов сокращения)
Теоретические границы sn≤pn проверены прямыми вычислениями для всех n=1,…,200.
Вычисление p7=17 (начиная с n=6):
- Использование s=2p6=26: h(26)≈16.941817904, получаем p7=⌈h(26)⌉=17
- Фактическое минимальное значение s6=8: h(8)≈16.5189076
- Теорема 10: sn≤2pn для всех n≥1
- Теорема 12: sn≤pn для всех n≥1 (через теорему Нагуры)
- Предложение 13: liminfn→∞σn=0 (безусловно)
- Теорема 14: C=limsupn→∞σn≤0.4332 (безусловно)
- Теорема 15: При гипотезе о близнецах C>logψ≈0.38225
Анализ данных для n=1 до 200 показывает:
- После удаления пяти выбросов распределение σn подобно бета-распределению
- Среднее ≈ 0.291, медиана ≈ 0.277, стандартное отклонение ≈ 0.087
- Параметры подгонки: Beta(α ≈ 7.64, β ≈ 18.62)
- Теоретическая мода ≈ 0.274, согласуется с эмпирической медианой
Теорема 18: Для любого c>c0≈0.5956 существует N0(c) такое, что для всех n≥N0(c) формула с использованием s=cpn корректно вычисляет pn+1.
- Ганди (1971): Первая рекуррентность с использованием примориала и функции Мёбиуса
- Голомб (1976): Введение методов аналитической теории чисел
- Келлер (2007): Независимое переоткрытие с альтернативным выводом
- Данная работа (2025): Первое преобразование формулы в конструктивную
- Формула Виллана: Прямая, но вычислительно неосуществима
- Формула Миллса: Основана на константе, но требует предварительного знания простых чисел
- Методы решета: Практичны, но не дают рекуррентной структуры
- Данный метод: Теоретически интересен, но вычислительно сложен
- Конструктивный прорыв: Впервые преобразована формула Голомба-Келлера из асимптотической в эффективно вычислимую
- Глубокие связи: Раскрыта внутренняя связь параметров рекуррентности простых чисел с созвездиями простых чисел и распределением промежутков
- Теоретические границы: Установлены точные границы и асимптотическое поведение параметра sn
- Супер-золотое сечение: Обнаружено новое применение в теории простых чисел
- Вычислительная сложность: O(npn3logpn) битовых операций; хотя полиномиальна, практически неосуществима
- Требования к точности: С ростом n требуется экспоненциальный рост рабочей точности
- Зависимость: Некоторые результаты зависят от недоказанных гипотез (например, гипотезы о близнецах)
- Оптимизация вычислений: Поиск методов снижения требований к точности
- Теоретическое совершенствование: Устранение зависимости от недоказанных гипотез
- Обобщение и применение: Расширение на другие теоретико-числовые функции
- Численные исследования: Вычисление sn для больших диапазонов с целью верификации гипотез
- Теоретическая инновация: Успешное решение долгостоящей конструктивной проблемы
- Методологическая строгость: Использованы глубокие методы аналитической теории чисел
- Глубокие связи: Установлена связь между рекуррентной формулой и теорией простых созвездий
- Богатые данные: Предоставлены подробная численная верификация и статистический анализ
- Ограниченная практичность: Хотя теоретически конструктивна, вычислительные затраты чрезмерны
- Вызовы точности: Требования высокой точности ограничивают масштабируемость метода
- Условность некоторых результатов: Важные результаты зависят от недоказанных гипотез
- Теоретический вклад: Предоставляет новую перспективу теории рекуррентностей для простых чисел
- Методологическая ценность: Демонстрирует, как преобразовать асимптотическую формулу в эффективный алгоритм
- Междисциплинарные связи: Соединяет аналитическую теорию чисел, вычислительную теорию чисел и геометрию простых чисел
- Теоретические исследования: Распределение простых чисел, теория промежутков, исследование созвездий
- Численные эксперименты: Верификация в малых масштабах и исследование закономерностей
- Педагогическое применение: Классическое применение методов аналитической теории чисел
Ключевые источники включают:
- 5 S. W. Golomb, Formulas for the next prime, Pacific J. Math. 63 (1976), 401–404
- 10 J. B. Keller, A recursion equation for prime numbers, arXiv:0711.3940, 2007
- 4 J. M. Gandhi, Formulae for the nth prime, Proc. Washington State Univ. Conf. Number Theory (1971), 96–107
- 12 J. Nagura, On the interval containing at least one prime number, Proc. Japan Acad. 28 (1952), 177–181
Данная статья имеет значительное теоретическое значение в области теории чисел. Хотя её практическое применение в вычислениях ограничено, теоретические идеи и методологические инновации открывают новые направления в исследовании теории простых чисел. Особенно примечательны обнаруженные связи с супер-золотым сечением и простыми созвездиями, которые демонстрируют неожиданные глубокие структуры в теории чисел.