A233698 Expansion of b(q^2) * c(q^2) / (3 * b(q)^2) in powers of q where b(), c() are cubic AGM functions.
1, 6, 25, 84, 248, 666, 1662, 3912, 8774, 18894, 39289, 79248, 155612, 298338, 559812, 1030224, 1862647, 3313494, 5807096, 10037796, 17129888, 28886052, 48170178, 79492824, 129900206, 210314976, 337545438, 537278124, 848509124, 1330069554, 2070183912
Offset: 0
Keywords
Examples
G.f. = 1 + 6*x + 25*x^2 + 84*x^3 + 248*x^4 + 666*x^5 + 1662*x^6 + 3912*x^7 + ... G.f. = q^2 + 6*q^5 + 25*q^8 + 84*q^11 + 248*q^14 + 666*q^17 + 1662*q^20 + ...
Links
- G. C. Greubel, Table of n, a(n) for n = 0..2500
- Vaclav Kotesovec, A method of finding the asymptotics of q-series based on the convolution of generating functions, arXiv:1509.08708 [math.CO], 2015-2016.
Programs
-
Mathematica
nmax=60; CoefficientList[Series[Product[((1-x^(2*k)) * (1-x^(3*k)) * (1-x^(6*k)) / (1-x^k)^3)^2,{k,1,nmax}],{x,0,nmax}],x] (* Vaclav Kotesovec, Oct 13 2015 *) eta[q_]:= q^(1/24)*QPochhammer[q]; a:= CoefficientList[Series[q^(-2/3) *(eta[q^2]*eta[q^3]*eta[q^6]/eta[q]^3)^2, {q, 0, 50}], q] (* G. C. Greubel, Aug 07 2018 *)
-
PARI
{a(n) = local(A); if( n<0, 0, A = x * O(x^n); polcoeff( (eta(x^2 + A) * eta(x^3 + A) * eta(x^6 + A) / eta(x + A)^3)^2, n))}
Formula
Expansion of (eta(q^2) * eta(q^3) * eta(q^6) / eta(q)^3)^2 in powers of q.
Euler transform of period 6 sequence [ 6, 4, 4, 4, 6, 0, ...].
a(n) ~ exp(2*Pi*sqrt(2*n/3)) / (2^(11/4) * 3^(9/4) * n^(3/4)). - Vaclav Kotesovec, Oct 13 2015
Comments