OFFSET
0,2
COMMENTS
Self-convolution of A035382.
In general, if a > 0, b > 0, GCD(a,b) = 1 and g.f. = Product_{k>=0} 1/(1 - x^(a*k+b))^2, then a(n) ~ Gamma(b/a)^2 * a^(b/a - 3/4) * exp(2*Pi*sqrt(n/(3*a))) * Pi^(2*b/a - 2) / (4 * 3^(b/a - 1/4) * n^(b/a + 1/4)).
FORMULA
a(n) ~ exp(2*Pi*sqrt(n)/3) * Gamma(1/3)^2 / (4 * sqrt(3) * Pi^(4/3) * n^(7/12)).
MATHEMATICA
nmax = 60; CoefficientList[Series[Product[1/(1 - x^(3*k+1))^2, {k, 0, nmax}], {x, 0, nmax}], x]
CROSSREFS
KEYWORD
nonn
AUTHOR
Vaclav Kotesovec, Aug 26 2015
STATUS
approved