|
|
A108851
|
|
a(n) = 4*a(n-1) + 3*a(n-2), a(0) = 1, a(1) = 2.
|
|
8
|
|
|
1, 2, 11, 50, 233, 1082, 5027, 23354, 108497, 504050, 2341691, 10878914, 50540729, 234799658, 1090820819, 5067682250, 23543191457, 109375812578, 508132824683, 2360658736466, 10967033419913, 50950109889050
(list;
graph;
refs;
listen;
history;
text;
internal format)
|
|
|
OFFSET
|
0,2
|
|
COMMENTS
|
Binomial transform of A083098, second binomial transform of (1, 0, 7, 0, 49, 0, 243, 0, ...).
|
|
LINKS
|
|
|
FORMULA
|
a(n) = ((2 + sqrt(7))^n + (2 - sqrt(7))^n) / 2.
G.f.: (1 - 2*x) / (1 - 4*x - 3*x^2).
E.g.f.: exp(2*x)*cosh(sqrt(7)*x).
a(n+1)/a(n) converges to 2 + sqrt(7) = 4.645751311064...
G.f.: G(0)/2, where G(k) = 1 + 1/(1 - x*(7*k-4)/(x*(7*k+3) - 2/G(k+1))); (continued fraction). - Sergei N. Gladkovskii, May 27 2013
|
|
MATHEMATICA
|
LinearRecurrence[{4, 3}, {1, 2}, 30] (* Harvey P. Dale, Jan 02 2022 *)
|
|
PROG
|
(Sage) [lucas_number2(n, 4, -3)/2 for n in range(0, 22)] # Zerinvary Lajos, May 14 2009
(Magma) [Floor(((2 + Sqrt(7))^n + (2 - Sqrt(7))^n) / 2): n in [0..30]]; // Vincenzo Librandi, Jul 18 2011
|
|
CROSSREFS
|
|
|
KEYWORD
|
easy,nonn
|
|
AUTHOR
|
|
|
STATUS
|
approved
|
|
|
|