В Е С Т

ВЕСТНИК ПЕРМСКОГО УНИВЕРСИТЕТА
2010
Математика. Механика. Информатика
Вып. 2(2)
УДК 511.2
Возможности вычислительных методов
в проблемах теории чисел
Е. Л. Тарунин
Пермский государственный университет, Россия, 614990, Пермь, ул. Букирева, 15
tarunin@psu.ru; (342) 2-237-10-31
Показано применение персональных компьютеров к решению двух проблем: выявлению
закономерности распределения простых чисел и подтверждению гипотезы Эйлера–
Гольдбаха.
Ключевые слова: теория простых чисел; численные методы; распределение простых чисел;
проблема Гольбаха.

нить, что теорией чисел занимались такие математики, как Ферми, Гаусс, Эйлер, Лагранж,
Риман, Дирихле, Чебышев, Виноградов и другие. Отметим, что И.М.Виноградов работал в
Пермском университете в 1918–1920 гг. вначале доцентом, а затем профессором.
В Интернете есть предупреждение для
тех, кто думает, что проблему можно просто
решить. Это предупреждение звучит так‫" ׃‬Не
пытайтесь решать великие проблемы, не поняв теории, которые их окружают. Сэкономьте нервы и себе, и окружающим". Невозможно не согласиться с этим предупреждением. И
все же мы полагаем, что у прикладников
больше возможностей подступиться к штурму
проблем теории чисел, несмотря на то, что им
(в отличие от математиков) теорию чисел не
излагают. Для ознакомления с теорией чисел
существует много учебников [1–5].
В статье содержится 8 разделов. В первых трех разделах приводятся в основном
сведения из теории простых чисел, которые
будут анализироваться в следующих разделах.
В теории чисел применяются тонкие
аналитические методы [1–3]. Несмотря на малость практических приложений, выдающиеся
математики называют эту теорию "жемчужиной математики". Со второй половины XX в.
для изучения теории чисел стали применять
вычислительную технику. И все же возможности вычислительной техники в этой области
используются недостаточно. Перечисление 18
задач теории, которые еще ждут своего решения, содержится в последней главе книги [2].
С помощью ПК удалось подкорректировать коэффициенты в оценках распределения простых чисел и найти новые закономерности. Что касается бинарной гипотезы Гольдбаха (любое четное число может быть представлено суммой двух простых чисел), то следует отметить, что новым результатом в этой
области будет доказательство выполнимости
гипотезы для чисел более 1010. Такая проверка
потребует значительных усилий от опытного
программиста. Но игра стоит свеч. За решение
проблемы обещан миллион долларов. Правда,
в Интернете есть непроверенная информация
о том, что условие о награде отменено в 2002
г. Кроме того, там же сказано, что будто бы
минский математик доказал гипотезу, но пока
его доказательство не проверено. Но даже если премии нет, решение этой проблемы (ей
более 3 столетий) даст большое моральное
удовлетворение и славу. Достаточно вспом-
1. Распределение простых чисел
Существование бесконечного множества простых чисел доказывается просто. Однако вопрос о том, как часто среди натуральных чисел встречаются простые и как простые числа распределены среди натуральных,
оказывается весьма сложным. Отметим, что в
данной статье мы касаемся лишь элементар-
© Е. Л. Тарунин, 2010
15
Е. Л. Тарунин
ных методов исследования. Так называются
[6] все методы, которые не опираются на
применение принципа Римана о нулях дзетафункции или на положения теории функций
комплексных переменных.
Число простых чисел, меньших или
равных х, как это принято, будем обозначать
функцией  ( x) . Этой функции приписывают
гладкость, что не соответствует действительности. Для того чтобы подчеркнуть целочисленные свойства соответствующей функции, в
нашей работе будем использовать следующие
обозначения‫ ׃‬i – номер простого числа, а N[i]
– значение этого числа. При этом будем считать, что первое простое число равно двум,
т.е. N(1)=2, N(2)=3, N(3)=5 и т.д. Функции
распределения  ( x) при этом соответствует
ступенчатая функция i(N). Производная этой
ступенчатой функции изменяется в широких
пределах‫ ׃‬от 1/2 на близнецах до 1/ d m ( d m –
максимальное расстояние между соседними
простыми числами).
Поиск формул, порождающих простые
числа, предпринимал еще Ферма (1601–1665).
Он высказал предположение, что все числа
n
вида 2 2  1 (n – натуральное число) являются
простыми [1]. Его предположение оказалось
неверным. В 1732 г. Эйлер показал, что при
n=5 число является составным. В 1875 г.
И.М.Первушин обнаружил, что при n=12 и
n=23 также получаются составные числа. В
1837 г. Лежен–Дирихле доказал, что среди
членов арифметической прогрессии вида
a  x  b с взаимно простыми числами a, b
имеется бесконечное множество простых чисел [1]. Доказывается теорема, что в качестве
формулы, порождающей все простые числа,
не может быть полином с целыми коэффициентами.
Оценку роста функции распределения
получали и уточняли многие ученые. Коснемся лишь тех оценок, которые будут проанализированы. К простейшим оценкам относятся
неравенства [5]
потеза о существовании определенного числа
простых чисел на интервале ( x,  x ) при
Первое неравенство является одним из
доказательств бесконечности простых чисел.
Последнее неравенство вытекает как следствие гипотезы Бертрана о существовании
простого числа на интервале x, 2 x . Справедливость этой гипотезы доказал Чебышев.
Расширением гипотезы Бертрана является ги-
видимо, в предварительных оценках. Доказательство П.Л.Чебышева о существовании коэффициентов a, A , обеспечивающих выполнение неравенств (4), содержится, например, в
[2]. Значения множителей позднее уточнялись. Уточняются они и в нашей работе. Из
теоремы Чебышева, в частности, следует, что

2.
В 1808 г. Лежандр опубликовал эмпирическую формулу [1, 2], предложенную для
больших значений x (от 104 до 106)‫׃‬
 x   f L x  
x
.
ln x   1.08366
(1)
Позднее Чебышевым было указано, что
более близкое значение к реальному распределению простых чисел дают функции
f1 x  
x
,
ln x   1
(2)
dx
 li  x   li 2 .
ln s 
2
(3)
x
FG  x   
Обрабатывая имеющиеся таблицы,
Гаусс еще в юношеские годы нашел, что
функция (3) хорошо аппроксимирует зависимость  ( x) , но опубликовал этот результатов
значительно позднее. Функция li x  в (3)
называется интегральным логарифмом‫׃‬
x
 1 s ds
