кортежи последовательных простых. ключ к 19-252

[Yadryara]

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

Эх, замедлилось продвижение к c19, всего лишь c27.

Ожидаемо замедлилось. Это регулярное замедление я видел на других паттернах. Поэтому и такой прогноз:

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

Устаканивание произойдёт, когда значение для $valids=19, len = 19$ станет отличным от нуля. Ориентировочно при 83# - 89#.

От c29 до c19 10 одиночных шагов и от 43# до 89# их тоже 10. И была надежда ещё хотя бы на один двойной шаг, например, как между 41# и 43#. Вроде топтания на месте не будет.

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

Это уже со счётом лишь половины (но числа в таблице сразу удвоил).

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

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

[Dmitriy40]

Посчиталось (c19-c50):

59#:  0,  0,  0,  0,  0,  0,  0, 2342, 169500, 4949320, 86498046, 1072162912, 9922991574, 69015165402, 361555473588, 1436573713066, 4372104902192, 10287665248786, 18821195521110, 26789470532208, 29566950432394, 25162208168768, 16402877618224, 8134573376362, 3050946438594, 862382999288, 183204299616, 28827888442, 3171808930, 210338416, 5822520,  0, sum=145544026521600

К сожалению надежды не оправдались, шаг одинарный, только c26>0.
Считалось 23ч (если в один поток).
В принципе можно посчитать и следующее значение 61#, за несколько дней, но это видимо пока максимум возможного, но будет похоже лишь с25. Смысла в этом не вижу (слишком много усилий надо чтобы запустить).
Прогноз на c19 остаётся 89#, досчитать нереально, даже ускорив ещё в разы.

PARI будет досчитывать 47# ещё часа два, если цифры вдруг разойдутся, то сообщу.

-- 19.03.2024, 13:12 --

Пожалуй покажу х64 код внутреннего цикла, обрабатывающего предпросчитанную таблицу 23# длиной 2304 с 128-битными масками a23 (по коду PARI выше) и входной маской (от больших простых) в RDX,RCX:

Используется синтаксис ASM

        mov     RAX,0
@@:     mov     R8,RCX
        and     R8,[tabl2304+RAX]
        popcnt  R8,R8
        mov     R9,RDX
        and     R9,[tabl2304+RAX+8]
        popcnt  R9,R9
        add     R8,R9
        add     qword [cnt+R8*8],2
        add     RAX,16
        cmp     RAX,2304*16/2
        jc      @b

Укорочение сделано по тупому (обработка лишь первой половины таблицы), но работает.
На х32 потребуется вдвое больше popcnt, которые и занимают всё время, каждая по такту, потому будет ровно вдвое медленнее. На процессорах старее Haswell (2013г) будет ещё на треть медленнее (три порта выполнения команд вместо 4-х).

[Yadryara]

Спасибо. Я пока в программерскую часть вникать не буду, мне жутко интересно можно ли предсказать значения $c$ на более ранних этапах. И уже есть успешное предсказание. Но это долго рассказывать. Ещё перепроверю идеи.

Я обсчитывал паттерны-тройки 6-6, 12-12, 18-18, 24-24 и 30-30. Вот здесь Ваша Асмовская прога могла бы пригодиться. Для обсчёта следующего: 36-36.

[Dmitriy40]

PARI досчитал 47# за 12.5ч, все цифры совпали, продолжаю верить результатам асма.

Пока sum<2e13 (т.е. мне хватит часа три посчитать) мне проще посчитать самому, чем переделывать программу под х32. Ещё и потому что программа тупая, использует готовые таблицы, получаемые из PARI.
Кстати если #a<65 (т.е. v[#v]<128+2*#v), то маски влезают в один х64 регистр и потому хватит одной popcnt (это точный аналог hammingweight) и считаться будет вдвое быстрее.
Так что показывайте паттерны, посмотрю как их считает PARI и досчитаю пару-тройку строк на асме. А, Вы похоже их уже и показали, непривычно мне в интервалах ... Сейчас посмотрю.

[Dmitriy40]

Yadryara
Держите все перечисленные паттерны (разворачивается):

Используется синтаксис Text

v=[0, 6, 12]
13#: 338, 750, 192, sum=1280
17#: 5482, 10134, 2304, sum=17920
19#: 97846, 156618, 32256, sum=286720
23#: 2113538, 3040254, 580608, sum=5734400
29#: 57992242, 77167566, 13934592, sum=149094400
31#: 1700950342, 2111393466, 362299392, sum=4174643200
37#: 59943705094, 70400583162, 11593580544, sum=141937868800

