Пентадекатлон мечты

[Dmitriy40]

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

Это понятно было, но сейчас-то очень похоже, что Вы именно реально посчитали за 10 часов паттерн b127.

Сколько раз вам надо повторить что результаты этого реального запуска нигде не используются? Разве я где-то сказал что паттерн b127 решений не даёт? Мною используется только оценка затраченного времени и её погрешность (из-за 8e34<8.021561347e34) меня полностью устраивает. Хоть реального запуска, хоть какого.

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

Во-первых. Даже если превысит, стартовое число всё равно может быть подходящим.

Я сразу сказал что буду считать только варианты с шагом меньше предела. Потому что меня не интересуют точные значения! Мне достаточно приблизительных оценок снизу (и/или сверху)! Сколько ещё раз и это тоже повторить?

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

Во-вторых, Ваша формулировка эквивалентна этой?

Ответьте сами. Или возводить в квадрат тоже не умеете и символов "sqrt()" в комментарии не заметили? Что ж за мания то задавать вопросы, ответ на которые можно самостоятельно посчитать на калькуляторе быстрее чем набирать вопрос ...

[Yadryara]

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

Сколько раз вам надо повторить что результаты этого реального запуска нигде не используются?

Не знаю сколько раз. Ну то есть я правильно понял, что это был реальный запуск.

Результаты этого реального запуска пока нигде не используются. Если это был корректый запуск, вполне можно использовать.

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

Разве я где-то сказал что паттерн b127 решений не даёт?

Вроде нигде. А что, даёт?

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

Потому что меня не интересуют точные значения!

А меня интересуют.

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

Ответьте сами. Или возводить в квадрат тоже не умеете и символов "sqrt()" в комментарии не заметили?

Меня интересовало Ваше мнение независимо от того, что я посчитал сам и умею ли я возводить в квадрат. sqrt не просто заметил, но и проверил численно.

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

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

Всё-таки я не называл какие-то Ваши вопросы манией.

А так, да сторонние(и не только) наблюдатели запросто могут посчитать манией всё, чем мы здесь занимаемся. 11 месяцев и две сотни страниц цепочки считаем.

Правильно ли я понимаю, что в Пари

i++ эквивалентно i=i+1

и

i+=x эквивалентно i=i+x ?

[Dmitriy40]

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

Правильно ли я понимаю, что в Пари
i++ эквивалентно i=i+1
и
i+=x эквивалентно i=i+x ?

Правильно.
И это действует не только на плюс, но и на минус, а второй вариант и на умножение и деление (и логические операции).

С операциями x++ и x-- есть тонкость: они фактически выполняются как ++x и --x (в синтаксисе С/С++), т.е. сначала изменение значения, а потом его использование. Например вполне корректно n=0;v[n++]=1;, присвоение будет в v[1] (а в языках типа C в v[0]).

[Dmitriy40]

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

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

Потому что меня не интересуют точные значения!

А меня интересуют.

Я так думаю это бессмысленно: все такие паттерны (у которых шаг больше предела или даже немногим меньше) проще и быстрее проверить принципиально другим способом, не перебором всех 6-ти простых, а лишь одного.

Собственно грубо прикинуть количество вариантов с шагом больше предела тоже несложно: возьмём к примеру паттерн b446, в нём простые подставляются на места с [13,3,5,1,3,11], соответственно на них надо подставлять простые до квадратного корня из 8e34/17/[13,3,5,19,3,11] или примерно [5e14,1e15,8.1e14,4.2e14,1e15,5.5e14] штук. Совершенно очевидно что нет никакой необходимости перебирать все возможные варианты (которых примерно 1e89), достаточно проверить лишь 4.2e14 вариантов простых на месте с 1 (подставляя для каждого из них все возможные значения qr до 7.2e9) и всё, причём до 90% из них точно не дают решений и их можно проверить очень быстро (без numdiv и даже isprime). Конечно перебирать 4.2e14 простых до 1.57e16, да ещё и для кучи разных qr, тоже весьма не быстро, однако это совершенно точно на много порядков быстрее вычисления всех 1e89 начальных чисел и проверки их на превышение порога. И потому реально знать сколько там неделимых паттернов (1e89 или несколько меньше) абсолютно незачем.

[Yadryara]

И чего все молчат...

Ну ладно впс, я не горю желанием разговаривать в таком тоне.

Наименьшее Общее Кратное(LCM).

Какие тут могут быть сложности... И Пари считает сразу для многих чисел и Альфа. Берёшь кэфы из паттерна и подставляешь.

Насчёт параметра $p$ в pcoul. Даже мне понятно, что он ограничивает подставляемые простые в квадратах(и бо́льших степенях) сверху.