ds 
.
li  x   lim s 0  
 

 0 ln s  1 s ln s  
Так как li 2  1.04 [5] и в теории инте-
ресуются большими значениями х, то обычно
не делают различия между интегральным логарифмом и FG x  , при этом функцию FG x 
также называют интегральным логарифмом.
В 1848–50 гг. были опубликованы работы П.Л.Чебышева, касающиеся роста функции  ( x) . Чебышев доказал, что при достаточно больших x справедливы неравенства
a  f x    x   A  f x ,
f x  
x
. (4)
ln x 
Для множителей в этих оценках он
предложил значения a  0.92129 , A  1.0555
[5]. В [14] указаны более грубые значения для
коэффициентов
a  ln2 / 2  0.347,
A  2 ln2  1.386 , полученные Чебышевым,
 x   ln x   1, pi  e i 1 , pi 1  2 pi .
16
Возможности вычислительных методов в проблемах теории чисел
при увеличении x отношение  x  / x  0 .
Чебышевым было высказано утверждение о
существовании асимптотической формулы
при x  
 x   lim
x
.
ln x 
(3). Из него следует, что относительная ошибка вычисления  x
 1 x  
1
ln  x 
становится менее 5% при x  10 9.
(5)
В 1914 г. Литллвуд доказал, что величина 1 x принимает бесконечно много
значений, как положительных, так и отрицательных. Анализ этой величины на РС (см.
ниже) вызывает сомнение в правильности
утверждения Литллвуда.
Кроме неравенств Чебышева (4) в теории чисел предлагались следующие неравенства:
Считается, что Чебышевым был решен
вопрос о существовании простой функции,
предложенной еще Гауссом, которая служит
наилучшим приближением для  x ‫׃‬
k  1! x  ...
ds
x
x

 2  ... k
ln s  ln x  ln x 
ln x 
2
x
li x   
Существование предела (5) или, точнее,
его следствия при x   дает утверждение
lim   lim
 x   ln x 
x
x
x
.
  x  
ln x   2
ln x   4
1,
 рассмат  x с па-
одно простое число. В качестве
ривались, например, значения
раметром  . Для этого параметра была получена оценка   48 / 77    0.6237   , а
позднее:   38 / 61    0.6230   . Существует недоказанное предположение Ле-

2
p  1  180 2127  1 [1], невозможно, так как
жандра [5] о том, что   0.5 (   x ). Ответы на указанные вопросы могут быть получены с помощью анализа большой таблицы
простых чисел на вычислительной машине.
В качестве примера рассмотрим вопрос
о справедливости предположения Дебова
(Desboves) [1, с.20]: на интервале от n 2 до
неизвестен его номер.
Доказательства Адамара и Вале–
Пуссена были основаны на рассмотрении дзета-функции Римана. Идея такого подхода к
проблеме распределения простых чисел была
дана Б. Риманом в 1859 г. Работы этих авторов показали, что функция (3), предложенная
Гауссом, дает более точное приближение к
 x , чем функция
f x  
x
.
ln x 
(10)
При решении проблемы распределения
простых чисел часто выясняют, на каком интервале от x до x   существует хотя бы
которое было доказано в 1896 г. независимо
друг от друга Адамаром и Вале–Пуссеном.
Представление о скорости приближения  к
единице дают полученные значения для номеров i  15001 и i  30001 ‫ ׃‬1.0992 и
1.09314 соответственно. К сожалению, вычислить значение  для известного большого
простого
числа,
например

(9)
n  12
находится не менее двух простых чисел. Имея таблицу простых чисел, легко убедиться в справедливости предположения и,
кроме того, выяснить, как часто появляется
число простых чисел 2 и более, находящихся
на указанном интервале. При числе простых
чисел im  30001 максимально возможное
для проверки значение n  588
Программа выяснила, что два простых
числа находятся лишь на двух интервалах с
n  3, 4 . Далее встречаются числа 3,4,5 с постепенным (не монотонным, как обычно) увеличением числа простых чисел. Расширение
интервала, равного 2n  1 , и увеличение простых чисел на них не оставляют сомнений в
справедливости предположения. В этих вы-
(7)
Отличие функции (3) от функции (7)
x
 1  x   
ds
x

 li  x   f  x  , (8)
ln s  ln  x 
2
анализировалось многими математиками. Последователями И.М. Виноградова методом
тригонометрических сумм получены оценки
модуля 1 x , содержащие функции с параметрами. Для существующих таблиц простых чисел обычно величина 1 x положительна. Этот факт следует из разложения
17
Е. Л. Тарунин
числительных экспериментах обнаружился
любопытный факт – чаще всего (14 раз, что
составляет 2.37% от общего числа рассмотренных интервалов) на интервалах Дебова
размещается 74 простых числа. То есть интервалы с меньшим числом простых чисел
встречались реже. Например, меньшее число
простых чисел равное двум, встретилось лишь
три раза при n =2, 3, 5. Максимальное число
простых чисел на интервалах Дебова при
n  590 равнялось 98.
ждение называют тернарной проблемой Гольдбаха, а второе – бинарной проблемой Гольдбаха или Гольдбаха–Эйлера. Из справедливости бинарной проблемы автоматически
следует справедливость тернарной проблемы.
Действительно, если четное число есть сумма
двух простых чисел (m=p1+p2), то, добавляя к
каждому четному числу 3, можно получить
все нечетные числа, равные сумме трех простых чисел (m+3=p1+p2+3).
Коротко осветим историю вопроса. В
1923 г. Харди и Литлвуд показали, что из
справедливости обобщенной гипотезы Римана
[3] следует справедливость гипотезы Гольдбаха для всех достаточно больших значений
нечетных чисел m  m . В 1937 г. Виноградов, не используя гипотезу Римана, а используя метод тригонометрических сумм, доказал
существование m . Последователь Виногра-
2. Разность между простыми числами
Существует много вопросов, касающихся разности между простыми числами [3].
Среди них выделяются вопросы поведения
d i  N i  1  N i  – разности между соседними простыми числами; количество и распределение близнецов, для которых d i  2 ,
или, более общо, пар с разностью равной
2k k  1, 2, 3,.... С помощью гипотезы Римана доказано, что
дова дал оценку для m , содержащую 6 млн
цифр. В 1989 г. Ванг и Чен опустили нижнюю
грань до числа, содержащего 43 тысячи цифр.
В 1997 г. четверо ученых показали, что справедливость гипотезы Римана обусловливает
справедливость тернарной проблемы Гольдбаха для m  m*  10 20 . Если предыдущие
оценки практически не позволяли убедиться в
справедливости гипотезы до указанных значений m , то последняя оценка дает такую
надежду при использовании современной вычислительной техники. Как следует из сказанного, последняя оценка получена из предположения справедливости гипотезы Римана,
которая не доказана. Мы не будем останавливаться на обсуждении этой гипотезы, сформулированной в 1859 г. Подтвердим лишь ее
сложность и значимость историческими фактами. Однажды Гильберта спросили, чем бы
он поинтересовался в первую очередь, если
бы проспал 500 лет. Гильберт сказал‫" ׃‬Я бы
спросил‫" ׃‬Доказана ли гипотеза Римана?" Гипотеза входит в список 7 проблем тысячелетия, за решение каждой из которых Математический институт Кембриджа обещал приз в
1 млн долларов США.
Обещанный миллион многим не дает
покоя. Поэтому не случайно в Интернете
можно найти объявления о якобы доказательстве гипотезы Римана. Так в Интернете (Компьюлента, 11 июня 2004 г.) помещена заметка
В.Парамонова с интригующим заголовком
"Найдено доказательство гипотезы Римана".
d i  1 i   N i  ln N i .
(11)
Эвристические рассуждения показывают, что, вероятно, справедлива оценка
d i   2 i   ln 2  N i .
(12)
Согласно [3, 6] лучшей к 1983 г. являлась оценка, полученная M.N.Huxley (1973)
по методу большого решета:
d i   3 i   N  i ,  
7
 .
