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