2014 dxdy logo

Научный форум dxdy

Математика, Физика, Computer Science, Machine Learning, LaTeX, Механика и Техника, Химия,
Биология и Медицина, Экономика и Финансовая Математика, Гуманитарные науки




Начать новую тему Ответить на тему На страницу Пред.  1 ... 19, 20, 21, 22, 23, 24, 25 ... 43  След.
 
 Re: Симметричные кортежи из последовательных простых чисел
Сообщение07.12.2023, 19:50 
Аватара пользователя


29/04/13
8067
Богородский
Ага. Кэфы $94.8$ и $96.8$. Вроде так и должно быть. Чем выше счёт, тем выше кэф. Хотя думал здесь будет $98-100$.

Интересно узнать кэфы на высоте 8е18.

 Профиль  
                  
 
 Re: Симметричные кортежи из последовательных простых чисел
Сообщение11.12.2023, 10:18 
Аватара пользователя


29/04/13
8067
Богородский
Yadryara в сообщении #1621374 писал(а):
Короткие-то можно и на PARI запрогать и сравнить.

Для нечётных сделано. Довольно долго игрался на разных высотах, пока тенденции кэфов не стали ясны:

1. Чем выше счёт, тем выше кэф.

2. Чем длиннее кортеж, тем выше кэф.

3. Чем длиннее кортеж, тем ближе соседний кэф на той же высоте счёта.

Пример для наиболее удачной высоты:

$\tikz[scale=.08]{
\fill[green!90!blue!50] (40,200) rectangle (160,210);
\draw[step=20cm] (0,200) grid +(160,40);
\draw (0,230) -- (160,230);
\draw (0,210) -- (160,210);
\node at (10,235) {\text{Length}};
\node at (10,215){\text{Total}};
\node at (10,225){\text{Height}};
\node at (10,205){\text{Coef}};
\node at (30,235){\text{1}};
\node at (30,215){\text{553288914}};
\node at (30,225){\text{$\approx$ 10 e15}};
\node at (50,235){\text{3}};
\node at (50,215){\text{8478142}};
\node at (50,225){\text{$\approx$ 10 e15}};
\node at (50,205){\text{65.3}};
\node at (70,235){\text{5}};
\node at (70,215){\text{108984}};
\node at (70,225){\text{$\approx$ 10 e15}};
\node at (70,205){\text{77.8}};
\node at (90,235){\text{7}};
\node at (90,215){\text{1305}};
\node at (90,225){\text{$\approx$ 10 e15}};
\node at (90,205){\text{83.5}};
\node at (110,235){\text{9}};
\node at (110,215){\text{7400220}};
\node at (110,225){\text{10 e15}};
\node at (110,205){\text{88.2}};
\node at (130,235){\text{11}};
\node at (130,215){\text{78024}};
\node at (130,225){\text{10 e15}};
\node at (130,205){\text{94.8}};
\node at (150,235){\text{13}};
\node at (150,215){\text{806}};
\node at (150,225){\text{10 e15}};
\node at (150,205){\text{96.8}};
}$

Для подсчёта количеств коротких цепочек до 7-ки включительно, брал интервал в 20 ярдов на высоте $4999 e12$. Для кэфа при 9-ке взял $\dfrac{7400220}{500000}\approx 14.80$ кортежей.

Прога ищет все кортежи длиной 3, 5, 7 в заданном интервале, записывает их количества, все удесятерённые кэфы, а также все количества различных диаметров. Например, самый большой найденный для 3-ки диаметр — 768.

(PARI)

Код:
{print();
l=70;kk=vector(l);skk=0;st=12;
po =precprime(4999*10^st);
poo=precprime(po-1);

forprime(p=4999*10^st,4999*10^st+20*10^9,kp++;
g2=p-po;
if(p-poo==g2*2, k3++;
if(poo-precprime(poo-1)==nextprime(p+1)-p, k5++;
po3=precprime(poo-1);p1=nextprime(p+1);
if(po3-precprime(po3-1)==nextprime(p1+1)-p1, k7++;
po4=precprime(po3-1);p2=nextprime(p1+1);di=p2-po4;
kk[di/12]++;
if(di>=216,
print(k7,"   ",po4,"    ",di)))));

poo=po;po=p);
write("10-7.txt",st,"     ",kp,"     ",k3,"     ",k5,"     ",k7);
write("10-7.txt",st,"                      ",round(kp*10/k3),"       ",round(k3*10/k5),"     ",round(k5*10/k7));
for(i=1,l,skk=skk+kk[i];
write("10-7.txt",i*12,"     ",kk[i]);
);
print();print(kp);print();
print();print(k3);print(k5);print(k7);print(skk);print();
print(kk);print();
}quit;

 Профиль  
                  
 
 Re: Симметричные кортежи из последовательных простых чисел
