Charles R Greathouse IV, <a href="/A020229/b020229_1.txt">Table of n, a(n) for n = 1..24767</a> (first 752 terms from R. J. Mathar)
Charles R Greathouse IV, <a href="/A020229/b020229_1.txt">Table of n, a(n) for n = 1..24767</a> (first 752 terms from R. J. Mathar)
editing
approved
proposed
approved
editing
proposed
<a href="/index/Ps#pseudoprimes">Index entries for sequences related to pseudoprimes</a>
<a href="/index/Ps#pseudoprimes">Index entries for sequences related to pseudoprimes</a>
reviewed
editing
proposed
reviewed
editing
proposed
sppQ[n_?EvenQ, _] := False; sppQ[n_?PrimeQ, _] := False; sppQ[n_, b_] := (s = IntegerExponent[n-1, 2]; d = (n-1)/2^s; If[PowerMod[b, d, n] == 1, Return[True], Do[If[PowerMod[b, d*2^r, n] == n-1, Return[True]], {r, 0, s-1}]]); A020229 = {}; lst = {}; k = 3; While[k < 500000, If[sppQ[k, 3], Print[k]; AppendTo[lst, k]]; k += 2]; lst (* Jean-François Alcover, Oct 20 2011, after _R. J. Mathar _ *)
approved
editing
editing
approved