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