|
| |
|
|
A090017
|
|
a(n) = 4a(n-1) + 2a(n-2), a(0) = 0, a(1) = 1.
|
|
16
| |
|
|
0, 1, 4, 18, 80, 356, 1584, 7048, 31360, 139536, 620864, 2762528, 12291840, 54692416, 243353344, 1082798208, 4817899520, 21437194496, 95384577024, 424412697088, 1888419942400, 8402505163776, 37386860539904, 166352452487168
(list; graph; refs; listen; history; internal format)
|
|
|
|
OFFSET
| 0,3
|
|
|
COMMENTS
| Starting with "1" = INVERT transform of A007482: (1, 3, 11, 39, 139,...). [From Gary W. Adamson (qntmpkt(AT)yahoo.com), Aug 06 2010]
|
|
|
LINKS
| Vincenzo Librandi, Table of n, a(n) for n = 0..1000
Index entries for sequences related to Chebyshev polynomials.
Index to sequences with linear recurrences with constant coefficients, signature (4,2).
|
|
|
FORMULA
| G.f.: x/(1-4*x-2*x^2).
a(n) = (-i*sqrt(2))^(n-1) U(n-1, i*sqrt(2)) where U is the Chebyshev polynomial of the second kind and i^2 = -1.
a(n) = ((2+sqrt(6))^n - (2-sqrt(6))^n)/(2 sqrt(6)). [From Al Hakanson (hawkuu(AT)gmail.com), Jan 05 2009, Jan 07 2009]
a(n+1) = Sum_{k, 0<=k<=n} A099089(n,k)*2^k. - Philippe DELEHAM, Nov 21 2011
|
|
|
MATHEMATICA
| a[n_Integer] := (-I Sqrt[2])^(n - 1) ChebyshevU[ n - 1, I Sqrt[2] ]
a[n_]:=(MatrixPower[{{1, 5}, {1, 3}}, n].{{1}, {1}})[[2, 1]]; Table[Abs[a[n]], {n, -1, 40}] [From Vladimir Orlovsky (4vladimir(AT)gmail.com), Feb 19 2010]
|
|
|
PROG
| (Sage) [lucas_number1(n, 4, -2) for n in xrange(0, 23)]# [From Zerinvary Lajos (zerinvarylajos(AT)yahoo.com), Apr 23 2009]
(MAGMA) I:=[0, 1]; [n le 2 select I[n] else 4*Self(n-1)+2*Self(n-2): n in [1..30]]; // Vincenzo Librandi, Oct 12 2011
(PARI) Vec(x/(1-4*x-2*x^2)+O(x^99)) \\ Charles R Greathouse IV, Oct 12 2011
|
|
|
CROSSREFS
| Cf. A007070, A084059, A007482.
Sequence in context: A177755 A037965 A045902 * A104631 A106391 A063881
Adjacent sequences: A090014 A090015 A090016 * A090018 A090019 A090020
|
|
|
KEYWORD
| easy,nonn
|
|
|
AUTHOR
| Paul Barry (pbarry(AT)wit.ie), Nov 19 2003
|
|
|
EXTENSIONS
| Edited by Stuart Clary, Oct 25 2009
|
| |
|
|
