OFFSET
1,11
LINKS
Amiram Eldar, Table of n, a(n) for n = 1..64
Jan Feitsma and William F. Galway, Tables of pseudoprimes and related data.
EXAMPLE
Below 2^10 = 1024 there is only one Euler-Jacobi pseudoprime, 561. Therefore a(10) = 1.
MATHEMATICA
ejpspQ[n_] := CompositeQ[n] && PowerMod[2, (n - 1)/2, n] == Mod[JacobiSymbol[2, n], n]; s = {}; c = 0; p = 2; n = 1; Do[If[ejpspQ[n], c++]; If[n > p, AppendTo[s, c]; p *= 2], {n, 1, 2^20 + 1, 2}]; s
CROSSREFS
KEYWORD
nonn
AUTHOR
Amiram Eldar, Nov 08 2019
STATUS
approved