Сообщение19.12.2023, 06:46 
Аватара пользователя


29/04/13
8067
Богородский
Dmitriy40 в сообщении #1621385 писал(а):
Томаш Брада выкладывал отдельно базу нечётных до 1e16, там было:
n9: 7400220шт

Именно про 9-ки не знаю где он выложил. А сколько из них с паттерном
$0,18,24,48,54,60,84,90,108$ ?

 Профиль  
                  
 
 Re: Симметричные кортежи из последовательных простых чисел
Сообщение19.12.2023, 14:19 
Заслуженный участник


20/08/14
11710
Россия, Москва
Yadryara в сообщении #1622968 писал(а):
Именно про 9-ки не знаю где он выложил. А сколько из них с паттерном
$0,18,24,48,54,60,84,90,108$ ?
Где выложил не помню, не сохранилось ссылки, а поиском быстро не находится. Файл был около 23 марта 2017 с именем odds.zip или odds_1e16.zip, мегов на 100 (или больше, тогда это уже я перепаковал в меньший). Могу перевыложить в своё облако.
Указанных паттернов в нём 87327шт.

 Профиль  
                  
 
 Re: Симметричные кортежи из последовательных простых чисел
Сообщение19.12.2023, 16:59 
Аватара пользователя


29/04/13
8067
Богородский
Благодарю.

Ой, здорово, что у Вас сохранились 9-ки! Именно они — недостающее звено для статистики, ибо 7-ки и ниже в подходящем количестве находятся достаточно легко, а 11-ки и выше в общедоступной Базе имеются.

А ведь если оставить только сами начальные числа и убрать паттерны, будет гораздо меньше места занимать. Паттерны-то всегда можно определить самому, разобраться будет недолго. В таком укороченном виде был бы рад получить по почте. Ну или в облаке, как Вы решите.

 Профиль  
                  
 
 Re: Симметричные кортежи из последовательных простых чисел
Сообщение19.12.2023, 18:18 
Заслуженный участник


20/08/14
11710
Россия, Москва
Yadryara
В папке https://cloud.mail.ru/public/MwdT/asL1TZgbf файл odds_1e16.zip. Содержимое точно исходное, сам архив возможно перепакованный и переименованный, не уверен.
В паттерне вместо левого нуля стоит длина цепочки. Ну вот так у Томаша хранились цепочки в БД (а это похоже выгрузка прямо из БД).

-- 19.12.2023, 18:28 --

7.4млн цепочек, даже только начальное число, по 8 байт на число, это уже под 60МБ бинарных данных. А если число текстом, то по 18 байтов на число, 133МБ текста. Всего втрое меньше полного размера (345МБ текста с паттернами без кавычек), не вижу смысла ужимать, с текстом работать сильно проще. Хотя, если только начальное число, то можно всё сразу загрузить в вектор в PARI ...
Впрочем получить только начальное число из имеющегося csv можно простой программой на PARI:
Код:
f=fileopen("9.csv","r"); while(s=filereadstr(f), write("n9.txt", eval(strsplit(s,",")[1]));); fileclose(f);

 Профиль  
                  
 
 Re: Симметричные кортежи из последовательных простых чисел
Сообщение22.12.2023, 06:13 
Аватара пользователя


29/04/13
8067
Богородский
Благодарю. В данном случае анализировал те самые центральные кортежи в легендарной $19-252$. Для тех, кто не очень в теме, попытался изобразить их левые половинки наглядно.

$ Central\;\;\; 3,\; diameter =  \;\,12  \hspace{6.8cm}                 6$

$ Central\;\;\; 5,\; diameter = \;\, 60        \hspace{6.cm}               24-6$

$ Central\;\;\; 7,\; diameter = \;\, 72        \hspace{5.4cm}               6-24-6$

$ Central\;\;\; 9,\; diameter =  108        \hspace{4.6cm}   18-6-24-6$

