A020306 Strong pseudoprimes to base 80.

9, 27, 49, 81, 169, 301, 481, 2107, 3439, 5461, 6253, 6401, 7107, 8261, 10349, 16609, 18361, 23281, 25351, 29089, 46531, 49051, 55231, 58969, 75851, 79501, 113401, 115231, 117569, 119341, 123139, 138881, 146611, 157753, 158701, 188501, 200911

Offset

1,1

Links

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

Index entries for sequences related to pseudoprimes

Mathematica

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 *)

Crossrefs

Sequence in context: A133636 A369564 A020208 * A330680 A069068 A051412

Adjacent sequences: A020303 A020304 A020305 * A020307 A020308 A020309

Keyword

nonn

Author

David W. Wilson

Status

approved