OFFSET
1,8
LINKS
Vaclav Kotesovec, Table of n, a(n) for n = 1..10000
FORMULA
G.f.: Sum_{k>=1} x^(4*k)/Product_{j>=4*k} (1-x^j).
a(n) ~ Pi^3 * exp(Pi*sqrt(2*n/3)) / (3 * 2^(5/2) * n^(5/2)) * (1 - (5*sqrt(6)/Pi + 169*Pi*sqrt(6)/144)/sqrt(n)). - Vaclav Kotesovec, May 21 2023
MATHEMATICA
nmax = 60; Rest[CoefficientList[Series[Sum[x^(4*k)/QPochhammer[x^(4*k), x], {k, 1, nmax/4}], {x, 0, nmax}], x]] (* Vaclav Kotesovec, May 20 2023 *)
PROG
(PARI) my(N=70, x='x+O('x^N)); concat([0, 0, 0], Vec(sum(k=1, N, x^(4*k)/prod(j=4*k, N, 1-x^j))))
CROSSREFS
KEYWORD
nonn
AUTHOR
Seiichi Manyama, May 19 2023
STATUS
approved