OFFSET
0,2
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.
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
EXAMPLE
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 + ...
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 *)
PROG
(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))}
CROSSREFS
KEYWORD
nonn
AUTHOR
Michael Somos, Dec 14 2013
STATUS
approved