%I #76 Jan 06 2021 02:07:44
%S 1,3,7,21,31,93,127,217,381,651,889,2667,3937,8191,11811,24573,27559,
%T 57337,82677,131071,172011,253921,393213,524287,761763,917497,1040257,
%U 1572861,1777447,2752491,3120771,3670009,4063201,5332341,7281799,11010027,12189603
%N Numbers that are a product of distinct Mersenne primes (elements of A000668).
%C Or, numbers n such that the sum of the divisors of n is a power of 2, see A094502.
%C Or, numbers n such that the number of divisors of n and the sum of the divisors of n are both powers of 2.
%C n is a product of distinct Mersenne primes iff sigma(n) is a power of 2: see exercise in Sivaramakrishnan, or Shallit.
%C Sequence gives n > 1 such that sigma(n) = 2*phi(sigma(n)). - _Benoit Cloitre_, Feb 22 2002
%C The graph of this sequence shows a discontinuity at the 4097th number because there is a large relative gap between the 12th and 13th Mersenne primes, A000043. Other discontinuities can be predicted using A078426. - _T. D. Noe_, Oct 12 2006
%C Supersequence of A051281 (numbers n such that sigma(n) is a power of tau(n)). Conjecture: numbers n such that sigma(n) = tau(n)^(a/b), where a, b are integers >= 1. Example: sigma(93) = 128 = tau(93)^(7/2) = 4^(7/2). - _Jaroslav Krizek_, May 04 2013
%D J.-M. De Koninck & A. Mercier, 1001 Problèmes en Théorie Classique des Nombres, Problem 264 pp. 188, Ellipses Paris 2004.
%D R. Sivaramakrishnan, Classical Theory of Arithmetic Functions. Dekker, 1989.
%H Amiram Eldar, <a href="/A046528/b046528.txt">Table of n, a(n) for n = 1..10000</a> (terms 1..5000 from T. D. Noe)
%H Kevin S. Brown, <a href="https://www.mathpages.com/home/kmath048/kmath048.htm">Sum of Divisors Equals a Power of 2</a>.
%H C. D. H. Cooper, <a href="http://www.jstor.org/stable/2319455">Problem E 2493</a>, The American Mathematical Monthly, Vol. 81, No. 8 (1974), p. 902; W. J. Dodge, <a href="http://www.jstor.org/stable/2319819">solution</a>, ibid., Vol. 82, No. 8 (1975), pp. 855-856.
%H Jeffrey Shallit, <a href="http://www.jstor.org/stable/2691073">Problem 1319, Diophantine Equation, sigma(n) = 2^m</a>, Math. Magazine, 63 (1990), 129.
%H Eric Weisstein's World of Mathematics, <a href="http://mathworld.wolfram.com/DivisorFunction.html">Divisor Function</a>.
%F Sum_{n>=1} 1/a(n) = Product_{p in A000668} (1 + 1/p) = 1.5855588879... (A306204) - _Amiram Eldar_, Jan 06 2021
%e a(20) = 82677 = 3*7*31*127, whose sum of divisors is 131072 = 2^17;
%e a(27) = 1040257 = 127*8191, whose sum of divisors is 1048576 = 2^20.
%p mersennes:= [seq(numtheory:-mersenne([i]),i=1..10)]:
%p sort(select(`<`,map(convert,combinat:-powerset(mersennes),`*`),numtheory:-mersenne([11]))); # _Robert Israel_, May 01 2016
%t {1}~Join~TakeWhile[Times @@@ Rest@ Subsets@ # // Sort, Function[k, k <= Last@ #]] &@ Select[2^Range[0, 31] - 1, PrimeQ] (* _Michael De Vlieger_, May 01 2016 *)
%o (PARI) isok(n) = (n==1) || (ispower(sigma(n), , &r) && (r==2)); \\ _Michel Marcus_, Dec 10 2013
%Y Cf. A000668, A000043, A056652 (product of Mersenne primes), A306204.
%K nonn
%O 1,2
%A _Labos Elemer_
%E More terms from _Benoit Cloitre_, Feb 22 2002
%E Further terms from _Jon Hart_, Sep 22 2006
%E Entry revised by _N. J. A. Sloane_, Mar 20 2007
%E Three more terms from _Michel Marcus_, Dec 10 2013