OFFSET
0,2
LINKS
G. C. Greubel, Table of n, a(n) for n = 0..1000
Jun Wang and Zhizheng Zhang, On extensions of Calkin's binomial identities, Discrete Math., 274 (2004), 331-342.
Index entries for linear recurrences with constant coefficients, signature (16,-112,448,-1120,1792,-1792,1024,-256).
FORMULA
a(n) = (1/21)*n*(n+1)*(381*n^5 + 921*n^4 + 381*n^3 + 39*n^2 - 746*n + 368) * 2^(n-7). (See Wang and Zhang, p. 334.)
From Chai Wah Wu, Jun 21 2016: (Start)
a(n) = 16*a(n-1) - 112*a(n-2) + 448*a(n-3) - 1120*a(n-4) + 1792*a(n-5) - 1792*a(n-6) + 1024*a(n-7) - 256*a(n-8) for n > 7.
G.f.: x*(2 + 227*x + 2364*x^2 + 4748*x^3 + 2096*x^4 - 72*x^5)/(1 - 2*x)^8. (End)
a(n) = 2^(n-7)*n*(381*n^6 + 1302*n^5 + 1302*n^4 + 420*n^3 - 707*n^2 - 378*n + 368)/21. - Ilya Gutkovskiy, Jun 21 2016
E.g.f.: (1/42)*x*(84 + 5271*x + 33278*x^2 + 57855*x^3 + 37086*x^4 + 9303*x^5 +
762*x^6)*exp(2*x). - G. C. Greubel, May 24 2022
MATHEMATICA
LinearRecurrence[{16, -112, 448, -1120, 1792, -1792, 1024, -256}, {0, 2, 259, 6284, 77180, 646960, 4235864, 23313408}, 40] (* G. C. Greubel, Jun 22 2016 *)
PROG
(Magma) [2^(n-7)*n*(381*n^6+1302*n^5+1302*n^4+420*n^3-707*n^2-378*n+368)/21: n in [0..40]]; // G. C. Greubel, May 24 2022
(SageMath) [n*(n+1)*(381*n^5 +921*n^4 +381*n^3 +39*n^2 -746*n +368)*2^(n-7)/21 for n in (0..40)] # G. C. Greubel, May 24 2022
CROSSREFS
KEYWORD
nonn,easy
AUTHOR
N. J. A. Sloane, Jan 04 2004
STATUS
approved