OFFSET
0,2
COMMENTS
(x,y)=(a(n),a(n+1)) are solutions of (x+y)^2/(1+xy)=9, the other solutions are in A033890. - Floor van Lamoen, Dec 10 2001
Sequence A033888 provides half of the solutions to the equation 5*x^2 + 4 is a square. The other half are found in A033890. Lim_{n->infinity} a(n)/a(n-1) = phi^4 = (7+3*sqrt(5))/2. - Gregory V. Richardson, Oct 13 2002
Fibonacci numbers divisible by 3. - Reinhard Zumkeller, Aug 20 2011
LINKS
Nathaniel Johnston, Table of n, a(n) for n = 0..300
Piero Filipponi and Marco Bucci, On the Integrity of Certain Fibonacci Sums, The Fibonacci Quarterly, Vol. 32, No. 3 (1994), pp. 245-252.
Tanya Khovanova, Recursive Sequences
Index entries for linear recurrences with constant coefficients, signature (7,-1).
FORMULA
a(n) = 7*a(n-1) - a(n-2).
a(n) = ((7+3*sqrt(5))^(n-1) - (7-3*sqrt(5))^(n-1)) / ((2^(n-1))*sqrt(5)). - Gregory V. Richardson, Oct 13 2002
a(n) = Sum_{k=0..n} F(3n-k)*binomial(n, k). - Benoit Cloitre, Jun 07 2004
a(n) = Lucas(2n) * Lucas(n) * Fibonacci(n). - Ralf Stephan, Sep 25 2004
G.f.: 3*x/(1-7*x+x^2). - Philippe Deléham, Nov 17 2008
a(n) = 3*A004187(n). - R. J. Mathar, Sep 03 2010
a(n) = Fibonacci[(4*n + 2)] modulo Fibonacci[(4*n + 1)]. - Artur Jasinski, Nov 15 2011 (corrected by Iain Fox, Dec 18 2017)
E.g.f.: 2*exp(7*x/2)*sinh(3*sqrt(5)*x/2)/sqrt(5). - Stefano Spezia, Feb 07 2021
a(n) = Sum_{k>=0} Fibonacci(2*n*k)/Lucas(2*n)^k (Filipponi and Bucci, 1994). - Amiram Eldar, Jan 17 2022
EXAMPLE
G.f. = 3*x + 21*x^2 + 144*x^3 + 987*x^4 + 6765*x^5 + 46368*x^6 + ...
MAPLE
MATHEMATICA
Table[Fibonacci[4*n], {n, 0, 14}] (* Vladimir Joseph Stephan Orlovsky, Jul 21 2008 *)
Table[Mod[Fibonacci[(4 n + 2)] , Fibonacci[(4 n + 1)]], {n, 1, 10}] (* Artur Jasinski, Nov 15 2011 (corrected by Iain Fox, Dec 18 2017) *)
PROG
(MuPAD) numlib::fibonacci(n*4) $ n = 0..30; // Zerinvary Lajos, May 08 2008
(Sage) [lucas_number1(n, 3, 1)*lucas_number2(n, 3, 1) for n in range(0, 21)] # Zerinvary Lajos, Jun 28 2008
(Sage) [fibonacci(4*n) for n in range(0, 20)] # Zerinvary Lajos, May 15 2009
(Magma) [ Fibonacci(4*n): n in [0..100]]; // Vincenzo Librandi, Apr 15 2011
(PARI) a(n)=fibonacci(4*n) \\ Charles R Greathouse IV, Feb 03 2014
(PARI) first(n) = Vec(3*x/(1 - 7*x + x^2) + O(x^n), -n) \\ Iain Fox, Dec 18 2017
(PARI) a(n) = fibonacci(4*n + 2) % fibonacci(4*n + 1) \\ Iain Fox, Dec 18 2017
CROSSREFS
KEYWORD
nonn,easy
AUTHOR
STATUS
approved