OFFSET
0,5
LINKS
Colin Barker, Table of n, a(n) for n = 0..1000
Index entries for linear recurrences with constant coefficients, signature (3,1,-5,-1,1)
FORMULA
a(n) = +3*a(n-1) + 1*a(n-2) - 5*a(n-3) - 1*a(n-4) + 1*a(n-5).
G.f.: x^3/(1-3*x-x^2+5*x^3+x^4-x^5) = x^3/((1+x)*(1-x-x^2)*(1-3*x+x^2)).
a(n) + a(n+1) = A056014(n+1). - R. J. Mathar, Jun 24 2011
a(n) = (2*F(n)^2 - F(n+4) + 3*F(n+1))/4, F(n) = A000045(n). - Gary Detlefs, Jan 05 2013
a(n) = Sum_{k=1..n-2} A122931(k). - J. M. Bergot, Apr 05 2013
a(n) = (2^(-1-n)*(-(-1)^n*2^(1+n) + sqrt(5)*(1-sqrt(5))^n + (3-sqrt(5))^n - sqrt(5)*(1+sqrt(5))^n + (3+sqrt(5))^n)) / 5. - Colin Barker, Mar 26 2017
EXAMPLE
a(7) = binomial(13,2) = 78.
MAPLE
with(combinat): seq(binomial(fibonacci(n), 2), n = 0 .. 30);
MATHEMATICA
Table[Binomial[Fibonacci[n], 2], {n, 0, 39}] (* Alonso del Arte, Apr 04 2013 *)
PROG
(PARI) a(n) = binomial(fibonacci(n), 2); \\ Michel Marcus, Sep 07 2015
(PARI) concat(vector(3), Vec(x^3 / ((1+x)*(1-x-x^2)*(1-3*x+x^2)) + O(x^40))) \\ Colin Barker, Mar 26 2017
(Python)
from sympy import binomial, fibonacci
def a(n): return binomial(fibonacci(n), 2) # Indranil Ghosh, Mar 26 2017
CROSSREFS
KEYWORD
nonn,easy
AUTHOR
Emeric Deutsch, Jun 21 2011
STATUS
approved