login
The OEIS is supported by the many generous donors to the OEIS Foundation.

 

Logo
Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
A163306 a(n) = 12*a(n-1) - 31*a(n-2) for n > 1; a(0) = 1, a(1) = 7. 2

%I #8 Sep 08 2022 08:45:46

%S 1,7,53,419,3385,27631,226637,1863083,15331249,126219415,1039364261,

%T 8559569267,70494539113,580587822079,4781723152445,39382455344891,

%U 324356046412897,2671416441263143,22001959856357909,181209608597137475

%N a(n) = 12*a(n-1) - 31*a(n-2) for n > 1; a(0) = 1, a(1) = 7.

%C Binomial transform of A090041. Inverse binomial transform of A163307.

%H G. C. Greubel, <a href="/A163306/b163306.txt">Table of n, a(n) for n = 0..1000</a>

%H <a href="/index/Rec#order_02">Index entries for linear recurrences with constant coefficients</a>, signature (12,-31).

%F a(n) = ((5+sqrt(5))*(6+sqrt(5))^n + (5-sqrt(5))*(6-sqrt(5))^n)/10.

%F G.f.: (1-5*x)/(1-12*x+31*x^2).

%F E.g.f.: (1/5)*exp(6*x)*(5*cosh(sqrt(5)*x) + sqrt(5)*sinh(sqrt(5)*x)). - _G. C. Greubel_, Dec 18 2016

%t LinearRecurrence[{12,-31}, {1,7}, 50] (* _G. C. Greubel_, Dec 18 2016 *)

%o (Magma) [ n le 2 select 6*n-5 else 12*Self(n-1)-31*Self(n-2): n in [1..20] ];

%o (PARI) Vec((1-5*x)/(1-12*x+31*x^2) + O(x^50)) \\ _G. C. Greubel_, Dec 18 2016

%Y Cf. A090041, A163307.

%K nonn,easy

%O 0,2

%A _Klaus Brockhaus_, Jul 24 2009

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

License Agreements, Terms of Use, Privacy Policy. .

Last modified April 23 07:34 EDT 2024. Contains 371905 sequences. (Running on oeis4.)