OFFSET
1,1
COMMENTS
No further terms up to 10^6.
No further terms up to 10^11. - Luke March, Jul 16 2014
EXAMPLE
a(1) = 997, since prime(997) = 7883 and 7883^996 % 997^2 = 1.
PROG
(PARI) forprime(p=2, 10^6, q=prime(p); if(Mod(q, p^2)^(p-1)==1, print1(p, ", ")))
CROSSREFS
KEYWORD
nonn,hard,more
AUTHOR
Felix Fröhlich, Jun 29 2014
STATUS
approved