The OEIS mourns the passing of Jim Simons and is grateful to the Simons Foundation for its support of research in many branches of science, including the OEIS.
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!)
A164598 a(n) = 12*a(n-1) - 34*a(n-2), for n > 1, with a(0) = 1, a(1) = 14. 8
1, 14, 134, 1132, 9028, 69848, 531224, 3999856, 29936656, 223244768, 1661090912, 12342768832, 91636134976, 679979479424, 5044125163904, 37410199666432, 277422140424448, 2057118896434688, 15253073982785024, 113094845314640896 (list; graph; refs; listen; history; text; internal format)
OFFSET
0,2
COMMENTS
Binomial transform of A164301. Sixth binomial transform of A164587. Inverse binomial transform of A164599.
This sequence is part of a class of sequences defined by the recurrence a(n,m) = 2*(m+1)*a(n-1,m) - ((m+1)^2 -2)*a(n-2,m) with a(0) = 1 and a(1) = m+9. The generating function is Sum_{n>=0} a(n,m)*x^n = (1 - (m-7)*x)/(1 - 2*(m+1)*x + ((m+1)^2 -2)*x^2) and have a series expansion in terms of Pell-Lucas numbers defined by a(n, m) = (1/2)*Sum_{k=0..n} binomial(n,k)*m^(n-k)*(5*Q(k) + 4*Q(k-1)). - G. C. Greubel, Mar 11 2021
LINKS
G. C. Greubel, Table of n, a(n) for n = 0..1000 (terms 0..100 from Vincenzo Librandi)
FORMULA
a(n) = ((1+4*sqrt(2))*(6+sqrt(2))^n + (1-4*sqrt(2))*(6-sqrt(2))^n)/2.
G.f.: (1+2*x)/(1-12*x+34*x^2).
E.g.f.: exp(6*x)*(cosh(sqrt(2)*x) + 4*sqrt(2)*sinh(sqrt(2)*x)). - G. C. Greubel, Aug 11 2017
From G. C. Greubel, Mar 11 2021: (Start)
a(n) = A147957(n) + 8*A081183(n).
a(n) = (1/2)*Sum_{k=0..n} binomial(n,k)*5^(n-k)*(5*Q(k) + 4*Q(k-1)), where Q(n) = Pell-Lucas(n) = A002203(n). (End)
MAPLE
m:=30; S:=series( (1+2*x)/(1-12*x+34*x^2), x, m+1):
seq(coeff(S, x, j), j=0..m); # G. C. Greubel, Mar 11 2021
MATHEMATICA
LinearRecurrence[{12, -34}, {1, 14}, 30] (* G. C. Greubel, Aug 11 2017 *)
PROG
(Magma) [ n le 2 select 13*n-12 else 12*Self(n-1)-34*Self(n-2): n in [1..30] ];
(PARI) my(x='x+O('x^30)); Vec((1+2*x)/(1-12*x+34*x^2)) \\ G. C. Greubel, Aug 11 2017
(Sage) [( (1+2*x)/(1-12*x+34*x^2) ).series(x, n+1).list()[n] for n in (0..30)] # G. C. Greubel, Mar 11 2021
CROSSREFS
Sequences in the class a(n, m): A164298 (m=1), A164299 (m=2), A164300 (m=3), A164301 (m=4), this sequence (m=5), A164599 (m=6), A081185 (m=7), A164600 (m=8).
Sequence in context: A021079 A174563 A233467 * A073554 A272652 A016801
KEYWORD
nonn,easy
AUTHOR
Klaus Brockhaus, Aug 17 2009
STATUS
approved

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 May 23 01:37 EDT 2024. Contains 372758 sequences. (Running on oeis4.)