Message boards : Number crunching : Модифицированная рабочая программа для Приложения 4
Message board moderation
| Author | Message |
|---|---|
Natalia MakarovaSend message Joined: 24 Dec 24 Posts: 2187 Credit: 16,991,036 RAC: 1,267 |
Публикую модифицированную рабочую программу для Приложения 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: 2187 Credit: 16,991,036 RAC: 1,267 |
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: 2187 Credit: 16,991,036 RAC: 1,267 |
Для всех! Вы можете протестировать опубликованную рабочую программу, записав в файл 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: 2187 Credit: 16,991,036 RAC: 1,267 |
На форуме 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: 2187 Credit: 16,991,036 RAC: 1,267 |
Здравствуйте, termit! Прошло две недели с опубликования стартового сообщения о замене рабочей программы в Приложении 4. Я не получила ответа ни здесь, ни в личной переписке. Я прекращаю загрузку заданий для Приложения 4. Предлагаю вам изменить рабочую программу на новую, которая здесь опубликована. Вы можете использовать для этого раздел, подготовленный для поиска 23-ки. Только желательно, конечно, переименовать подпроект. А действующее Приложение 4 пусть досчитывается до конца. Как я уже говорила, задания будут использоваться из того же пакета, который был сгенерирован раньше; буду просто продолжать загрузку следующих заданий. Мой домашний ящик вроде в порядке, письма в него поступают. Но от вас писем я не получаю уже почти два месяца. Причина мне не известна. Поэтому пишу на форуме. Настоятельно прошу заменить рабочую программу в Приложении 4. Нет никакого смысла искать 15-ки в том диапазоне, в котором они все уже найдены другой командой. The TBEG BOINC project maintains the database https://boinc.tbrada.eu/spt/explore.php |
Natalia MakarovaSend message Joined: 24 Dec 24 Posts: 2187 Credit: 16,991,036 RAC: 1,267 |
Всем привет и пожелание всего самого лучшего! Я рискнула разморозить Приложение 4 - Calc15Tuples. Оно было остановлено мной, я перестала загружать в него задания. Я просила termit перезапустить это Приложение с новой рабочей программой, которую модифицировала. Смотрите сообщения выше. Прошёл почти год. По неведомым мне причинам termit не стал перезапускать это Приложение. Я решила попробовать возобновить работу Приложения, просто загрузив новую порцию заданий. Как отмечено выше, все задания, сгенерированные ранее, годятся и для модифицированной рабочей программы. Я просто изменила диапазон поиска. Более того, БД всех заданий для Приложения у меня хранится. Это к тому, если БД заданий на сервере пропала или случайно удалена. Я могу постоянно загружать задания из своего архива. БД для этого Приложения огромная! termit, пожалуйста, обратите внимание на это Приложение. Я загрузила сейчас из БД заданий пакет 568.7z. Если задания загрузились нормально, пусть Приложение работает. Только центральные 15-ки в текущем диапазоне все уже найдены командой на форуме dxdy.ru. Однако не найдены центральные 11-ки и центральные 13-ки, а их в разы больше, чем центральных 15-к. И их никто не искал и не ищет! Кроме того, я запустила модифицированную программу в ручном проекте, обрабатываю пакет заданий 567.7z. Но у меня работает всего один поток The TBEG BOINC project maintains the database https://boinc.tbrada.eu/spt/explore.php |
Natalia MakarovaSend message Joined: 24 Dec 24 Posts: 2187 Credit: 16,991,036 RAC: 1,267 |
Приложение 4 работает! Результаты получены и проверены. И опубликованы https://boinc.mak.termit.me/odlk2025/forum_thread.php?id=47&postid=2311 Спасибо большое всем, кто подключился к этому Приложению! Это самое эффективное Приложение из трёх активных Приложений. The TBEG BOINC project maintains the database https://boinc.tbrada.eu/spt/explore.php |
©2026 Natalia Makarova and Co