OFFSET
0,2
COMMENTS
The a(n+1) are the numerators of A178381(4*n+3)/A178381(4*n+2). For the denominators see A179133(n). - Johannes W. Meijer, Jul 01 2010
LINKS
Vincenzo Librandi, Table of n, a(n) for n = 0..1000
Index entries for linear recurrences with constant coefficients, signature (0,0,18,0,0,-1).
FORMULA
a(n) = Lucas(2*n)*(Fibonacci(n) mod 2 + 1)/2, Lucas(n)=A000032, Fibonacci(n)=A000045. - Gary Detlefs, Jan 19 2001
From Colin Barker, Jun 27 2013: (Start)
a(n) = 18*a(n-3) - a(n-6).
G.f: -(3*x^5 + 7*x^4 + 9*x^3 - 7*x^2 - 3*x - 1) / ((x^2 - 3*x + 1)*(x^4 + 3*x^3 + 8*x^2 + 3*x + 1)). (End)
With L(n) the Lucas number A000032, a(n) = L(2*n)/2 or L(2*n) according as n is, or is not, divisible by 3. - David Callan, Jul 17 2019
MAPLE
with(combinat): nmax:=25; for n from 0 to nmax do a(n):= (fibonacci(2*n-1)+fibonacci(2*n+1))*(5/6-cos(2*Pi*n/3)/3) od: seq(a(n), n=0..nmax); # Johannes W. Meijer, Jul 01 2010
MATHEMATICA
LinearRecurrence[{0, 0, 18, 0, 0, -1}, {1, 3, 7, 9, 47, 123}, 40] (* Vincenzo Librandi, Jul 17 2019 *)
PROG
(Magma) I:=[1, 3, 7, 9, 47, 123]; [n le 6 select I[n] else 18*Self(n-3)-Self(n-6): n in [1..30]]; // Vincenzo Librandi, Jul 17 2019
CROSSREFS
KEYWORD
easy,nonn
AUTHOR
Paul Barry, Feb 13 2007
STATUS
approved