OFFSET
1,1
LINKS
Amiram Eldar, Table of n, a(n) for n = 1..366 (terms below 10^11; terms 1..146 from Robert G. Wilson v)
Eric Weisstein's World of Mathematics, Fermat Pseudoprime.
MATHEMATICA
Do[ f=PowerMod[ 41, 2n-1, 2n ]; If[ f==1, Print[ 2n ] ], {n, 2, 500000} ]
lst = {}; Do[ If[ PowerMod[41, 2n - 1, 2n] == 1, AppendTo[lst, 2n]], {n, 2, 2^31}]; lst (* Robert G. Wilson v, Jun 01 2007 *)
PROG
(PARI) is(k) = k > 2 && !(k % 2) && Mod(41, k)^(k-1) == 1; \\ Amiram Eldar, Sep 29 2024
CROSSREFS
Cf. A020169 (pseudoprimes to base 41).
Cf. A006935 (even pseudoprimes (or primes) to base 2: n divides 2^n - 2, n even).
Cf. A130433 (even pseudoprimes to base 3).
Cf. A090082 (even pseudoprimes to base 5).
KEYWORD
nonn
AUTHOR
Alexander Adamchuk, May 26 2007
EXTENSIONS
More terms from Robert G. Wilson v, Jun 01 2007
STATUS
approved