Message boards :
Cafe :
К штурму 19-ки
Message board moderation
Previous · 1 . . . 15 · 16 · 17 · 18 · 19 · 20 · 21 . . . 55 · Next
Author | Message |
---|---|
Send message Joined: 11 Jul 17 Posts: 174 Credit: 4,964,801 RAC: 1 |
Примерная статистика по наполненности паттернов: for tuple diameter=360 number of pretuples=2543 after (%5,7,11) good:226 for tuple diameter=468 number of pretuples=27452 after (%5,7,11) good:4272 for tuple diameter=480 number of pretuples=98019 after (%5,7,11) good:18700 for tuple diameter=492 number of pretuples=67494 after (%5,7,11) good:12798 for tuple diameter=504 number of pretuples=103675 after (%5,7,11) good:19624 for tuple diameter=636 number of pretuples=968315 after (%5,7,11) good:202563 for tuple diameter=756 number of pretuples=7503945 after (%5,7,11) good:1833597 for tuple diameter=768 number of pretuples=4196984 after (%5,7,11) good:1044080 for tuple diameter=828 number of pretuples=9303442 after (%5,7,11) good:2361304 for tuple diameter=972 number of pretuples=54527964 after (%5,7,11) good:14552230 for tuple diameter=1056 number of pretuples=136861807 after (%5,7,11) good:38905195 |
Send message Joined: 6 Apr 17 Posts: 14381 Credit: 0 RAC: 0 |
Demis спасибо! Очень интересно! |
Send message Joined: 6 Apr 17 Posts: 14381 Credit: 0 RAC: 0 |
Новая центральная 9-ка в терьтьем диапазоне (09:57) gp > \r gris9_3.txt logfile = "res_gris9_3.txt" 32686985961925324679543, 32686985961925324679561, 32686985961925324679567, 32686985961925324679591, 32686985961925324679597, 32686985961925324679603, 32686985961925324679627, 32686985961925324679633, 32686985961925324679651, 32686985961925324679471, 32686985961925324679543, 32686985961925324679561, 32686985961925324679567, 32686985961925324679591, 32686985961925324679597, 32686985961925324679603, 32686985961925324679627, 32686985961925324679633, 32686985961925324679651, 32686985961925324679471, 32686985961925324679477, 32686985961925324679483, 32686985961925324679501, 32686985961925324679513, 32686985961925324679543, 32686985961925324679561, 32686985961925324679567, 32686985961925324679591, 32686985961925324679597, 32686985961925324679603, 32686985961925324679627, 32686985961925324679633, 32686985961925324679651, 32686985961925324679681, 32686985961925324679693, 32686985961925324679711, 32686985961925324679717, 32686985961925324679723, 9-ка не продолжается. И очень мало вписанных последовательных простых чисел - всего 10. |
Send message Joined: 6 Apr 17 Posts: 14381 Credit: 0 RAC: 0 |
Две центральные 9-ки для паттернов с диаметром 396 4852494900776740487, 4852494900776740499, 4852494900776740517, 4852494900776740541, 4852494900776740559, 4852494900776740577, 4852494900776740601, 4852494900776740619, 4852494900776740631, 4852494900776740361, 4852494900776740363, 4852494900776740379, 4852494900776740411, 4852494900776740417, 4852494900776740421, 4852494900776740429, 4852494900776740487, 4852494900776740499, 4852494900776740517, 4852494900776740541, 4852494900776740559, 4852494900776740577, 4852494900776740601, 4852494900776740619, 4852494900776740631, 4852494900776740643, 4852494900776740687, 4852494900776740697, 4852494900776740361, 4852494900776740367, 4852494900776740391, 4852494900776740409, 4852494900776740421, 4852494900776740487, 4852494900776740499, 4852494900776740517, 4852494900776740541, 4852494900776740559, 4852494900776740577, 4852494900776740601, 4852494900776740619, 4852494900776740631, 4852494900776740697, 4852494900776740709, 4852494900776740727, 4852494900776740751, 4852494900776740757, 4852501434234981167, 4852501434234981179, 4852501434234981209, 4852501434234981221, 4852501434234981239, 4852501434234981257, 4852501434234981269, 4852501434234981299, 4852501434234981311, 4852501434234981041, 4852501434234981077, 4852501434234981143, 4852501434234981151, 4852501434234981167, 4852501434234981179, 4852501434234981209, 4852501434234981221, 4852501434234981239, 4852501434234981257, 4852501434234981269, 4852501434234981299, 4852501434234981311, 4852501434234981313, 4852501434234981413, 4852501434234981431, 4852501434234981041, 4852501434234981047, 4852501434234981071, 4852501434234981089, 4852501434234981101, 4852501434234981167, 4852501434234981179, 4852501434234981209, 4852501434234981221, 4852501434234981239, 4852501434234981257, 4852501434234981269, 4852501434234981299, 4852501434234981311, 4852501434234981377, 4852501434234981389, 4852501434234981407, 4852501434234981431, 4852501434234981437, 9-ки не продолжаются. Вписанных последовательных простых много, в первом решении аж 19 штук. То есть имеем не симметричный кортеж длины 19 из последовательных простых чисел с диаметром 336. 9-ка в нём расположена не по центру. Кортеж начинается "сладкой парочкой". А дальше ещё и кузены :) |
Send message Joined: 6 Apr 17 Posts: 14381 Credit: 0 RAC: 0 |
Сформировала программу на все 360 паттернов с диаметром 372. Запустила на тестирование в интервале на 100000. На старте программа ошибок не выдала. На очереди у меня паттерны с диаметром 384, их 190 штук. Сейчас запущу генерацию паттернов и соответствующих формул. Готово! Ахиллес трудится. |
Send message Joined: 6 Apr 17 Posts: 14381 Credit: 0 RAC: 0 |
Кстати, о "сладких парочках". В BOINC-проекте TBEG искались расстояния между рядом стоящими близнецами. Смотрите страницу https://boinc.tbrada.eu/spt/tpt_gaps.php?t=0&f=0 Показываю начало и конец этой страницы # Copyright Tomas Brada, ask on forum about reuse or citation. # tp_gaps where k=0 5: 4 4 [1 2189806] 137: 10 [2 2189806] 1931: 16 [3 2189806] 2969: 28 [4 2189806] 20441: 34 [5 2189806] 48677: 52 [6 2189806] 173357: 70 [7 2189806] 838247: 100 [8 2189806] 4297091: 28 106 [9 2189806] 14982551: 124 [10 2189806] 30781187: 130 [11 2189806] 34570661: 136 46 [12 2189806] 43891037: 148 [13 2189806] 79167731: 184 [14 2189806] 875971469: 196 [15 2189806] 1209266801: 208 [16 2189806] 2505898931: 238 [17 2189806] 3399081821: 244 [18 2189806] 5002002407: 250 [19 2189806] . . . . . . . . . 641445536620451: 736 [48 2503770] 1500286960671221: 748 4 [49 2504207] 1835020825207679: 808 [50 2508524] 3545774410393511: 826 [51 2509394] 4127074165753079: 886 [52 2509690] 11370805696999121: 898 [53 2535957] 22155036026187017: 910 [54 2648183] 24255834268843691: 988 [55 2649340] 81264833949895859: 1018 [56 2728390] 369466778811719969: 1036 [57 2918951] 584899147844503181: 1078 [58 2422743] 930954769609130597: 1120 [59 3469090] 1168126825571331887: 1150 [60 3659635] 1406616941413893059: 1300 [61 4013231] # count = 61 Что стоит в квадратных скобках, я без понятия, Томаш ничего не расшифровал. Что в результатах, более-менее понятно, кроме первой строки. Например, две пары рядом стоящих близняшек: 137, 139, 149, 151 Расстояние между ними (gap) считается. как понимаю, так: 149-139=10. Не понимаю эту строку 4297091: 28 106 [9 2189806] Тут два числа указаны. Что они означают? И главный вопрос: какая это последовательность OEIS? PS. Ну вот две пары рядом стоящих близняшек 4297091, 4297093, 4297199, 4297201, gap=106 А что означает число 28? Аналогичная строка далее 34570661: 136 46 [12 2189806] тоже два числа указаны. А, ну здесь понятно: три пары рядом стоящих близнецов 34570661, 34570663, 34570799, 34570801, 34570847, 34570849, Между первыми парами gap=136, а между второй и третьей парой gap=46. Тогда для первой строки будет так: 4297061, 4297063, 4297091, 4297093, 4297199, 4297201, первая пара близняшек неправильно указана. Между первыми парами близнецов gap=28. а между второй и третьей парой gap=106. Первую строчку тоже поняла. Три рябом стоящие пары близнецов: 5, 7, 11, 13, 17, 19. gap между этими парами близнецов будут 4 и 4. |
Send message Joined: 6 Apr 17 Posts: 14381 Credit: 0 RAC: 0 |
Программа для всех 226 паттернов с диаметром 360 выдала первое решение ? \r d360_9ki_all.txt log = 1 (on) [logfile is "rez_d360_9ki_all.txt"] *** Warning: new stack size = 536870912 (512.000 Mbytes). range of search 500272000001 (p=4852483315689699690 ) 500272300000 (p=4852486225587000000 ) 4852483766425491919, 4852483766425491949, 4852483766425491961, 4852483766425491967, 4852483766425491979, 4852483766425491991, 4852483766425491997, 4852483766425492009, 4852483766425492039, 4852483766425491799, 4852483766425491893, 4852483766425491913, 4852483766425491919, 4852483766425491949, 4852483766425491961, 4852483766425491967, 4852483766425491979, 4852483766425491991, 4852483766425491997, 4852483766425492009, 4852483766425492039, 4852483766425492063, 4852483766425492133, 4852483766425492159, 4852483766425491799, 4852483766425491811, 4852483766425491829, 4852483766425491877, 4852483766425491907, 4852483766425491919, 4852483766425491949, 4852483766425491961, 4852483766425491967, 4852483766425491979, 4852483766425491991, 4852483766425491997, 4852483766425492009, 4852483766425492039, 4852483766425492051, 4852483766425492081, 4852483766425492129, 4852483766425492147, 4852483766425492159, Вписаны 15 последовательных простых чисел. Хорошо вписались! Вполне себе 15-ка с 4 "дырками" {4852483766425491799, 4852483766425491893, 4852483766425491913, 4852483766425491919, 4852483766425491949, 4852483766425491961, 4852483766425491967, 4852483766425491979, 4852483766425491991, 4852483766425491997, 4852483766425492009, 4852483766425492039, 4852483766425492063, 4852483766425492133, 4852483766425492159} И ещё одна 4852485511978012697, 4852485511978012703, 4852485511978012709, 4852485511978012739, 4852485511978012763, 4852485511978012787, 4852485511978012817, 4852485511978012823, 4852485511978012829, 4852485511978012583, 4852485511978012609, 4852485511978012621, 4852485511978012651, 4852485511978012697, 4852485511978012703, 4852485511978012709, 4852485511978012739, 4852485511978012763, 4852485511978012787, 4852485511978012817, 4852485511978012823, 4852485511978012829, 4852485511978012859, 4852485511978012873, 4852485511978012919, 4852485511978012939, 4852485511978012583, 4852485511978012607, 4852485511978012613, 4852485511978012619, 4852485511978012637, 4852485511978012697, 4852485511978012703, 4852485511978012709, 4852485511978012739, 4852485511978012763, 4852485511978012787, 4852485511978012817, 4852485511978012823, 4852485511978012829, 4852485511978012889, 4852485511978012907, 4852485511978012913, 4852485511978012919, 4852485511978012943, 17 вписанных последовательных простых чисел, вписались симметрично, но числа почти все не совпадают, кроме, конечно, чисел центральной 9-ки и первого элемента кортежа. Можно посмотреть как на очень дырявую 17-ку. |
Send message Joined: 6 Apr 17 Posts: 14381 Credit: 0 RAC: 0 |
Нашла последовательность в OEIS о расстояниях между парами рядом стоящих близнецов https://oeis.org/A329165 A329165 Let P1, P2, P3, P4 be consecutive primes with P2-P1=P4-P3=2. a(n)=(P3-P1)/6 when the length of the gap with no primes between the two pairs of twin primes sets a record. 1, 2, 3, 5, 6, 9, 12, 17, 18, 21, 22, 23, 25, 31, 33, 35, 40, 41, 42, 47, 48, 49, 51, 52, 53, 57, 58, 62, 63, 66, 71, 75, 77, 78, 81, 83, 85, 90, 91, 93, 98, 100, 105, 108, 111, 115, 119, 123, 125, 135, 138, 148, 150, 152, 165, 170, 173, 180 a(29)-a(58) found by Tomáš Brada, Natalia Makarova, May 12 2020 Зачем-то в этой последовательности расстояние между парами близнецовделится на 6. И вторая последовательность - по сути та же самая, только в ней почему-то a(n)=(P1+P2)/12 https://oeis.org/A329164 A329164 Let P1, P2, P3, P4 be consecutive primes, with P2-P1=P4-P3=2. a(n)=(P1+P2)/12 when P3-P2 sets a new record. 1, 23, 322, 495, 3407, 8113, 28893, 139708, 716182, 2497092, 5130198, 5761777, 7315173, 13194622, 145995245, 201544467, 417649822, 566513637, 833667068, 2266818768, 4710228962, 5186737183, 5192311957, 7454170028, 9853412390, 11817808908 |
Send message Joined: 6 Apr 17 Posts: 14381 Credit: 0 RAC: 0 |
Новая центральная 9-ка в пятом диапазоне (паттерн с минимальным диаметром) 32687470674588723523529, 32687470674588723523547, 32687470674588723523553, 32687470674588723523577, 32687470674588723523583, 32687470674588723523589, 32687470674588723523613, 32687470674588723523619, 32687470674588723523637, 32687470674588723523457, 32687470674588723523529, 32687470674588723523547, 32687470674588723523553, 32687470674588723523577, 32687470674588723523583, 32687470674588723523589, 32687470674588723523613, 32687470674588723523619, 32687470674588723523637, 32687470674588723523457, 32687470674588723523463, 32687470674588723523469, 32687470674588723523487, 32687470674588723523499, 32687470674588723523529, 32687470674588723523547, 32687470674588723523553, 32687470674588723523577, 32687470674588723523583, 32687470674588723523589, 32687470674588723523613, 32687470674588723523619, 32687470674588723523637, 32687470674588723523667, 32687470674588723523679, 32687470674588723523697, 32687470674588723523703, 32687470674588723523709, Всего 10 вписанных последовательных простых чисел. Если добавить 11-ое последовательное простое, которое не вписалось, получим 11-ку с одной "дыркой" {32687470674588723523457, 32687470674588723523529, 32687470674588723523547, 32687470674588723523553, 32687470674588723523577, 32687470674588723523583, 32687470674588723523589, 32687470674588723523613, 32687470674588723523619, 32687470674588723523637, 32687470674588723523723} Симпатичненькая почти 11-ка с диаметром 252. |
Send message Joined: 6 Apr 17 Posts: 14381 Credit: 0 RAC: 0 |
Я снова о процедуре. Должна же быть в PARI/GP процедура! В теме о PARI/GP читала, что можно создать свою функцию. И процедуру тоже, наверное, можно создать. Процедура - вот она, я уже неоднократно её показывала for( i=i1,i2, bp=i*period; for( n=1,lena, bpt=bp+a[n]; if( ispseudoprime (bpt) && ispseudoprime (bpt+ 72) && ispseudoprime (bpt+ 90) && ispseudoprime (bpt+ 96) && ispseudoprime (bpt+120) && ispseudoprime (bpt+126) && ispseudoprime (bpt+132) && ispseudoprime (bpt+156) && ispseudoprime (bpt+162) && ispseudoprime (bpt+180) && nextprime(bpt+ 73) == (bpt+ 90) && nextprime(bpt+ 91) == (bpt+ 96) && nextprime(bpt+ 97) == (bpt+120) && nextprime(bpt+121) == (bpt+126) && nextprime(bpt+127) == (bpt+132) && nextprime(bpt+133) == (bpt+156) && nextprime(bpt+157) == (bpt+162) && nextprime(bpt+163) == (bpt+180) , print1(bpt+72,", ");print1(bpt+90,", ");print1(bpt+96,", ");print1(bpt+120,", ");print1(bpt+126,", ");print1(bpt+132,", ");print1(bpt+156,", "); print1(bpt+162,", ");print1(bpt+180,", ");print(); forprime(k=bpt, bpt+252, print(k,", "););print(); print1(bpt,", ");print1(bpt+6,", ");print1(bpt+12,", ");print1(bpt+30,", ");print1(bpt+42,", ");print1(bpt+72,", ");print1(bpt+90,", "); print1(bpt+96,", ");print1(bpt+120,", ");print1(bpt+126,", ");print1(bpt+132,", ");print1(bpt+156,", ");print1(bpt+162,", ");print1(bpt+180,", "); print1(bpt+210,", ");print1(bpt+222,", ");print1(bpt+240,", ");print1(bpt+246,", ");print1(bpt+252,", ");print(); ); ); ); Сейчас я вставляю эту процедуру в текст программы сотни раз в Ворде. Но она большая и за один раз не вставляется, пришлось разделить на 7 вставок. А если бы это был просто вызов процедуры, тогда вставлялось бы за один раз. Сама процедура создаётся (описывается) один раз, в начале программы. Ну вот же помню отлично, что в Паскале процедура создавалась и вызывалась. Неужели в PARI/GP такой возможности нет? PS. Ну конечно, здесь исчезла lena. Просто скопировала блок не оттуда. Исправила. Ах! И не только это, конечно. Блок скопирован для конкретного паттерна, наверное, с минимальным диаметром. |
Send message Joined: 6 Apr 17 Posts: 14381 Credit: 0 RAC: 0 |
Для паттернов с диаметром 456 нашлись две новые 9-ки, показываю только первую 4852019466266118829, 4852019466266118889, 4852019466266118901, 4852019466266118919, 4852019466266118931, 4852019466266118943, 4852019466266118961, 4852019466266118973, 4852019466266119033, 4852019466266118703, 4852019466266118719, 4852019466266118829, 4852019466266118889, 4852019466266118901, 4852019466266118919, 4852019466266118931, 4852019466266118943, 4852019466266118961, 4852019466266118973, 4852019466266119033, 4852019466266119039, 4852019466266119057, 4852019466266119093, 4852019466266119097, 4852019466266119109, 4852019466266118703, 4852019466266118709, 4852019466266118733, 4852019466266118739, 4852019466266118823, 4852019466266118829, 4852019466266118889, 4852019466266118901, 4852019466266118919, 4852019466266118931, 4852019466266118943, 4852019466266118961, 4852019466266118973, 4852019466266119033, 4852019466266119039, 4852019466266119123, 4852019466266119129, 4852019466266119153, 4852019466266119159, 9-ка продолжается до 11-ки в одну сторону, вот 11-ка с одной "дыркой" {4852019466266118719, 4852019466266118829, 4852019466266118889, 4852019466266118901, 4852019466266118919, 4852019466266118931, 4852019466266118943, 4852019466266118961, 4852019466266118973, 4852019466266119033, 4852019466266119039} Вторая 9-ка не продолжается даже до почти 11-ки. Поскольку 9-ок находится довольно много, дальше буду показывать только 9-ки с продолжением хотя бы до почти 11-ки. |
Send message Joined: 11 Jul 17 Posts: 174 Credit: 4,964,801 RAC: 1 |
Примерная статистика по наполненности паттернов:Посчитал сегодня полную статистику в пределах первой тысячи: for tuple diameter=276 number of pretuples=7 after (%5,7,11) good:0 for tuple diameter=288 number of pretuples=18 after (%5,7,11) good:0 for tuple diameter=300 number of pretuples=260 after (%5,7,11) good:11 for tuple diameter=312 number of pretuples=305 after (%5,7,11) good:10 for tuple diameter=324 number of pretuples=445 after (%5,7,11) good:27 for tuple diameter=336 number of pretuples=1570 after (%5,7,11) good:116 for tuple diameter=348 number of pretuples=342 after (%5,7,11) good:34 for tuple diameter=360 number of pretuples=2543 after (%5,7,11) good:226 for tuple diameter=372 number of pretuples=3638 after (%5,7,11) good:360 for tuple diameter=384 number of pretuples=2484 after (%5,7,11) good:190 for tuple diameter=396 number of pretuples=4312 after (%5,7,11) good:655 for tuple diameter=408 number of pretuples=2961 after (%5,7,11) good:621 for tuple diameter=420 number of pretuples=29579 after (%5,7,11) good:4302 for tuple diameter=432 number of pretuples=17985 after (%5,7,11) good:3549 for tuple diameter=444 number of pretuples=17281 after (%5,7,11) good:3401 for tuple diameter=456 number of pretuples=29252 after (%5,7,11) good:4435 for tuple diameter=468 number of pretuples=27452 after (%5,7,11) good:4272 for tuple diameter=480 number of pretuples=98019 after (%5,7,11) good:18700 for tuple diameter=492 number of pretuples=67494 after (%5,7,11) good:12798 for tuple diameter=504 number of pretuples=103675 after (%5,7,11) good:19624 for tuple diameter=516 number of pretuples=121738 after (%5,7,11) good:20625 for tuple diameter=528 number of pretuples=136693 after (%5,7,11) good:24300 for tuple diameter=540 number of pretuples=326599 after (%5,7,11) good:61067 for tuple diameter=552 number of pretuples=284912 after (%5,7,11) good:55180 for tuple diameter=564 number of pretuples=222397 after (%5,7,11) good:50565 for tuple diameter=576 number of pretuples=380000 after (%5,7,11) good:73998 for tuple diameter=588 number of pretuples=496087 after (%5,7,11) good:98872 for tuple diameter=600 number of pretuples=1060518 after (%5,7,11) good:227311 for tuple diameter=612 number of pretuples=702580 after (%5,7,11) good:145913 for tuple diameter=624 number of pretuples=611350 after (%5,7,11) good:135515 for tuple diameter=636 number of pretuples=968315 after (%5,7,11) good:202563 for tuple diameter=648 number of pretuples=958590 after (%5,7,11) good:198053 for tuple diameter=660 number of pretuples=3119054 after (%5,7,11) good:705132 for tuple diameter=672 number of pretuples=2676251 after (%5,7,11) good:615370 for tuple diameter=684 number of pretuples=2233796 after (%5,7,11) good:483349 for tuple diameter=696 number of pretuples=1768372 after (%5,7,11) good:457143 for tuple diameter=708 number of pretuples=2243099 after (%5,7,11) good:514174 for tuple diameter=720 number of pretuples=6403183 after (%5,7,11) good:1455012 for tuple diameter=732 number of pretuples=4511022 after (%5,7,11) good:1074343 for tuple diameter=744 number of pretuples=4698089 after (%5,7,11) good:1058599 for tuple diameter=756 number of pretuples=7503945 after (%5,7,11) good:1833597 for tuple diameter=768 number of pretuples=4196984 after (%5,7,11) good:1044080 for tuple diameter=780 number of pretuples=13616134 after (%5,7,11) good:3491411 for tuple diameter=792 number of pretuples=10197425 after (%5,7,11) good:2559623 for tuple diameter=804 number of pretuples=8415809 after (%5,7,11) good:1999905 for tuple diameter=816 number of pretuples=9152351 after (%5,7,11) good:2626749 for tuple diameter=828 number of pretuples=9303442 after (%5,7,11) good:2361304 for tuple diameter=840 number of pretuples=40008365 after (%5,7,11) good:10289762 for tuple diameter=852 number of pretuples=17059304 after (%5,7,11) good:4612534 for tuple diameter=864 number of pretuples=16330430 after (%5,7,11) good:4606165 for tuple diameter=876 number of pretuples=20718852 after (%5,7,11) good:5359783 for tuple diameter=888 number of pretuples=21206925 after (%5,7,11) good:5352510 for tuple diameter=900 number of pretuples=50643396 after (%5,7,11) good:13689263 for tuple diameter=912 number of pretuples=30694975 after (%5,7,11) good:8443329 for tuple diameter=924 number of pretuples=56202731 after (%5,7,11) good:15349179 for tuple diameter=936 number of pretuples=43311591 after (%5,7,11) good:11137615 for tuple diameter=948 number of pretuples=38660329 after (%5,7,11) good:9902080 for tuple diameter=960 number of pretuples=92529681 after (%5,7,11) good:24904241 for tuple diameter=972 number of pretuples=54527964 after (%5,7,11) good:14552230 for tuple diameter=984 number of pretuples=59766381 after (%5,7,11) good:16779239 for tuple diameter=996 number of pretuples=71114143 after (%5,7,11) good:18921126 |
Send message Joined: 6 Apr 17 Posts: 14381 Credit: 0 RAC: 0 |
Посчитал сегодня полную статистику в пределах первой тысячи: Замечательно! Теперь можно организовать штурм ;) Думаю, что диаметров в пределах тысячи для штурма вполне достаточно. Demis как насчёт штурма? :) |
Send message Joined: 6 Apr 17 Posts: 14381 Credit: 0 RAC: 0 |
gris пытается сочинить функцию вместо постоянного блока, который я вставляю сейчас в текст формируемой программы. Он говорит, что это может быть функция. Хорошо, пусть будет функция, дело не в названии - процедура или функция. Главное, чтобы работало. Итак, вот функция, созданная gris {ffff(lena)=lena; for( i=i1,i2, bp=i*period; for( n=1,384, bpt=bp+a[n]; if( ispseudoprime (bpt) && ispseudoprime (bpt+ 72) && ispseudoprime (bpt+ 90) && ispseudoprime (bpt+ 96) && ispseudoprime (bpt+120) && ispseudoprime (bpt+126) && ispseudoprime (bpt+132) && ispseudoprime (bpt+156) && ispseudoprime (bpt+162) && ispseudoprime (bpt+180) && nextprime(bpt+ 73) == (bpt+ 90) && nextprime(bpt+ 91) == (bpt+ 96) && nextprime(bpt+ 97) == (bpt+120) && nextprime(bpt+121) == (bpt+126) && nextprime(bpt+127) == (bpt+132) && nextprime(bpt+133) == (bpt+156) && nextprime(bpt+157) == (bpt+162) && nextprime(bpt+163) == (bpt+180) , print1(bpt+72,", ");print1(bpt+90,", ");print1(bpt+96,", ");print1(bpt+120,", ");print1(bpt+126,", ");print1(bpt+132,", ");print1(bpt+156,", "); print1(bpt+162,", ");print1(bpt+180,", ");print(); forprime(k=bpt, bpt+252, print(k,", "););print(); print1(bpt,", ");print1(bpt+6,", ");print1(bpt+12,", ");print1(bpt+30,", ");print1(bpt+42,", ");print1(bpt+72,", ");print1(bpt+90,", "); print1(bpt+96,", ");print1(bpt+120,", ");print1(bpt+126,", ");print1(bpt+132,", ");print1(bpt+156,", ");print1(bpt+162,", ");print1(bpt+180,", "); print1(bpt+210,", ");print1(bpt+222,", ");print1(bpt+240,", ");print1(bpt+246,", ");print1(bpt+252,", ");print(); ); ); ); return(lena)} Фактически показанная выше процедура описана как функция. Это мне и требовалось. Осталось грамотно оформить весь текст программы. По шагам. gris, пожалуйста, не торопитесь. Я не могу быстро, у меня бо-бо :) |
Send message Joined: 6 Apr 17 Posts: 14381 Credit: 0 RAC: 0 |
Цитирую gris после переменных блоков будет вставляться только ffff(lena); Так, сейчас... соображаю. Пишу сначала функцию, потом текст программы. Потом всё это вместе запускаю. Правильно? Сейчас напишу для двух паттернов с диаметром 264. |
Send message Joined: 6 Apr 17 Posts: 14381 Credit: 0 RAC: 0 |
Ну вот, gris прислал письмо, что функция не такая. gris, пожалуйста, разберитесь. Я выложила свою процедуру здесь https://boinc.progger.info/odlk/forum_thread.php?id=237&postid=11249 Вот её и надо записать в вашу функцию. Или я что-то не понимаю?? gris Ну вот, с этого и надо было начинать! Выложенная по указанной ссылке процедура не та. Просто скопировала не оттуда. [Ошибку, понятное дело, никто не указал, некому. Вот только сейчас и обнаружили. Я тоже могу ошибаться!] Конечно, она заточена под конкретный паттерн с минимальным диаметром 252. Правильно будет такая процедура for( i=i1,i2, bp=i*period; for( n=1,lena, bpt=bp+a[n]; if( ispseudoprime (bpt) && ispseudoprime (bpt+p[6]) && ispseudoprime (bpt+p[7] ) && ispseudoprime (bpt+p[8]) && ispseudoprime (bpt+p[9]) && ispseudoprime (bpt+p[10]) && ispseudoprime (bpt+p[11]) && ispseudoprime (bpt+p[12]) && ispseudoprime (bpt+p[13]) && ispseudoprime (bpt+p[14]) && nextprime((bpt+p[6])+1) == (bpt+p[7]) && nextprime((bpt+p[7])+1) == (bpt+p[8]) && nextprime((bpt+p[8])+1) == (bpt+p[9]) && nextprime((bpt+p[9])+1) == (bpt+p[10]) && nextprime((bpt+p[10])+1) == (bpt+p[11]) && nextprime((bpt+p[11])+1) == (bpt+p[12]) && nextprime((bpt+p[12])+1) == (bpt+p[13]) && nextprime((bpt+p[13])+1) == (bpt+p[14]) , print1(bpt+p[6],", ");print1(bpt+p[7],", ");print1(bpt+p[8],", ");print1(bpt+p[9],", ");print1(bpt+p[10],", ");print1(bpt+p[11],", ");print1(bpt+p[12],", "); print1(bpt+p[13],", ");print1(bpt+p[14],", ");print(); forprime(k=bpt, bpt+p[19], print(k,", "););print(); print1(bpt,", ");print1(bpt+p[2],", ");print1(bpt+p[3],", ");print1(bpt+p[4],", ");print1(bpt+p[5],", ");print1(bpt+p[6],", ");print1(bpt+p[7],", "); print1(bpt+p[8],", ");print1(bpt+p[9],", ");print1(bpt+p[10],", ");print1(bpt+p[11],", ");print1(bpt+p[12],", ");print1(bpt+p[13],", ");print1(bpt+p[14],", "); print1(bpt+p[15],", ");print1(bpt+p[16],", ");print1(bpt+p[17],", ");print1(bpt+p[18],", ");print1(bpt+p[19],", ");print(); ); ); ); |
Send message Joined: 6 Apr 17 Posts: 14381 Credit: 0 RAC: 0 |
Теперь функция будет так записываться {ffff(lena)=lena; for( i=i1,i2, bp=i*period; for( n=1,lena, bpt=bp+a[n]; if( ispseudoprime (bpt) && ispseudoprime (bpt+p[6]) && ispseudoprime (bpt+p[7] ) && ispseudoprime (bpt+p[8]) && ispseudoprime (bpt+p[9]) && ispseudoprime (bpt+p[10]) && ispseudoprime (bpt+p[11]) && ispseudoprime (bpt+p[12]) && ispseudoprime (bpt+p[13]) && ispseudoprime (bpt+p[14]) && nextprime((bpt+p[6])+1) == (bpt+p[7]) && nextprime((bpt+p[7])+1) == (bpt+p[8]) && nextprime((bpt+p[8])+1) == (bpt+p[9]) && nextprime((bpt+p[9])+1) == (bpt+p[10]) && nextprime((bpt+p[10])+1) == (bpt+p[11]) && nextprime((bpt+p[11])+1) == (bpt+p[12]) && nextprime((bpt+p[12])+1) == (bpt+p[13]) && nextprime((bpt+p[13])+1) == (bpt+p[14]) , print1(bpt+p[6],", ");print1(bpt+p[7],", ");print1(bpt+p[8],", ");print1(bpt+p[9],", ");print1(bpt+p[10],", ");print1(bpt+p[11],", ");print1(bpt+p[12],", "); print1(bpt+p[13],", ");print1(bpt+p[14],", ");print(); forprime(k=bpt, bpt+p[19], print(k,", "););print(); print1(bpt,", ");print1(bpt+p[2],", ");print1(bpt+p[3],", ");print1(bpt+p[4],", ");print1(bpt+p[5],", ");print1(bpt+p[6],", ");print1(bpt+p[7],", "); print1(bpt+p[8],", ");print1(bpt+p[9],", ");print1(bpt+p[10],", ");print1(bpt+p[11],", ");print1(bpt+p[12],", ");print1(bpt+p[13],", ");print1(bpt+p[14],", "); print1(bpt+p[15],", ");print1(bpt+p[16],", ");print1(bpt+p[17],", ");print1(bpt+p[18],", ");print1(bpt+p[19],", ");print(); ); ); ); return(lena)} gris, правильно? |
Send message Joined: 6 Apr 17 Posts: 14381 Credit: 0 RAC: 0 |
Цитата Цитирую gris Начинаю сначала. 1) Пишу функцию {ffff(lena)=lena; for( i=i1,i2, bp=i*period; for( n=1,lena, bpt=bp+a[n]; if( ispseudoprime (bpt) && ispseudoprime (bpt+p[6]) && ispseudoprime (bpt+p[7] ) && ispseudoprime (bpt+p[8]) && ispseudoprime (bpt+p[9]) && . . . . . . . . . nextprime((bpt+p[13])+1) == (bpt+p[14]) , print1(bpt+p[6],", ");print1(bpt+p[7],", ");print1(bpt+p[8],", ");print1(bpt+p[9],", ");print1(bpt+p[10],", ");print1(bpt+p[11],", ");print1(bpt+p[12],", "); print1(bpt+p[13],", ");print1(bpt+p[14],", ");print(); forprime(k=bpt, bpt+p[19], print(k,", "););print(); print1(bpt,", ");print1(bpt+p[2],", ");print1(bpt+p[3],", ");print1(bpt+p[4],", ");print1(bpt+p[5],", ");print1(bpt+p[6],", ");print1(bpt+p[7],", "); print1(bpt+p[8],", ");print1(bpt+p[9],", ");print1(bpt+p[10],", ");print1(bpt+p[11],", ");print1(bpt+p[12],", ");print1(bpt+p[13],", ");print1(bpt+p[14],", "); print1(bpt+p[15],", ");print1(bpt+p[16],", ");print1(bpt+p[17],", ");print1(bpt+p[18],", ");print1(bpt+p[19],", ");print(); ); ); ); return(lena)} дальше пишу текст программы \l rez_d264_9ki.txt; {period = 9699690; pl=19; dm=264; N=period; \\N=#19 print("range of search"); i1=500226000001; print(i1," (p=", i1*9699690," )"); i2=500226030000; print(i2," (p=", i2*9699690," )"); p=[0, 12, 24, 30, 42, 54, 84, 90, 114, 132, 150, 174, 180, 210, 222, 234, 240, 252, 264]; a=[19739,30827,103757,179339,234317,290009,307247,322339,373169,446099,473269,501077,528247,576659,601177,643439,649589,660677,676759,716369,771347,781577,846929,870589,919859,943519,954607,985067,1003019,1030189,1103119,1113689,1140859,1158097,1168327,1213789,1273289,1278997,1296949,1338497,1411427,1438597,1500439,1567219,1632427,1705357,1770709,1780939,1821539,1835917,1894469,1908847,1959677,2119277,2229947,2247899,2262277,2313107,2386037,2413207,2451389,2461619,2516597,2572289,2583377,2589527,2616697,2656307,2731889,2745319,2755549,2786867,2810527,2859797,2883457,2925719,2942957,2959039,2970127,2998649,3025819,3043057,3053627,3080797,3129209,3153727,3202139,3213227,3229309,3236887,3312469,3385399,3423139,3440377,3496069,3507157,3555569,3582739,3655669,3706499,3710647,3720877,3761477,3834407,3849499,3909989,4052989,4103819,4176749,4187837,4203919,4263419,4336349,4374089,4391327,4401557,4447019,4512227,4530179,4557349,4630279,4671827,4685257,4695487,4733669,4760839,4800449,4833769,4865657,4898977,4938587,4965757,5003939,5014169,5027599,5069147,5142077,5169247,5187199,5252407,5297869,5308099,5325337,5363077,5436007,5495507,5511589,5522677,5595607,5646437,5789437,5849927,5865019,5937949,5978549,5988779,5992927,6043757,6116687,6143857,6192269,6203357,6259049,6276287,6314027,6386957,6462539,6470117,6486199,6497287,6545699,6570217,6618629,6645799,6656369,6673607,6700777,6729299,6740387,6756469,6773707,6815969,6839629,6888899,6912559,6943877,6954107,6967537,7043119,7082729,7109899,7116049,7127137,7182829,7237807,7248037,7286219,7313389,7386319,7437149,7451527,7469479,7580149,7739749,7790579,7804957,7863509,7877887,7918487,7928717,7994069,8066999,8132207,8198987,8260829,8287999,8360929,8402477,8420429,8426137,8485637,8531099,8541329,8558567,8585737,8596307,8669237,8696407,8714359,8744819,8755907,8779567,8828837,8852497,8917849,8928079,8983057,9022667,9038749,9049837,9055987,9098249,9122767,9171179,9198349,9226157,9253327,9326257,9377087,9392179,9409417,9465109,9520087,9595669,9668599,9679687]; lena=#a; стоп! Что здесь надо написать? Функцию? Как? |
Send message Joined: 6 Apr 17 Posts: 14381 Credit: 0 RAC: 0 |
Пришла подсказка от gris В функции перед закрывающей фигурной скобкой поставьте точку с запятой. Пишу снова {ffff(lena)=lena; for( i=i1,i2, bp=i*period; for( n=1,lena, bpt=bp+a[n]; if( ispseudoprime (bpt) && ispseudoprime (bpt+p[6]) && ispseudoprime (bpt+p[7] ) && ispseudoprime (bpt+p[8]) && ispseudoprime (bpt+p[9]) && . . . . . . . . . nextprime((bpt+p[13])+1) == (bpt+p[14]) , print1(bpt+p[6],", ");print1(bpt+p[7],", ");print1(bpt+p[8],", ");print1(bpt+p[9],", ");print1(bpt+p[10],", ");print1(bpt+p[11],", ");print1(bpt+p[12],", "); print1(bpt+p[13],", ");print1(bpt+p[14],", ");print(); forprime(k=bpt, bpt+p[19], print(k,", "););print(); print1(bpt,", ");print1(bpt+p[2],", ");print1(bpt+p[3],", ");print1(bpt+p[4],", ");print1(bpt+p[5],", ");print1(bpt+p[6],", ");print1(bpt+p[7],", "); print1(bpt+p[8],", ");print1(bpt+p[9],", ");print1(bpt+p[10],", ");print1(bpt+p[11],", ");print1(bpt+p[12],", ");print1(bpt+p[13],", ");print1(bpt+p[14],", "); print1(bpt+p[15],", ");print1(bpt+p[16],", ");print1(bpt+p[17],", ");print1(bpt+p[18],", ");print1(bpt+p[19],", ");print(); ); ); ); return(lena);} \l rez_d264_9ki.txt; {period = 9699690; pl=19; dm=264; N=period; \\N=#19 print("range of search"); i1=500226000001; print(i1," (p=", i1*9699690," )"); i2=500226030000; print(i2," (p=", i2*9699690," )"); p=[0, 12, 24, 30, 42, 54, 84, 90, 114, 132, 150, 174, 180, 210, 222, 234, 240, 252, 264]; a=[19739,30827,103757,179339,234317,290009,307247,322339,373169,446099,473269,501077,528247,576659,601177,643439,649589,660677,676759,716369,771347,781577,846929,870589,919859,943519,954607,985067,1003019,1030189,1103119,1113689,1140859,1158097,1168327,1213789,1273289,1278997,1296949,1338497,1411427,1438597,1500439,1567219,1632427,1705357,1770709,1780939,1821539,1835917,1894469,1908847,1959677,2119277,2229947,2247899,2262277,2313107,2386037,2413207,2451389,2461619,2516597,2572289,2583377,2589527,2616697,2656307,2731889,2745319,2755549,2786867,2810527,2859797,2883457,2925719,2942957,2959039,2970127,2998649,3025819,3043057,3053627,3080797,3129209,3153727,3202139,3213227,3229309,3236887,3312469,3385399,3423139,3440377,3496069,3507157,3555569,3582739,3655669,3706499,3710647,3720877,3761477,3834407,3849499,3909989,4052989,4103819,4176749,4187837,4203919,4263419,4336349,4374089,4391327,4401557,4447019,4512227,4530179,4557349,4630279,4671827,4685257,4695487,4733669,4760839,4800449,4833769,4865657,4898977,4938587,4965757,5003939,5014169,5027599,5069147,5142077,5169247,5187199,5252407,5297869,5308099,5325337,5363077,5436007,5495507,5511589,5522677,5595607,5646437,5789437,5849927,5865019,5937949,5978549,5988779,5992927,6043757,6116687,6143857,6192269,6203357,6259049,6276287,6314027,6386957,6462539,6470117,6486199,6497287,6545699,6570217,6618629,6645799,6656369,6673607,6700777,6729299,6740387,6756469,6773707,6815969,6839629,6888899,6912559,6943877,6954107,6967537,7043119,7082729,7109899,7116049,7127137,7182829,7237807,7248037,7286219,7313389,7386319,7437149,7451527,7469479,7580149,7739749,7790579,7804957,7863509,7877887,7918487,7928717,7994069,8066999,8132207,8198987,8260829,8287999,8360929,8402477,8420429,8426137,8485637,8531099,8541329,8558567,8585737,8596307,8669237,8696407,8714359,8744819,8755907,8779567,8828837,8852497,8917849,8928079,8983057,9022667,9038749,9049837,9055987,9098249,9122767,9171179,9198349,9226157,9253327,9326257,9377087,9392179,9409417,9465109,9520087,9595669,9668599,9679687]; lena=#a; ffff(lena); p=[0, 12, 24, 42, 54, 72, 84, 90, 114, 132, 150, 174, 180, 192, 210, 222, 240, 252, 264]; a=[1889,10937,72817,99949,113039,129797,156967,165157,178247,313057,331999,340189,353279,380449,397207,415159,499309,508357,583327,591517,610459,675667,823567,837647,850699,925669,939749,1102027,1116107,1179989,1264007,1291139,1348157,1356347,1366109,1431317,1450259,1458449,1542467,1606349,1626617,1633519,1690499,1761427,1774517,1782707,1801649,1809877,1863529,1866857,1876619,1884809,1911979,1928737,1941827,1968959,2039887,2052977,2103769,2116859,2144029,2152219,2271937,2280127,2293217,2320387,2355097,2374039,2382229,2395319,2422489,2439247,2550397,2614279,2662729,2790637,2865607,2879687,2892739,3054719,3103169,3119927,3222029,3231077,3306047,3333179,3390197,3398387,3546287,3565229,3573419,3573457,3600589,3613679,3630437,3648389,3675559,3732539,3741587,3803467,3816557,3824747,3843689,3851917,3908897,3915799,4026949,4043707,4056797,4083929,4083967,4092157,4111099,4145809,4158899,4186069,4194259,4313977,4322167,4335257,4362427,4426309,4537459,4554217,4602667,4656319,4704769,4832677,4866749,4994657,5043107,5096759,5145209,5161967,5273117,5336999,5364169,5377259,5385449,5505167,5513357,5540527,5553617,5588327,5607269,5615459,5615497,5642629,5655719,5672477,5783627,5790529,5847509,5855737,5874679,5882869,5895959,5957839,5966887,6023867,6051037,6068989,6085747,6098837,6125969,6126007,6134197,6153139,6301039,6309229,6366247,6393379,6468349,6477397,6579499,6596257,6644707,6806687,6819739,6833819,6908789,7036697,7085147,7149029,7260179,7276937,7304107,7317197,7325387,7344329,7379039,7406209,7419299,7427489,7547207,7555397,7582567,7595657,7646449,7659539,7730467,7757599,7770689,7787447,7814617,7822807,7832569,7835897,7889549,7897777,7916719,7924909,7937999,8008927,8065907,8072809,8093077,8156959,8240977,8249167,8268109,8333317,8343079,8351269,8408287,8435419,8519437,8583319,8597399,8759677,8773757,8848727,8861779,8875859,9023759,9088967,9107909,9116099,9191069,9200117,9284267,9302219,9318977,9346147,9359237,9367427,9386369,9521179,9534269,9542459,9569629,9586387,9599477,9626609,9688489,9697537]; lena=#a; ffff(lena); } Изобразила :)) Надеюсь, что теперь всё правильно. Будет ли работать? Как я понимаю, это всё мне надо записать в текстовый файл и запустить его. Щас... всё исполню... :) Нет, не работает, ошибки выдаёт (20:32) gp > \r test3_gris_d264.txt *** syntax error, unexpected $end, expecting )-> or ',' or ')': *** ...rint(););););return(lena); *** ^- log = 1 (on) [logfile is "rez_d264_9ki.txt"] range of search 500226000001 (p=4852037129949699690 ) 500226030000 (p=4852037420930700000 ) *** at top-level: ...9668599,9679687];lena=#a;ffff(lena);p=[0,12,24 *** ^--------------------- *** not a function in function call (20:33) gp > gris, что я сделала не так? Вроде всё делала точно следуя вашим инструкциям. Не хочет работать :( Кстати, перевод ошибки "not a function in function call" "не функция в вызове функции" Таким образом, функцию программа не понимает. |
Send message Joined: 6 Apr 17 Posts: 14381 Credit: 0 RAC: 0 |
Так и хочется воскликнуть "Б-л-и-и-и-н!" На dxdy.ru полно спецов по PARI. Многие читают мои блоги. Неужели никто не может подсказать?! gris сам недавно начал осваивать PARI и пока знает его примерно так же, как и я. Вот с функциями он вопрос задавал на dxdy.ru, я это видела. Ему там wrest ответил. Но дальше вопроса gris вряд ли сильно продвинулся. В теории знать что-то - это одно, а на практике уметь применить - это совсем другое. В теории всегда всё просто. |
©2025 (C) Progger