OFFSET
0,1
LINKS
Vincenzo Librandi, Table of n, a(n) for n = 0..1000
Index entries for linear recurrences with constant coefficients, signature (1,8,6,-1,-1).
FORMULA
a(n) = a(n-1) + 8*a(n-2) + 6*a(n-3) - a(n-4) - a(n-5).
G.f.: (5-4*x-24*x^2-12*x^3+x^4)/((1+x)*(1-2*x-6*x^2+x^4)). - Maksym Voznyy (voznyy(AT)mail.ru), Aug 12 2009
MAPLE
seq(coeff(series((5-4*x-24*x^2-12*x^3+x^4)/((1+x)*(1-2*x-6*x^2+x^4)), x, n+1), x, n), n = 0..30); # G. C. Greubel, Oct 30 2019
MATHEMATICA
CoefficientList[Series[(5-4*x-24*x^2-12*x^3+x^4)/((1+x)*(1-2*x-6*x^2+ x^4)), {x, 0, 30}], x] (* Vincenzo Librandi, May 11 2017 *)
PROG
(Magma) I:=[5, 1, 17, 43, 181]; [n le 5 select I[n] else Self(n-1) + 8*Self(n-2) + 6*Self(n-3) - Self(n-4) - Self(n-5): n in [1..30]]; // Vincenzo Librandi, May 11 2017
(PARI) my(x='x+O('x^30)); Vec((5-4*x-24*x^2-12*x^3+x^4)/((1+x)*(1-2*x-6*x^2+x^4))) \\ G. C. Greubel, Oct 30 2019
(Sage)
def A077952_list(prec):
P.<x> = PowerSeriesRing(ZZ, prec)
return P((5-4*x-24*x^2-12*x^3+x^4)/((1+x)*(1-2*x-6*x^2+x^4))).list()
A077952_list(30) # G. C. Greubel, Oct 30 2019
(GAP) a:=[5, 1, 17, 43, 181];; for n in [6..30] do a[n]:=a[n-1]+8*a[n-2] +6*a[n-3] -a[n-4]-a[n-5]; od; a; # G. C. Greubel, Oct 30 2019
CROSSREFS
KEYWORD
easy,nonn
AUTHOR
Stephen G Penrice, Dec 21 1999
EXTENSIONS
More terms from Michael Lugo (mlugo(AT)thelabelguy.com), Dec 22 1999
STATUS
approved