Thread 'Модифицированная рабочая программа для Приложения 4'

Message boards : Number crunching : Модифицированная рабочая программа для Приложения 4
Message board moderation

To post messages, you must log in.

AuthorMessage
ProfileNatalia Makarova
Volunteer moderator
Project scientist
Avatar

Send message
Joined: 24 Dec 24
Posts: 2187
Credit: 16,991,036
RAC: 1,267
Message 731 - Posted: 11 Oct 2025, 8:09:14 UTC
Last modified: 11 Oct 2025, 8:54:12 UTC

Публикую модифицированную рабочую программу для Приложения 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
ID: 731 · Rating: 0 · rate: Rate + / Rate - Report as offensive     Reply Quote
ProfileNatalia Makarova
Volunteer moderator
Project scientist
Avatar

Send message
Joined: 24 Dec 24
Posts: 2187
Credit: 16,991,036
RAC: 1,267
Message 732 - Posted: 11 Oct 2025, 8:29:42 UTC
Last modified: 11 Oct 2025, 8:32:13 UTC

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
ID: 732 · Rating: 0 · rate: Rate + / Rate - Report as offensive     Reply Quote
ProfileNatalia Makarova
Volunteer moderator
Project scientist
Avatar

Send message
Joined: 24 Dec 24
Posts: 2187
Credit: 16,991,036
RAC: 1,267
Message 733 - Posted: 11 Oct 2025, 8:37:53 UTC

Для всех!

Вы можете протестировать опубликованную рабочую программу, записав в файл 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
ID: 733 · Rating: 0 · rate: Rate + / Rate - Report as offensive     Reply Quote
ProfileNatalia Makarova
Volunteer moderator
Project scientist
Avatar

Send message
Joined: 24 Dec 24
Posts: 2187
Credit: 16,991,036
RAC: 1,267
Message 734 - Posted: 11 Oct 2025, 8:43:07 UTC
Last modified: 11 Oct 2025, 8:45:14 UTC

На форуме dxdy.ru утверждают, что опубликованная рабочая программа может быть ещё оптимизирована с достижением ускорения в 3-4 раза.

Если среди участников проекта есть программисты, хорошо знающие PARI/GP, прошу высказать предложения по оптимизации программы с целью ускорения.
Необходимо учитывать, что оптимизация не должна нарушать алгоритм программы.

Сделать это надо быстро, пока рабочая программа в Приложении 4 не заменена на новую.

Если у вас есть вопросы по программе, пожалуйста, задавайте их здесь или в ЛС.
The TBEG BOINC project maintains the database
https://boinc.tbrada.eu/spt/explore.php
ID: 734 · Rating: 0 · rate: Rate + / Rate - Report as offensive     Reply Quote
ProfileNatalia Makarova
Volunteer moderator
Project scientist
Avatar

Send message
Joined: 24 Dec 24
Posts: 2187
Credit: 16,991,036
RAC: 1,267
Message 760 - Posted: 26 Oct 2025, 9:04:49 UTC
Last modified: 26 Oct 2025, 9:08:13 UTC

Здравствуйте, termit!

Прошло две недели с опубликования стартового сообщения о замене рабочей программы в Приложении 4.
Я не получила ответа ни здесь, ни в личной переписке.

Я прекращаю загрузку заданий для Приложения 4.
Предлагаю вам изменить рабочую программу на новую, которая здесь опубликована.
Вы можете использовать для этого раздел, подготовленный для поиска 23-ки.
Только желательно, конечно, переименовать подпроект.
А действующее Приложение 4 пусть досчитывается до конца.

Как я уже говорила, задания будут использоваться из того же пакета, который был сгенерирован раньше; буду просто продолжать загрузку следующих заданий.

Мой домашний ящик вроде в порядке, письма в него поступают.
Но от вас писем я не получаю уже почти два месяца.
Причина мне не известна.
Поэтому пишу на форуме.

Настоятельно прошу заменить рабочую программу в Приложении 4.
Нет никакого смысла искать 15-ки в том диапазоне, в котором они все уже найдены другой командой.
The TBEG BOINC project maintains the database
https://boinc.tbrada.eu/spt/explore.php
ID: 760 · Rating: 0 · rate: Rate + / Rate - Report as offensive     Reply Quote
ProfileNatalia Makarova
Volunteer moderator
Project scientist
Avatar

Send message
Joined: 24 Dec 24
Posts: 2187
Credit: 16,991,036
RAC: 1,267
Message 2297 - Posted: 25 Aug 2026, 13:22:03 UTC
Last modified: 26 Aug 2026, 6:37:42 UTC

Всем привет и пожелание всего самого лучшего!

Я рискнула разморозить Приложение 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
ID: 2297 · Rating: 0 · rate: Rate + / Rate - Report as offensive     Reply Quote
ProfileNatalia Makarova
Volunteer moderator
Project scientist
Avatar

Send message
Joined: 24 Dec 24
Posts: 2187
Credit: 16,991,036
RAC: 1,267
Message 2310 - Posted: 27 Aug 2026, 2:41:24 UTC
Last modified: 27 Aug 2026, 3:03:15 UTC

Приложение 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
ID: 2310 · Rating: 0 · rate: Rate + / Rate - Report as offensive     Reply Quote

Message boards : Number crunching : Модифицированная рабочая программа для Приложения 4

©2026 Natalia Makarova and Co