OFFSET
1,3
LINKS
Vincenzo Librandi, Table of n, a(n) for n = 1..1000
Index entries for linear recurrences with constant coefficients, signature (1,0,0,1,-1).
FORMULA
a(n) = a(n-1)+a(n-4)-a(n-5) for n>5. G.f.: x^2*(8*x^3+7*x^2+8*x+1) / ((x-1)^2*(x+1)*(x^2+1)). - Colin Barker, Nov 29 2012
a(n) = (12*n+3*i^(n*(n-1))-2*(-1)^n-17)/2, where i=sqrt(-1). - Bruno Berselli, Nov 29 2012
MAPLE
A151973:=n->(12*n+3*I^(n*(n-1))-2*I^(2*n)-17)/2: seq(A151973(n), n=1..100); # Wesley Ivan Hurt, Jun 07 2016
MATHEMATICA
CoefficientList[Series[x (8 x^3 + 7 x^2 + 8 x + 1) / ((x - 1)^2 (x + 1) (x^2 + 1)), {x, 0, 40}], x] (* Vincenzo Librandi, Jun 19 2013 *)
PROG
(Magma) [n: n in [0..350] | IsZero(n*(n-1) mod 24)]; // Bruno Berselli, Nov 29 2012
(PARI) is(n)=(n^2-n)%24==0 \\ Charles R Greathouse IV, Oct 16 2015
CROSSREFS
KEYWORD
nonn,easy
AUTHOR
N. J. A. Sloane, Aug 23 2009
STATUS
approved