cp's OEIS Frontend

This is a front-end for the Online Encyclopedia of Integer Sequences, made by Christian Perfect. The idea is to provide OEIS entries in non-ancient HTML, and then to think about how they're presented visually. The source code is on GitHub.

A308051 Decimal expansion of lim_{m->oo} (sqrt(log(m))/m^2) Sum_{k=1..m} sigma(k)/d(k), where d(k) is the number of divisors of k (A000005) and sigma(k) is their sum (A000203).

Original entry on oeis.org

3, 5, 6, 9, 0, 4, 9, 6, 5, 2, 4, 9, 9, 5, 7, 0, 7, 6, 1, 2, 2, 0, 0, 5, 3, 0, 2, 0, 1, 3, 9, 9, 6, 4, 5, 9, 1, 3, 6, 0, 6, 6, 6, 8, 2, 6, 2, 5, 7, 3, 8, 4, 4, 2, 9, 6, 8, 7, 8, 8, 0, 2, 0, 1, 2, 7, 7, 4, 3, 4, 4, 2, 1, 4, 1, 8, 7, 2, 1, 3, 8, 5, 5, 3, 2, 1, 5
Offset: 0

Views

Author

Amiram Eldar, May 10 2019

Keywords

Examples

			0.35690496524995707612200530201399645913606668262573...
		

References

  • V. I. Arnold, Dynamics, Statistics, and Projective Geometry of Galois Fields, Cambridge University Press, Cambridge, 2011, p. 78.

Crossrefs

Programs

  • Mathematica
    $MaxExtraPrecision = 1000; m = 1000; f[x_] := Log[1 + x]/x/Sqrt[1 - x]; c = Rest[CoefficientList[Series[Log[f[x]], {x, 0, m}], x]]; RealDigits[(1/2/ Sqrt[Pi])*Exp[NSum[Indexed[c, k]*PrimeZetaP[k], {k, 2, m}, NSumTerms -> m, WorkingPrecision -> m]], 10, 100][[1]]

Formula

Equals (1/(2*sqrt(Pi))) * Product_{p prime} p^(3/2) * log(1 + 1/p) / sqrt(p-1).