Huz в сообщении #1574556 писал(а):

Here are the highest primes that appear as a square or higher power for current upper bounds of $D(12,k)$: $(2, 7, 7, 17, 19, 13, 13, 29, 29, 23, 67, 113, 37, 47, 101)$.

А в других цепочках встречалась когда-нибудь простое больше $113$? О-о-очень вряд ли больше $200$. Это я к тому, что ограничение в $1000$ явно завышенное.

[EUgeneUS]

EUgeneUS в сообщении #1575619 писал(а):

Dmitriy40 в сообщении #1575531

писал(а):
Полтора месяца назад запускал pcoul -p100 -x6e26 -g9 -f13 12 13, заняло два дня, нашлась известная 13-ка, и больше ничего.
Спасибо, что написали/напомнили. Это, конечно, резко снижает вероятность найти более "низкую" цепочку.

У меня уже около недели считается такой запуск:
pcoul -p100 -x5e26:12e26 -X -g9 -f13 12 13

Правда опцию -X включил только после того, как был найден первый (он же известный кандидат). Надеюсь, это не повлияло на корректность работы.
Смысл этого эксперимента - посмотреть, а нет ли решений (пусть и с таким ограничением на простые) в диапазоне в два раза большем, чем известная минимальная цепочка.

Также у меня помаленьку (в одном потоке) считаются паттерны с $LCM=7214407200$ (два в пятой, остальные обязательные простые в квадратах), с ограничениями на простые: -p100:1000. А ограничение по диапазону - до 10e26.

И если ничего не найдется в этих экспериментах (пусть даже и цепочки большие, чем известная минимальная), то, видимо, нужно будет признать, что известная минимальная 13-ка - это некий "статистический выброс вниз", то есть цепочка, появившаяся сильно раньше, чем её следовало ожидать.

[EUgeneUS]

EUgeneUS в сообщении #1576525 писал(а):

У меня уже около недели считается такой запуск:
pcoul -p100 -x5e26:12e26 -X -g9 -f13 12 13

Посчитался примерно за 431 тысячу секунд. Других находок, кроме известной - нет.
Запустил аналогичный расчет для диапазона -x12e26:18e26.

[Dmitriy40]

Заморочился тут исключением вариантов $qr$ при $p^2$, до $qr<100$ кроме встречающихся только вместе $10p^2,14p^2,15p^2,21p^2$ и проверяемого отдельно $22p^2$ остались не запрещёнными лишь два варианта, $34p^2,93p^2$, причём $34p^2$ только для длин 13-, все остальные варианты $qr<100$ запретил. $34p^2$ до 6e26 проверяется быстро, а вот с $93p^2$ оказалась засада, там один подвариант тоже только 14-, один подвариант с большим шагом и проверяется быстро, а вот один подвариант с малым шагом и до 8.1e34 проверялся дней пять и только что завершился.
Нигде ничего не найдено.
Т.е. никакие варианты $qr<100$ не дают никаких решений до известных минимальных 13+ (включая и пятнашку). Это позволяет несколько снизить порог -p и для 13, и для 14, и для 15.

Конечно на PARI проверять простые до 3e16 всё равно весьма долго, но по этому поводу есть ещё несколько разных мыслей, постепенно их проверяю, может какая и сработает и позволит уменьшить -p до адекватных значений.

[Yadryara]

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

до 8.1e34 проверялся

Другое дело

Что касается НОК, то разумеется моя прога на Пари, которая генерит паттерны, заодно считает и LCM:

shag=lcm(pro[l..l+12]);

Это уже обсуждалось. Спасибо Дмитрию, что подсказал эту краткую форму. Это для 13-к. Для 15-к будет l+14.

А ежели считать LCM перемножая простые в степенях, то не забываем, что в число основных 3408 паттернов не входит кэф 243. Так что 3-ка в этом произведении всегда только в квадрате.

[EUgeneUS]

У меня завершился еще один эксперимент с интересным результатом.
$D(12,9)$ как видно, даже меньше, чем $D(12,8)$.
Поэтому для построения трендов, я рассматривал значение $D(12,9)$, как "правильное" значение для $D(12,8)$ (но цепочка по реализации некоторых небольших шансов удлинилась до 9).
Тогда в графике на девятом месте оказывается дырка. И мне нужно было следующее наименьшее значение для $U(12,9)$.
Для этого был сделан такой запуск:

pcoul -x:1e18 -X -f7 -g10 -p200 -rb12-09.log 12 9

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

(Оффтоп)

