OFFSET
1,2
COMMENTS
Subsequence of odd terms of A277285.
Sequence is infinite. In particular, k = (21p)^6 for primes p = 5 or p > 7 are such: d(k) = 7^3 divides 2^147 - 1, which in turn divides 2^k - 1. - Max Alekseyev, 30 Sep 2016
Note that there are also other forms. For example, (A002110(6)/2)^6 is a term.
LINKS
Amiram Eldar, Table of n, a(n) for n = 1..38
EXAMPLE
Odd number 729 is a term because 2^729-1 is divisible by d(729) = 7.
MATHEMATICA
Select[Range[1000], Mod[ 2^# - 1, DivisorSigma[0, # ]] == 0 && OddQ[#] &] (* G. C. Greubel, Oct 18 2016 *)
Select[Range[1, 1.2*10^6, 2]^2, # == 1 || PowerMod[2, #, DivisorSigma[0, #]] == 1 &] (* Amiram Eldar, May 12 2024 *)
PROG
(PARI) is(n) = (2^n-1) % numdiv(n) == 0 && n % 2 == 1;
(PARI) lista(kmax) = {my(d); forstep(k = 1, kmax, 2, d = vecprod(apply(x -> 2*x+1, factor(k)[, 2])); if(Mod(2, d)^(k^2) == 1, print1(k^2, ", "))); } \\ Amiram Eldar, May 12 2024
CROSSREFS
KEYWORD
nonn
AUTHOR
Altug Alkan, Oct 03 2016
EXTENSIONS
a(4)-a(9) from Giovanni Resta, Oct 03 2016
a(10)-a(16) from Amiram Eldar, May 12 2024
STATUS
approved