# Greetings from The On-Line Encyclopedia of Integer Sequences! http://oeis.org/

Search: id:a361013
Showing 1-1 of 1

%I A361013 #10 Feb 28 2023 07:46:22
%S A361013 2,9,6,0,0,8,0,3,0,2,0,2,4,9,4,1,4,1,0,4,8,1,8,2,0,4,7,8,1,1,0,8,9,4,
%T A361013 6,9,3,9,2,8,4,3,9,0,9,5,9,2,5,1,6,3,4,1,1,9,6,7,5,0,4,4,8,0,8,6,6,3,
%U A361013 3,9,3,5,7,8,7,3,7,3,8,2,4,9,5,8,4,6,2,6,7,3,8,5,0,1,0,8,0,5,1,7,8,6,0,6,6
%N A361013 Decimal expansion of a constant related to the asymptotics of A361012.
%F A361013 Equals limit_{n->oo} A361012(n) / n.
%F A361013 Equals Product_{p prime} (1 + Sum_{e>=2} (sigma(e) - sigma(e-1)) / p^e), where sigma = A000203.
%e A361013 2.960080302024941410481820478110894693928439095925163411967504480866339...
%t A361013 $MaxExtraPrecision = 1000; smax = 500; Do[Clear[f]; f[p_] := 1 + Sum[(DivisorSigma[1, e] - DivisorSigma[1, e-1])/p^e, {e, 2, emax}]; cc = Rest[CoefficientList[Series[Log[f[1/x]], {x, 0, smax}], x, smax + 1]]; Print[f[2] * f[3] * f[5] * f[7] * Exp[N[Sum[cc[[n]]*(PrimeZetaP[n] - 1/2^n - 1/3^n - 1/5^n - 1/7^n), {n, 2, smax}], 120]]], {emax, 100, 1000, 100}]
%Y A361013 Cf. A361012, A327837, A327838.
%K A361013 nonn,cons
%O A361013 1,1
%A A361013 _Vaclav Kotesovec_, Feb 28 2023

# Content is available under The OEIS End-User License Agreement: http://oeis.org/LICENSE