OFFSET
0,4
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,0,2,-2,0,-1,1).
FORMULA
G.f.: (1+4*x^3+x^6)/((1-x)*(1-x^3)^2).
a(n) = a(n-1)+2*a(n-3)-2*a(n-4)-a(n-6)+a(n-7), n>7. - Wesley Ivan Hurt, Jun 23 2015
A003215 with each term repeated three times: a(n) = A003215(floor(n/3)). - Robert Israel, Jul 14 2015
MAPLE
f:= gfun:-rectoproc({q(n+3)-3*q(n+2)+3*q(n+1)-q(n), q(0) = 1, q(1) = 7, q(2) = 19}, q(n), remember):
seq(f(i)$3, i=0..30); # Robert Israel, Jul 14 2015
MATHEMATICA
CoefficientList[Series[(1 + 4*x^3 + x^6)/((1 - x)*(1 - x^3)^2), {x, 0, 50}], x] (* Wesley Ivan Hurt, Jun 23 2015 *)
LinearRecurrence[{1, 0, 2, -2, 0, -1, 1}, {1, 1, 1, 7, 7, 7, 19}, 60] (* Vincenzo Librandi, Jul 13 2015 *)
With[{c=LinearRecurrence[{3, -3, 1}, {1, 7, 19}, 20]}, {c, c, c}]//Flatten//Sort (* Harvey P. Dale, Aug 03 2019 *)
PROG
(Magma) I:=[1, 1, 1, 7, 7, 7, 19]; [n le 7 select I[n] else Self(n-1)+2*Self(n-3)-2*Self(n-4)-Self(n-6)+Self(n-7): n in [1..70]]; // Vincenzo Librandi, Jul 13 2015
CROSSREFS
KEYWORD
nonn,easy
AUTHOR
N. J. A. Sloane, Mar 29 2004
STATUS
approved