OFFSET
0,4
LINKS
Vincenzo Librandi, Table of n, a(n) for n = 0..1000
M. Baake, J. Hermisson, P. Pleasants, The torus parametrization of quasiperiodic LI-classes, J. Phys. A 30 (1997), no. 9, 3029-3056. See (31).
Index entries for linear recurrences with constant coefficients, signature (1,6,-1,-10,-1,6,1,-1).
FORMULA
G.f.: (1-x-x^2)^2*(1+x-x^2)^2/((1-3*x+x^2)*(1-x)^2*(1+x)^4).
From Colin Barker, Mar 13 2019: (Start)
a(n) = 2^(-4-n)*(-768*sqrt(5)*((3-sqrt(5))^n - (3+sqrt(5))^n) - 5*(517*(-2)^n + 75*2^n)*n + 25*(-1)^n*2^(1+n)*n^3) / 375 for n>1.
a(n) = a(n-1) + 6*a(n-2) - a(n-3) - 10*a(n-4) - a(n-5) + 6*a(n-6) + a(n-7) - a(n-8) for n>8.
(End)
MATHEMATICA
CoefficientList[Series[(1 - x - x^2)^2 (1 + x - x^2)^2 / ((1 - 3 x + x^2) (1 - x)^2 (1 + x)^4), {x, 0, 33}], x] (* Vincenzo Librandi, Mar 13 2019 *)
LinearRecurrence[{1, 6, -1, -10, -1, 6, 1, -1}, {1, 1, 1, 6, 12, 36, 91, 241, 632}, 30] (* Harvey P. Dale, Sep 06 2019 *)
PROG
(Magma) m:=30; R<x>:=PowerSeriesRing(Integers(), m); Coefficients(R!(1-x-x^2)^2*(1+x-x^2)^2/((1-3*x+x^2)*(1-x)^2*(1+x)^4)); // Vincenzo Librandi, Mar 13 2019
(PARI) Vec((1 + x - x^2)^2*(1 - x - x^2)^2 / ((1 - x)^2*(1 + x)^4*(1 - 3*x + x^2)) + O(x^30)) \\ Colin Barker, Mar 13 2019
CROSSREFS
KEYWORD
nonn,easy
AUTHOR
N. J. A. Sloane, Mar 12 2019
STATUS
approved