Message boards : Number crunching : Модифицированная рабочая программа для Приложения 4
Message board moderation
| Author | Message |
|---|---|
Natalia MakarovaSend message Joined: 24 Dec 24 Posts: 569 Credit: 8,342,279 RAC: 32,490 |
Публикую модифицированную рабочую программу для Приложения 4 default(timer,1);
{pt=[0, 18, 30, 60, 78, 84, 108, 114, 120, 144, 150, 168, 198, 210, 228];
tin = fileopen("in.txt");
rost17=eval(filereadstr(tin));
rost19=eval(filereadstr(tin));
rost23=eval(filereadstr(tin));
rost29=eval(filereadstr(tin));
rost31=eval(filereadstr(tin));
rost37=eval(filereadstr(tin));
rost41=eval(filereadstr(tin));
vivod=eval(filereadstr(tin));
fileclose(tin);
mod=117288381359406970983270;
vmy=vector(50);
pat9=vector(9); pat11=vector(11); pat13=vector(13); pat15=vector(15);
pt9=[0, 18, 24, 48, 54, 60, 84, 90, 108];
pt11=[0, 30, 48, 54, 78, 84, 90, 114, 120, 138, 168];
pt13=[0,12,42,60,66,90,96,102,126,132,150,180,192];
v2=[1]; \\1
v3=[1, 2]; \\2
v5=[3, 4]; \\2
v7=[1, 2]; \\2
v11=[5, 9]; \\2
v13=[2, 4]; \\2
v17=[12, 13, 14, 15]; \\4
v19=[4, 5, 7, 9, 10, 12, 14, 15]; \\8
v23=[3, 4, 6, 10, 12, 13, 15, 19, 21, 22]; \\10
v29=[7, 10, 12, 13, 14, 15, 16, 17, 18, 19, 20, 21, 23, 26]; \\14
v31=[3, 6, 8, 12, 14, 17, 21, 22, 23, 24, 25, 26, 27, 28, 29, 30]; \\16
v37=[1, 2, 5, 6, 8, 9, 10, 11, 13, 15, 16, 18, 20, 21, 22, 23, 25,
26, 29, 30, 32, 36]; \\22
v41=[1, 2, 5, 6, 8, 10, 12, 13, 16, 17, 19, 21, 24, 25, 26, 27, 28, 29, 30,
31, 32, 33, 34, 35, 38, 40]; \\26
v43=[1, 3, 6, 7, 10, 11, 12, 14, 16, 18, 19, 20, 23, 24, 27, 29,
31, 32, 33, 34, 35, 36, 37, 38, 39, 40, 41, 42]; \\28
v47=[1, 2, 3, 4, 5, 6, 8, 9, 11, 12, 13, 14, 15, 18, 19, 22, 23, 24,
26, 28, 30, 31, 32, 35, 36, 39, 40, 41, 42, 43, 45, 46]; \\32
v53=[1, 3, 4, 5, 6, 7, 8, 10, 11, 12, 13, 16, 17, 18, 19, 20, 21, 24,
25, 26, 27, 29, 30, 31, 32, 33, 34, 36, 38, 40, 41, 42, 43, 47,
48, 49, 50, 52]; \\38
v59=[1, 2, 3, 5, 6, 7, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20, 21,
22, 23, 24, 25, 28, 30, 31, 32, 35, 36, 37, 39, 42, 43, 44, 45,
46, 47, 48, 49, 50, 51, 52, 53, 54, 55, 56]; \\44
v61=[3, 4, 5, 6, 7, 9, 10, 11, 12, 13, 17, 18, 19, 20,
21, 22, 23, 24, 25, 26, 27, 28, 29, 30, 32, 35, 36, 37,
40, 41, 42, 45, 47, 48, 49, 50, 51, 52,
53, 54, 55, 56, 57, 58, 59, 60]; \\46
fout=fileopen(vivod".txt","a");
a=[2, 3, 5, 7, 11, 13, 17, 19, 23, 29, 31, 37, 41, 43, 47, 53, 59, 61];
r=[1, 1, 3, 1, 5, 2, rost17, rost19, rost23, rost29, rost31, rost37, rost41, 0, 0, 0, 0, 0];
M=vector(18); MMO=vector(18);
for(i=1,18, M[i]=mod/a[i]; );
for(i=1,18, MMO[i]=M[i]*gcdext(M[i],a[i])[1]; );
x1 = sum(i = 1, 13, r[i] * MMO[i]);
for(i43=1,28,r[14]=v43[i43];
for(i47=1,32,r[15]=v47[i47];
for(i53=1,38,r[16]=v53[i53];
for(i59=1,44,r[17]=v59[i59];
for(i61=1,46,r[18]=v61[i61];
x2 = sum(i = 14, 18, r[i] * MMO[i]);
bpt = (x1 + x2)%mod+mod;
\\кортеж 9
if(ispseudoprime(bpt+pt[4])
&& ispseudoprime(bpt+pt[5]) && ispseudoprime(bpt+pt[6])
&& ispseudoprime(bpt+pt[7]) && ispseudoprime(bpt+pt[8])
&& ispseudoprime(bpt+pt[9]) && ispseudoprime(bpt+pt[10])
&& ispseudoprime(bpt+pt[11])&& ispseudoprime(bpt+pt[12]),
l=0;
forprime(p=bpt+pt[4],bpt+pt[12], l++; vmy[l]=p; );
if(l==9,
for(m=2,9, pat9[m]=vmy[m]-vmy[1]; );
if(pat9==pt9,
w1=strprintf("%d: %d\n",vmy[1],pat9);
print(w1);
filewrite(fout,w1); );
\\кортеж 11
if(ispseudoprime(bpt+pt[3]) && ispseudoprime(bpt+pt[13]),
l=0;
forprime(p=bpt+pt[3],bpt+pt[13], l++; vmy[l]=p; );
if(l==11,
for(m=2,11, pat11[m]=vmy[m]-vmy[1]; );
if(pat11==pt11,
w1=strprintf("%d: %d\n",vmy[1],pat11);
print(w1);
filewrite(fout,w1); );
\\кортеж 13
if(ispseudoprime(bpt+pt[2]) && ispseudoprime(bpt+pt[14]),
l=0;
forprime(p=bpt+pt[2],bpt+pt[14], l++; vmy[l]=p; );
if(l==13,
for(m=2,13, pat13[m]=vmy[m]-vmy[1]; );
if(pat13==pt13,
w1=strprintf("%d: %d\n",vmy[1],pat13);
print(w1);
filewrite(fout,w1); );
\\кортеж 15
if(ispseudoprime(bpt) && ispseudoprime(bpt+pt[15]),
l=0;
forprime(p=bpt,bpt+pt[15], l++; vmy[l]=p; );
if(l==15,
for(m=2,15, pat15[m]=vmy[m]-vmy[1]; );
if(pat15==pt,
w1=strprintf("%d: %d\n",vmy[1],pat15);
print(w1);
filewrite(fout,w1); );
))))))); \\ if l==9
); \\ if ispseudoprime
);););););
w1=strprintf("end ");
print(w1);
filewrite(fout,w1);
fileclose(fout)
}The TBEG BOINC project maintains the database https://boinc.tbrada.eu/spt/explore.php |
Natalia MakarovaSend message Joined: 24 Dec 24 Posts: 569 Credit: 8,342,279 RAC: 32,490 |
termit прошу вас заменить работающую в Приложении 4 рабочую программу на новую (см. предыдущее сообщение). Изменения очень важные. Работающая сейчас программа ищет центральные 15-ки в диапазоне (0, 61#). Попутно ищет также центральные 9-ки, 11-ки и 13-ки в этом диапазоне. Понятно, что главная цель найти все центральные 15-ки в указанном диапазоне. Но команда, работающая на форуме dxdy.ru, уже нашла все центральные 15-ки в указанном диапазоне. Поэтому необходимо перевести поиск в другой диапазон, что и сделано в новой программе. Кроме того, мне удалось с помощью коллеги ускорить программу в два раза. В связи с этим такой вопрос: оставить ли диапазон таким со временем в два разе меньшим или увеличить диапазон вдвое, время тоже увеличится вдвое. Сейчас программа выполняется на моём ПК 14 минут. Очень важно! Для новой программы используются те же задания, который вы уже сгенерировали, и они используются в настоящий момент в Приложении 4. Убедительно прошу заняться этим вопросом как можно быстрее. The TBEG BOINC project maintains the database https://boinc.tbrada.eu/spt/explore.php |
Natalia MakarovaSend message Joined: 24 Dec 24 Posts: 569 Credit: 8,342,279 RAC: 32,490 |
Для всех! Вы можете протестировать опубликованную рабочую программу, записав в файл in.txt следующие значения 12 4 4 18 3 13 28 out225 В выходном файле out225.txt вы увидите результаты. Время выполнения программы вы увидите в консоли. The TBEG BOINC project maintains the database https://boinc.tbrada.eu/spt/explore.php |
Natalia MakarovaSend message Joined: 24 Dec 24 Posts: 569 Credit: 8,342,279 RAC: 32,490 |
На форуме dxdy.ru утверждают, что опубликованная рабочая программа может быть ещё оптимизирована с достижением ускорения в 3-4 раза. Если среди участников проекта есть программисты, хорошо знающие PARI/GP, прошу высказать предложения по оптимизации программы с целью ускорения. Необходимо учитывать, что оптимизация не должна нарушать алгоритм программы. Сделать это надо быстро, пока рабочая программа в Приложении 4 не заменена на новую. Если у вас есть вопросы по программе, пожалуйста, задавайте их здесь или в ЛС. The TBEG BOINC project maintains the database https://boinc.tbrada.eu/spt/explore.php |
Natalia MakarovaSend message Joined: 24 Dec 24 Posts: 569 Credit: 8,342,279 RAC: 32,490 |
Здравствуйте, termit! Прошло две недели с опубликования стартового сообщения о замене рабочей программы в Приложении 4. Я не получила ответа ни здесь, ни в личной переписке. Я прекращаю загрузку заданий для Приложения 4. Предлагаю вам изменить рабочую программу на новую, которая здесь опубликована. Вы можете использовать для этого раздел, подготовленный для поиска 23-ки. Только желательно, конечно, переименовать подпроект. А действующее Приложение 4 пусть досчитывается до конца. Как я уже говорила, задания будут использоваться из того же пакета, который был сгенерирован раньше; буду просто продолжать загрузку следующих заданий. Мой домашний ящик вроде в порядке, письма в него поступают. Но от вас писем я не получаю уже почти два месяца. Причина мне не известна. Поэтому пишу на форуме. Настоятельно прошу заменить рабочую программу в Приложении 4. Нет никакого смысла искать 15-ки в том диапазоне, в котором они все уже найдены другой командой. The TBEG BOINC project maintains the database https://boinc.tbrada.eu/spt/explore.php |
©2025 Natalia Makarova and Co