12
Вычислительные эксперименты дают
поправочные коэффициенты к функциям
(11)–(13) и позволяют расставить соответствующие функции в следующем порядке
(при достаточно больших N ):
1 (i )  3 (i )  2 (i ) ,
из которого следует, что лучшей является
оценка (12), а не (13).
3. Проблема Гольдбаха
В 1742 г. прусский математик Кристиан
Гольдбах в письме Л. Эйлеру высказал предположение‫" ׃‬Каждое нечетное число больше 5
можно представить в виде суммы трех простых чисел". По этому поводу Эйлер выдвинул более сильную гипотезу‫" ׃‬Каждое четное
число больше двух можно представить в виде
суммы двух простых чисел". Первое утвер-
18
Возможности вычислительных методов в проблемах теории чисел
Профессор математики из университета Пердью Луи де Бранж де Бурсил утверждает, что
нашел доказательство гипотезы Римана. Его
изыскания на 23 листах выложены в свободном доступе [12] и ждут опровержения.
Бинарная проблема Гольдбаха далека от
решения. Есть работы [1], в которых доказывается, что если есть доля чисел, не представимых в сумме двух простых, то она мала.
Есть работы, в которых увеличено число слагаемых до 6 и более. Утверждается, что на
июль 2008 г. бинарная проблема проверена до
n  1.2  1018 [8], а на сайте "Элементы‫ ׃‬Проблема Гольдбаха" указано значительно меньшее значение m  100000. .
В сентябре 2004 г. журналистка Татьяна
Нечапайко разместила в Интернете заметку о
том, что якобы 75-летний минский математик,
кандидат наук Виктор Карпов решил "задачу"
Гольдбаха–Эйлера [8]. Карпов отнес свою работу в Белорусский государственный университет, педагогический университет и АН. По
его словам, от него "отмахиваются, как от
назойливой мухи". В заметке описана запутанная история с обещанием в 2000 г. выплатить 1 млн долларов за решение проблемы (по
непроверенным сведениям обещание действовало лишь до 2002 г.). Там же отмечено,
что греческий ученый и писатель Апостолос
Доксиядис проблему не решил, но, написав
роман об этой проблеме [9], приобрел всемирную славу. После того как книга стала
бестселлером, английское издательство "Faber
and Faber" и американское "Bloomsbery" и пообещали приз за решение проблемы. Журналисты "Комсомольской правды в Беларуси"
пытаются выяснить судьбу обещанного миллиона и ждут подтверждения правильности
доказательства В.Карпова.
for i1:=2 to im do begin if N[i1] > 0.5*nc
then goto 12;
if nc mod N[i1] =0 then begin j:=j+2; goto
10; end; end;
12: N[i+1]:=nc; end; {for i}
Использованный алгоритм на ПК Toshiba требует около 4 сек для получения 30 тысяч простых чисел при программе в Delphi.
Одна и та же программа на Паскале в Delphi
работает примерно в 2.5 раза быстрее. Для
выявления зависимости времени счета от последнего номера простого числа im может
быть использована формула:
(14)
t  5  10 9  im 2 сек.
Из нее, в частности, следует, что для
получения 1 млн простых чисел требуется
около 1.5 часа счета.
Отметим особенности программы. Целые
переменные и массив для целых чисел предназначены для хранения longint. Первый кандидат
на простое число получается добавлением к
предыдущему числу j=2. Если полученное число делится на младшие предыдущие простые
числа, добавка к новому кандидату на испытание увеличивается еще на 2 оператором j:=j+2.
Процесс проверки повторяется до тех пор, пока
не выяснятся, что нет делителей у числа nc. И
если делителей нет, происходит запись нового
простого числа N[i+1]:=nc. Оператор в третьей
строке служит для сокращения излишних проверок. Если его убрать, время счета увеличивается примерно вдвое.
Выяснено, что простое число становится более чем в 10 раз больше своего номера
при i  6473 , так как
N6472  64717, N6473  64747.
Среди последних цифр простых чисел
(1, 3, 7, 9) нет выделенных, их доли при
i  200 составляют  (25  1)%.
4. Вычисление таблицы простых чисел
5. Анализ аппроксимаций
распределения простых чисел
Наиболее древним способом получения
простых чисел является использование алгоритма Эратосфена. Работа с алгоритмом на
вычислительных машинах требует первоначального массива гораздо большего размера,
чем размер окончательного массива N i  , поэтому вызывает определенные трудности.
В использованном алгоритме после задания первых четырех простых чисел следующие
числа вычислялись по программе (Paskal):
for i:=4 to im do begin j:=2;
10: nc:=N[i]+j; {candidate for next number}
Опишем вначале отличие функций
FG ( x), f ( x) , описанных в (3) и (4) соответственно. Определенный интеграл в функции
Гаусса вычислялся по трем приближенным
формулам прямоугольников:
S1 i  
S i  
N i 1
1
 lnk  1 ,
