OFFSET
1,1
COMMENTS
The subsequence of primes in this sequence is A128667.
LINKS
Giovanni Resta, Table of n, a(n) for n = 1..140 (terms < 3*10^11)
MATHEMATICA
Select[Range@ 1000000, Mod[13^EulerPhi[#], #^2] == 1 &] (* Michael De Vlieger, Jul 27 2015 *)
PROG
(PARI) for(n=2, 1e9, if(Mod(13, n^2)^(eulerphi(n))==1, print1(n, ", ")))
CROSSREFS
KEYWORD
nonn
AUTHOR
Felix Fröhlich, Jul 26 2015
STATUS
approved