A195734 G.f.: exp( Sum_{n>=1} (2*sigma(n^2) - sigma(n)^2) * x^n/n ).

1, 1, 3, 6, 11, 22, 40, 72, 123, 215, 363, 605, 991, 1618, 2598, 4139, 6537, 10229, 15871, 24476, 37487, 56995, 86177, 129531, 193662, 287992, 426254, 627841, 920708, 1344331, 1954987, 2831688, 4086168, 5875087, 8417724, 12020250, 17108958, 24275947, 34340966

Offset

0,3

Links

Vaclav Kotesovec, Table of n, a(n) for n = 0..10000 (terms 0..1000 from Paul D. Hanna)

Formula

Logarithmic derivative equals A195735.

log(a(n)) ~ 3*(5*zeta(3)*(12 - Pi^2))^(1/3) * n^(2/3) / (2*Pi^(2/3)). - Vaclav Kotesovec, Oct 31 2024

Example

G.f.: A(x) = 1 + x + 3*x^2 + 6*x^3 + 11*x^4 + 22*x^5 + 40*x^6 + 72*x^7 +...
where
log(A(x)) = x + 5*x^2/2 + 10*x^3/3 + 13*x^4/4 + 26*x^5/5 + 38*x^6/6 + 50*x^7/7 + 29*x^8/8 +...+ A195735(n)*x^n/n +...

Mathematica

nmax = 40; CoefficientList[Series[Exp[Sum[(2*DivisorSigma[1, k^2] - DivisorSigma[1, k]^2) * x^k / k, {k, 1, nmax}]], {x, 0, nmax}], x] (* Vaclav Kotesovec, Oct 31 2024 *)

Prog

(PARI) {a(n)=polcoeff(exp(sum(k=1, n, (2*sigma(k^2)-sigma(k)^2)*x^k/k)+x*O(x^n)), n)}

Crossrefs

Cf. A195735 (log), A156302, A156303, A000203.

Sequence in context: A192896 A115030 A228206 * A018177 A103322 A284474

Adjacent sequences: A195731 A195732 A195733 * A195735 A195736 A195737

Keyword

nonn

Author

Paul D. Hanna, Sep 22 2011

Status

approved