OFFSET
0,1
REFERENCES
Steven R. Finch, Mathematical Constants, Cambridge University Press, 2003, Section 2.11 Abundant numbers density constant p. 126.
LINKS
Eric Weisstein's World of Mathematics, e-Divisor
FORMULA
B = lim_{N->inf} (1/N^2) * Sum_{n=1..N} sigma^(e)(n), where sigma^(e)(n) is the sum of all exponential divisors of n.
B = (1/2) * Product_{p prime} (1 + 1/(p*(p^2 - 1)) - 1/(p^2 - 1) + (1 - 1/p)*Sum_{k>=2} p^k/(p^(2k)-1)).
EXAMPLE
0.5682854937...
MATHEMATICA
digits = 10; maxPi = 10^5;
B = (1/2)*Product[1 + 1/(p*(p^2-1)) - 1/(p^2-1) + (1-1/p)*((Log[-(1/p^2)] - Log[1/p^2] + QPolyGamma[0, -(Log[-(1/p^2)]/Log[p]), p] - QPolyGamma[0, -(Log[1/p^2]/Log[p]), p])/(2*Log[p])), {p, Prime[Range[maxPi]]}];
RealDigits[N[B] // Chop, 10, digits][[1]]
$MaxExtraPrecision = 2000; Do[m = 2000; Clear[f]; f[p_] := (1 + 1/(p*(p^2 - 1)) - 1/(p^2 - 1) + (1 - 1/p)*Sum[p^k/(p^(2 k) - 1), {k, 2, kmax}]); c = Rest[CoefficientList[Series[Log[f[1/x]], {x, 0, m}], x]]; Print[f[2]/2 * Exp[N[Sum[Indexed[c, n]*(PrimeZetaP[n] - 1/2^n), {n, 2, m}], 112]]], {kmax, 100, 1000, 100}] (* Vaclav Kotesovec, Jun 19 2020 *)
CROSSREFS
KEYWORD
nonn,cons
AUTHOR
Jean-François Alcover, Jul 29 2016
EXTENSIONS
More digits from Robert G. Wilson v, Feb 25 2019
More digits from Vaclav Kotesovec, Jun 19 2020
STATUS
approved