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: 569
Credit: 8,342,279
RAC: 32,490
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: 569
Credit: 8,342,279
RAC: 32,490
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: 569
Credit: 8,342,279
RAC: 32,490
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: 569
Credit: 8,342,279
RAC: 32,490
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: 569
Credit: 8,342,279
RAC: 32,490
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

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

©2025 Natalia Makarova and Co