OFFSET
1,20
COMMENTS
In general, if m>=1 and g.f. = x^m * Sum_{k>=0} x^(2*m*k) / Product_{j=1..2*k} (1-x^j), then a(n) ~ Pi^(m-1) * (m-1)! * exp(Pi*sqrt(2*n/3)) / (2^((m+5)/2) * 3^(m/2) * n^((m+1)/2)). - Vaclav Kotesovec, Jun 20 2025
LINKS
Vaclav Kotesovec, Table of n, a(n) for n = 1..5000
FORMULA
G.f.: x^6 * Sum_{k>=0} x^(12*k)/Product_{j=1..2*k} (1-x^j). - Seiichi Manyama, May 15 2023
a(n) ~ 5 * Pi^5 * exp(Pi*sqrt(2*n/3)) / (9 * 2^(5/2) * n^(7/2)). - Vaclav Kotesovec, Jun 20 2025
MATHEMATICA
nmax = 100; p = 1; s = 1; Do[p = Expand[p*(1 - x^(2*k))*(1 - x^(2*k - 1))]; p = Take[p, Min[nmax + 1, Exponent[p, x] + 1, Length[p]]]; s += x^(12*k)/p; , {k, 1, nmax}]; Join[{0, 0, 0, 0, 0}, CoefficientList[Series[s, {x, 0, nmax}], x]] (* Vaclav Kotesovec, Jun 20 2025 *)
CROSSREFS
KEYWORD
nonn
AUTHOR
EXTENSIONS
More terms from Vladeta Jovovic, Aug 01 2009
STATUS
approved
