OFFSET
1,2
LINKS
Index entries for linear recurrences with constant coefficients, signature (1,398,-398,-1,1).
FORMULA
a(1)=0, a(2)=39, a(3)=398*a(1)+198-a(2), a(4)=398*a(2)+198-a(1), a(n) = 398*a(n-2) + 198 - a(n-4) for n>4.
From Bruno Berselli, Feb 20 2013: (Start)
G.f.: 3*x*(13+40*x+13*x^2)/((1-x)*(1-20*x+x^2)*(1+20*x+x^2)).
a(n) = a(-n+1) = -1/2+((11+2*t*(-1)^n)*(10-3*t)^(2*floor(n/2))+(11-2*t*(-1)^n)*(10+3*t)^(2*floor(n/2)))/44, where t=sqrt(11). (End)
MATHEMATICA
LinearRecurrence[{1, 398, -398, -1, 1}, {0, 39, 159, 15720, 63480}, 18] (* Bruno Berselli, Feb 20 2013 *)
PROG
(Magma) m:=17; R<x>:=PowerSeriesRing(Integers(), m); [0] cat Coefficients(R!(3*(13+40*x+13*x^2)/((1-x)*(1-20*x+x^2)*(1+20*x+x^2)))); // Bruno Berselli, Feb 20 2013
(Maxima) makelist(expand(-1/2+((11+2*sqrt(11)*(-1)^n)*(10-3*sqrt(11))^(2*floor(n/2))+(11-2*sqrt(11)*(-1)^n)*(10+3*sqrt(11))^(2*floor(n/2)))/44), n, 1, 18); /* Bruno Berselli, Feb 20 2013 */
CROSSREFS
KEYWORD
nonn,easy
AUTHOR
Pierre CAMI, Apr 22 2005
EXTENSIONS
More terms from Bruno Berselli, Feb 20 2013
STATUS
approved