v=[0, 12, 24]
13#: 2, 58, 304, 568, 348, sum=1280
17#: 86, 1362, 5352, 7640, 3480, sum=17920
19#: 2738, 31134, 97848, 113240, 41760, sum=286720
23#: 85894, 787242, 2100984, 2092120, 668160, sum=5734400
29#: 3020486, 23883018, 56699976, 50791400, 14699520, sum=149094400
31#: 108456626, 758241438, 1626573576, 1328583080, 352788480, sum=4174643200
37#: 4445766722, 28275114606, 56036103672, 42597229400, 10583654400, sum=141937868800

v=[0, 18, 36]
13#: 0, 0, 6, 66, 330, 548, 258, 64, 8, sum=1280
17#: 0, 10, 264, 1978, 5998, 6484, 2574, 556, 56, sum=17920
19#: 10, 678, 9630, 49706, 104396, 86768, 29632, 5452, 448, sum=286720
23#: 878, 32142, 322458, 1262586, 2104176, 1479312, 453012, 74460, 5376, sum=5734400
29#: 54970, 1448466, 11526750, 37456182, 53688432, 33783624, 9581460, 1457748, 96768, sum=149094400
31#: 2987626, 62162082, 412064046, 1151158278, 1457440488, 834512112, 220996356, 31386852, 1935360, sum=4174643200
37#: 163741366, 2875476798, 16639524306, 41515668570, 47895775200, 25526829384, 6407605932, 862927884, 50319360, sum=141937868800

v=[0, 24, 48]
13#: 0, 0, 0, 0, 2, 44, 286, 454, 366, 120, 8, sum=1280
17#: 0, 0, 0, 12, 246, 1952, 5258, 6112, 3478, 818, 44, sum=17920
19#: 0, 0, 36, 1242, 12344, 50126, 93184, 84626, 37882, 6972, 308, sum=286720
23#: 0, 72, 4332, 69716, 443326, 1297482, 1884968, 1412596, 533592, 84928, 3388, sum=5734400
29#: 72, 10464, 313116, 3376772, 16240582, 38556930, 47587532, 31108060, 10369008, 1477656, 54208, sum=149094400
31#: 12480, 908760, 18271332, 149381628, 582558618, 1172334582, 1264682124, 736221324, 220679064, 28617544, 975744, sum=4174643200
37#: 1333080, 66531744, 1032827508, 6961064940, 23338431450, 41585795622, 40564648740, 21643408260, 5995213560, 725196040, 23417856, sum=141937868800

v=[0, 30, 60]
13#: 0, 0, 0, 0, 0, 0, 0, 26, 232, 800, 1030, 412, 58, 2, sum=2560
17#: 0, 0, 0, 0, 0, 6, 316, 2240, 8344, 13116, 9118, 2434, 258, 8, sum=35840
19#: 0, 0, 0, 2, 136, 2658, 20716, 84484, 176118, 182530, 87412, 17818, 1526, 40, sum=573440
23#: 0, 0, 20, 758, 17380, 168292, 868190, 2444364, 3700444, 2927118, 1131744, 195824, 14346, 320, sum=11468800
29#: 0, 38, 2598, 82832, 1177908, 8487904, 33786472, 75277014, 93091582, 61986798, 20841306, 3234106, 215762, 4480, sum=298188800
31#: 38, 6222, 316044, 6782432, 70709312, 397940624, 1270241482, 2325549950, 2419712822, 1386162222, 410718674, 57568946, 3510432, 67200, sum=8349286400
37#: 7514, 837414, 30460704, 493092640, 4110982480, 19161726988, 51845611146, 82147551226, 75387993370, 38761242290, 10490506582, 1366210942, 78103104, 1411200, sum=283875737600

