OFFSET
1,1
COMMENTS
a(24) > 7*10^10.
The corresponding pseudoprimes a(n)^2 + 1 are 46657, 435601, 2202257, 2433601, 67371265, ...
a(25) > 7.5*10^11. - Giovanni Resta, Mar 15 2020
EXAMPLE
216 is a term since 216^2 + 1 = 46657 is a Fermat pseudoprime to base 3.
MATHEMATICA
Select[Range[10^3], CompositeQ[#^2 + 1] && PowerMod[3, #^2, #^2 + 1] == 1 &]
CROSSREFS
KEYWORD
nonn,more
AUTHOR
Amiram Eldar, Mar 14 2020
EXTENSIONS
a(24) from Giovanni Resta, Mar 15 2020
STATUS
approved