1,1

R. J. Mathar, Table of n, a(n) for n = 1..300

strongPseudoprimeQ[b_, n_] := Module[{rems = Table[PowerMod[b, (n - 1)/2^expo, n], {expo, 0, IntegerExponent[n - 1, 2]}]}, (rems[[-1]] == 1 || MemberQ[rems, n - 1]) && PowerMod[b, n - 1, n] == 1]; max = 5000; Select[Complement[Range[2, max], Prime[Range[PrimePi[max]]]], strongPseudoprimeQ[80, #] &] (* Alonso del Arte, Aug 10 2018 *)

Sequence in context: A133636 A369564 A020208 * A330680 A069068 A051412

Adjacent sequences: A020303 A020304 A020305 * A020307 A020308 A020309

approved