OFFSET
0,2
LINKS
Vincenzo Librandi, Table of n, a(n) for n = 0..1000
C. Krattenthaler, Advanced determinant calculus Séminaire Lotharingien de Combinatoire, B42q (1999), 67 pp, (see p. 54).
Index entries for linear recurrences with constant coefficients, signature (7,-21,35,-35,21,-7,1).
FORMULA
G.f.: (1 + 6*x + 21*x^2 + 56*x^3 + 126*x^4 + 252*x^5 + 462*x^6) / (1-x)^7.
a(n) = 7*a(n-1) - 21*a(n-2) + 35*a(n-3) - 35*a(n-4) + 21*a(n-5) - 7*a(n-6) + a(n-7).
a(n) = (120 - 8042*n + 20581*n^2 - 17380*n^3 + 7645*n^4 -1518*n^5 + 154*n^6)/120.
MATHEMATICA
Table[(120 - 8042 n + 20581 n^2 - 17380 n^3 + 7645 n^4 - 1518 n^5 + 154 n^6)/120, {n, 0, 40}] (* or *) CoefficientList[Series[(1 + 6 x + 21 x^2 + 56 x^3 + 126 x^4 + 252 x^5 + 462 x^6)/(1 - x)^7, {x, 0, 40}], x]
PROG
(Magma) [(120-8042*n+20581*n^2-17380*n^3+7645*n^4-1518*n^5+ 154*n^6)/120: n in [0..40]];
(Magma) I:=[1, 13, 91, 455, 1820, 6188, 18564]; [n le 7 select I[n] else 7*Self(n-1)-21*Self(n-2)+35*Self(n-3)-35*Self(n-4)+21*Self(n-5)-7*Self(n-6)+Self(n-7): n in [1..40]];
(Sage) m=6; [sum((binomial(2*m, k)*binomial(n, k)) for k in (0..m)) for n in (0..40)] # Bruno Berselli, Sep 22 2014
CROSSREFS
KEYWORD
nonn,easy
AUTHOR
Vincenzo Librandi, Sep 22 2014
STATUS
approved