k 2
 lnk  0.5 ,
k 2
19
1
N i 1
S 2 i  
N i 1
1
 lnk .
k 2
Е. Л. Тарунин
В силу монотонного убывания подынтегральной функции для этих сумм выполняются неравенства S1  S  S 2 .
Уже при N =1000 относительная погрешность интегралов менее 0.08% и с ростом
N она убывает. Приближенные значения интегралов с избытком и недостатком удобно
использовать для получения гарантированных
оценок сверху и снизу.
Согласно теории значения рассматриваемых функций при x   должны сближаться,
а их отношения к номеру простого числа стремиться к 1. При конечных x обычно уверенно
выполняется
неравенство
FG x   f x  .
Представление о поведении этих функций
дает табл. 1. Напомним, что номер простого
числа обозначаем буквой i , а соответствующее значение простого числа – буквой N i  .
Таблица 1
i
N i 
1000
4000
10000
20000
30000
7919
37813
104729
224737
350377
1.00259
1.001837
0.91188
0.91482
FG N i  / i
1.0163
1.00625
f N i  / i
0.8821
0.89686
Как видно, функция FG  N  дает завышенные значения номера простого числа, а
функции f N  – заниженные. Их сближение
происходит медленно. Относительное превышение FG  N  над номером простого числа
может быть охарактеризовано также отношением  i   100  S i  / i  1 .
Эта величина с ростом номера числа в
целом монотонно убывает (нарушения монотонности характерны для значений i  100 ).
При i  1000   1.63 ,  становится
меньше 1% при i  2100 , а при i  15000
  0.33%.
Указанные зависимости характерны и
для бо'льших номеров простых чисел. В [4]
приведена теорема Е.Мейсселя, позволяющая
вычислять  x далеко за пределам обычных
таблиц простых чисел, и в качестве примера
приведено значение  10 9  i  50847478 .
Аргумент x в этом примере может отличаться от соответствующего значения N на расстояние между простыми числами, но при достаточно больших x это слабо сказывается на
точности оценок. Использование этого примера позволяет оценить коэффициенты
табл.1, указанные в третье и четвертой строках‫ ׃‬1.000035, 0.94901 соответственно.
Заметим, что вычисление определенно9
го интеграла в формуле Гаусса при N  10
требует немалых затрат машинного времени
(при N  10 9 около 10 мин). Для сокращения
затрат можно использовать идеи распараллеливания алгоритма.
1.0038
6
0.9069
3
Перейдем к нахождению коэффициентов a, A , которые удовлетворяют неравенствам Чебышева (4) на реальных таблицах
N i  .
В этом случае алгоритм заключался в
циклическом переборе всех табличных значений и нахождении коэффициентов, которые
удовлетворяют неравенствам( 10  i  im  1 )‫׃‬
a  N i / lnN i  i  A  N i / lnN i. (15)
Перед началом проверки задавались
значения коэффициентов, которые заведомо
изменятся, и отслеживался номер простого
числа, при котором произошла последняя
корректировка. При im  30000 оба коэффициента оказались бо'льшими единицы‫׃‬
a  1.0926 i  26963, A  1.25506 i  30. .
В скобках указан номер последней корректировки. Заметим, что если цикл проверки
начать с i  31 , то величина А уменьшится до
А=1.22358 ( i =46). Значение a  1 подтверждает факт заниженной величины f ( x) при
малых и умеренных значениях x . Для значительно бо'льших x коэффициент a , возможно, станет меньше 1. Для проверки полученных значений в цикле по i вычислялись относительные отклонения функции‫׃‬
 i  k  N i / ln N i  / i  1
(16)
и по ним вычислялись следующие интегральные
характеристики относительного отклонения‫׃‬
 
 im 
E1  max  i , E11  min  i , E2     i  / im  i0,
 i0 
(17)
20
Возможности вычислительных методов в проблемах теории чисел
N i  338473.
 im

 im

E 21     i  / im  i 0, E3     i2  / im  i 0.
 i0

 i0

