Вернулся на шаг назад.

,


,

Программа для этого алгоритмма:
Код:
b=10^1; c=10^9;
v=[1,3,7,9]; e=0; n=0;
{forprime(pr = b, c, qr=pr+2;
if(isprime(qr)==1,
d=digits(pr); rp=fromdigits(Vecrev(d)); l=#d;
d1=digits(rp); w=d1[l];
d2=digits(qr); rq=fromdigits(Vecrev(d2));
d3=digits(rq); w1=d3[l]; fl=0;
for(k=1,4, if(w==v[k], next));
for(i=1,10, a=i-1;
x=pr*10^(l+1)+a*10^l+rp; x1=rp*10^(l+1)+a*10^l+pr;
if(isprime(x)==1, fl=1; e+=1; print1(x, ", "));
if(isprime(x1)==1, fl=1; e+=1; print1(x1, ", ")));
for(k=1,4, if(w1==v[k], next));
for(i=1,10, a=i-1;
y=qr*10^(l+1)+a*10^l+rq; y1=rq*10^(l+1)+a*10^l+qr;
if(isprime(y)==1, fl=1; e+=1; print1(x, ", "));
if(isprime(y1)==1, fl=1; e+=1; print1(y1, ", ")));
if(fl==0, n+=1; print("*** (", pr, ", ", qr, ")")); if(fl==1, print(" (", pr, ", ", qr, ")"));
));
print; print(" es=", e, ", no=", n, " ");
}
Обозначение «***» - для простых близнецов, которым не найдено решения с помощь этого алгоритма.

- количество найденнных решений,

-количесво не найденных решений.
До

До

.
До

.
До

7.
До

.
До

.
До миллиона: es=6379728, no=581199, es/ no= 11.
Отношение

уменьшается, но функця не линейная. Интересно станет ли

меньше 1.