$ Central\; 11,\; diameter =  168        \hspace{3.8cm}  30-18-6-24-6$

$ Central\; 13,\; diameter =  192        \hspace{3.cm}  12-30-18-6-24-6$

$ Central\; 15,\; diameter =  228        \hspace{2.2cm}  18-12-30-18-6-24-6$

$ Central\; 17,\; diameter =  240        \hspace{1.6cm}  6-18-12-30-18-6-24-6$

$ Central\; 19,\; diameter =  252        \hspace{1cm}  6-6-18-12-30-18-6-24-6$

Здесь указаны гэпы между простыми числами кортежей. Кортежи достраиваются вправо до полных, исходя из симметрии относительно центра. Пример для 9-ки:

$18-6-24-6\qquad \to \qquad 18-6-24-6-6-24-6-18$

Ну а диаметр равен сумме гэпов полного кортежа.

Таким образом, в таблице ниже — статистика именно по таким центральным кортежам.

$\tikz[scale=.08]{
\fill[green!90!blue!50] (40,200) rectangle (120,210);
\draw[step=20cm] (0,200) grid +(140,40);
\draw (0,230) -- (140,230);
\draw (0,210) -- (140,210);
\node at (10,235) {\text{Length}};
\node at (10,215){\text{Total}};
\node at (10,225){\text{Height}};
\node at (10,205){\text{Coef}};
\node at (30,235){\text{1}};
\node at (30,215){\text{27921439}};
\node at (30,225){\text{$\approx 10^{16}$}};
\node at (50,235){\text{3}};
\node at (50,215){\text{109620}};
\node at (50,225){$\approx 10^{16}$};
\node at (50,205){\text{255}};
\node at (70,235){\text{5}};
\node at (70,215){\text{221.7}};
\node at (70,225){\text{$\approx 10^{16}$}};
\node at (70,205){\text{494}};
\node at (90,235){\text{7}};
\node at (90,215){\text{2.68}};
\node at (90,225){\text{$\approx 10^{16}$}};
\node at (90,205){\text{82.7}};
\node at (110,235){\text{9}};
\node at (110,215){\text{87327}};
\node at (110,225){$10^{16}$};
\node at (110,205){\text{307}};
\node at (130,235){\text{11}};
\node at (130,215){\text{175}};
\node at (130,225){\text{$10^{16}$}};
\node at (130,205){\text{499}};
}$


В тех столбцах, где стоит символ $\approx$ , указано среднее количество соответствующих кортежей на один миллиард чисел натурального ряда. И это среднее, если брать его на нужной высоте, примерно равно среднему по всему диапазону.

В прошлый раз я неправильно выбрал высоту для таких подсчётов:

Yadryara в сообщении #1621893 писал(а):
Для подсчёта количеств коротких цепочек до 7-ки включительно, брал интервал в 20 ярдов на высоте $4999 e12$.

В этот раз брал на высоте $3572 e12$ и поблизости. Здесь учёл в общей сложности 31 ярд. И итоговые первые цифры усреднённого значения, то есть $279214...$ гораздо лучше соответствуют $\pi(10^{16})=279238...$ А на указанной высоте количество простых совпадает ещё лучше: $279240...$

Если количество кортежей, взятое для расчёта кэфа, меньше тысячи, то я считаю этот кэф не шибко надёжным и зелёным его не закрашиваю. Именно так получилось в последнем случае, где центральных 11-к в указанном диапазоне всего лишь $175$.

В дальнейшем попробую забраться повыше и оценить частотность для более длинных кортежей, включая и $19-252$. Для этой оценки необязательно брать именно центральные кортежи меньшей длины.

 Профиль  
                  
 
 Re: Симметричные кортежи из последовательных простых чисел
Сообщение22.12.2023, 15:07 
Заслуженный участник


20/08/14
11710
Россия, Москва
В прошлый раз зелёненькие цифирки выглядели понадёжнее - меньше различались.

-- 22.12.2023, 15:23 --

Для большей точности лучше взять много небольших кусочков (но не менее десятков/сотни миллионов) равномерно по всему диапазону 1e16.

 Профиль  
                  
 
 Re: Симметричные кортежи из последовательных простых чисел
Сообщение23.12.2023, 11:20 
Аватара пользователя


29/04/13
8067
Богородский
Dmitriy40 в сообщении #1623430 писал(а):
В прошлый раз зелёненькие цифирки выглядели понадёжнее - меньше различались.

