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