[go: up one dir, main page]

login
A155144
a(n) = 9*a(n-1) + 9*a(n-2), n>2; a(0)=1, a(1)=8, a(2)=80.
11
1, 8, 80, 792, 7848, 77760, 770472, 7634088, 75641040, 749476152, 7426054728, 73579777920, 729052493832, 7223690445768, 71574686456400, 709185392119512, 7026840707183208, 69624234893724480, 689859680408169192
OFFSET
0,2
FORMULA
G.f.: (1-x-x^2)/(1-9*x-9*x^2).
From G. C. Greubel, Mar 24 2021: (Start)
a(n) = (1/9)*[n=0] - 8*3^(n-2)*ChebyshevU(n, -3*i/2).
a(n) = (1/9)*[n=0] + 8*3^(n-2)*Fibonacci(n+1, 3). (End)
MAPLE
m:=30; S:=series( (1-x-x^2)/(1-9*x-9*x^2), x, m+1):
seq(coeff(S, x, j), j=0..m); # G. C. Greubel, Mar 24 2021
MATHEMATICA
Join[{1}, LinearRecurrence[{9, 9}, {8, 80}, 20]] (* or *) CoefficientList[ Series[ (1-x-x^2)/(1-9x-9x^2), {x, 0, 20}], x] (* Harvey P. Dale, Jun 19 2011 *)
PROG
(Magma) [1] cat [n le 2 select 8*(9*n-8) else 9*(Self(n-1) + Self(n-2)): n in [1..30]]; // G. C. Greubel, Mar 24 2021
(Sage) [1]+[-8*(3*i)^(n-2)*chebyshev_U(n, -3*i/2) for n in (1..30)] # G. C. Greubel, Mar 24 2021
CROSSREFS
Sequences of the form a(n) = m*(a(n-1) + a(n-2)) with a(0)=1, a(1) = m-1, a(2) = m^2 -1: A155020 (m=2), A155116 (m=3), A155117 (m=4), A155119 (m=5), A155127 (m=6), A155130 (m=7), A155132 (m=8), this sequence (m=9), A155157 (m=10).
Sequence in context: A024101 A368920 A291181 * A299871 A136949 A346178
KEYWORD
nonn
AUTHOR
Philippe Deléham, Jan 21 2009
STATUS
approved