login
A246645
Expansion of 1/(1 - 22*x + 81*x^2), used in A246643.
5
1, 22, 403, 7084, 123205, 2136706, 37027927, 641541208, 11114644489, 192557340910, 3335975296411, 57794311907332, 1001260862952013, 17346399720450394, 300518663950795615, 5206352229561021616, 90197737270328030737, 1562635689352773925318, 27071968446864455867299
OFFSET
0,2
COMMENTS
This sequence is used in the formula for the curvature in a touching circle problem considered in A247512 and A246643.
FORMULA
O.g.f.: 1/(1 - 22*x + 81*x^2).
a(n) = 9^n*S(n, 22/9) with Chebyshev's S-polynomials (see A049310).
a(n) = 22*a(n-1) - 81*a(n-2), n >= 1, a(-1) = 0 and a(0) = 1.
a(n) = 9^n*(ap^(n+1) - am^(n+1))/(ap - am)), n >= 1, with ap:= (11 + 2*sqrt(10))/9 and am = 1/ap = (11 - 2*sqrt(10))/9 (Binet - de Moivre formula). a(0) = 1 (via L'Hopital's rule).
a(n) = 9^(n+1)*sinh(2*(n + 1)*arccsch(3))/(2*sqrt(10)). - Federico Provvedi, Feb 02 2021
MATHEMATICA
CoefficientList[Series[1/(1 - 22*x + 81*x^2), {x, 0, 50}], x] (* or *) LinearRecurrence[{22, -81}, {1, 22}, 50] (* G. C. Greubel, Dec 20 2017 *)
PROG
(PARI) Vec(1/(1 - 22*x + 81*x^2) + O(x^40)) \\ Michel Marcus, Sep 30 2014
(Magma) I:=[1, 22]; [n le 2 select I[n] else 22*Self(n-1) - 81*Self(n-2): n in [1..30]]; // G. C. Greubel, Dec 20 2017
CROSSREFS
KEYWORD
nonn,easy
AUTHOR
Wolfdieter Lang, Sep 30 2014
STATUS
approved