По этим характеристикам можно судить
об отклонении функции f ( x) , умноженной
меру i  29080
на коэффициент k , от реальной зависимости
N i  (точнее, от обратной функции iN  ).
При k  1 анализируется отклонение функции f x   x / ln x  . В соответствии с расче-
При k  0.8 значения функции интегрального логарифма лежат ниже значений i
от 20 до 30000.
Кроме характеристик относительных отклонений рассматривались и соответствующие
характеристики абсолютных отклонений. В
табл. 2 приведены характеристики отклонений
для трех значений множителя k .
среднеквадратичной ошибки E3  0.00491 
достигается
при
E3  0.00587 ).
тами при im  30000 и k  1 E1  0 (все
относительные отклонения отрицательны).
Положительные отклонения  i появляются
при значениях k  a . При значениях k  A
исчезают все отрицательные отклонения  i
i
k  0.995
(при
k 1
Таблица 2
( E11  0 ). При значении множителя k в интервале от a до A есть отклонения
Минимум
разных
знаков. При k  1.104 достигается минимум
среднеквадратичной ошибки относительного
отклонения E3  0.64  10 2 % (при k  1 эта
величина в 8.5 раз больше); при этом максимальное отклонение E1  1.04% соответствует i  26963 , а минимальное значение
E2  12.04% соответствует i  30 .
Перейдем теперь к анализу функции,
определенной через интегральный логарифм
(3). В этом случае коэффициенты a , A определялись из неравенств
k
Е1
Е11
E2
E21
E3
1
0.998
0.995
94.4(26217)
57.6(26217)
23.4
5.32(45)
2.38(29080)
-102.5
48.6
27.6
48.6
27.6
31.8
0.295
0.167
0.237
-26.7
Экстремальные значения отклонений соответствуют другим номерам простых чисел
для относительных отклонений. По величине
E11 заключаем, что смена знака отклонения
действительно происходит при множителе
k  0.998 . Минимум среднеквадратичной
ошибки не соответствует значению k  0.95 ,
которое определено для относительной ошибки.
Все это указывает на то, что следует использовать различные критерии близости функций к
реальным табличным данным.
Рассмотрим аппроксимационные свойства функции с параметром b ‫׃‬
a  k  S 2 i   iN i, A  k  S1 i   iN i. (18)
f  x, b  
В этих формулах S1 – приближенное
значение интеграла в функции Гаусса с недостатком, а S 2 – значение с избытком. Неравенства (18) используются для поиска ‫״‬оптимальных‫ ״‬коэффициентов для функции Гаусса, которая считается наиболее приближенной
к  x . Характеристиками отклонения были
x
.
ln x   b
(19)
Очевидно, что положительные значения
параметра b  0 увеличивают значения
функции, а отрицательные – уменьшают. Значение b  1 рекомендовал Чебышев, значение
b  1.08366 – Лагранж. Для построения
функций, ограничивающих  x сверху и
снизу, предложены [4] значения 2 и –4. Функция (19) дает существенно завышенные значения при b  1 и малых x . Так, при b  1
относительное превышение менее 1% характерно лишь для простых чисел с номерами
те же величины (17), что и для функции f x  .
При значении коэф-фициента k  1 все относительные откло-нения положительны (интегральный лога-рифм больше всех iN  в интервале от i  10 до i  30000 ); максимальное относительное отклонение достигается
при i  11 N 11  31 . Уменьшение коэффициента ( k  0.9986 ) приводит к появлению
отрицательных значений отклонений. Первое
отрицательное отклонение соответствует но-
i  116 ( N i  641).
Результаты анализа зависимости (19)
для различных значений параметра b в интервале от 2 до –4 представлены в табл. 3 ( i
изменялось от 30 до 30001).
21
Е. Л. Тарунин
Таблица 3
b
E1 i1 %
2.0
1.5
1.115
1.10
44.04(31)
22.50(31)
9.857(31)
9.417(31)
1.08633
9.019(31)
1.0755
8.706(31)
1.05
7.975(31)
1.00
0.5
0
-0.5
-1.0
-2.0
-4.0
6.571(31)
0
0
0
0
0
0
E11 i 2 %
0
0
0
0.1330(2688)
0.2646(2688)
0.4044(2688)
0.6853(2688)
-1.3807(258)
-10.90(30)
-20.42(30)
-27.94(30)
-34.23(30)
-44.01930)
-56.84(30)
В табл. 3 указаны характеристики относительных ошибок в %. Определение этих
величин приведены в (17). При значениях параметра b  1.115 E11  0 и, следовательно,
значения функции (19) лежат выше значений
iN  в таблице (при этом, естественно,
E2 %
E21 %
9.703
4.211
0.3505
0.2060
9.703
4.211
0.3505
0.2063
0.05708
0.02481
0.002763
0.002142
0.07466
0.09194
0.001748
-0.02913
0.10775
0.001636
-0.2726
0.3114
0.002162
-0.7466
-5.245
-9.3489
-13.107
-16.56
-22.70
-32.60
0.7630
5.245
9.3489
13.107
16.56
22.70
32.60
0.004489
0.03045
0.005430
0.07611
0.09615
0.1317
0.1890
E3 %
ровки формулы типа Лагранжа согласно анализу Чебышева логично переписать ее в виде
f L x  
x
.
(20)
ln x   1  dbx 
В формуле Чебышева db  0 , а в формуле Лагранжа db  const  0.08633 (уточненное нами значение db =0.0755). Если потребовать от функции dbx  стремления к нулю при
достаточно больших значениях x , то не будет
E2  E21 ).
При b  0.5 E1  0 и, следовательно, значения функции (19) меньше значений в таблице простых чисел.
В довольно узком интервале b от 1.0 до
1.1 функция (19) дает как завышающие, так и
занижающие значения. Напомним, что в этом
интервале находится значение b  1 , указанное Чебышевым для больших значений N , а
также значение b  1.08633 , предложенное
Лагранжем. Для рассмотренного интервала
простых чисел с номерами от i0  30 до
im  30000 формула Лагранжа точнее формулы Чебышева (отношение соответствующих среднеквадратичных относительных
ошибок E3  2.56). Это не означает, что Чебышев неправ, так как согласно его утверждениям преимущества его формулы скажутся за пределами рассмотренных значений
простых чисел. Вызывает восхищение то, что
Лагранж сумел найти значение параметра b ,
близкое к оптимальному. Расчеты показали,
что
оптимальное
значение
параметра
b  1.0755 (при этом величина E3 уменьшается по сравнению с величиной в варианте
Лагранжа примерно на 6.8%). Для корректи-
противоречия, о котором упоминал Чебышев.
Однако подобрать хорошую зависимость dbx 
с упомянутым свойством непросто.
Вычисления
показали,
что
при
db  db  0.0755 формула (19) дает завышенные значения при i  5000 , а при
i  6000 – заниженные; в районе i  5060 появляются как завышенные, так и заниженные
значения. Этот анализ подсказал уменьшить
db при малых i и слегка увеличить при больших
i . Вычислительный эксперимент со значениями
dbt   0.0755
i
при i  5060 ,
5060
(21)
0.0135i  5060
dbi   0.0755 
при i  5060
im  5060
уменьшил среднеквадратичное отклонение
E3 на 18%. При этом максимальные отклонения (не относительные) укладываются в интервал от –22.85 до 37.56 (среднее отклонение по
модулю менее 7.04); при постоянном значении
22
Возможности вычислительных методов в проблемах теории чисел
статье i0  30 ) и с критериями погрешности
(относительными или абсолютными).
db этот интервал шире (от –43.91 до 43.91 при
среднем отклонении по модулю 12.23).
В предложенной формуле (21) величина
db увеличивалась с ростом номера простого
числа. Казалось бы, это противоречит требованию согласования с замечанием Чебышева.
Противоречия нет. Просто значения, при которых, по-видимому, следует требовать
асимптотического стремления
(22)
lim x db  0 ,
велики. Сам Чебышев указывал [13], что отличия формулы Лагранжа от его формулы
будут значимыми лишь при значениях x  10 7
(примерное
равенство
при
6
x  1.247646  10 ). Указанным значениям x
соответствуют номера i 90 000. Отметим,
что способ получения оценок [13] требует
проверки, так как простой анализ показывает,
что разность между аппроксимацией Лагранжа и Чебышева ( db  0 ) монотонно убывает,
достигая относительной разности менее 1%
при x  16650 .
Подведем итоги выполненных оценок
функций, полагая, что  x – реальная зависимость для номера простого числа, которой соответствует iN  в таблице. Полученные оценки трех функций сверху и снизу имеют вид
6. Расстояние между простыми числами
Расстояние между простыми числами
d i  N i  1  N i  ведет себя не монотонно,
но в целом расстояние d i , усредненное по
большому интервалу i , и максимальное d m
имеют тенденцию к увеличению. Известно
несколько формул оценки сверху максимального расстояния:
(27)
d m  F1  N i  ln N i.
d m  F2  ln 2  N i .
7
d m  F3  N  i ,  
 .
