login
Strong pseudoprimes to base 80.
1

%I #20 Aug 12 2024 13:32:16

%S 9,27,49,81,169,301,481,2107,3439,5461,6253,6401,7107,8261,10349,

%T 16609,18361,23281,25351,29089,46531,49051,55231,58969,75851,79501,

%U 113401,115231,117569,119341,123139,138881,146611,157753,158701,188501,200911,203941,205927

%N Strong pseudoprimes to base 80.

%H R. J. Mathar, <a href="/A020306/b020306.txt">Table of n, a(n) for n = 1..300</a>

%H <a href="/index/Ps#pseudoprimes">Index entries for sequences related to pseudoprimes</a>

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

%K nonn

%O 1,1

%A _David W. Wilson_