login
Expansion of Product_{k>=1} (1 - x^k)^(-k^2).
28

%I #61 Sep 08 2022 08:44:47

%S 1,1,5,14,40,101,266,649,1593,3765,8813,20168,45649,101591,223654,

%T 486046,1045541,2225167,4692421,9804734,20318249,41766843,85218989,

%U 172628766,347338117,694330731,1379437080,2724353422,5350185097,10449901555,20304465729,39254599832

%N Expansion of Product_{k>=1} (1 - x^k)^(-k^2).

%C In general, if g.f. = Product_{k>=1} 1/(1 - x^k)^(c2*k^2 + c1*k + c0) and c2 > 0, then a(n) ~ exp(4*Pi * c2^(1/4) * n^(3/4) / (3*15^(1/4)) + c1*Zeta(3) / Pi^2 * sqrt(15*n/c2) + (Pi * 5^(1/4) * c0 / (2*3^(3/4) * c2^(1/4)) - 15^(5/4) * c1^2 * Zeta(3)^2 / (2*c2^(5/4) * Pi^5)) * n^(1/4) + c1/12 + 75 * c1^3 * Zeta(3)^3 / (c2^2 * Pi^8) - 5*c0 * c1 * Zeta(3) / (4*c2 * Pi^2) - c2*Zeta(3) / (4*Pi^2)) * Pi^(c1/12) * (c2/15)^(1/8 + c0/8 + c1/48) / (A^c1 * 2^((c0 + 3)/2) * n^(5/8 + c0/8 + c1/48)), where A is the Glaisher-Kinkelin constant A074962. - _Vaclav Kotesovec_, Nov 08 2017

%H Seiichi Manyama, <a href="/A023871/b023871.txt">Table of n, a(n) for n = 0..10000</a> (first 1001 terms from Alois P. Heinz)

%H G. Almkvist, <a href="https://projecteuclid.org/euclid.em/1047674152">Asymptotic formulas and generalized Dedekind sums</a>, Exper. Math., 7 (No. 4, 1998), pp. 343-359.

%H Vaclav Kotesovec, <a href="http://arxiv.org/abs/1509.08708">A method of finding the asymptotics of q-series based on the convolution of generating functions</a>, arXiv:1509.08708 [math.CO], Sep 30 2015, p. 21.

%F a(n) = 1/n*Sum_{k=1..n} a(n-k)*sigma_3(k), n > 0, a(0)=1, where sigma_3(n) = A001158(n) = sum of cubes of divisors of n. - _Vladeta Jovovic_, Jan 20 2002

%F G.f.: Prod_{n>=1} exp(sigma_3(n)*x^n/n), where sigma_3(n) is the sum of cubes of divisors of n (=A001158(n)). - _N-E. Fahssi_, Mar 28 2010

%F G.f. (conjectured): 1/prod(n>=1, E(x^n)^J2(n))) where E(x) = prod(n>=1,1-x^n) and J2(n) = A007434(n). - _Joerg Arndt_, Jan 25 2011

%F a(n) ~ exp(4 * Pi * n^(3/4) / (3^(5/4) * 5^(1/4)) - Zeta(3) / (4*Pi^2)) / (2^(3/2) * 15^(1/8) * n^(5/8)), where Zeta(3) = A002117 = 1.2020569031595942853997... . - _Vaclav Kotesovec_, Feb 27 2015

%p with(numtheory):

%p a:= proc(n) option remember; `if`(n=0, 1,

%p add(add(d*d^2, d=divisors(j)) *a(n-j), j=1..n)/n)

%p end:

%p seq(a(n), n=0..35); # _Alois P. Heinz_, Nov 02 2012

%t max = 31; Series[ Product[ 1/(1-x^k)^k^2, {k, 1, max}], {x, 0, max}] // CoefficientList[#, x]& (* _Jean-François Alcover_, Mar 05 2013 *)

%o (PARI) m=40; x='x+O('x^m); Vec(prod(k=1, m, 1/(1-x^k)^k^2)) \\ _G. C. Greubel_, Oct 29 2018

%o (Magma) m:=40; R<x>:=PowerSeriesRing(Rationals(), m); Coefficients(R! ( (&*[1/(1-x^k)^k^2: k in [1..m]]) )); // _G. C. Greubel_, Oct 29 2018

%o (SageMath) # uses[EulerTransform from A166861]

%o b = EulerTransform(lambda n: n^2)

%o print([b(n) for n in range(32)]) # _Peter Luschny_, Nov 11 2020

%Y Euler transform of squares (A000290).

%Y Cf. A000219, A023872-A023878, A294530.

%Y Column k=2 of A144048. - _Alois P. Heinz_, Nov 02 2012

%K nonn

%O 0,3

%A _Olivier Gérard_

%E Definition corrected by _Franklin T. Adams-Watters_ and _R. J. Mathar_, Dec 04 2006