v=[0, 36, 72]
13#: 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 48, 262, 474, 380, 106, 10, sum=1280
17#: 0, 0, 0, 0, 0, 0, 0, 0, 62, 624, 2694, 5954, 5736, 2422, 408, 20, sum=17920
19#: 0, 0, 0, 0, 0, 0, 56, 922, 7560, 32706, 77202, 94432, 56046, 15740, 1976, 80, sum=286720
23#: 0, 0, 0, 0, 10, 444, 9046, 80152, 390398, 1095858, 1755202, 1532930, 696038, 157296, 16442, 584, sum=5734400
29#: 0, 0, 0, 92, 3976, 79740, 817076, 4753016, 16691298, 35495044, 44604362, 31789682, 12234252, 2400058, 218908, 6896, sum=149094400
31#: 0, 0, 260, 17878, 487640, 6641892, 50493968, 229799540, 644261040, 1110622436, 1153892762, 698314552, 235091702, 41448078, 3469356, 102096, sum=4174643200
37#: 0, 520, 61954, 2504778, 47838660, 495578676, 3022427884, 11358675900, 26746388964, 39304488520, 35374886360, 18882335120, 5710842458, 918980622, 70918560, 1939824, sum=141937868800
41#: 520, 143148, 9744678, 279021870, 4104407820, 34488663612, 176228423588, 566090064612, 1156132065600, 1493579030680, 1198202504400, 578353157736, 160428651994, 23967375390, 1731142800, 44615952, sum=5393639014400

Показаны начиная с c3 и по последний ненулевой элемент.
Все 37# считались по минуте, 41# считалось 40 минут. Думаю считать 41# для остальных не нужно, там уже младший c3 ненулевой.

[Yadryara]

Благодарю.

Я, правда, только 36-36 просил посчитать. Похоже, что Вы так настроили, что для каждого паттерна считается сразу вся таблица. Круто!

(Статистика)

6-6      3                    5

 5#                           2          5-5 = 0

 7#                           4          7-5 = 2          4/2/2 = 1



12-12    3                    7

7#                           12          7-7 = 0

11#                          58         11-7 = 4        58/12/4 = 1 5/24

13#                         348         13-7 = 6       348/58/6 = 1



18-18    3                   11

11#                           2        11-11 = 0      

13#                           8        13-11 = 2         8/ 2/2 = 2 

17#                          56        17-11 = 6        56/ 8/6 = 1 1/6 

19#                         448        19-11 = 8       448/56/8 = 1 
 


24-24    3                   13

13#                           8        13-13 = 0       

17#                          44        17-13 = 4       44/   8/ 4 = 1 3/8 

19#                         308        19-13 = 6      308/  44/ 6 = 1 1/6 

23#                        3388        23-13 =10     3388/ 308/10 = 1 1/10 

29#                       54208        29-13 =16    54208/3388/16 = 1 



30-30    3                   16

13#                           2        13-16 = -3        

17#                           8        17-16 = 1        8/   2/ 1 = 4  

19#                          40        19-16 = 3       40/   8/ 3 = 1 2/3 

23#                         320        23-16 = 7      320/  40/ 7 = 1 1/7 

29#                        4480        29-16 =13     4480/ 320/13 = 1 1/13 

31#                       67200        31-16 =15    67200/4480/15 = 1   



36-36    3                   18

17#                          20        17-18 = -1         

19#                          80        19-18 = 1        80/     20/ 1 = 4

23#                         584        23-18 = 5       584/     80/ 5 = 1 23/50 

29#                        6896        29-18 =11      6896/    584/11 = 1 59/803 

31#                      102096        31-18 =13    102096/   6896/13 = 1 778/5603  

37#                     1939824        37-18 =19   1939824/ 102096/19 = 1

41#                    44615952        41-18 =23  44615952/1939824/23 = 1

Это статистика по последнему $c$ для всех паттернов. И ведь как всё здорово шло: кэфы(крайнее число справа) предсказуемо сходились к 1-це. Число 4480 я предсказал. Но кэф для 36-36 сначала болтануло вниз, потом вверх, а потом он раньше времени сравнялся с 1-цей...

[Dmitriy40]

(Оффтоп)

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

Я, правда, только 36-36 просил посчитать. Похоже, что Вы так настроили, что для каждого паттерна считается сразу вся таблица. Круто!

