Эвристический метод нахождения некоторых простых близнецов

[wrest]

Rak so dna в сообщении #1723463 писал(а):

Грубо говоря, шанс наугад взятого $160$-значного числа оказаться простым близнецом
$\approx\frac{1}{160^2}\approx 0,00004$

Примерно в 5-6 раз меньше :) Логарифм же натуральный надо, а не десятичный.

[Cantata]

Dmitriy40
Благодарю Вас за помощь!
На самом деле без ваших объяснений я бы не смогла разобраться!

Так бы и разглядывала диагональ в треугольном расположении чисел, где от 11 сразу много пар чисел простых близнецов идут, и не знала бы как к ним подойти)

[Yadryara]

Cantata в сообщении #1723464 писал(а):

У меня обычный домашний компьютер, который не потянет такие вычисления.

У меня вроде тоже обычный комп. Однако мировые рекорды устанавливать удавалось.

А нет ли желания поискать цепочки посложнее чем близнецы?

Cantata в сообщении #1723464 писал(а):

Я даже Pari GP всё считаю онлайн :)

Интересно почему. Коротенькие проги в интерпретаторе gp обычно очень легко считаются. Вставил код, да нажал на "Enter".

Или Вы пока PARI/gp себе не установили?

[Cantata]

Yadryara в сообщении #1723486 писал(а):

Однако мировые рекорды устанавливать удавалось.

Поздравляю!

Yadryara в сообщении #1723486 писал(а):

А нет ли желания поискать цепочки посложнее чем близнецы?

К кортежам? Вот эта задача с n звездочками. Сейчас точно не готова - недавно тему о кортежах читала и поняла, что у меня пока нет идей.

Yadryara в сообщении #1723486 писал(а):

Или Вы пока PARI/gp себе не установили?

Да, так и есть)

[Rak so dna]

wrest в сообщении #1723479 писал(а):

Примерно в 5-6 раз меньше :) Логарифм же натуральный надо, а не десятичный.

Я делал оценку по порядку величины (там ведь ещё множитель вроде был какой-то). Сейчас пересчитал по первой гипотезе Харди-Литтлвуда: $\approx 0,000009781$

Dmitriy40 в сообщении #1723477 писал(а):

Как видите для 300000 последовательных $k$ после $7\cdot10^{80}$ выдаёт сотню простых близнецов.

Ну, выходит что ТС нашла способ получения простых близнецов примерно в $30$ раз эффективнее обычного угадывания (если эти расчёты подтвердятся и на других выборках). Неужели и для бОльших $k$ преимущество сохранится?

[Cantata]

Rak so dna в сообщении #1723490 писал(а):

Ну, выходит что ТС нашла способ получения простых близнецов примерно в $30$ раз эффективнее обычного угадывания

Сейчас посмотрела, оказывается результаты даны для формулы лидера, которую нашел tolstopuz.
Она по праву называется лидером!

Моя формула скромнее дает результаты:

? {
ss = 10^5;
step = 10^5;
n = 0;
for(kk = 1, 3*10^5,
    k = kk + 7*10^80;
    x = 33*k^2 + 23595*k;
    if(ispseudoprime(x-1) && ispseudoprime(x+1), n++);
    if(kk == ss,
        print(kk, ":", n);
        ss += step;
    )
);
}
100000:25
200000:49
300000:86

-- добавлено через 1 минуту --

Dmitriy40
Код изумительный! Такие числа посчитал у меня онлайн меньше чем за минуту!

[wrest]

Rak so dna в сообщении #1723490 писал(а):

Ну, выходит что ТС нашла способ получения простых близнецов примерно в $30$ раз эффективнее обычного угадывания

Ну так даже $30n \pm 1$ даёт это преимущество, ну а бОльшие праймориалы - ещё больше.

-- добавлено через 6 минут --

Rak so dna в сообщении #1723490 писал(а):

Неужели и для бОльших $k$ преимущество сохранится?

Решето же, отсеивает сразу для всех k

[Cantata]

Оказывается эффективность зависит от локальной модулярной структуры интервала.
Выбор лучшей формулы может определяться конкретным диапазоном поиска.

Вот моя формула на другом интервале:

? {
    ss = 10^5;
    step = 10^5;
    n = 0;
    for(kk = 1, 3*10^5,
        k = kk + 3*10^80;
        x = 33*k^2 + 23595*k;
        if(ispseudoprime(x-1) && ispseudoprime(x+1), n++);
        if(kk == ss,
            print(kk, ":", n);
            ss += step;
        )
    );
}
100000:32
200000:58
300000:106

[Rak so dna]

wrest в сообщении #1723493 писал(а):

Ну так даже $30n \pm 1$ даёт это преимущество

Не даёт. Проверил $5$ случайных диапазонов длинной $100000$ на числах порядка $10^{160}$, везде получилось в $3$-$4$ раза хуже.

[wrest]

Rak so dna в сообщении #1723495 писал(а):

везде получилось в $3$-$4$ раза хуже.

Хуже чем что?

[Yadryara]

Cantata в сообщении #1723489 писал(а):

Сейчас точно не готова - недавно тему о кортежах читала и поняла, что у меня пока нет идей.

На всякий случай. О кортежах речь идёт не только в тех темах, где в названии есть слово "кортеж", но и например в той, где есть слово "мечта"...