OFFSET
0,3
COMMENTS
LINKS
G. C. Greubel, Table of n, a(n) for n = 0..850
Index entries for linear recurrences with constant coefficients, signature (9,72).
FORMULA
a(n+1) = Sum_{k=0..n} A154929(n,k)*8^(n-k).
G.f.: (1 - 8*x - 64*x^2)/(1 - 9*x - 72*x^2). - G. C. Greubel, Apr 20 2021
MAPLE
a[0] := 1: a[1] := 1: a[2] := 17: for n from 3 to 25 do a[n] := 9*a[n-1]+72*a[n-2] end do: seq(a[n], n = 0 .. 17); # Emeric Deutsch, Jan 21 2009
MATHEMATICA
LinearRecurrence[{9, 72}, {1, 1, 17}, 20] (* Harvey P. Dale, Apr 26 2016 *)
PROG
(Magma) I:=[1, 17]; [1] cat [n le 2 select I[n] else 9*(Self(n-1) +8*Self(n-2)): n in [1..30]]; // G. C. Greubel, Apr 20 2021
(Sage)
def A155001_list(prec):
P.<x> = PowerSeriesRing(ZZ, prec)
return P( (1-8*x-64*x^2)/(1-9*x-72*x^2) ).list()
A155001_list(30) # G. C. Greubel, Apr 20 2021
CROSSREFS
KEYWORD
nonn
AUTHOR
Philippe Deléham, Jan 18 2009
EXTENSIONS
Corrected by Philippe Deléham, Jan 21 2009
Corrected and extended by Emeric Deutsch and R. J. Mathar, Jan 21 2009
STATUS
approved