OFFSET
1,8
FORMULA
G.f.: Sum_{k>=1} x^(5*k+1)/Product_{j=k..4*k+1} (1-x^j).
MATHEMATICA
nmax = 100; p = 1; s = 0; Do[p = Simplify[p*(1 - x^(4*k - 3))*(1 - x^(4*k - 2))*(1 - x^(4*k - 1))*(1 - x^(4*k))/(1 - x^k)]; p = Normal[p + O[x]^(nmax + 1)]; s += x^(5*k + 1)/(1 - x^k)/(1 - x^(4*k + 1))/p; , {k, 1, nmax}]; Rest[CoefficientList[Series[s, {x, 0, nmax}], x]] (* Vaclav Kotesovec, Jun 19 2025 *)
PROG
(PARI) my(N=70, x='x+O('x^N)); concat([0, 0, 0, 0, 0], Vec(sum(k=1, N, x^(5*k+1)/prod(j=k, 4*k+1, 1-x^j))))
CROSSREFS
KEYWORD
nonn
AUTHOR
Seiichi Manyama, May 17 2023
STATUS
approved
