OFFSET
1,1
LINKS
Felix Fröhlich, Table of n, a(n) for n = 1..3119 (all terms up to 10^6)
EXAMPLE
5 is a term because 5 is prime and (2^5 mod 5^2) = 32 mod 25 = 7 also prime.
MATHEMATICA
Select[Prime[Range[500]], PrimeQ[PowerMod[2, #, #^2]] &] (* Alonso del Arte, Mar 28 2011 *)
PROG
(PARI) forprime(p=2, 10^3, if(isprime(2^p%p^2), print1(p, ", "))) \\ Felix Fröhlich, Jun 28 2014
CROSSREFS
KEYWORD
nonn
AUTHOR
Juri-Stepan Gerasimov, Mar 28 2011
STATUS
approved