Вот программа подана ИИ на анализ.
Код:
{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);
vmy=vector(50);
pat9=vector(9); pat11=vector(11); pat13=vector(13); pat15=vector(15);
pat1=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");
k=0;
for(i43=1,28,
for(i47=1,32,
for(i53=1,38,
for(i59=1,44,
for(i61=1,46,
bpt=lift(chinese([Mod(1,2),Mod(2,3),Mod(4,5),Mod(2,7),Mod(9,11),Mod(4,13),
Mod(rost17,17),Mod(rost19,19),Mod(rost23,23),Mod(rost29,29),Mod(rost31,31),
Mod(rost37,37),Mod(rost41,41),Mod(v43[i43],43),Mod(v47[i47],47),
Mod(v53[i53],53),Mod(v59[i59],59),Mod(v61[i61],61)]));
k++;
\\кортеж 9
if(bpt>=2079914861571286679 && ispseudoprime(bpt+pt[4]) && 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); );
pat1=vector(15,i,(pat15[i]==pt[i]));
vlds=vecsum(pat1);
if(vlds>9,
code=fromdigits(pat1[2..14],2);
w1=strprintf("%d: %d\n %d\n",vmy[1],pat15,code);
filewrite(fout,w1);
print(w1); );
))))))); \\ if l==9
);\\ if ispseudoprime
);););););
print("sgenerirovano dobavok ",k);
w1=strprintf("end ");
print(w1);
filewrite(fout,w1);
fileclose(fout)
}
ИИ сказал:
Я посмотрел на вашу программу на PARI/GP (файл 15_61_0period_boinc.txt ). Она
довольно тяжёлая: вложенные циклы по ~28×32×38×44×46 ≈ 67 миллионов
итераций, и внутри каждой проверяются условия на простоту, формируются кортежи
и сравниваются с шаблонами.
Вот что можно реально ускорить/оптимизировать:
1 Работа с китайской системой сравнений
Сейчас каждый раз считается
gp
Copy code
bpt = lift(chinese([...]));
Это дорогая операция.
Можно вынести все константные модули ( 2,3,5,7,11,13,17,19,... ) в один общий
модуль
gp
Copy code
M = lcm([2,3,5,7,11,13,17,19,23,29,31,37,41,43,47,53,59,61]);
и заранее посчитать обратные элементы для каждого модуля (CRT через Garner или
precomputed inverses). Тогда вместо chinese вы делаете быстрый ручной CRT на
базе предвычислений.
2 Использование primepi вместо forprime
Сейчас у вас:
forprime(p=bpt+pt[4], bpt+pt[1S2k]i,p tlo+ c+o;n tvemnyt[l]=p; );
Но forprime — это перебор всех простых в интервале. Если интервал узкий и
нужно только число простых и их позиции, лучше:
Сначала проверить, что количество простых равно нужному
( primepi(bpt+pt[12]) - primepi(bpt+pt[4]-1) == 9 ),
И только если совпало — тогда уже реально получать список простых через
vecPrimes(bpt+pt[4], bpt+pt[12]) .
Это уберёт 99% ненужных проходов.
3 Избавиться от повторных ispseudoprime
Много проверок вида:
gp
if(ispseudoprime(bpt+pt[4]) && ispseudoprime(bpt+pt[12]), ...)
После этого снова запускается forprime , который всё равно проверяет простоту.
Оптимизация:
Сначала считать vpr = vecPrimes(start,end) ,
Потом уже проверять границы и нужные позиции.
4 Векторы шаблонов ( pt9 , pt11 , …)
Сейчас для каждого кандидата вы делаете:
gp
for(m=2,9, pat9[m]=vmy[m]-vmy[1]);
if(pat9==pt9, ...)
Это можно свернуть в одну проверку через vector_equal или напрямую через
сравнение разностей (чтобы избежать выделения новых массивов каждый раз).
Например:
Copy code
gp
Skip to content
ok = 1;
for(m=2,9, if(vmy[m]-vmy[1] != pt9[m], ok=0; break));
if(ok, ...)
Это уменьшит расход памяти и GC.
5 Вынести файлы из цикла
filewrite на каждой итерации — очень медленно. Лучше собирать результаты в
строковый буфер (например, вектор строк), и писать в файл блоками (раз в 1000
записей или в конце цикла).
6 Параллелизация
PARI/GP умеет параллельные циклы через parfor . У вас есть 5 внешних циклов
( i43,i47,i53,i59,i61 ).
Можно заменить на что-то вроде:
gp
Copy code
parfor(i43=1,28,
/* внутри вложенные циклы */
)
или ещё лучше — распараллелить по самому крупному измерению ( i61 ).
7 Микрооптимизации
ispseudoprime можно заменить на isprime — он тоже использует
детерминированный тест, но может быть чуть быстрее для чисел в вашем
диапазоне.
vecsum(pat1) можно заменить на sum(i=1,#pat1, pat1[i]) . Иногда быстрее,
так как не создаёт новый объект.
fromdigits(pat1[2..14],2) — лучше вручную посчитать сдвигами.
В целом: самый большой выигрыш даст (1) предвычисленный CRT вместо
вызова chinese , (2) использование primepi для отсечения неверных кандидатов
до вызова forprime , и (3) параллелизация внешних циклов. Это может ускорить
программу в десятки раз.
Хотите, я могу переписать ваш кусок с проверкой кортежей (например, для 9 и 11) в
оптимизированном стиле, чтобы вы сразу видели прирост?
Буду признателен за высказывание мнений и советов.
Или, наоборот, отвратные 