Не, каждая строка считается отдельным запуском, но между ними всего лишь запятую переставить, и перекомпилить, вот это сделал сразу, да.
А для смены паттерна запускать PARI и брать из него новые таблицы. Это минута, тут важнее внимательность/аккуратность, не перепутать ничего, где что раскомментировать. Ну и потом сравнение пары-тройки строк от PARI (в котором специально без удвоения для большей надёжности) и от асма и если ОК, то запуск, копирование результата, перестановка запятой, компиляция, запуск, ...
С одним вон паттерном, где 13#=2560 вместо 1280, сначала удивился чего это цифры вдруг не совпали, не сразу сообразил что таблица длиннее, для неё поленился делать автоопределение размера и налетел, хорошо что проверял начальные строки.
В общем куча ручной работы, но для малого количества паттернов и быстрого счёта (секунды/минуты) проще так чем долго писать универсальную программу.

[Yadryara]

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

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

Тогда тем более жаль, что Вы считали то, что я ранее уже обсчитал.

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

41# считалось 40 минут. Думаю считать 41# для остальных не нужно, там уже младший c3 ненулевой.

Я уже настолько уверен в этих красивых формулах, что полагаю, не надо считать и этот шаг, который делает младший $c$ ненулевым. Назову его последним. Так вот гораздо интереснее предпоследний(-1) и предпредпоследний(-2) шаги. Есть гипотеза, что тот кэф не проваливается ниже 1. Либо ровно 1 досрочно (уже два таких случая), либо чуть больше.

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

Пока что вывод вроде тривиальный: сильно загрязнённые паттерны плохо предсказуемы. Ведь тот паттерн 36-36 содержит $18-3=15$ лишних простых.

[Yadryara]

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

Я уже настолько уверен в этих красивых формулах, что полагаю, не надо считать и этот шаг, который делает младший $c$ ненулевым. Назову его последним.

И однако же я обнаружил несимметричный паттерн 6-6-6-12-6 для которого на одном шаге младший $c$ становится ненулевым, но последний шагом здесь надо называть следующий, ведь именно на нём происходит окончательное устаканивание:

p#     6    7    8    9   10  11      Sum

11#              1    4    4   1       10

13#         2   17   26   21   4       70

17#    2   54  226  282  178  28      770

19#   80 1100 3332 3532 1742 224    10010



17#            11*0+1*2   10*2+2*17    9*17+3*26    8*26+4*21   7*21+5*4  6*4

Otkl. ot tabl.        0           0            5           10        -11   -4


19#           13*2+1*54 12*54+2*226 11*226+3*282 10*282+4*178 9*178+5*28 8*28

Otkl. ot tabl.        0           0            0            0          0    0

[Yadryara]

И поэтому я теперь добавил проверку того самого окончательного устаканивания:

r=[0,0,0,0,0,0,2,17,26,21,4,0];sot=0;

for(i=cmin,cmax,
ot=(f-i)*r[i] + (i-cmin+1)*r[i+1] - vc[i];sot=sot+ot;
print1(ot, "    ");
);
print(, "    ",sot);

Добавил предыдущую строку r, к которой специально приписал ноль справа, чтобы легко посчитать по своим формулам. И теперь строку считаю последней, только если все отклонения (ot) для неё равны нулю.

[Dmitriy40]

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

Тогда тем более жаль, что Вы считали то, что я ранее уже обсчитал.

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

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

я обнаружил несимметричный паттерн 6-6-6-12-6

Не знаю нужно ли, но вот его данные:

17#: 2, 54, 226, 282, 178, 28, sum=770
19#: 80, 1100, 3332, 3532, 1742, 224, sum=10010
23#: 2460, 24264, 60576, 56416, 23766, 2688, sum=170170
29#: 80844, 654960, 1441344, 1223384, 464994, 48384, sum=3913910
31#: 2676060, 18601728, 36821064, 28774424, 10006794, 967680, sum=97847750
37#: 101559588, 631693968, 1154134128, 845711048, 275021838, 25159680, sum=3033280250
41#: 4186279548, 23785863168, 40623559368, 28162840888, 8651475378, 754790400, sum=106164808750

[Yadryara]

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

Впрочем вот сейчас сел и таки сделал выдачу всех строк подряд одним запуском.

Я тоже сделал такую штуку на ПАРИ. И сделал подсчёт не только абсолютных, но и относительных погрешностей. Вот анатомия сходимости для паттернов 24-24 и 30-30:

24-24                                ... 11   12   13

17#                   30   -40   -9   0   9   23   38

19#             65   -12   -13   -4   3  10   14   17

23#        -6   -4    -2    -2   -1   1   4    7   10

29#    0    0    0     0     0    0   0   0    0    0



30-30                                                ... 14   15   16

17#                                 111  -18  -30   11   41  100  300

