OFFSET
0,5
LINKS
G. Nebe, E. M. Rains and N. J. A. Sloane, Self-Dual Codes and Invariant Theory, Springer, Berlin, 2006.
Index entries for linear recurrences with constant coefficients, signature (1,1,-1,1,-1,-1,1,1,-1,-1,1,-1,1,1,-1).
FORMULA
G.f.: ( 1+x^4+x^8-x^2-x^6 ) / ( (1+x^4) * (1+x^2)^2 * (1+x)^3 * (1-x)^4 ). - R. J. Mathar, Dec 18 2014
a(n) = (2*n^3 + 21*n^2 + 175*n + 441 + 3*(n^2 + 7*n + 29)*(-1)^n + 30*(2*n + 7)*(-1)^((2*n - 1 + (-1)^n)/4) + 30*(-1)^((6*n - 1 + (-1)^n)/4) + 48*((-1)^((2*n - 3 + (-1)^n + 2*(-1)^((2*n - 1 + (-1)^n)/4))/8) - (-1)^((6*n - 5 + 3*(-1)^n + 2*(-1)^((2*n - 1 + (-1)^n)/4))/8)))/768. - Luce ETIENNE, Mar 31 2015
Shorter version of above:
a(n) = (2*n^3 + 21*n^2 + 175*n + 441 + 3*(n^2 + 7*n + 29)*(-1)^n + 30*(2*n + 7)*sign(1-n%4+n%2) + 30*sign(2-(n+1)%4-n%2) + 48*((5*n-n^2)%4)*sign(5-n%8))/768, where sign(x) = x/abs(x), and a%b = a (mod b). - Derek Orr, Apr 05 2015
MATHEMATICA
CoefficientList[Series[(1 + x^4 + x^8 - x^2 - x^6) / ((x^4 + 1) (1 + x^2)^2 (1 + x)^3 (x-1)^4), {x, 0, 70}], x] (* Vincenzo Librandi, Apr 05 2015 *)
PROG
(PARI) Vec((1+x^10)/((1-x)*(1-x^4)^2*(1-x^8)) + O(x^80)) \\ Michel Marcus, Apr 05 2015
(PARI)
a(n) = (2*n^3 + 21*n^2 + 175*n + 441 + 3*(n^2 + 7*n + 29)*(-1)^n + 30*(2*n + 7)*sign(1-n%4+n%2) + 30*sign(2-(n+1)%4-n%2) + 48*((5*n-n^2)%4)*sign(5-n%8))/768
vector(100, n, a(n-1)) \\ Derek Orr, Apr 05 2015
(Magma) m:=70; R<x>:=PowerSeriesRing(Integers(), m); Coefficients(R!((1+x^4+x^8-x^2-x^6)/((1+x^4)*(1+x^2)^2*(1+x)^3*(1-x)^4))); // Bruno Berselli, Apr 07 2015
CROSSREFS
KEYWORD
nonn,easy
AUTHOR
N. J. A. Sloane, Apr 08 2004
STATUS
approved
