OFFSET
1,1
COMMENTS
a(12) > 23053. - Robert Israel, Jun 10 2014
FORMULA
a(n) = prime(A117587(n)). - Daniel Suteu, Jun 25 2022
EXAMPLE
5 is in this sequence because 5 and 2^5 - 5^2 = 7 are both prime.
MAPLE
select(p -> isprime(p) and isprime(2^p - p^2), {2} union {seq(2*i+1, i=1..2000)}); # Robert Israel, Jun 10 2014
MATHEMATICA
Select[Prime[Range[300]], PrimeQ[2^# - #^2] &] (* Alonso del Arte, May 27 2014 *)
PROG
(Magma) [p: p in PrimeUpTo(2200) | IsPrime(2^p - p^2)];
(PARI) isok(p) = isprime(p) && ispseudoprime(2^p - p^2); \\ Daniel Suteu, Jun 25 2022
CROSSREFS
KEYWORD
nonn,more
AUTHOR
Juri-Stepan Gerasimov, May 27 2014
EXTENSIONS
a(9) from Alonso del Arte, May 27 2014
a(10) from Alois P. Heinz, May 28 2014
a(11) from Robert Israel, Jun 10 2014
a(12) added by Daniel Suteu, Jun 25 2022
STATUS
approved