OFFSET
1,7
LINKS
Vaclav Kotesovec, Table of n, a(n) for n = 1..5000
FORMULA
G.f.: Sum_{k>=1} x^(4*k+1)/Product_{j=k..3*k+1} (1-x^j).
a(n) ~ c * A376815^sqrt(n) / sqrt(n), where c = 0.33761... - Vaclav Kotesovec, Jun 20 2025
MATHEMATICA
nmax = 100; p = 1; s = 0; Do[p = Simplify[p*(1 - x^(3*k - 2))*(1 - x^(3*k - 1))*(1 - x^(3*k))/(1 - x^k)]; p = Normal[p + O[x]^(nmax + 1)]; s += x^(4*k + 1)/(1 - x^k)/(1 - x^(3*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], Vec(sum(k=1, N, x^(4*k+1)/prod(j=k, 3*k+1, 1-x^j))))
CROSSREFS
KEYWORD
nonn
AUTHOR
Seiichi Manyama, May 17 2023
STATUS
approved
