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

[Dmitriy40]

А, так да, для далёкого асм+AVX+дельфи (правда пока она считает только vc[] длиной до 12, но как нарастить понятно, просто будет немного медленнее), хотя до 1e10-1e11 вполне реально и на PARI, это всего час или часов 12, для близкого асм или PARI (разные алгоритмы и считают по разному).

Первая программа на PARI, которая вполне может посчитать vc[] до 1e11 за разумное время (где-то выше уже показывал):

? vc=[3075976, 172443036, 2982765864, 22582804840, 87976549320, 195008570124, 262408889876, 220882584188, 117236333144, 40358558480, 9629086752, 1478049120, 90478080]; cm=9;
? for(k=1,8, c=vc*1.0; forprime(p=47,10^k, for(i=1,#c-1, c[i]=c[i]*(p-i-cm+1)+c[i+1]*i); c[#c]*=p-#c-cm+1); print("10^",k,": c=",c[1]);)
10^1: c=3075976.0000000000000000000000000000000
10^2: c=9700723135938319098329396524.0000000000
10^3: c=1.0376559211320620230155745256956034679 E407
10^4: c=9.0253029272976793965590757856321081266 E4288
10^5: c=2.6851992474995041923249653671105627727 E43283
10^6: c=1.8136762603118047528812914570381931279 E433626
10^7: c=6.9304824582518476787558725783564086233 E4340840
10^8: c=3.9260010206257956280639458393354973103 E43424108

[Yadryara]

Эх, а я-то радовался, что у нас нет путаницы с обозначениями.

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

хотя до 1e10-1e11 вполне реально и на PARI, это всего час или часов 12

Да нет же, это 1e20-1e22.

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

10^6: c=1.8136762603118047528812914570381931279 E433626
10^7: c=6.9304824582518476787558725783564086233 E4340840
10^8: c=3.9260010206257956280639458393354973103 E43424108

И она скачет через порядок: 1е12, 1е14, 1е16.

А надо бы на каждом порядке сравнение делать в надежде увидеть ту самую линейность. Тогда прогноз может быть жутко точным.

Вот уже немало 9-к и 11-к ждут обсчёта по формулам. И не так важно, какое превышение, лишь бы регулярность сохранялась. Разделим на кэф превышения и вуаля

[Dmitriy40]

9-ки есть до 1e16, это на PARI считать до 1e8, хватает минуты. Ну пусть пары-тройки минут если комп не быстрый. Думаю Вы и сами легко справитесь.
11-ки есть всего две, 9-132, их надо считать до $\sqrt{2.148\cdot10^{18}}=1.46\cdot10^9$, хватит получаса, тоже справитесь ведь.
Что там у Вас считалось сутки до 1e15 я не понял, PARI до 3.2e7# (т.е. 1e15) считает секунды, я для этого и показал снова реальный код, у меня он 1e8# (т.е. 1e16) считает 35с, у Вас пусть даже пару минут, это не сутки.
Ну вот как считает оба 11-132 сразу до конца:

? vc=[27819218530, 3205709529406, 132518657338014, 2726047695649162, 32359996014192550, 242025952481691690, 1206902329923127624, 4154854508647781376, 10048123918875325956, 17137545399442822324, 20490408641203428008, 16943429732449871952, 9484928332352457234, 3481588777782070874, 799760667043630092, 108168644430424952, 8222129621109552, 377531049332304, 11794935398400];
? cm=11; c=vc*1.0; k=1; forprime(p=68,sqrt(2.148e18), while(p^2>10^k, print("1e",k,": c=",c[1]); k++); for(i=1,#c-1, c[i]=c[i]*(p-i-cm+1)+c[i+1]*i); c[#c]*=p-#c-cm+1); print("2.148e18: c=",c[1]);
1e1: c=27819218530.000000000000000000000000000
1e2: c=27819218530.000000000000000000000000000
1e3: c=27819218530.000000000000000000000000000
1e4: c=306276338750203471406050.00000000000000
1e5: c=4.1090018069442649589861921642232028517 E116
1e6: c=2.6731494395382713524633468425181707088 E404
1e7: c=2.4336655210292254413082792310556668517 E1318
1e8: c=4.9583091736421711248612059433751748828 E4286
1e9: c=2.1211219820805541819435967885456180889 E13609
1e10: c=1.7999028457677614457796673455706212189 E43281
1e11: c=6.7951639493874492336320029076899529076 E136975
1e12: c=1.2404711191261473064124537454839465774 E433624
1e13: c=8.2854214085890781014101543313386788662 E1372327
1e14: c=4.4999599935406569012035752366277858300 E4340838
1e15: c=6.2626560085586352255806996285256878171 E13731275
1e16: c=2.3441380791500246537126873145322601837 E43424106
1e17: c=1.1392201268430852366222174630489810356 E137328721
1e18: c=1.0094068176926317770949690894110851561 E434280996
2.148e18: c=1.7934002489110223173819298073458894157 E636486311
time = 10min, 45,891 ms.

? vc=[7501811456, 1208970989984, 65169226165182, 1669289597490466, 23733434231630646, 204487221464295138, 1128253605829605168, 4134296731958743400, 10309800229759011216, 17779901928236428416, 21415991524970872398, 18104049561984185730, 10734599774458772586, 4425573782555262318, 1239633885066283956, 224417971568292884, 23494744611658872, 1038094551860184];
? cm=11; c=vc*1.0; k=1; forprime(p=68,sqrt(2.148e18), while(p^2>10^k, print("1e",k,": c=",c[1]); k++); for(i=1,#c-1, c[i]=c[i]*(p-i-cm+1)+c[i+1]*i); c[#c]*=p-#c-cm+1); print("2.148e18: c=",c[1]);
1e1: c=7501811456.0000000000000000000000000000
1e2: c=7501811456.0000000000000000000000000000
1e3: c=7501811456.0000000000000000000000000000
1e4: c=177213447296352897165548.00000000000000
1e5: c=3.7431137777060540843747418853695741096 E116
1e6: c=2.6054019299995507438134781223868701018 E404
1e7: c=2.4316133322517102971134910164903371443 E1318
1e8: c=5.0190993074846403736739386397600644464 E4286
1e9: c=2.1644430424380707115051866438516420445 E13609
1e10: c=1.8470518288493454101527931216534526674 E43281
1e11: c=7.0025906225103140909108507455943920748 E136975
1e12: c=1.2825721089396853198899446922808888365 E433624
1e13: c=8.5896272834846896816901172268658130869 E1372327
1e14: c=4.6755828187754578990983109106040234671 E4340838
1e15: c=6.5193587123589648305540178133655981722 E13731275
1e16: c=2.4441820381085603974137516487900334177 E43424106
1e17: c=1.1895170944470854939804130066522257885 E137328721
1e18: c=1.0552804581046507644559826027776230802 E434280996
2.148e18: c=1.8756150153451242307891129894634751348 E636486311
time = 10min, 12,507 ms.

До 1e14 оно считает 6с, 18с до 1e15, 53с до 1e16с, 160с до 1e17.

Отдельно отмечу что отношение sum второго к первому, равное 16/15=6.7%, для конечного результата уменьшилось до 1.8756/1.7934=1.046 или 4.6%. Если же второй результат поделить на отношение sum равное 16/15 получим 1.758389, что отличается от точного значения первого всего на 0.3%.

-- 11.04.2024, 15:06 --

Если нужны знаменатели, то они считаются ещё сильно быстрее и только один раз:

? s=1.0; k=1; forprime(p=1,sqrt(2.148e18), while(p^2>10^k, print("1e",k,": √p#=",s); k++); s*=p; ); print("2.148e18: √p#=",s);
1e1: √p#=6.0000000000000000000000000000000000000
1e2: √p#=210.00000000000000000000000000000000000
1e3: √p#=200560490130.00000000000000000000000000
1e4: √p#=2305567963945518424753102147331756070.0
1e5: √p#=6.1076929465933196099278943388997855150 E127
1e6: √p#=1.9590340644999083431262508198206381046 E415
1e7: √p#=1.9522882315128211808496464634974874762 E1329
1e8: √p#=5.9490679579998635868564022242503445659 E4297
1e9: √p#=4.3052353855950690494878129398668630681 E13620
1e10: √p#=6.5235761466597667813712621729206200503 E43292
1e11: √p#=4.4587461565944708330737038375609650291 E136987
1e12: √p#=1.4705552760839308761649229845696490594 E433636
1e13: √p#=1.7536647295156189735811651573073679657 E1372340
1e14: √p#=1.6725476714034538506990478976102144813 E4340851
1e15: √p#=4.0129671196962765407337488720246638814 E13731288
1e16: √p#=2.5401483842215342462788890842604146196 E43424119
1e17: √p#=2.0482350313774922452835522308588480309 E137328734
1e18: √p#=2.9557521529413355924858770497467279652 E434281009
2.148e18: √p#=6.1501327818293311097988392787319372235 E636486324
time = 26,802 ms.

[Yadryara]

Благодарю.

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

До 1e14 оно считает 6с, 18с до 1e15, 53с до 1e16с, 160с до 1e17.

Хммм... Действительно, эта прога и у меня за минуты считает!

А почему же в тот раз так долго было... Буду разбираться.

[Yadryara]

Я занимался постройкой

gris в сообщении #1197339 писал(а):

арифметического моста

Но одного моста мне было недостаточно и я построил два.

На берегу речки стоят два одинаковых числа: 1485. Это количество св-близнецов до 17#. Как я их считал, показано вот здесь. Платой за переход по мосту является существенное уменьшение чисел до 343 (для паттерна 3-36) и до 132 (для 5-36). Причина уменьшения — запреты по модулям.

      3-36 13#-->17#                                 5-36 13#-->17#

      1485                                           1485

  u% 5 == 3                                     u% 5 == 3 &&
                                                u% 7 <> 3 && u%7 <> 5 && 
                                                u%11 <> 5 &&
                                                u%13 <> 9

       495                                            240
-             u%7==2 && (u%11==4 || u%11==8)  + 3    
-             u%7==5 && (u%11==4 || u%11==8)
-             u%7==3 && (u%11==2 || u%11==6)
-             u%7==6 && (u%11==2 || u%11==6)  + 7
      4*22                                           2*30
-             u%7==2 && (u%13==3 || u%13==8)  + 7
-             u%7==5 && (u%13==3 || u%13==8)
-             u%7==3 && (u%13==2 || u%13==10)
-             u%7==6 && (u%13==2 || u%13==10) +11
      4*14                                           2*15
-            u%11==2 && (u%13==3 || u%13==8)  + 7
-            u%11==6 && (u%13==3 || u%13==8)  + 7
             u%11==7 && (u%13==3 || u%13==8)        -
             u%11==3 && (u%13==2 || u%13==10)       -
-            u%11==4 && (u%13==2 || u%13==10) +11
-            u%11==8 && (u%13==2 || u%13==10) +11
      4* 2                                           6* 2

       343           (7*7*7)     (11*12)              132

        35                                              8

c1g1   378                                            140

Затем 35 и 8 св-близнецов соответственно выбивают по c1g1 сразу по два кортежа каждый.

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

Не знаю, надо ли подробно объяснять схему, может это никому не интересно. А вопрос такой:

Можно ли другие, более простые мосты построить для тех же чисел?

[Yadryara]

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

Действительно, эта прога и у меня за минуты считает!

Но, правда, после 1е17 переполнение показывает и дальше не считает:

overflow in expo()

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

А почему же в тот раз так долго было... Буду разбираться.

Кстати, не сутки, а примерно 14 часов было, судя по постам.

Вот, нашёл ту прогу.

(PARI)

{print();
f = [0, 0, 0, 999, 999, 999, 999, 999, 999, 8, 20, 73, 403, 2315, 13514, 81589];

mor0=2*3*5*7*11*13*17*19*23*29*31*37*41;

stfin=16;

c0=[0,0,0,0,0,0,0,0,25928, 2194424, 50013522, 460777720, 2074673880, 5147266584, 7622973198, 6995091476, 4001889292, 1465245728, 372741528, 62117280, 4112640, 0];

s=vecsum(c0);

cm=9;

for(k=4,stfin,

c=c0;mor=mor0;

print();

forprime(p=43,sqrt(10^k),

mor=mor*p;

for(i=cm,#c-1,c[i]=c[i]*(p-i)+c[i+1]*(i+1-cm);););

print("10^",k,": c9/sum = ",round(c[cm]/vecsum(c[cm..#c-1])*10^20),"   ",round(c[cm]*10^k*10^3/mor/f[k])););
print();

}quit;

Неужто это из-за c=vc*1.0 которого здесь нет?

[Yadryara]

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

Неужто это из-за c=vc*1.0 которого здесь нет?

А ведь у меня ещё и подсчёт огромного праймориала. Понял, надо сделать mor=mor0*1.0; и тогда считается в сотни раз быстрее.

[gris]

"и немедленно сэкспериментировал"
(08:45) gp > m=3; for( i=1,10^8,m0=m*1.0 );
time = 1min, 4,678 ms.
(08:46) gp > m=3; for( i=1,10^8,m0=m+0.0 );
time = 32,465 ms.
Ускорение, однако!

[Yadryara]

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

Ну вот как считает оба 11-132 сразу до конца:

Дык, а какой из них какой?

Я записываю их по возрастанию (лексикографически) :

1. [0, 6, 12, 36, 42, 66, 90, 96, 120, 126, 132]

2. [0, 6, 30, 42, 60, 66, 72, 90, 102, 126, 132]

Подозреваю, что у Вас наоборот: тот, который 0, 6, 12, — внизу.

Пришлось в разных вариантах считать. Вот, для 0, 6, 12, из Базы Томаша выудил:

10^   13  14  15   16    17   18   end
       1, 16, 79, 491, 2750

gris, на Ерофеевское намёк?

Вы лучше скажите, у Вас 32 разряда? И для той же проги тоже после 1е17 переполнение?

Арифмосты через речки строить будем?

[Dmitriy40]

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

Неужто это из-за c=vc*1.0 которого здесь нет?

Именно! Без этого PARI пытается считать офигенно длинные числа, они же изначально целые, а это разумеется страшно медленно.

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

Дык, а какой из них какой?

Ой, я разве их не показывал раньше? Забыл. Вот:

v=[0, 6, 12, 36, 42, 66, 90, 96, 120, 126, 132]+
67#: 27819218530, 3205709529406, 132518657338014, 2726047695649162, 32359996014192550, 242025952481691690, 1206902329923127624, 4154854508647781376, 10048123918875325956, 17137545399442822324, 20490408641203428008, 16943429732449871952, 9484928332352457234, 3481588777782070874, 799760667043630092, 108168644430424952, 8222129621109552, 377531049332304, 11794935398400, sum=84141570156134400000, OK

v=[0, 6, 30, 42, 60, 66, 72, 90, 102, 126, 132]+
67#: 7501811456, 1208970989984, 65169226165182, 1669289597490466, 23733434231630646, 204487221464295138, 1128253605829605168, 4134296731958743400, 10309800229759011216, 17779901928236428416, 21415991524970872398, 18104049561984185730, 10734599774458772586, 4425573782555262318, 1239633885066283956, 224417971568292884, 23494744611658872, 1038094551860184, sum=89751008166543360000, OK

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

Но, правда, после 1е17 переполнение показывает и дальше не считает:

Беда. У меня тоже gp32 дальше не считает. Даже просто p#. Значит gp32 не умеет работать с показателем степени больше примерно 161.6млн. Надо в программу добавлять проверку на 1e100000000 и при превышении делить на это число, подсчитывая реальную степень отдельно, как Вы уже и делали.

-- 12.04.2024, 19:58 --

gris в сообщении #1636136 писал(а):

"и немедленно сэкспериментировал"

Это вообще не то, тут числа в цикле не удлиняются на каждой итерации. Вы попробуйте вот это запустить:

? m=1;for(i=1,10^5,m=m*i);
time = 2,544 ms.
? #m
%2 = 47398
? sizebyte(m)
%3 = 189600
? m=1;for(i=1,3*10^5,m=m*i);
time = 26,832 ms.
? #m
%5 = 157050
? sizebyte(m)
%6 = 628208
? m=1.0;for(i=1,3*10^5,m=m*i);
time = 63 ms.
? #m
%8 = 3
? sizebyte(m)
%9 = 20
? m=1.0;for(i=1,10^7,m=m*i);
time = 1,888 ms.

Yadryara
Видите разницу между целыми числами и плавающими и в размерах и во времени счёта ...
Так что да, *1.0 у меня стоит совсем не просто так.

[Dmitriy40]

За примерно день посчитал все 7-120 (9-120 кажется уже были посчитаны ранее), по 1.5ч-2ч на каждый, на PARI.

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

v=[0, 6, 30, 36, 60, 84, 90, 114, 120]+
41#: 0, 0, 0, 3884, 254400, 6057410, 70286424, 460499340, 1882956650, 5133980902, 9667853048, 12651291142, 11312882740, 6676339958, 2452840186, 503916712, 46356184, 902780, sum=50866421760
43#: 0, 0, 12048, 1166020, 38303884, 595065412, 5110017650, 26822994836, 92403044258, 218864609474, 364252402918, 425312060232, 341463110756, 181789244604, 60501937274, 11328007710, 958751484, 17611280, sum=1729458339840
47#: 0, 24024, 3891012, 191538778, 4224377970, 49738744532, 347778633408, 1555640307808, 4695575203268, 9896428412632, 14775538416876, 15561221886192, 11330459135934, 5504523556562, 1682901464378, 291687570712, 23093708634, 410041200, sum=65719416913920
53#: 21878, 8206116, 695877266, 23527477292, 399589017576, 3875879257306, 23350698443590, 92759472392704, 253465877770904, 488882866535924, 671657372179712, 652886336603722, 439857461378006, 198253194355232, 56410991791374, 9138699965712, 680209513846, 11463424320, sum=2891654344212480
59#: 9227860, 1771770278, 102512426302, 2666235467552, 37283898115152, 311194324377394, 1662793260040910, 5979245255978260, 14981647069132772, 26694474128938812, 34025387071165894, 30785176645106836, 19369453875911290, 8182848434444088, 2190820283942890, 335475074846556, 23757603304274, 389756426880, sum=144582717210624000
61#: 2251618998, 295385136782, 13124327717756, 279781130370656, 3345598731414266, 24602892805982978, 118343206753729680, 388919213072083876, 899442738202291276, 1488116258256027856, 1767703200085142744, 1494625688960315756, 881155184684224744, 349802572918519892, 88283296912528160, 12792699422190956, 861899578210824, 13641474940800, sum=7518301294952448000, OK

v=[0, 12, 18, 42, 60, 78, 102, 108, 120]+
37#: 0, 0, 0, 66, 3588, 93996, 1504854, 14479464, 82446024, 286013266, 624069190, 872144340, 784311990, 450661934, 160383356, 32417722, 3023630, 62580, sum=3311616000
41#: 0, 0, 340, 22540, 667838, 12157472, 139037474, 986209720, 4375270630, 12396867958, 22848209306, 27643785338, 21942763190, 11283818970, 3621989620, 663479038, 56357954, 1074612, sum=105971712000
43#: 0, 736, 84938, 3530404, 81170880, 1159959532, 10451254426, 59700837592, 219447114530, 528332165198, 844399067102, 900916327726, 639359751126, 296913607198, 86623129574, 14496555514, 1133274888, 20376636, sum=3603038208000
47#: 656, 179698, 12629476, 418688534, 8048308732, 95572325952, 717593718640, 3471957281466, 11030948100038, 23381676985236, 33397939449154, 32229492653028, 20879135039090, 8906971639390, 2396691854026, 371457038556, 27073881540, 462130788, sum=136915451904000
53#: 204796, 32942850, 1782146640, 48815415516, 785215611960, 7865176863802, 50568721119268, 213351286901698, 601235609131834, 1146467603840024, 1489860684379960, 1319771510762752, 790324149834276, 313228065487638, 78582013098166, 11394835900260, 781492636748, 12887497812, sum=6024279883776000
59#: 42275066, 5083935484, 228287770066, 5360191034408, 74564459057386, 651035277418524, 3690441004904794, 13907174614166278, 35446771441946486, 61809929329458732, 74133245328272790, 61070338346307116, 34217081884825192, 12748168896387958, 3017752629847492, 414415335608292, 27049868862488, 432166721448, sum=301213994188800000
61#: 7282238916, 715856249816, 27494961606524, 560907196915536, 6834270421847148, 52741304068099392, 267110508524784470, 909397029173054398, 2115947307410773972, 3399159414449453376, 3785370025596835456, 2914488854816494060, 1534409471046051134, 539427123776995250, 120890829968329076, 15766165319306612, 981142113314184, 15125835250680, sum=15663127697817600000, OK

v=[0, 18, 30, 42, 60, 78, 90, 102, 120]+
37#: 0, 0, 0, 64, 3908, 110060, 1505060, 11456436, 53377420, 160584166, 320635374, 428538820, 379325906, 214780112, 72150104, 12444046, 878028, 18496, sum=1655808000
41#: 0, 0, 118, 12908, 531824, 10285452, 106998388, 660965484, 2586474990, 6663805902, 11542581682, 13526305410, 10608265610, 5374936218, 1630928884, 256824538, 16620000, 318592, sum=52985856000
43#: 0, 236, 40966, 2429572, 65179418, 913133642, 7413869668, 37539341444, 124354560304, 277343156152, 422323427842, 439433556522, 307883643902, 139980481914, 38375000520, 5546613310, 338371628, 6296960, sum=1801519104000
47#: 222, 84496, 8223952, 328649622, 6522998704, 72408717518, 488138008794, 2115961461968, 6126581922454, 12126115239624, 16578329298670, 15621566393526, 9973046234478, 4151612509828, 1048069366930, 140813530314, 8079903540, 143407360, sum=68457725952000
53#: 103270, 22191802, 1439048796, 41675215678, 642538205222, 5844650144770, 33551193832942, 127227672531914, 328433205824906, 587477640407504, 733119699518624, 635476848698924, 375656759801064, 145645880673554, 34429109688232, 4352328951502, 235373750736, 3903298560, sum=3012139941888000
59#: 26835248, 3871511328, 189152355250, 4414911295942, 57407474665710, 454824042577394, 2327235117500684, 7996972694162854, 18927869286331240, 31304094934634908, 36339142651896954, 29432989746180414, 16314684551875344, 5950319591794448, 1328681845524962, 159741360474360, 8293122631920, 132712151040, sum=150606997094400000
61#: 5266944224, 575751788228, 22702351650326, 445960552163998, 5029678996841050, 35340140706141622, 163031624264171442, 511286725527978350, 1114563103010288732, 1709467508708270584, 1850006878587656622, 1402528794215901102, 729941536768341584, 250664009917332940, 52886030537063956, 6043120299662040, 300808521316800, 4644925286400, sum=7831563848908800000, OK

v=[0, 18, 30, 48, 60, 72, 90, 102, 120]+
37#: 0, 0, 0, 512, 16392, 221438, 1778534, 9732164, 38240914, 107025546, 208294674, 279555518, 254085694, 147075988, 49068860, 8246490, 529276, sum=1103872000
41#: 0, 0, 1174, 65858, 1356586, 15007918, 106137900, 520593618, 1808930164, 4420353194, 7531271328, 8893356094, 7115057822, 3643796622, 1090985076, 167142552, 9848094, sum=35323904000
43#: 0, 3228, 294838, 8561830, 123545500, 1089850266, 6533950358, 27830964764, 84708735548, 182555551352, 276267711816, 290730075402, 207210854244, 94718186304, 25497749390, 3545477366, 191223794, sum=1201012736000
47#: 2904, 648398, 33715948, 748216892, 9213943396, 72557868848, 393226057150, 1509312339904, 4120071238706, 7970303887422, 10877280168014, 10365791579278, 6716076613312, 2809387651936, 699086996148, 90783611202, 4609428542, sum=45638483968000
53#: 776174, 95388688, 3669301552, 67816454222, 735299586202, 5214774454328, 25575073924746, 88750175212280, 219295051560218, 385917846533000, 482153292141830, 422854282090886, 253225501177800, 98427845656204, 22926816196028, 2810198902112, 135555235730, sum=2008093294592000
59#: 134021822, 11972337064, 377160139504, 6067962804618, 59122367148216, 382533669592984, 1722767204790150, 5509596757933500, 12592297523101770, 20579228399763000, 23955957058144870, 19616211166691478, 10987278335090442, 4006466320011648, 879774953483212, 102267850763972, 4705893781750, sum=100404664729600000
61#: 18941471808, 1364909469272, 37061895389054, 533819646019146, 4750541971079288, 28315685699611148, 117814468725881400, 348670234291821660, 739274146614344880, 1124466391771257700, 1221928519275690798, 936111997855435902, 491575195563769104, 168569035829219240, 34965465993821636, 3859204778774964, 169412176143000, sum=5221042565939200000, OK

v=[0, 18, 42, 48, 60, 72, 78, 102, 120]+
41#: 0, 0, 0, 758, 89326, 3091248, 48566340, 414105864, 2088594974, 6471677666, 12460760032, 14849009836, 10769914192, 4612709288, 1116014794, 142693268, 8472894, 155520, sum=52985856000
43#: 0, 0, 4390, 530378, 21487376, 412806252, 4378469698, 27836251608, 110824126266, 281537034032, 457439437096, 471174259286, 301884510368, 116685618830, 26021766254, 3123992950, 175854336, 2954880, sum=1801519104000
47#: 16, 17344, 2430596, 122662020, 3007686390, 41061297390, 336840342136, 1738592109294, 5794508019544, 12606936033052, 17892583690224, 16393250487312, 9499388887826, 3370438448122, 698245584264, 78526737802, 4156511308, 65007360, sum=68457725952000
53#: 19926, 6126122, 497825984, 17628699880, 331000689394, 3644135419834, 24957444254544, 110219224886608, 320268297765740, 617225466011100, 787239816688522, 656799233642490, 350663643068320, 115776320110616, 22490860515166, 2384931265750, 119614691924, 1820206080, sum=3012139941888000
59#: 6912060, 1260891928, 74923996110, 2107029856718, 32849389481952, 309187357781636, 1848802297041018, 7243464579416510, 18913017041862990, 33111816649867794, 38741363913501246, 29917884731045294, 14910186933085454, 4632395455819022, 853254783553430, 86366684673002, 4162168607116, 61887006720, sum=150606997094400000
61#: 1620319048, 214153480548, 10067289375654, 234642020906990, 3122707484041876, 25624619597983000, 135749157719802398, 477260042408646870, 1130179099690781706, 1811221755079327602, 1956234016408550566, 1405555517169882502, 656628618249065446, 192608989746689878, 33719182045125370, 3262162030614930, 150890148970416, 2166045235200, sum=7831563848908800000, OK

v=[0, 24, 36, 54, 60, 66, 84, 96, 120]+
37#: 0, 0, 0, 4, 1484, 71126, 1228168, 10799116, 55788012, 179271254, 364914810, 474541240, 394432356, 206648150, 65855392, 11697802, 926630, 19656, sum=1766195200
41#: 0, 0, 20, 8902, 462082, 9752738, 107472574, 699234748, 2859951790, 7554333962, 13019802034, 14686682898, 10788687944, 5069711064, 1465814676, 238462420, 17514740, 353808, sum=56518246400
43#: 0, 28, 24996, 2036226, 61149276, 909361638, 7733022398, 40773359164, 139134689790, 313456751586, 469889453820, 468873711022, 308889250120, 131583913098, 34755224388, 5197097370, 354610328, 6722352, sum=1921620377600
47#: 28, 48674, 6687922, 302867184, 6391586192, 74119028208, 519377685556, 2327177947328, 6882357098496, 13649139287324, 18262619830792, 16460074226588, 9889071256520, 3869354858662, 942776181322, 130367730740, 8290135520, 147891744, sum=73021574348800
53#: 44980, 14200602, 1114459386, 36176814280, 603856557120, 5844075723778, 35279784837864, 138953830411368, 366524708879368, 656279349511564, 801427828304988, 665523489561984, 371350796365312, 135792162598620, 31045962274826, 4039034350280, 242945482048, 4140968832, sum=3212949271347200
59#: 16115632, 2870142824, 159359377744, 4058097297904, 56342293412212, 470368035896256, 2508014476785566, 8866473177300458, 21244559769505060, 34890227812538544, 39401422876920202, 30468324974667602, 15922187485195112, 5478321618848040, 1182733300139724, 145777668236576, 8353083009088, 136651971456, sum=160647463567360000
61#: 3708155688, 465096039512, 20142260780912, 424215941246144, 5056270263267456, 37155384547837428, 177433978173239242, 568947771134561090, 1248772680171069536, 1894294024708359412, 1990011335551992106, 1440267573783713026, 708105680452829000, 230212809337029696, 47130530428858152, 5527423052898720, 303034071841920, 4782819000960, sum=8353668105502720000, OK

v=[0, 30, 42, 48, 60, 72, 78, 90, 120]+
41#: 0, 0, 0, 196, 24338, 1018420, 19344982, 193530418, 1105342826, 3794936678, 8123881914, 11119446278, 9867280834, 5658712102, 2034982700, 425719676, 43205238, 1258200, sum=42388684800
43#: 0, 0, 820, 126340, 6468272, 152183580, 1914631976, 13962609542, 62011376688, 173502693398, 313390394548, 371845153058, 292435647196, 151423625302, 49927904564, 9693451376, 923457152, 25559388, sum=1441215283200
47#: 0, 3502, 671638, 41146492, 1181151638, 18421948736, 168778060154, 954290778594, 3444573607304, 8130754189922, 12767003270718, 13480524349240, 9594781689226, 4548350253956, 1384670399886, 249925044630, 22299658204, 584537760, sum=54766180761600
53#: 4194, 1879094, 185124298, 7426631322, 153489214112, 1831830345904, 13462447815800, 63434157187262, 196887233904736, 410040527882562, 579892941380892, 559972045832950, 368036673427898, 161864212022404, 45789223802102, 7684132119182, 639541875728, 15883059960, sum=2409711953510400
59#: 2035538, 446441684, 30035710876, 929621134884, 15715069236782, 158891334857612, 1014395132471302, 4234767977889966, 11834714116113248, 22490524646947290, 29341066099921510, 26360696169185082, 16222351094330848, 6711521532260766, 1792203878201920, 284725420808444, 22527186330664, 537911641584, sum=120485597675520000
61#: 552289660, 82839947636, 4290648948452, 108411712556444, 1548779997653596, 13554263533135576, 76305551938909654, 285242271933954454, 723142142931508522, 1260503220817948570, 1522292434057739322, 1275456756068558538, 736143823692623878, 286840194052996754, 72374628683799620, 10895275551203052, 820123205810832, 18826907455440, sum=6265251079127040000, OK

А вот 5-120 на PARI не считаются, памяти 7ГБ мало, надо порядка 30ГБ (очень навскидку). На асме же они считаются по несколько суток каждый (в 4 потока!). Пока вот один из них будет считаться 53# часов 14, значит 59# с неделю, и ещё недели полторы-две 61# (если его считать). Паттернов 9шт, это больше 2 месяцев даже без 61#, долговато.

[Yadryara]

Ну вот наконец-то статистика по трём 9-кам с самыми маленькими диаметрами.

Столбцы:

Степень десятки, кэф превышения, доля чистых среди всех.

[0, 12, 24, 30, 42, 54, 60, 72, 84]

12   1.689    15.2 %

13   1.749    17.8 %

14   1.790    20.4 %

15   1.850    22.9 %

16   1.894    25.3 %


[0, 6, 18, 36, 48, 60, 78, 90, 96]

12   1.829     9.0 %

13   1.705    11.1 %

14   1.710    13.2 %

15   1.754    15.3 %

16   1.826    17.5 %


[0, 18, 24, 48, 54, 60, 84, 90, 108] Центральная 9-ка

12   1.559     5.8 %

13   1.614     7.5 %

14   1.709     9.2 %

15   1.752    11.0 %

16   1.788    12.8 %

Как видим, прекрасная линейность кэфа, наблюдающаяся для паттерна 9-84, куда-то подевалась, но тенденция его роста с ростом диапазона сохранилась. А вот с ростом диаметра, наоборот, кэф снижается, но растёт с ростом длины кортежей. Для 11-132 он превысил 2.

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

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

9-120 кажется уже были посчитаны ранее

Было объявлено об этом. Но ни стартовых, ни окейных строк предъявлено не было. А жаль. Может я бы их уже обсчитал и сравнил. И таких 9-к в таблице было бы не 3, а 10.

[Dmitriy40]

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

Было объявлено об этом. Но ни стартовых, ни окейных строк предъявлено не было.

Извиняюсь: постоянно правлю свой файлик с полным списком, не помню что показал а что нет, выкладывать его в прошлое сообщение возможности нет, выкладывать каждый день по новой ... лень. Думал в облаке обновлять, видимо так удобнее всего, вот текущий: https://cloud.mail.ru/public/JEZx/DFCKzMBwi Обновляться будет по этой же ссылке, через пару часов будет 53# для одного из 5-120.

[Yadryara]

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

Обновляться будет по этой же ссылке, через пару часов будет 53# для одного из 5-120.

Ура! И всё чётко ведь, начиная с 3#. Вдруг пригодится.

Что делать с 3-ми, 5-ми и 7-ми всё равно пока непонятно.

Самые ценные пока — 9-ки (7.4 млн) и 11-ки (9.0 млн). 13-ки (даже если будут обсчитаны) вроде пока не добрались до 0.4 млн.

С Базой 9-к еле разобрался.

А вот Базу Томаша не знаю как проверять, кроме как копипастить каждую из 46 страниц 11-к.

[Dmitriy40]

53# посчиталось, файл обновил.

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

А вот Базу Томаша не знаю как проверять, кроме как копипастить каждую из 46 страниц 11-к.

Ну да, а как ещё. Но можно сделать один раз, потом пользоваться локальной.
Чтобы не делать всё руками я использую cURL: качаем архив, берём из него один curl.exe, запускаем в виде
curl -k "https://boinc.tbrada.eu/spt/tuples.php?spt=11&p=[1-46]" -o spt11.p#1
Она сама скачает все указанные в [] страницы и запишет их в файлы. Потом достаточно просто объединить все 46 файлов в один, чем угодно, хотя бы даже командой DOS type spt11.p* 2>nul 1>spt11 и потом спокойно им пользоваться. В нём правда переводы строк будут Unix like (только 0x0A вместо 0x0D,0x0A под виндой), но и findstr и PARI кажется его вполне понимают, хотя в блокноте винды будет страх и ужас (одна длинная строка) (как под Win10+ я не проверял).
Впрочем, выложил Вам n11 (SPT11) от Томаша в облако (уже с правильными переносами строк): https://cloud.mail.ru/public/Yx9B/8ZVH5wCVY - 517МБ текста! Скачаете, как вспомню - уберу.