A020161 - OEIS

The OEIS is supported by the many generous donors to the OEIS Foundation.

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

(list; graph; refs; listen; history; text; internal format)

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