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
Links
- Vincenzo Librandi, Table of n, a(n) for n = 0..1000
- Index entries for linear recurrences with constant coefficients, signature (5,1).
Programs
-
Magma
[ n le 2 select (n) else 5*Self(n-1)+Self(n-2): n in [1..25] ]; // Vincenzo Librandi, Sep 12 2013
-
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 *)
-
PARI
Vec((1-3*x)/(1-5*x-x^2) + O(x^40)) \\ Colin Barker, Oct 13 2015
Formula
G.f.: (1-3*x)/(1-5*x-x^2).
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
From Peter Bala, Jul 08 2025: (Start)
The following series telescope:
Sum_{n >= 1} 1/(a(n) - 7*(-1)^n/a(n)) = 3/10, since 1/(a(n) - 7*(-1)^n/a(n)) = b(n) - b(n+1), where b(n) = (1/5) * (a(n) + a(n-1)) / (a(n)*a(n-1)).
Sum_{n >= 1} (-1)^(n+1)/(a(n) - 7*(-1)^n/a(n)) = 1/10, since 1/(a(n) - 7*(-1)^n/a(n)) = c(n) + c(n+1), where c(n) = (1/5) * (a(n) - a(n-1)) / (a(n)*a(n-1)). (End)