A020161 Pseudoprimes to base 33.

4, 8, 16, 32, 85, 496, 545, 703, 1057, 1105, 1417, 1649, 1729, 1853, 2465, 2501, 2821, 3368, 4033, 4492, 4681, 5461, 5713, 5833, 6533, 6601, 7861, 8911, 9061, 9073, 9265, 10585, 11305, 11359, 12209, 12403, 13741, 15841, 16589, 16745, 17968, 18103, 19909

Offset

1,1

Comments

Composite numbers n such that 33^(n-1) == 1 (mod n).

Links

T. D. Noe, Table of n, a(n) for n = 1..1000

Index entries for sequences related to pseudoprimes

Mathematica

base = 33; t = {}; n = 1; While[Length[t] < 100, n++; If[! PrimeQ[n] && PowerMod[base, n-1, n] == 1, AppendTo[t, n]]]; t (* T. D. Noe, Feb 21 2012 *)
Select[Range[20000], CompositeQ[#]&&PowerMod[33, #-1, #]==1&] (* Harvey P. Dale, Oct 05 2015 *)

Crossrefs

Cf. A001567 (pseudoprimes to base 2).

Sequence in context: A338313 A049932 A272712 * A068890 A073856 A019458

Adjacent sequences: A020158 A020159 A020160 * A020162 A020163 A020164

Keyword

nonn

Author

David W. Wilson

Status

approved