OFFSET
1,1
PROG
(PARI) is(n, k=primepi(n))=ispseudoprime((2^k-n)*2^k-1) && isprime(n)
c=0; forprime(p=2, 1e9, if(is(p, c++), print1(p", "))) \\ Charles R Greathouse IV, Nov 25 2014
CROSSREFS
KEYWORD
nonn
AUTHOR
Pierre CAMI, Jul 09 2011
EXTENSIONS
Incorrect a(34) removed by Charles R Greathouse IV, Nov 25 2014
a(34) from Charles R Greathouse IV, Dec 22 2014
STATUS
approved