В прошлый раз была куча мала — статистика сразу по всем находкам по всем паттернам. Сейчас — по конкретным. И нынешняя, особенно по 3-кам и 5-кам, надёжнее — больше выборка.

Dmitriy40 в сообщении #1623430 писал(а):
Для большей точности лучше взять много небольших кусочков (но не менее десятков/сотни миллионов) равномерно по всему диапазону 1e16.

Для большей точности чего?

Вот, поднялся на порядок выше:

$\tikz[scale=.08]{
\fill[green!90!blue!50] (40,200) rectangle (100,210);
\draw[step=20cm] (0,200) grid +(100,40);
\draw (0,230) -- (100,230);
\draw (0,210) -- (100,210);
\node at (10,235) {\text{Length}};
\node at (10,215){\text{Total}};
\node at (10,225){\text{Height}};
\node at (10,205){\text{Coef}};
\node at (30,235){\text{1}};
\node at (30,215){\text{26235078}};
\node at (30,225){\text{$\approx 10^{17}$}};
\node at (50,235){\text{3}};
\node at (50,215){\text{91699}};
\node at (50,225){$\approx 10^{17}$};
\node at (50,205){\text{286}};
\node at (70,235){\text{5}};
\node at (70,215){\text{171.5}};
\node at (70,225){\text{$\approx 10^{17}$}};
\node at (70,205){\text{535}};
\node at (90,235){\text{7}};
\node at (90,215){\text{1.72}};
\node at (90,225){\text{$\approx 10^{17}$}};
\node at (90,205){\text{99.7}};
}$

Видимо, теперь надо и вниз спуститься, на 9-ки посмотреть.

 Профиль  
                  
 
 Re: Симметричные кортежи из последовательных простых чисел
Сообщение23.12.2023, 13:10 
Заслуженный участник


20/08/14
11710
Россия, Москва
Yadryara в сообщении #1623510 писал(а):
Для большей точности чего?
Вы пытаетесь оценить площадь под неизвестной кривой (общее количество цепочек заданной длины) по одной точке (вероятность/частота в одной точке или мельчайшем интервальчике), а лучше взять несколько точек по всему интересующему диапазону. Или зафиксировать вид кривой (как например для $\pi(x)=\frac{x}{\ln x}$). Иначе вообще непонятно какую точность Вы получаете. Вот к примеру насколько точна оценка 9-ок в ваших интервальчиках? Насколько совпадает с 87327 или 7400220 после пересчёта на 1e16?

 Профиль  
                  
 
 Re: Симметричные кортежи из последовательных простых чисел
Сообщение23.12.2023, 15:00 
Аватара пользователя


29/04/13
8067
Богородский
Dmitriy40 в сообщении #1623516 писал(а):
Вот к примеру насколько точна оценка 9-ок в ваших интервальчиках?

Конкретика. Вот я обсчитал интервал в 500 ярдов:
$$35720\cdot10^{11}-35725 \cdot10^{11} $$
Там нашлись 1340 тех самых 7-к. А для 9-к это именно интервальчик, а не интервал, там должно быть в среднем 4 тех самых 9-ки, а нашлись 6.

Всего лишь 6, а нужна хотя бы тысяча. Нечего и говорить о точности для 9-к.

А о точности для 7-к говорить можно.

А для 5-к вполне достаточно интервала в 30 ярдов.

 Профиль  
                  
 
 Re: Симметричные кортежи из последовательных простых чисел
Сообщение23.12.2023, 16:04 
Аватара пользователя


29/04/13
8067
Богородский
Yadryara в сообщении #1623527 писал(а):
А для 9-к это именно интервальчик, а не интервал,

А какой нужно взять интервал для 9-к? Согласно таблице, одна центральная 9-ка приходится в среднем на 307 центральных 7-к. Значит 200 трлн. должно хватить за глаза, чтобы набрать тысячу 9-к.

Для начала выбрать область вокруг 357, например такую:

$$3470\cdot10^{12}-3670 \cdot10^{12} $$

Скорее всего надо перекашивать, а не просто отнять и прибавить по 100 трлн. Умножить количество найденных 9-к на 50, затем сравнить с 87327. В зависимости от результата поискать ещё 9-ки пониже или повыше. Подгонка под ответ, которая позволяет лучше понять как ведут себя цепочки (да, мне тоже удобней говорить цепочки, а не кортежи).