19#                      1055   59  -12  -13   -1   11   21   38   67

23#             529   88   10   -5   -6   -3    2    7   11   13   14

29#       -32   -24  -12   -7   -4   -2    0    2    3    4    5    8

31#    0    0     0    0    0    0    0    0    0    0    0    0    0

Здесь приведены относительные погрешности в %. Считались по формуле

for(i=cmin+0,cmax,
if(r[i]>0,ot = round((vc[i] - (f-i)*r[i] - (i-cmin+1)*r[i+1])/r[i]/(f-i)*100);

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

[Yadryara]

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

Впрочем вот сейчас сел и таки сделал выдачу всех строк подряд одним запуском.

Я так понимаю, не просто так, а чтобы что-нибудь посчитать.

У меня хорошие новости: паттерн 36-36, на правый край которого я ругался, всё равно демонстрирует хорошую регулярность. Покажу позже.

Так что примите заказ Обсчитать следующий: 42-42.

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

Все 37# считались по минуте, 41# считалось 40 минут.

Здесь для 41# скорость будет примерно такая же? Тогда 43# обсчитывать не надо: верю что сойдётся.

[Yadryara]

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

Но кэф для 36-36 сначала болтануло вниз, потом вверх, а потом он раньше времени сравнялся с 1-цей...

Так это не только правый край болтануло вниз-вверх, а потом он сошёлся, это вся строка на шаг раньше сошлась:

36-36                                                       ... 16    17    18

19#                                   292   32  -28  -14    7   38   105   300

23#                          231   64   6  -10   -7    1   10   18    29    46

29#               698   173   42    7  -3   -4   -1    2    5    6     7     7

31#          -14   -7    -5   -4   -4  -3   -1    0    2    4    7    10    14

37#       0    0    0     0    0    0   0    0    0    0    0    0     0     0

41#   0   0    0    0     0    0    0   0    0    0    0    0    0     0     0

[Dmitriy40]

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

Здесь для 41# скорость будет примерно такая же?

Нет, так как суммы в 1.5 раза больше, то и время выполнения в 1.5 раза больше, практически час.
Держите:

Используется синтаксис Text

v=[0, 42, 84]
13#: 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 4, 98, 354, 674, 572, 198, 18, 2, sum=1920
17#: 0, 0, 0, 0, 0, 0, 0, 0, 0, 8, 260, 1468, 4746, 8240, 7870, 3626, 606, 52, 4, sum=26880
19#: 0, 0, 0, 0, 0, 0, 0, 0, 364, 3636, 20162, 62440, 114794, 126874, 76448, 22242, 2844, 260, 16, sum=430080
23#: 0, 0, 0, 0, 0, 0, 164, 4942, 48572, 268492, 904698, 1890088, 2468138, 1941704, 858204, 192868, 21780, 1854, 96, sum=8601600
29#: 0, 0, 0, 0, 4, 1424, 43518, 525434, 3446946, 13788040, 35010632, 57067894, 59207540, 37700430, 13841974, 2706348, 279322, 21134, 960, sum=223641600
31#: 0, 0, 0, 104, 10584, 336882, 4870920, 39082760, 191844980, 603369458, 1239355838, 1666266362, 1444275908, 778699190, 246690058, 42795174, 4073424, 281638, 11520, sum=6261964800
37#: 0, 0, 462, 52702, 2102242, 39430922, 408369236, 2563767024, 10290807434, 27167143640, 47713932606, 55600328394, 42220299484, 20163731988, 5737487326, 912331278, 81594398, 5228224, 195840, sum=212906803200
41#: 0, 922, 173122, 10114884, 264451498, 3690470360, 30525389506, 159502826194, 546738842952, 1254433960054, 1939829127142, 2010825208458, 1370908155696, 593640303464, 155014292414, 22997963896, 1954975470, 118152928, 4112640, sum=8090458521600

v=[0, 48, 96]
13#: 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 2, 40, 194, 456, 376, 160, 48, 4, sum=1280
17#: 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 4, 156, 874, 3132, 5506, 5000, 2520, 648, 76, 4, sum=17920
19#: 0, 0, 0, 0, 0, 0, 0, 0, 0, 8, 296, 2614, 14266, 43902, 78566, 80344, 47578, 16020, 2860, 256, 10, sum=286720
23#: 0, 0, 0, 0, 0, 0, 0, 4, 256, 4992, 42496, 215418, 676678, 1323616, 1596232, 1179664, 534372, 139596, 19638, 1392, 46, sum=5734400
29#: 0, 0, 0, 0, 0, 4, 120, 4208, 66384, 589706, 3263600, 11480502, 26131526, 38534928, 36530574, 22101296, 8319096, 1835558, 222532, 13952, 414, sum=149094400
31#: 0, 0, 0, 4, 60, 2004, 50680, 801230, 7668400, 46975200, 188854498, 505358980, 904315338, 1076322346, 844430892, 430996368, 138999620, 26808736, 2892062, 162498, 4284, sum=4174643200
37#: 0, 0, 16, 554, 15890, 444100, 7923546, 88915728, 644016702, 3098941600, 10075178142, 22312103678, 33652260394, 34346414458, 23482489766, 10585788910, 3053186642, 534424614, 52960628, 2737588, 65844, sum=141937868800
41#: 0, 32, 2244, 83302, 2773756, 62621704, 882704766, 7953709458, 47306902950, 190299946364, 525224314390, 1000940896542, 1316531453644, 1188026343384, 727593074588, 296905718420, 78125698774, 12570810480, 1154941526, 55767040, 1251036, sum=5393639014400

v=[0, 54, 108]
13#: 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 14, 36, 182, 382, 432, 194, 40, sum=1280
17#: 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 44, 226, 1058, 2790, 5378, 5100, 2704, 592, 28, sum=17920
19#: 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 64, 776, 5152, 18654, 48128, 79516, 76276, 42876, 13622, 1604, 52, sum=286720
23#: 0, 0, 0, 0, 0, 0, 0, 0, 0, 54, 1302, 14194, 84998, 319536, 823744, 1421412, 1532780, 1021014, 416522, 90914, 7730, 200, sum=5734400
29#: 0, 0, 0, 0, 0, 0, 0, 10, 944, 20392, 220088, 1377662, 5621650, 15777688, 30391356, 39085790, 32552426, 17328896, 5655910, 989160, 70866, 1562, sum=149094400
31#: 0, 0, 0, 0, 0, 4, 606, 20192, 342682, 3415458, 22034594, 96954030, 300013282, 654759532, 990695942, 1015196012, 690608264, 304873962, 82768206, 12181790, 763376, 15268, sum=4174643200
37#: 0, 0, 0, 0, 56, 6204, 219006, 4056022, 45097598, 328919636, 1652902232, 5884950794, 15008681762, 27251522652, 34665601096, 30329057478, 17910754396, 6948229422, 1673075024, 221717364, 12837526, 240532, sum=141937868800

v=[0, 60, 120]
13#: 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 2, 2, 32, 232, 580, 894, 612, 192, 14, sum=2560
17#: 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 8, 32, 208, 1574, 4796, 9314, 10604, 6858, 2210, 230, 6, sum=35840
19#: 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 28, 168, 1108, 6984, 32010, 84474, 144866, 157638, 102764, 36704, 6298, 388, 10, sum=573440
23#: 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 54, 532, 3742, 24324, 147744, 591646, 1546426, 2662502, 3034322, 2203656, 978518, 245040, 29106, 1172, 16, sum=11468800
29#: 0, 0, 0, 0, 0, 0, 0, 0, 0, 518, 7096, 64538, 461904, 2694512, 11083696, 30885334, 58486846, 75566800, 65660220, 37209254, 13129424, 2660844, 267486, 10188, 140, sum=298188800
31#: 0, 0, 0, 0, 0, 0, 0, 582, 11318, 129972, 1158066, 8530542, 49288390, 209325666, 628036846, 1320060384, 1940132792, 1977332136, 1371759468, 628824150, 181590134, 30455796, 2565850, 83356, 952, sum=8349286400
37#: 0, 0, 0, 0, 0, 98, 6234, 133122, 1731852, 17440694, 140181554, 868027882, 3998950130, 13377997688, 32148384894, 55276476042, 67618730748, 58169237868, 34514725416, 13709237126, 3474995856, 518683614, 39593124, 1191282, 12376, sum=283875737600

Повезло что у 48-48 тоже c5>0 уже при 37#, потому для него досчитал и 41# за 40 минут.
Считать 43# за 30-40 часов влом.
Как и остальные, там нет надежды получить даже c5>0.