001 pcoul(12 9) -p200 -f7 -g10 -x1000000000000000000 *RT*
202 Candidate 530991107334472543 (404.83s)
202 Candidate 867736356029018143 (406.23s)
305 97^2 2^5 3^5 2 5.53^2 2^2.3 7.23^2 2 3.179^5 (542.64s)
202 Candidate 22353166734014175 (589.62s)
305 113^2 2^5 3.19^2 2.5^2 7.83^2 2^2.3 17^2 2 3^2.5: 0 / 0 (1083.38s)
305 41^2 2^3 3^2 2.89^2 5^2 2^2.3 7.59^2 2.101^2 3.13^2 (1626.52s)
305 89^2 2^3 3.43^2 2.101^2 5^2 2^2.3 7.11^2 2.13^2 3^2 (2170.45s)
202 Candidate 647527014904007646 (2535.33s)
305 2.61^5 3.5^2 2^5 139^2 2.3^2 7.109^2 2^2.5 3.13^2 2 (2707.22s)
202 Candidate 322337479870315421 (2900.70s)
305 41^2 2.29^2 3.83^2 2^5 5^2 2.3^2 7^2 2^2 3.17^2 (3248.33s)
202 Candidate 274283185690329820 (3621.34s)
202 Candidate 209536626496829020 (3625.11s)
202 Candidate 140820895491341020 (3639.41s)
202 Candidate 144544816376861020 (3652.80s)
305 2^2 5.103^2 2.3^2 7.199^2 2^5 3 2.5^2 137^2 2^2.3: 1 / 1 (3787.64s)
202 Candidate 165901919263463644 (3797.12s)
202 Candidate 237291359466618844 (3807.12s)
202 Candidate 859143285735282844 (3808.01s)
202 Candidate 885750266052335644 (3845.56s)
202 Candidate 587233688475048348 (3930.55s)
202 Candidate 25128134623824348 (3977.67s)
202 Candidate 968268739084569948 (3983.12s)
202 Candidate 586251962819865948 (3993.66s)
305 3.29^2 2^2 179^2 2.3^2 5^2 2^5 3.7^2 2 11^5: 0 / 0 (4327.92s)
202 Candidate 15724736975643 (4638.84s)
202 Candidate 378116262316752347 (4736.05s)
305 197^2 2^2.3 157^2 2.7 3.5^2 2^5 43^2 2.3^2 181^5 (4872.31s)
202 Candidate 869022436864632346 (5409.98s)
305 2.13^2 43^2 2^2.3 11^2 2.5^2 3.7^2 2^5 103^2 2.3^2: 1 / 3 (5413.95s)
305 3^2.5 2 13^2 2^2.3 7.113^2 2.5^2 3.79^2 2^5 149^2 (5960.83s)
305 29^2 2.131^5 3.7^2 2^2.5 59^2 2.3^2 149^2 2^5 3.5^2 (6503.73s)
202 Candidate 635713945023705817 (6513.59s)
202 Candidate 209107569431177049 (6624.59s)
202 Candidate 135671566001719449 (6661.03s)
305 3.47^2 2.31^5 7^2 2^2.3 5.13^2 2 3^2 2^5 11^5 (7044.17s)
305 43^2 2.47^2 3.31^2 2^2 5.101^2 2.3^2 7^2 2^3 3.11^2 (7590.53s)
367 coul(12, 9): recurse 469450034, walk 945629010, walkc 455912049 (7680.00s)
200 f(12, 9) = 15724736975643 (7680.00s)

После чего решил посмотреть, сколько цепочек будет с более большими простыми числами в квадрате.
Был сделан такой запуск:

pcoul -x:1e18 -X -f7 -g10 -p200:1000 -rb12-09-p200-1000.log 12 9

Считалось 57,5 тысяч секунд. И не нашлось ни одной цепочки.

[Huz]

EUgeneUS в сообщении #1576624 писал(а):

После чего решил посмотреть, сколько цепочек будет с более большими простыми числами в квадрате.
Был сделан такой запуск:

pcoul -x:1e18 -X -f7 -g10 -p200:1000 -rb12-09-p200-1000.log 12 9

Считалось 57,5 тысяч секунд. И не нашлось ни одной цепочки.

I have mentioned before that I do not believe two runs with "-p200" and "-p200:1000" are equivalent to a single run with "-p1000", so I don't know if this test is reliable. It would be worth checking whether a run with "-p100:200" correctly finds all examples it should out of your first list - I think there are five.

If I checked correctly, highest squared primes from each of those 23 candidates, in order, go like: 23 31 43 23 127 149 61 157 41 53 19 23 151 29 17 41 43 29 29 31 41 103 71.

[EUgeneUS]

Huz в сообщении #1576628 писал(а):

I have mentioned before that I do not believe two runs with "-p200" and "-p200:1000" are equivalent to a single run with "-p1000", so I don't know if this test is reliable.

