login
Base-2 Fermat pseudoprimes k such that the multiplicative order of 2 modulo k is odd.
3

%I #11 Nov 12 2023 00:06:15

%S 2047,4681,15841,42799,52633,90751,220729,256999,271951,486737,514447,

%T 647089,741751,916327,1082401,1145257,1730977,1969417,2205967,2304167,

%U 2748023,2811271,2953711,2976487,3567481,4188889,4469471,4835209,4863127,5016191,5049001,5681809

%N Base-2 Fermat pseudoprimes k such that the multiplicative order of 2 modulo k is odd.

%C The corresponding sequence for primes is A014663.

%C These pseudoprimes seem to be relatively rare: among the 118968378 base-2 Fermat pseudoprimes below 2^64 only 6292535 are terms of this sequence.

%C These pseudoprimes appear in a theorem by Rotkiewicz and Makowski (1966) about pseudoprimes that are products of two Mersenne numbers (see A367229).

%H Amiram Eldar, <a href="/A367230/b367230.txt">Table of n, a(n) for n = 1..10000</a>

%H Andrzej Rotkiewicz and Andrzej Makowski, <a href="https://doi.org/10.5169/seals-24655">On Pseudoprime Numbers of the Form M_p M_t</a>, Elemente der Mathematik, Vol. 21 (1966), pp. 133-134.

%t Select[2*Range[10^6] + 1, PowerMod[2, # - 1, #] == 1 && CompositeQ[#] && OddQ[MultiplicativeOrder[2, #]] &]

%o (PARI) is(n) = n > 1 && n % 2 && Mod(2, n)^(n-1) == 1 && !isprime(n) && znorder(Mod(2, n)) % 2;

%Y Intersection of A001567 and A036259.

%Y A367231 is a subsequence.

%Y Cf. A002326, A014663, A306413, A367229.

%K nonn

%O 1,1

%A _Amiram Eldar_, Nov 11 2023