12
(29)
Первая оценка (27) получена в работе
[6] с использованием гипотезы Римана. Вторая оценка (28) получена посредством эвристических (вероятностных) рассуждений.
Третья оценка (29) получена Huxley M.N. в
1973 г.
Эти формулы были проверены на интервале по номеру простого числа от 6 до
30 000. Оказалось, что все формулы дают завышенные значения (правые части неравенств
могут быть уменьшены). В результате анализа
определены поправочные коэффициенты для
первой и второй формул и величина  для
третьей формулы. Кроме того, указаны значения номеров i1 ,i2 , в которых производилась
последняя поправка для выполнения неравенств.
На
рассмотренном
интервале
d m  86  N 14358  N 14357 . В формуле
(27) правая часть может быть умножена на
k1  0.43253i1  6 , а в формуле (28) – на
x
x
, (23)
  x   1.25506
ln x 
ln x 
(24)
0.80  li x   x  0.9987  li x ,
x
x
.
(25)
  x  
ln x   0.5
ln x   1.115
1.0926
Эти функции ближе к табличным значениям iN  по величине относительной
ошибки при следующих "оптимальных" значениях коэффициентов‫׃‬
 x  
(28)
k 2  0.67155i2  3385 . С учетом указанного в работе [4] значения d  112 для
N  370261 и N  d коэффициент для вто-
x
x
 0.995  1ix   1.104
