OFFSET
1,1
EXAMPLE
a(1) = 2 because 2 and (2-1)*2^2 + 1 = 5 are both prime.
a(2) = 3 because 3 and (3-1)*2^3 + 1 = 17 are both prime.
a(3) = 7 because 7 and (7-1)*2^7 + 1 = 769 are both prime.
MATHEMATICA
Select[Prime[Range[9000]], PrimeQ[(#-1)2^#+1]&] (* Harvey P. Dale, Jan 19 2012 *)
PROG
(PARI) forprime(n=1, 1e4, if(ispseudoprime((n-1)<<n+1), print1(n", "))) \\ Charles R Greathouse IV, Oct 09 2011
CROSSREFS
KEYWORD
nonn
AUTHOR
Juri-Stepan Gerasimov, Oct 02 2011
EXTENSIONS
a(7) corrected by Michael S. Branicky, May 14 2025
STATUS
approved
