A164581 a(n) = 5*a(n - 1) + a(n - 2), with a(0)=1, a(1)=2.

1, 2, 11, 57, 296, 1537, 7981, 41442, 215191, 1117397, 5802176, 30128277, 156443561, 812346082, 4218173971, 21903215937, 113734253656, 590574484217, 3066606674741, 15923607857922, 82684645964351, 429346837679677, 2229418834362736, 11576441009493357

Offset

0,2

Links

Vincenzo Librandi, Table of n, a(n) for n = 0..1000

Index entries for linear recurrences with constant coefficients, signature (5,1).

Formula

a(n) = 5*a(n-1)+a(n-2) = A052918(n)-3*A052918(n-1).

G.f.: (1-3*x)/(1-5*x-x^2).

a(n) = A052918(n) + A015449(n). - R. J. Mathar, Jul 06 2012

a(n) = (2^(-1-n)*((5-sqrt(29))^n*(1+sqrt(29))+(-1+sqrt(29))*(5+sqrt(29))^n))/sqrt(29). - Colin Barker, Oct 13 2015

a(n) = Sum_{k=0..n-2} A168561(n-2,k)*5^k + 2 * Sum_{k=0..n-1} A168561(n-1,k)*5^k, n>0. - R. J. Mathar, Feb 14 2024

a(n) = A052918(n) -3*A052918(n-1). - R. J. Mathar, Feb 14 2024

Mathematica

LinearRecurrence[{5, 1}, {1, 2}, 40] (* or *) Rest[CoefficientList[Series [x (1 - 3 x) / (1 - 5 x - x^2), {x, 0, 40}], x]] (* Harvey P. Dale, May 02 2011 *)

Prog

(Magma) [ n le 2 select (n) else 5*Self(n-1)+Self(n-2): n in [1..25] ]; // Vincenzo Librandi, Sep 12 2013
(PARI) Vec((1-3*x)/(1-5*x-x^2) + O(x^40)) \\ Colin Barker, Oct 13 2015

Crossrefs

Sequence in context: A037570 A240888 A120464 * A054130 A037738 A037633

Adjacent sequences: A164578 A164579 A164580 * A164582 A164583 A164584

Keyword

nonn,easy

Author

Vincenzo Librandi, Aug 17 2009

Status

approved