. (26)
ln x   1.115
ln x 
Функции указаны в порядке возрастания относительной среднеквадратичной погрешности со значениями
рой формулы следует увеличить примерно на
1.5% до значения k 2  0.6813 . Третья формула остается справедливой при   0 . Кроме того, для третьей формулы был найден
меньший
показатель
вместо
‫׃‬
7 / 12  0.583(3)   0.57144404 (последняя
E3  0.0016, 0.0049, 0.0064%
соответственно.
Указанные коэффициенты ждут своего
уточнения за счет расширения таблицы простых чисел. При этом следует определиться с
начальным значением интервала i0 (в данной
поправка была выполнена для i  9 ).
23
Е. Л. Тарунин
После корректировок было выяснено,
что эти ограничения при x  130 удовлетворяют неравенству
(30)
F1  F3  F2 .
Например, между близнецами (71, 73) и
ближайшими к ним (101, 103) находится 4 простых числа (расстояние между этими близнецами N = 101–73 =28).
В то же время между близнецами (101,
103) и следующими (107, 109) нет ни одного
простого числа. В 1959 г. была опубликована
таблица близнецов в пределах 1.1 млн [2]. В
работе [6] есть пример для близнецов‫׃‬
x  8 004 119, x  2 , а в [1] – для еще более
далеких‫ ׃‬x  1000061081, x  2 . Есть основание предполагать, что близнецов бесчисленное множество, но это до сих пор не удалось доказать. На интервале до i  30001 доля близнецов составляет  11.4% , с ростом
интервала эта доля убывает.
Программа, вычисляющая число обобщенных близнецов (31), на выбранном интервале от i1 до i2 имела простой вид‫׃‬
Отсюда следует, что наиболее жесткие
ограничения для максимального расстояния
между простыми числами d m дает вторая
формула, а не третья, как утверждалось в [6].
При x  1000 , например, F3 / F1  0.457 , а
F2 / F1  0.057.
Гаусс обнаружил [1], что 26 379-я сотня
не содержит простых чисел. Понятно, что появление такой сотни требует выполнения неравенства d m  101 (в худшем случае 200). Расчет показывает, что такая сотня обязательно
найдется в интервале от 212 000 до 31 247 000.
Сотня,
найденная
Гауссом
(2 637 8012 637 899), попадает в указанный интервал.
Простая заниженная оценка для d m
следует из известного и легко доказываемого
факта [6] об отсутствии простых чисел на интервале n!2, n!n . Отсюда следует, что
for k:=1 to km do nn[k]:=0;
for i:=i1 to i2 do for k:=1 to km do
if (N[i+1]-N[i])=2*k then nn[k]:=nn[k]+1;
for k:=1 to km begin s:=100*nn[k]/(i2-i1);
writel n(k:8,s:9:6); end;
d m  n  2 . Уже при n  4 реальное d m  4 ,
В табл. 4 указан процент чисел с разницей равной 2k на первой тысяче простых чисел (i1=1, i2=1000).
а по оценке d m  2 .
7. Распределение близнецов
Близнецами (twins) называют соседние
простые числа, отличающиеся на 2. Обобщением близнецов являются соседние простые
числа, отличающиеся на число, кратное двум:
N i  1  N i  2  k , k  1, 2,... (31)
Расстояние между близнецами с ростом
номера простого числа в целом увеличивается,
но с нарушением монотонности.
Таблица 4
2k
nnk  
2
4
6
8
10
12
14
17.417
16.917
24.424
8.408
10.010
7.507
4.294
1
2
3
4
5
6
14.538
12.337
11.437
11.637
11.337
10.837
Таблица 5
j
nn1%
Заметно нарушение монотонного убывания числа близнецов. Возникает вопрос‫׃‬
какой выбрать интервал по номерам чисел
Любопытно, что оказалось очень много
пар для k  3 . Распределение близнецов
( k  1 ) на последовательных интервалах с
длиной i  3000 представлено в табл. 5 ( j
– номер интервала).
24
Возможности вычислительных методов в проблемах теории чисел
по причине повышенного числа таких пар в
табл. 4.
i  i 2  i1 , при котором число близнецов
будет монотонной функцией от номера интервала?
В
табл.
6
для
интервала
i  const  5000 представлены доли nnk 
для k  1 и k  3 (значение k  3 выбрано
Таблица 6
j
nn1%
nn3%
1
2
3
4
5
6
13.602
11.802
11.302
10.722
10.582
10.462
21.504
18.743
18.083
17.843
17.283
17.283
Как следует из табл. 6, число близнецов
монотонно убывает с ростом номера интервала j ; однако при k  3 значения пар на 5 и 6
интервалах совпали.
Таблица 7
В табл. 7 приведено среднее число пар в
процентном отношении к полному интервалу
простых чисел до i =300 000‫׃‬
k 1
k 2
k 3
k 4
k 5
k 6
11.410
11.323
18.453
7.530
9.366
10.240
Доля следующих пар ( k  7 ) менее 6%,
но при расширении интервала их доля будет
естественно возрастать.
Перейдем к исследованию расстояния
между ближайшими (соседними) близнецами.
Дадим пояснения к соответствующей программе. Пусть ( i11, i12 ) – номера простых
чисел для близнецов ( N i12  N i11  2 ).
Начальные значения первых номеров близнецов равнялись i11  5, i12  6 (это соответствовало паре 11,13). Далее программа в
цикле отыскивала номера простых чисел для
ближайших близнецов ( i 21, i 22 ) и вычисляла расстояние между ними по формуле
(32)
L  2k  N i21  N i12 .
При этом в элементах массива nnk 
суммировались случаи с соответствующим
значением k в формуле (31). Значения всех
элементов этого массива позволяло вычислить и среднее расстояние между близнецами:
L  SL 
2 k  nnk 
Sk
, Sk   nnk . (33)
Верхняя граница суммирования в (33)
определялась экспериментально. Признаком
правильной верхней границы служила неизменность суммы Sk при существенном расширении верхней границы. В нашем случае
достаточно было выбрать k  318 . Этому
значению соответствовало наибольшее за-
11
фиксированное расстояние между близнецами: L  2k  636.
Анализ массива nnk  позволил выяснить, что для проанализированного числа
простых чисел ( im  30001 ) найдено 3417
близнецов (без учета тех, для которых
x  13 ). Любопытно, что чаще всего встречались
расстояния,
соответствующие
k  14 L  28, nn14  278 . При этом
наименьшее расстояние, равное L  4 , встретилось лишь 76 раз. На втором и третьем месте по численности оказались значения
nn20  229 и nn5  222 для L  40 и
L  10 соответственно.
Число близнецов на интервалах с шагом
i  5000 простых чисел монотонно уменьшалось от 679 на первом интервале до 523 на
последнем, шестом интервале. При этом
среднее значение расстояния между близнецами L монотонно возрастало‫ ׃‬на первом
интервале ( i  5001) оно равнялось  69.6, а
на полном интервале ( i  30001)  99.7 .
Наиболее интересный результат заключается в том, что многие расстояния между
близнецами отсутствуют. Ненулевыми элементами массива nnk  оказались лишь те,
для которых k  2  3n n  0,1, 2,... и соответственно L  4  6n . Формула
k  2  3n n  0,1, 2,...
Е. Л. Тарунин
была справедлива до k  188 . Затем стали
появляться значения с L  6 j  j  2 . Таким образом, для определения расстояния
между близнецами применяется формула
L  4  6nj с целочисленными значениями
n, j .
for k:=0 to km do if i2=k then
nn[k]:=nn[k]+1;
goto 30; end;
end; {for i1}
20: end; {for i} writeln(‘ impssible for
M=’,M:9); goto 40;
30: end; {for j} writeln(‘ all was found’);
writeln(‘ jm,M=’,jm:9,M:9);
40: {end program}
8. Бинарная проблема
Гольдбаха–Эйлера
Дадим пояснения к программе. Задание
величины MM определяет наибольшее значение числа M, для которого подбираются простые числа, удовлетворяющие гипотезе Гольдбаха (34). Верхняя граница цикла jm по j вычисляется по формуле исходя из условия, что
образование числа M идет по формуле
M:=20+2j. В цикле по i вычисляется величина
разности dm=M–N[i]. Далее организуется
После создания таблицы простых чисел
анализировалось выполнение бинарной проблемы
M  N1  N 2  N i 2  N i3 (34)
для всех четных чисел M . Известно, что это
равенство может выполняться разными способами.
Например,
8=1+7=3+5,
10=3+7=5+5=7+3. Для однозначности представления полагалось, что N[i2] < N[i3] (или,
что то же, i2 < i3). Это означает, что первое
простое число меньше второго. Во всех проверенных случаях алгоритм обнаруживал выполнимость строгого неравенства N[i2] <
N[i3]. Это не означает, что не было случаев
равенства N[i2] = N[i3]. Просто алгоритм
начинал проверку равенства с постепенным
увеличением первого простого числа и прекращал работу при нахождении подходящей
пары чисел. Число вариантов с равенством
M  2 N i 2 подсчитывалось суммированием
ситуаций, для которых выполнено неравенство 2 N i  MM .
Отметим, что равенство Гольдбаха (34)
для некоторых M может выполняться тремя
и более способами даже при условии
N1  N 2 . Например, 22=3+19=5+7=11+11,
24=1+23=5+19= 7+17=11+13. Эти случаи заслуживают особого внимания. Они обеспечивают некоторый "запас" выполнимости гипотезы (34).
Программа нахождения номеров простых чисел i2, i3, удовлетворяющих равенству
(34), имела вид
MM:=400000; {maxM} jm:= MM shr 1-10;
for j:=1 to jm do begin M:=20+2*j;
for i:=0 to im do begin dm:=M-N[i];
i0:=i+Round(dm/sqr(ln(M))); if i0 > im
then goto 20;
for i1:=i0 to im do begin
if N[i1] > dm then goto 20;
if dm=N[i1] then begin i2:=I; i3:=i1;
26
Возможности вычислительных методов в проблемах теории чисел
цикл по i1 с начального значения i0 (эта величина сокращает перебор). Если выполняется
одно из неравенств‫ ׃‬i0 > im, N[i1] > dm происходит выход из поиска программы на печать
информации о том, что для равенства (34) не
найдена соответствующая пара простых чисел. При выполнении равенства dm=N[i1] (а
следовательно (34)) и нахождения номеров i2,
i3 может быть вставлена программа обработки полученных значений. При завершении
цикла по j печатается информация об удачном
завершении прог-раммы и выдаются параметры, для которых выполнена проверка равенства (34).
В качестве верхней границы проверяемых чисел ММ можно брать значение почти
вдвое большее N[im]. Так, при N[30001] =
350381 программа находит пары чисел до
значения ММ=699890  1.997  N[30001] (для
меньших im соответствующий коэффициент
при N im – меньше). Программа работает
Заметим, что если задано завышенное
число ММ, программа укажет то число, для
которого не хватает нужной пары чисел из
полученной таблицы простых чисел. После
этого следует задать число ММ, меньшее на 2,
и повторить счет. Время счета растет при увеличении ММ быстрее, чем ММ  lnMM  , и
при ММ=500000 составляет примерно 54 сек.
Алгоритм поиска чисел в равенстве (34) может быть оптимизирован. Можно, например,
использовать информацию о предыдущей паре простых чисел.
Для анализа выполнимости гипотезы
суммировались числа случаев с i2 =k для k,
меняющегося от 0 до назначаемого номера
km. Соответствующие суммы накапливались в
массиве nn[k]. Значение nn[0] указывает число
случаев, когда первое число в гипотезе равно
1. Значение nn[1] естественно оказалось равным нулю, так как N[1]=2 – единственное
четное простое число. Представление о распределении nn[k] при различных значения
ММ дает табл. 8.
быстро при M  N im и резко замедляется
при M  N im.
Таблица 8
MM
N1=N2
nn[0]%
nn[2]%
nn[3]%
nn[4]%
nn[5]%
1000
2000
3000
9.59%
8.99
8.665
17.98
15.00
10.00
15.82
13.29
8.86
13.69
11.59
7.73
10.23
8.82
5.88
8.50
7.39
4.92
В первом столбце таблицы указано число
рассмотренных чисел, во втором – процент случаев, когда простые числа в формуле (34) совпадают. Число таких случаев равно числу простых чисел, удовлетворяющих неравенству
2N[i] < MM. Заметим, что в программе всегда
находились простые числа с номерами i2 < i3.
Видно, что доля чисел с малыми значениями
номеров простых чисел монотонно убывает.
Обнаружена интересная особенность
для чисел nn[k] (k<14)‫ ׃‬при значениях ММ >
N[im] эти числа практически не изменялись,
так как равенства (34) удовлетворялись при
бо'льших номерах первого числа. Это свидетельствует о том, что для выполнения равенства (34) имеется большой запас, не говоря
уже о выполнимости его при одинаковых значениях простых чисел. Проиллюстрируем
разницу в выполнении равенства (34) для
двух случаев: М > N[im] и М < N[im]. Выберем одно и то же значение М=200 000 для
N[15001]=163847 и N[30001]=350381. В первом случае выполняется первое неравенство, а
во втором – второе. В первом случае для М
найдено представление
M=N[3845]+N[14995]=36229+163771,
а во втором случае
M=N[0]+N[17984]=1+199999
с наименьшим номером i2=0.
Видоизменение программы позволит выяснить число случаев, когда равенство (34) может быть выполнено более чем двумя способами.
Выводы
1. Приведены алгоритмы вычисления таблицы простых чисел и проверки бинарной
проблемы Гольдбаха.
2. Скорректированы оценки функций распределения простых чисел (логарифмический
интеграл, функции Чебышева и Лагранжа)
сверху и снизу.
3. Найдены приближенные значения коэффициентов, обеспечивающих минимум среднеквадратичного отклонения.
4. Предложена аппроксимация, согласующая формулу Лагранжа с асимптотикой Чебышева.
5. Уточнены формулы для оценки максимального расстояния между простыми числами, указана их относительная точность.
6. Проанализирована зависимость расстояния между близнецами и найдена формула,
27
Е. Л. Тарунин
исключающая появление "запрещенных" расстояний.
7. Проверена справедливость гипотезы
Гольдбаха до 700 тысяч и указаны возможности для анализа неоднозначного выполнения
этой гипотезы.
Напомним, что все оценки получены
для сравнительно малого (по современным
меркам) числа простых чисел im=30 001.
Предложенные алгоритмы могут быть использованы при параллельных вычислениях
на современной технике. А для начала можно
обойтись и без распараллеливания алгоритма.
Для этого следует лишь воспользоваться возможностями динамического массива для хранения расширенной таблицы простых чисел.
Следует считать результаты существенными
при числе простых чисел более 10 миллионов.
6. Математическая энциклопедия (Распределение простых чисел). М.‫ ׃‬Советская энциклопедия, 1984. Т.4
7. Интернет, Википедия, Проблема Гольдбаха, 19.11.2009.
8. Weisstein, EricW. Goldbach Conjecture (на
сайте Wolfram MathWorld).
9. Доксиадис А. Дядя Петрос и проблема
Гольдбаха. М.‫ ׃‬АСТ. 2002.
10. Петров С. Абсолютное программирование.
Рекурсия – пример типичной псевдоматематической попытки доказательства проблемы
Гольдбаха методом просеивания.
11.Интернет,
http://news.bbc.co.uk/hi/russian/sci/tech/newsi
d_3687000/3687852.stm
12. Интернет (доказательство гипотезы Римана на 23 листах)
http://riemann.narod.ru/index.html
13. Чебышев П.Л. Об определении числа
простых чисел, не превосходящих данной
величины // Изб. тр. М.‫ ׃‬Изд-во АН. 1955.
С.9–32.
14. Математический энциклопедический
словарь. М.‫ ׃‬Советская энциклопедия, 1988.
Список литературы
1. Сушкевич А.К. Теория чисел (Элементарный
курс). М‫ ׃‬Вузовская книга, 2007. 240 с.
2. Бухштаб А.А. Теория чисел: уч. пос. СПб.;
М.; Краснодар‫ ׃‬Лань, 2008. 384 с.
3. Виноградов И.М. Основы теории чисел. М.‫׃‬
Лань, 2004 (десятое издание). 176 с.
4. Трост Э. Простые числа. М.‫ ׃‬ИЛ, 1959. 136 с.
5. Прахар К. Распределение простых чисел.
М.‫ ׃‬Мир, 1967. 512 с.
Possibilities of numerical methods for problems
of number theory
E. L. Tarunin
Perm State University, Russia, 614990, Perm, Bukireva st., 15
tarunin@psu.ru; (342) 237-10-31
A theory of numbers is famed for using of analytical methods. Grate mathematicians call the theory as a pearl of mathematics. Numerical methods and computers were used in the theory only in
the second half of the 20 century. But I dare say that computers may be used more widely. The
main aim of the article is to give impulse in that direction. Classical analytical methods deal with
smooth functions but the real distribution of simple numbers is a discrete function. So it is more
suitable for computers. It was shown that computers permit to find "optimal" parameters of different functions (Chebishev, Lagrange, integral logarithm) that can describe the real distribution of
simple numbers more exactly. It was found also new properties of a distribution for so called wins.
Corrections of the results will be done by next investigators by expansion of the tested table of
simple numbers.
Key wоrds: simple numbers; numerical methods; displasment of simple numbers; goldbods
problem.
28