OFFSET
0,4
LINKS
Seiichi Manyama, Table of n, a(n) for n = 0..1000
Index entries for linear recurrences with constant coefficients, signature (0,0,9,0,0,-36,0,3,84,0,-18,-126,0,45,126,-3,-60,-84,9,45,36,-9,-18,-8,3,3,1).
FORMULA
G.f.: ((1-x^3)^6 + 7*x^8*(1-x^3)^3 + x^16) / ((1-x^3)^3 - x^8)^3.
a(n) = 9*a(n-3) - 36*a(n-6) + 3*a(n-8) + 84*a(n-9) - 18*a(n-11) - 126*a(n-12) + 45*a(n-14) + 126*a(n-15) - 3*a(n-16) - 60*a(n-17) - 84*a(n-18) + 9*a(n-19) + 45*a(n-20) + 36*a(n-21) - 9*a(n-22) - 18*a(n-23) - 8*a(n-24) + 3*a(n-25) + 3*a(n-26) + a(n-27).
MATHEMATICA
CoefficientList[Series[((1-x^3)^6+7*x^8*(1-x^3)^3+x^16)/((1-x^3)^3-x^8)^3, {x, 0, 60}], x] (* Vincenzo Librandi, Jan 16 2026 *)
PROG
(PARI) my(N=50, x='x+O('x^N)); Vec(((1-x^3)^6+7*x^8*(1-x^3)^3+x^16)/((1-x^3)^3-x^8)^3)
(Magma) m:=60; R<x>:=PowerSeriesRing(Integers(), m); Coefficients(R! ((1-x^3)^6 + 7*x^8*(1-x^3)^3 + x^16) / ((1-x^3)^3 - x^8)^3); // Vincenzo Librandi, Jan 16 2026
CROSSREFS
KEYWORD
nonn,easy
AUTHOR
Seiichi Manyama, Jan 06 2026
STATUS
approved
