A218483 Fermat pseudoprimes to base 2 which are congruent to 1 (mod 8).
561, 1105, 1729, 1905, 2465, 4033, 4369, 4681, 6601, 8321, 8481, 10585, 11305, 12801, 15841, 16705, 18705, 18721, 23001, 23377, 25761, 30121, 30889, 31417, 31609, 33153, 34945, 39865, 41041, 41665, 46657, 52633, 62745, 65281, 74665, 75361, 83665, 85489
Offset: 1
Keywords
Links
- Amiram Eldar, Table of n, a(n) for n = 1..10000
- Eric Weisstein's World of Mathematics, Fermat Pseudoprime.
- Eric Weisstein's World of Mathematics, Poulet Number.
Programs
-
Maple
select(t -> 2 &^ t mod t = 2 and not isprime(t), [seq(1+8*j,j=0..10^5)]); # Robert Israel, Dec 07 2014
-
Mathematica
Select[8 * Range[10^4] + 1, PowerMod[2, # - 1, #] == 1 && CompositeQ[#] &] (* Amiram Eldar, Mar 30 2021 *)
-
PARI
is(n)=n%8==1 && Mod(2,n)^n==2 && !isprime(n) \\ Charles R Greathouse IV, Dec 07 2014
Extensions
Corrected by Charles R Greathouse IV, Dec 07 2014
New name from Charles R Greathouse IV, Dec 07 2014
Comments