login
This site is supported by donations to The OEIS Foundation.

 

Logo


Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
A269028 a(n) = 40*a(n - 1) - a(n - 2) for n>1, a(0) = 1,  a(1) = 1. 0
1, 1, 39, 1559, 62321, 2491281, 99588919, 3981065479, 159143030241, 6361740144161, 254310462736199, 10166056769303799, 406387960309415761, 16245352355607326641, 649407706263983649879, 25960062898203738668519, 1037753108221885563090881 (list; graph; refs; listen; history; text; internal format)
OFFSET

0,3

COMMENTS

In general, the ordinary generating function for the recurrence relation b(n) = k*b(n - 1) - b(n - 2) with n>1 and b(0)=1, b(1)=1, is (1 - (k - 1)*x)/(1 - k*x +x^2). This recurrence gives the closed form b(n) = (2^( -n - 1)*((k - 2)*(k - sqrt(k^2 - 4))^n + sqrt(k^2 - 4)*(k - sqrt(k^2 - 4))^n - (k - 2)*(sqrt(k^2 - 4) + k)^n + sqrt(k^2 - 4)*(sqrt(k^2 - 4) + k)^n))/sqrt(k^2 - 4).

LINKS

Table of n, a(n) for n=0..16.

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

FORMULA

G.f.: (1 - 39*x)/(1 - 40*x + x^2).

a(n) = cosh(n*log(20 + sqrt(399))) - sqrt(19/21)*sinh(n*log(20 + sqrt(399))).

a(n) = (2^(-n - 2)*(38*(40 - 2*sqrt(399))^n + 2*sqrt(399)*(40 - 2*sqrt(399))^n - 38*(40 + 2*sqrt(399))^n + 2*sqrt(399)*(40 + 2*sqrt(399))^n))/sqrt(399).

Sum_{n>=0} 1/a(n) = 2.0262989201139499769986...

MATHEMATICA

Table[Cosh[n Log[20 + Sqrt[399]]] - Sqrt[19/21] Sinh[n Log[20 + Sqrt[399]]], {n, 0, 17}]

Table[(2^(-n - 2) (38 (40 - 2 Sqrt[399])^n + 2 Sqrt[399] (40 - 2 Sqrt[399])^n - 38 (40 + 2 Sqrt[399])^n + 2 Sqrt[399] (40 + 2 Sqrt[399])^n))/Sqrt[399], {n, 0, 17}]

LinearRecurrence[{40, -1}, {1, 1}, 17]

PROG

(MAGMA) [n le 2 select 1 else 40*Self(n-1)-Self(n-2): n in [1..20]]; // Vincenzo Librandi, Feb 19 2016

CROSSREFS

Cf. A001519, A001835, A001653, A049685, A070997, A070998, A072256, A078922, A160682, A007805, A075839, A157014, A159664, A159668, A157877, A238379, A097315.

Sequence in context: A218741 A112617 A009983 * A158768 A139191 A176073

Adjacent sequences:  A269025 A269026 A269027 * A269029 A269030 A269031

KEYWORD

nonn,easy

AUTHOR

Ilya Gutkovskiy, Feb 18 2016

STATUS

approved

Lookup | Welcome | Wiki | Register | Music | Plot 2 | Demos | Index | Browse | More | WebCam
Contribute new seq. or comment | Format | Style Sheet | Transforms | Superseeker | Recent | More pages
The OEIS Community | Maintained by The OEIS Foundation Inc.

License Agreements, Terms of Use, Privacy Policy .

Last modified February 21 12:59 EST 2018. Contains 299411 sequences. (Running on oeis4.)