OFFSET
1,13
LINKS
Vincenzo Librandi, Table of n, a(n) for n = 1..1000
FORMULA
G.f.: Sum_{j>=1} q^(4*j+1) * Product_{k=j+1..3*j} (1+q^k).
MATHEMATICA
Nmax=60; a=CoefficientList[Series[Sum[q^(4 j+1) Product[1+q^k, {k, j+1, 3 j}], {j, 1, Nmax}], {q, 0, Nmax}], q][[2;; ]] (* Vincenzo Librandi, Mar 07 2026 *)
PROG
(Magma) N := 60; R<q> := PowerSeriesRing(Integers(), N+1); gf := &+[ q^(4*j+1) * &*[1+q^k : k in [j+1..3*j]] : j in [1..N] ]; [Coefficient(gf, n) : n in [1..N]]; // Vincenzo Librandi, Mar 07 2026
CROSSREFS
KEYWORD
nonn
AUTHOR
Seiichi Manyama, Mar 07 2026
STATUS
approved