Ну и попробовать применить найденные закономерности на других высотах и для других цепочек.

 Профиль  
                  
 
 Re: Симметричные кортежи из последовательных простых чисел
Сообщение23.12.2023, 18:16 
Заслуженный участник


20/08/14
11710
Россия, Москва
Yadryara в сообщении #1623527 писал(а):
Вот я обсчитал интервал в 500 ярдов:
Там нашлись 1340 тех самых 7-к.
Т.е. в полном 1e16 по вашему их должно быть $1340\dfrac{10^{16}}{500\cdot10^9}=26.8\cdot10^6$. Я (тоже на PARI) посчитал за 36 минут 2241 интервал по 23#=223e6 (т.е. суммарно практически 5e11) равномерно по 1e16 и нашёл 1295 семёрок, т.е. в полном 1e16 их должно быть 25.9млн. Разницу в оценке, 26.8 или 25.9 миллионов можно считать несущественной погрешностью, хоть и порядка 3%.
Для сравнения посчитал и интервал 35720e11-35725e11 с шагом 5:1 и нашёл 291 семёрок, т.е. в полном 5e11 их было бы около 1455шт, по сравнению с вашими 1340шт, согласие считаю неплохое, хотя погрешность и почти 9%.
Так что да, вообще говоря вполне можно брать вероятность в одной точке, вопрос лишь одинакова ли она для всех длин, это стоило бы проверить. А проход по всему диапазону с большим шагом не требует такой проверки, именно в этом его большая точность.

 Профиль  
                  
 
 Re: Симметричные кортежи из последовательных простых чисел
Сообщение23.12.2023, 18:33 
Аватара пользователя


29/04/13
8067
Богородский
Dmitriy40 в сообщении #1623557 писал(а):
Для сравнения посчитал и интервал 35720e11-35725e11 с шагом 5:1 и нашёл 291 семёрок,

А сколько времени у Вас займёт полный обсчёт этого интервала? Подозреваю, что меньше часа. У меня чуть ли не 8 часов. Неплохо бы сверить, вдруг у кого-то из нас ошибка в проге.

Dmitriy40 в сообщении #1623557 писал(а):
Т.е. в полном 1e16 по вашему их должно быть $1340\dfrac{10^{16}}{500\cdot10^9}=26.8\cdot10^6$.

Да, и если разделить на 87327 как раз и получается табличный кэф 307.

 Профиль  
                  
 
 Re: Симметричные кортежи из последовательных простых чисел
Сообщение23.12.2023, 18:34 
Заслуженный участник


20/08/14
11710
Россия, Москва
Подсчитал по своему и меньшие длины:
3: 109485шт/1e9,
5: 232.9шт/1e9.
Согласие с выбранной Вами точкой ещё сильно лучше. Зато без выбора.

-- 23.12.2023, 18:44 --

Yadryara в сообщении #1623563 писал(а):
А сколько времени у Вас займёт полный обсчёт этого интервала? Подозреваю, что меньше часа. У меня чуть ли не 8 часов. Неплохо бы сверить, вдруг у кого-то из нас ошибка в проге.
Да те же минут 40. Сейчас запущу.
Я прохожу не по простым, а лишь по допустимым начальным числам цепочки, их 552960/23# (/223e6), потому и быстрее. Но ищется только один паттерн (и более длинные если надо, но это не делал). Вряд ли у меня ошибка, я сверял для начала интервала с проходом по всем простым. И думаю вся 9% разница между нами лишь из-за флуктуаций цепочек по диапазону, всё же три сотни не столь уж много для точного совпадения, вон на тройки и пятёрки посмотрите.

 Профиль  
                  
Показать сообщения за:  Поле сортировки  
Начать новую тему Ответить на тему  [ Сообщений: 641 ]  На страницу Пред.  1 ... 19, 20, 21, 22, 23, 24, 25 ... 43  След.

Модераторы: Karan, Toucan, PAV, maxal, Супермодераторы



Кто сейчас на конференции

Сейчас этот форум просматривают: нет зарегистрированных пользователей


Вы не можете начинать темы
Вы не можете отвечать на сообщения
Вы не можете редактировать свои сообщения
Вы не можете удалять свои сообщения
Вы не можете добавлять вложения

Найти:
Powered by phpBB © 2000, 2002, 2005, 2007 phpBB Group