OFFSET
1,1
COMMENTS
Must be divisible by the square of a Mirimanoff prime, A014127. - Charles R Greathouse IV, Jun 21 2014
LINKS
Felix Fröhlich and Charles R Greathouse IV, Table of n, a(n) for n = 1..10000 (first 196 terms from Fröhlich)
PROG
(PARI) for(n=2, 10^9, if(!isprime(n) && Mod(3, n)^(n-1)==1 && !issquarefree(n), print1(n, ", ")))
(PARI) list(lim)=my(M=[11, 1006003], v=List(), p2); for(i=1, #M, p2=M[i]^2; forstep(n=p2, lim, p2, if(Mod(3, n)^(n-1)==1, listput(v, n)))); Set(v) \\ Good for lim <= 9.4 * 10^29; Charles R Greathouse IV, Jun 21 2014
CROSSREFS
KEYWORD
nonn
AUTHOR
Felix Fröhlich, Jun 19 2014
STATUS
approved