Oh. Thanks for this tip.

Huz в сообщении #1576628 писал(а):

It would be worth checking whether a run with "-p100:200" correctly finds all examples it should out of your first list - I think there are five.

I ran this check. In the first hour, four of the five expected chains were found. The fifth should be found at the end of the second hour. UPD: The fifth chain was also found

This is surprising if we assume that the probability of finding a chain is approximately the same for each checking item, regardless of the size of the substituting primes in the power of two or higher.
Then we should expect that the quantity of chains found should be roughly proportional to the quantity of items checked. And in proportion to the time of the check in the end.
But in this case, it doesn't seem to work.
It seems that when substituting large primes into squares (and higher powers), a large overhead occurs, which reduces the likelihood of finding chains in each individual item check.
This is indirectly indicated by the fact that in this case none of the two acceleration methods proposed earlier by Dmitry was used.

[Yadryara]

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

Итого, системы паттернов исчерпываются такими количествами паттернов:

$D(12, 11)$ — 1163 (1044);
$D(12, 12)$ — 578 (504);
$D(12, 13)$ — 3704 (3408);
$D(12, 14)$ — 1116 (1044);
$D(12, 15)$ — 528 (488).

В скобках количество основных(бесквадратных) паттернов.

Для полноты картины посчитал и для

$D(12, 10)$ — 103 (78).

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

А слабо матрицу оформить как матрицу?

Не слабо:

(103)

