|
PROG
|
(PARI) twiprimesl(n, a) = { local(pr, x, y, j); for(j=1, n, pr=1; for(x=1, j, pr*=twinl(x); ); y=pr+a; if(ispseudoprime(y), print1(y", ") ) ) } 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)) }
|