A244065 Pseudoprimes to base 3 that are not squarefree.

121, 3751, 4961, 7381, 11011, 29161, 32791, 142901, 228811, 239701, 341341, 551881, 566401, 595441, 671671, 784201, 856801, 1016521, 1237951, 1335961, 1433971, 1804231

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

Cf. A005935, A158358, A014127.

Sequence in context: A206629 A014758 A069600 * A231649 A212802 A272775

Adjacent sequences: A244062 A244063 A244064 * A244066 A244067 A244068

Keyword

nonn

Author

Felix Fröhlich, Jun 19 2014

Status

approved