{[    9,  10,  1,  12,  1,  14,  15,  32,  1,  18;  
9,  10,  1,  12,  1,  98,  15,  32,  1,  18;  
9,  10,  1,  12,  1,  14,  15,  32,  1,  18;  
9,  10,  1,  12,  1,  14,  15,  32,  1,  18;  
45,  2,  1,  12,  7,  50,  3,  32,  1,  18;  
45,  2,  1,  12,  1,  50,  147,  32,  1,  18;  
45,  2,  1,  12,  49,  50,  3,  32,  1,  18;  
45,  98,  1,  12,  1,  50,  3,  32,  7,  18;  
45,  2,  1,  12,  7,  50,  3,  32,  1,  18;  
45,  2,  1,  12,  16807,  50,  3,  32,  1,  18;  
45,  2,  1,  12,  7,  50,  3,  32,  1,  18;  
243,  10,  1,  12,  1,  14,  15,  32,  1,  18;  
243,  10,  1,  12,  1,  98,  15,  32,  1,  18;  
243,  10,  1,  12,  1,  14,  15,  32,  1,  18;  
243,  10,  1,  12,  1,  14,  15,  32,  1,  18;  
10,  1,  12,  1,  14,  15,  32,  1,  18,  1;  
10,  1,  12,  1,  98,  15,  32,  1,  18,  1;  
10,  1,  12,  1,  14,  15,  32,  1,  18,  1;  
10,  1,  12,  1,  14,  15,  32,  1,  18,  1;  
2,  1,  12,  7,  50,  3,  32,  1,  18,  5;  
2,  1,  12,  1,  50,  147,  32,  1,  18,  5;  
2,  1,  12,  49,  50,  3,  32,  1,  18,  5;  
98,  1,  12,  1,  50,  3,  32,  7,  18,  5;  
2,  1,  12,  7,  50,  3,  32,  1,  18,  5;  
2,  1,  12,  16807,  50,  3,  32,  1,  18,  5;  
2,  1,  12,  7,  50,  3,  32,  1,  18,  5;  
1,  12,  1,  14,  15,  32,  1,  18,  1,  20;  
1,  12,  1,  98,  15,  32,  1,  18,  1,  20;  
1,  12,  1,  14,  15,  32,  1,  18,  1,  20;  
1,  12,  1,  14,  15,  32,  1,  18,  1,  20;  
1,  12,  1,  2,  75,  32,  7,  18,  1,  20;  
1,  12,  1,  14,  75,  32,  1,  18,  1,  20;  
1,  12,  1,  2,  75,  32,  49,  18,  1,  20;  
1,  12,  1,  98,  75,  32,  1,  18,  1,  20;  
1,  12,  1,  14,  75,  32,  1,  18,  1,  20;  
1,  12,  1,  2,  75,  32,  16807,  18,  1,  20;  
1,  12,  1,  14,  75,  32,  1,  18,  1,  20;  
1,  12,  7,  50,  3,  32,  1,  18,  5,  28;  
1,  12,  1,  50,  3,  32,  7,  18,  5,  4;  
1,  12,  1,  50,  3,  32,  49,  18,  5,  4;  
1,  12,  1,  50,  147,  32,  1,  18,  5,  4;  
1,  12,  49,  50,  3,  32,  1,  18,  5,  28;  
1,  12,  1,  50,  3,  32,  16807,  18,  5,  4;  
3,  28,  5,  18,  1,  32,  3,  50,  7,  12;  
3,  4,  5,  18,  7,  32,  3,  50,  1,  12;  
3,  28,  5,  18,  1,  32,  3,  50,  49,  12;  
3,  4,  5,  18,  1,  32,  147,  50,  1,  12;  
3,  4,  5,  18,  49,  32,  3,  50,  1,  12;  
3,  4,  5,  18,  7,  32,  3,  50,  1,  12;  
3,  28,  5,  18,  1,  32,  3,  50,  16807,  12;  
3,  4,  5,  18,  16807,  32,  3,  50,  1,  12;  
3,  4,  5,  18,  7,  32,  3,  50,  1,  12;  
3,  4,  5,  18,  343,  32,  3,  50,  1,  12;  
3,  20,  1,  18,  7,  32,  75,  2,  1,  12;  
3,  20,  1,  18,  49,  32,  75,  2,  1,  12;  
3,  20,  1,  18,  7,  32,  75,  2,  1,  12;  
3,  20,  1,  18,  16807,  32,  75,  2,  1,  12;  
3,  20,  1,  18,  7,  32,  75,  2,  1,  12;  
3,  20,  1,  18,  343,  32,  75,  2,  1,  12;  
12,  1,  14,  15,  32,  1,  18,  1,  20,  21;  
12,  1,  14,  75,  32,  1,  18,  1,  20,  21;  
12,  1,  2,  75,  32,  7,  18,  1,  20,  3;  
12,  1,  2,  75,  32,  49,  18,  1,  20,  3;  
12,  1,  2,  75,  32,  7,  18,  1,  20,  3;  
12,  1,  2,  75,  32,  16807,  18,  1,  20,  3;  
12,  1,  2,  75,  32,  7,  18,  1,  20,  3;  
12,  7,  50,  3,  32,  1,  18,  5,  28,  3;  
12,  1,  50,  3,  32,  7,  18,  5,  4,  3;  
12,  1,  50,  3,  32,  49,  18,  5,  4,  3;  
12,  1,  50,  147,  32,  1,  18,  5,  4,  3;  
12,  49,  50,  3,  32,  1,  18,  5,  28,  3;  
12,  1,  50,  3,  32,  7,  18,  5,  4,  3;  
12,  1,  50,  3,  32,  16807,  18,  5,  4,  3;  
12,  1,  50,  3,  32,  7,  18,  5,  4,  3;  
28,  5,  18,  1,  32,  3,  50,  7,  12,  1;  
4,  5,  18,  7,  32,  3,  50,  1,  12,  1;  
28,  5,  18,  1,  32,  3,  50,  49,  12,  1;  
4,  5,  18,  1,  32,  147,  50,  1,  12,  1;  
4,  5,  18,  49,  32,  3,  50,  1,  12,  1;  
4,  5,  18,  7,  32,  3,  50,  1,  12,  1;  
28,  5,  18,  1,  32,  3,  50,  16807,  12,  1;  
4,  5,  18,  16807,  32,  3,  50,  1,  12,  1;  
4,  5,  18,  7,  32,  3,  50,  1,  12,  1;  
4,  5,  18,  343,  32,  3,  50,  1,  12,  1;  
20,  1,  18,  7,  32,  75,  2,  1,  12,  1;  
20,  1,  18,  1,  32,  75,  98,  1,  12,  1;  
20,  1,  18,  49,  32,  75,  2,  1,  12,  1;  
20,  1,  18,  7,  32,  75,  2,  1,  12,  1;  
20,  1,  18,  16807,  32,  75,  2,  1,  12,  1;  
20,  1,  18,  7,  32,  75,  2,  1,  12,  1;  
20,  1,  18,  343,  32,  75,  2,  1,  12,  1;  
5,  18,  1,  32,  3,  50,  7,  12,  1,  2;  
5,  18,  7,  32,  3,  50,  1,  12,  1,  98;  
5,  18,  1,  32,  3,  50,  49,  12,  1,  2;  
5,  18,  1,  32,  147,  50,  1,  12,  1,  2;  
5,  18,  1,  32,  3,  50,  16807,  12,  1,  2;  
18,  1,  32,  3,  50,  7,  12,  1,  2,  45;  
18,  7,  32,  3,  50,  1,  12,  1,  98,  45;  
18,  1,  32,  3,  50,  49,  12,  1,  2,  45;  
18,  1,  32,  147,  50,  1,  12,  1,  2,  45;  
18,  1,  32,  3,  50,  7,  12,  1,  2,  45;  
18,  1,  32,  3,  50,  16807,  12,  1,  2,  45;  
18,  1,  32,  3,  50,  7,  12,  1,  2,  45      ]}

