OFFSET
8,1
LINKS
Andrew Howroyd, Table of n, a(n) for n = 8..1000
Index entries for linear recurrences with constant coefficients, signature (3,-1,-4,3,-1,3,0,-3,1,-3,4,1,-3,1).
FORMULA
G.f.: x^8*(2+5*x+9*x^2+8*x^3+5*x^4+x^6)/((1-x)^7*(1+x)^3*(1+x^2)*(1+x+x^2)) - Colin Barker, Nov 07 2012
a(n) = -[n=0] + (60*n^6 -1236*n^5 +9450*n^4 -33520*n^3 +59940*n^2 -66294*n +48065)/69120 +(-1)^n*(4*n^2 -38*n +89)/512 +(3/32)*(-1)^floor((n+1)/2) + ChebyshevU(n, -1/2)/27. - G. C. Greubel, Nov 09 2023
MATHEMATICA
Drop[CoefficientList[Series[x^8*(2+5*x+9*x^2+8*x^3+5*x^4+x^6)/((1-x)^3*(1-x^2)^2*(1-x^3)*(1-x^4)), {x, 0, 50}], x], 8] (* G. C. Greubel, Nov 09 2023 *)
PROG
(PARI) Vec((2 + 5*x + 9*x^2 + 8*x^3 + 5*x^4 + x^6)/((1 - x)^7*(1 + x)^3*(1 + x^2)*(1 + x + x^2)) + O(x^40)) \\ Andrew Howroyd, Aug 28 2018
(Magma) R<x>:=PowerSeriesRing(Integers(), 50); Coefficients(R!( x^8*(2+5*x+9*x^2+8*x^3+5*x^4+x^6)/((1-x)^3*(1-x^2)^2*(1-x^3)*(1-x^4)) )); // G. C. Greubel, Nov 09 2023
(SageMath)
def A055329_list(prec):
P.<x> = PowerSeriesRing(ZZ, prec)
return P( x^8*(2+5*x+9*x^2+8*x^3+5*x^4+x^6)/((1-x)^3*(1-x^2)^2*(1-x^3)*(1-x^4)) ).list()
a=A055329_list(50); a[8:] # G. C. Greubel, Nov 09 2023
CROSSREFS
KEYWORD
nonn,easy
AUTHOR
Christian G. Bower, May 12 2000
STATUS
approved