OFFSET
1,1
COMMENTS
The corresponding sequence for primes is A014663.
These pseudoprimes seem to be relatively rare: among the 118968378 base-2 Fermat pseudoprimes below 2^64 only 6292535 are terms of this sequence.
These pseudoprimes appear in a theorem by Rotkiewicz and Makowski (1966) about pseudoprimes that are products of two Mersenne numbers (see A367229).
LINKS
Amiram Eldar, Table of n, a(n) for n = 1..10000
Andrzej Rotkiewicz and Andrzej Makowski, On Pseudoprime Numbers of the Form M_p M_t, Elemente der Mathematik, Vol. 21 (1966), pp. 133-134.
MATHEMATICA
Select[2*Range[10^6] + 1, PowerMod[2, # - 1, #] == 1 && CompositeQ[#] && OddQ[MultiplicativeOrder[2, #]] &]
PROG
(PARI) is(n) = n > 1 && n % 2 && Mod(2, n)^(n-1) == 1 && !isprime(n) && znorder(Mod(2, n)) % 2;
CROSSREFS
KEYWORD
nonn
AUTHOR
Amiram Eldar, Nov 11 2023
STATUS
approved