Прогу готовлю к публикации. Даже для 10-к она весьма длинная — 260 строк. А для 14-15-к почти вдвое длиннее — 480 строк.

[Yadryara]

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

Прогу готовлю к публикации. Даже для 10-к она весьма длинная — 260 строк.

Подсократил. Некоторые пустые строки оставил для удобства восприятия.

(PARI, 245 строк)

kvar=0;kvold=0;kosn=0;
kv=[0,0,0,0,0,0];
{print();write1("xxx.txt","{[  ");
u=[10,14,15,21,343]; lu=5;
ul=vector(343);ul[10]=1;ul[14]=2;ul[15]=3;ul[21]=4;ul[343]=5;

print("Mnojiteli pri kvadratah: ",u);print();print();

print("Ostatki po modulu 64");print();
ost64=matrix(lu,13,i,j,-1);
for(i=1,lu,v=vector(64,n,-1);
for(j=1,64, x=(u[i]*j^2)%64; v[x+1]=x);
vres=setminus(Set(v),[-1]);nom=0;
print1(u[i],"      ");
for(j=1,#vres,
if(25<=vres[j] && vres[j]<=39 && vres[j]%4<>0,
nom++;ost64[i,nom]=vres[j];print1(vres[j],"  ")));
print();print());print();

print("Ostatki po modulu 3");print();
ost3=matrix(lu,3,i,j,-1);
for(i=1,lu,v=vector(3,n,-1);
for(j=1,3, x=(u[i]*j^2)%3; v[x+1]=x);
vres=setminus(Set(v),[-1]);
print1(u[i],"      ");
for(j=1,#vres, ost3[i,j]=vres[j];print1(vres[j],"  "));
print();print());print();

print("Ostatki po modulu 9");print();
ost9=matrix(lu,4,i,j,-1);
for(i=1,lu,v=vector(9,n,-1);
for(j=1,9, x=(u[i]*j^2)%9; v[x+1]=x);
vres=setminus(Set(v),[-1]);
print1(u[i],"      ");
for(j=1,#vres, ost9[i,j]=vres[j];print1(vres[j],"  "));
print();print());print();

print("Ostatki po modulu 243");print();
ost243=matrix(lu,13,i,j,-1);
for(i=1,lu,v=vector(243,n,-1);
for(j=1,243, x=(u[i]*j^2)%243; v[x+1]=x);
vres=setminus(Set(v),[-1]);nom=0;
print1(u[i],"      ");
for(j=1,#vres,
if(vres[j]<15 || vres[j]>243-15, nom++;
ost243[i,nom]=vres[j];print1(vres[j],"  ")));
print();print());print();

print("Ostatki po modulu 5");print();
ost5=matrix(lu,4,i,j,-1);
for(i=1,lu,v=vector(5,n,-1);
for(j=1,5, x=(u[i]*j^2)%5; v[x+1]=x);
vres=setminus(Set(v),[-1]);
print1(u[i],"      ");
for(j=1,#vres, ost5[i,j]=vres[j];print1(vres[j],"  "));
print();print());print();

print("Ostatki po modulu 25");print();
ost25=matrix(lu,11,i,j,-1);
for(i=1,lu,v=vector(25,n,-1);
for(j=1,25, x=(u[i]*j^2)%25; v[x+1]=x);
vres=setminus(Set(v),[-1]);
print1(u[i],"      ");
for(j=1,#vres, ost25[i,j]=vres[j];print1(vres[j],"  "));
print();print());print();

print("Ostatki po modulu 7");print();
ost7=matrix(lu,4,i,j,-1);
for(i=1,lu,v=vector(7,n,-1);
for(j=1,7, x=(u[i]*j^2)%7; v[x+1]=x);
vres=setminus(Set(v),[-1]);
print1(u[i],"      ");
for(j=1,#vres, ost7[i,j]=vres[j];print1(vres[j],"  "));
print();print());print();

print("Ostatki po modulu 49");print();
ost49=matrix(lu,13,i,j,-1);
for(i=1,lu,v=vector(49,n,-1);
for(j=1,49, x=(u[i]*j^2)%49; v[x+1]=x);
vres=setminus(Set(v),[-1]);nom=0;
print1(u[i],"      ");
for(j=1,#vres,
if(vres[j]<15 || vres[j]>49-15, nom++;
ost49[i,nom]=vres[j];print1(vres[j],"  ")));
print();print());print();

print("Ostatki po modulu 343");print();
ost343=matrix(lu,13,i,j,-1);
for(i=1,lu,v=vector(343,n,-1);
for(j=1,343, x=(u[i]*j^2)%343; v[x+1]=x);
vres=setminus(Set(v),[-1]);nom=0;
print1(u[i],"      ");
for(j=1,#vres,
if(vres[j]<15 || vres[j]>343-15, nom++;ost343[i,nom]=vres[j];print1(vres[j],"  ")));
print();print());print();

print("Ostatki po modulu 16807");print();
ost16807=matrix(lu,14,i,j,-1);
for(i=1,lu,v=vector(16807,n,-1);
for(j=1,16807, x=(u[i]*j^2)%16807; v[x+1]=x);
vres=setminus(Set(v),[-1]);nom=0;
print1(u[i],"      ");
for(j=1,#vres,
if(vres[j]<15 || vres[j]>16807-15, nom++;
ost16807[i,nom]=vres[j];print1(vres[j],"  ")));
print();print());print();print();print();

pro   =[0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,  1,1,1,1,1,1,1, 1,1,1,1,1,    1,1,  1];
v2    =[0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,  1,2,1,4,1,2,1,32,1,2,1,4,    1,2,  1];
vs31  =[0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,  9,1,1,3,1,1,3, 1,1,9,1,1,    3,1,  1];
vs32  =[0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,  1,1,3,1,1,9,1, 1,3,1,1,3,    1,1,  9];
vs33  =[0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,243,1,1,3,1,1,3, 1,1,9,1,1,    3,1,  1];
vs34  =[0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,  1,1,3,1,1,9,1, 1,3,1,1,3,    1,1,243];
vst5  =[0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,  1,1,1,1,5,1,1, 1,1,5,1,1,    1,1, 25,1, 1,1,1,5,1,1,1,1,5,1,1,1,1,5,1];
vst7  =[0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,  1,1,7,1,1,1,1, 1,1,7,1,1,    1,1,  1,1,49,1,1,1,1,1,1,7,1,1,1,1,1,1,7,1,1];
vstp7 =[0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,  1,1,1,1,1,7,1, 1,1,1,1,1,16807,1,  1,1, 1,1,1,7,1,1,1,1,1,1,7,1,1,1,1,1,1];
vcub7 =[0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,  1,1,1,1,1,7,1, 1,1,1,1,1,  343,1,  1,1, 1,1,1,7,1,1,1,1,1,1,7,1,1,1,1,1,1];
v3=vector(39);v5=vector(55);v7=vector(57);
for(l=25,30,r=l+9;
for(i=1,39,v3[i]=vs31[i]);
for(i3=1,3,
if(i3==2,for(i=1,39,v3[i]=vs32[i]));
if(i3*l>75 && i3*l<90,next);
if(i3*l==75,for(i=1,39,v3[i]=vs33[i]));
if(i3*l==90,for(i=1,39,v3[i]=vs34[i]));
for(i=1,55,v5[i]=vst5[i]);
for(i5=1,15,
for(i=1,57,v7[i]=vst7[i]);
for(i7=1,37,
if(i7==18,for(i=1,57,v7[i]=vstp7[i]));
if(i7==28,for(i=1,57,v7[i]=vcub7[i]));

zap=0;kd4=0;k12=0;k18=0;
k3=0;m3=0;k5=0;m5=0;k7=0;m7=0;
k9=0;m9=0;k25=0;m25=0;k49=0;m49=0;
k343=0;m343=0;
k243=0;m243=0;
k16807=0;m16807=0;
km=0;ksov=0;kpr=0;y=0;
m=vector(7);mm=vector(7);
for(i=l,r,pro[i]=v2[i]*v3[i]*v5[i]*v7[i];
if(v3[i]==3 && k3==0, k3++;m3=i);
if(v3[i]==9 && k9==0, k9++;m9=i);
if(v3[i]==243, k243++; m243=i);
if(v5[i]==5 && k5==0, k5++;m5=i);
if(v7[i]==7 && k7==0, k7++;m7=i);
if(v7[i]==49, k49++;m49=i);
if(v7[i]==343, m343=i);
if(v7[i]==16807, k16807++;m16807=i);
);
for(i=l,r,
if(numdiv(pro[i]) > 6, zap=zap+1;next);
if(numdiv(pro[i]) == 4, kd4++);
if(pro[i] == 35, zap=zap+1;next);
if(pro[i] == 12, k12=k12+1);
if(pro[i] == 18, k18=k18+1);
if(pro[i] == 50 || pro[i] == 75, k25=k25+1;m25=i);
if(pro[i] == 10, if(i == 26, km=km+1;m[km]=pro[i]; mm[km]=i, zap=zap+1;next));
if(pro[i] == 14, if(i == 30, km=km+1;m[km]=pro[i]; mm[km]=i, zap=zap+1;next));
if(pro[i] == 15, if(i == 31 || i == 39, km=km+1;m[km]=pro[i];mm[km]=i, zap=zap+1;next));
if(pro[i] == 21, if(i == 29 || i == 37, km=km+1;m[km]=pro[i];mm[km]=i, zap=zap+1;next));
if(pro[i] == 343, if(i == 31 || i == 39, k343=k343+1;km++;m[km]=pro[i];mm[km]=i, zap=zap+1));
);

if(zap==0 && k12==1 && k18==1,
if(kd4 > 0,
for(y=1,km,
if(ksov == kpr, 
if(m[y] <> 0, x=ul[m[y]];

if(k3 == 1, kpr=kpr+1;
d=mm[y]-m3;d=d%3;
for(j=1,3, if(ost3[x,j]==d, ksov=ksov+1));
if(ksov<kpr, next));

if(k5 == 1, kpr=kpr+1;
d=mm[y]-m5;d=d%5;
for(j=1,4, if(ost5[x,j]==d, ksov=ksov+1));
if(ksov<kpr, next));

if(k7 == 1, kpr=kpr+1;
d=mm[y]-m7;d=d%7;
for(j=1,4, if(ost7[x,j]==d, ksov=ksov+1));
if(ksov<kpr, next));

if(k9 == 1, kpr=kpr+1;
d=mm[y]-m9;d=d%9; 
for(j=1,4, if(ost9[x,j]==d, ksov=ksov+1));
if(ksov<kpr, next));

if(ksov==0, quit);

if(k25 == 1, kpr=kpr+1;
d=mm[y]-m25; d=d%25;
for(j=1,11, if(ost25[x,j]==d, ksov=ksov+1));
if(ksov<kpr, next));

if(k49 == 1, kpr=kpr+1;
d=mm[y]-m49; d=d%49;
for(j=1,13, if(ost49[x,j]==d, ksov=ksov+1));
if(ksov<kpr, next));

if(k243 == 1, kpr=kpr+1;
d=mm[y]-m243; d=d%243;
for(j=1,13, if(ost243[x,j]==d, ksov=ksov+1));
if(ksov<kpr, next));

if(k343 == 1, kpr=kpr+1;
d=mm[y]-m343; d=d%343;
for(j=1,13, if(ost343[x,j]==d, ksov=ksov+1));
if(ksov<kpr, next));

if(k16807 == 1, kpr=kpr+1;
d=mm[y]-m16807; d=d%16807;
for(j=1,14, if(ost16807[x,j]==d, ksov=ksov+1));
if(ksov<kpr, next));
))));

if(ksov==kpr,

kvar++;if(kd4==0,kosn++);
print();print1(kvar,"   ",kd4,"   ",km, "  ",kpr, "   ",ksov);
print1("     ");

for(i=l,r,print1("   ",pro[i]);write1("xxx.txt",pro[i]);
if(i<r, write1("xxx.txt",",  ")));
print();print();write("xxx.txt",";  ");

shag=lcm(pro[l..l+9]);
));

for(i=l,l+26,v7[i]=v7[i+1]));
for(i=l,l+24,v5[i]=v5[i+1]));
);

kv[l-24]=kvar-kvold;kvold=kvar;

forstep(i=55,l,-1,vst5[i]=vst5[i-1]);
forstep(i=57,l,-1,vst7[i]=vst7[i-1];
vstp7[i]=vstp7[i-1];vcub7[i]=vcub7[i-1]);
);
print();print();write("xxx.txt","  ]}");
print(kv[1],"   ", kv[2],"   ", kv[3],"   ", kv[4],"   ",kv[5],"   ",kv[6]);
print();print("Vsego patternov:            ",kvar," (",kosn,")");print();
} quit;

[Yadryara]

И опять тишина. И никто даже не сказал сколько паттернов для 10-к генерит pcoul. Я же не знаю и сравнить не могу.

Если всё в порядке, то с системами паттернов для 12-ти делителей закончу.

И можно двигаться дальше:

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

Так что если даже Вы откажетесь делать новые проги, может быть, всё равно сможем это сделать, хотя и уйдут месяцы.

Дмитрий не отказался. Но непонятно, есть ли какие-то подвижки в прогах Асм+Пари.

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

а вот один подвариант с малым шагом и до 8.1e34 проверялся дней пять и только что завершился.

И не сказано как именно проверялся. Асм+Пари, только Пари или pcoul. А может и 4-м способом...

Чем заниматься-то дальше. 13-ками? Что уже посчитано? От Ахиллеса и Черепашки есть регулярная инфа. Периодически что-то сообщает Евгений. Но что-то не видно таблицы, посмотрев на которую можно понять что обсчитано, а что нет.