login

Year-end appeal: Please make a donation to the OEIS Foundation to support ongoing development and maintenance of the OEIS. We are now in our 61st year, we have over 378,000 sequences, and we’ve reached 11,000 citations (which often say “discovered thanks to the OEIS”).

Decimal expansion of Product_{p>=3} (1+1/p) over the Mersenne primes.
3

%I #14 Jan 06 2021 02:05:03

%S 1,5,8,5,5,5,8,8,8,7,9,2,5,6,3,8,7,7,6,9,7,8,6,3,7,0,2,3,2,1,9,2,3,8,

%T 4,7,6,0,6,9,4,0,5,8,6,7,9,4,7,0,2,8,1,1,3,2,9,8,1,2,6,7,8,9,2,8,8,5,

%U 9,7,5,4,5,7,6,7,8,5,5,6,9,0,5,3,5,0,0,7,9,1,1,7,9,9,3,5,6,1,9,5

%N Decimal expansion of Product_{p>=3} (1+1/p) over the Mersenne primes.

%C This is equal to Product_{q>=1} (1-1/2^q)^(-1) over all q with 2^q - 1 a Mersenne prime.

%H Tomohiro Yamada, <a href="/A306204/b306204.txt">Table of n, a(n) for n = 1..99</a>

%H Tomohiro Yamada, <a href="http://mathematica-pannonica.ttk.pte.hu/index_elemei/mp19-1/MP19-1(2008)pp37-47.pdf">Unitary super perfect numbers</a>, Mathematica Pannonica, Volume 19, No. 1, 2008, pp. 37-47, using this constant with only a rough upper bound (4/3)*exp(4/21) < 1.6131008.

%F Equals Sum_{n>=1} 1/A046528(n). - _Amiram Eldar_, Jan 06 2021

%e Decimal expansion of (4/3) * (8/7) * (32/31) * (128/127) * (8192/8191) * (131072/131071) * (524288/524287) * ... = 1.585558887...

%o (PARI) t=1.0;for(i=1,500,p=2^i-1;if(isprime(p),t=t*(p+1)/p))

%Y Cf. A065446 (the corresponding product over all Mersenne numbers, prime or composite).

%Y Cf. A173898 (the sum of reciprocals of the Mersenne primes).

%Y Cf. A065442 (the sum of reciprocals of the Mersenne numbers, prime or composite).

%Y Cf. A046528.

%K nonn,cons

%O 1,2

%A _Tomohiro Yamada_, Jan 29 2019