|
PROGRAM
|
(PARI) twiprimesl2(n, a) = { local(pr, x, y, j); for(a=0, n, for(j=1, n, pr=1; for(x=1, j, pr*=twinl(x); ); y=pr+2^a; if(ispseudoprime(y), print1(y", "); break ) ) ) } twinl(n) = \The n-th lower twin prime { local(c, x); c=0; x=1; while(c<n, if(isprime(prime(x)+2), c++); x++; ); return(prime(x-1)) }
|