OFFSET
0,2
COMMENTS
A fast algorithm follows from Theorem 5 of Brent et al. article.
LINKS
Colin Barker, Table of n, a(n) for n = 0..800
Richard P. Brent, Hideyuki Ohtsuka, Judy-anne H. Osborn, Helmut Prodinger, Some binomial sums involving absolute values, arXiv:1411.1477v2 [math.CO], 2016.
Index entries for linear recurrences with constant coefficients, signature (80,-2560,40960,-327680,1048576).
FORMULA
a(n) = Sum_{k=-n..n} (Sum_{l=-n..n} binomial(2*n, n+k)*binomial(2*n, n+l)*abs(k^2 - l^2)^4).
From Colin Barker, Feb 11 2016: (Start)
a(n) = 4^(2*n-1)*n*(36*n^3-84*n^2+67*n-17).
a(n) = 80*a(n-1)-2560*a(n-2)+40960*a(n-3)-327680*a(n-4)+1048576*a(n-5) for n>4.
G.f.: 8*x*(1+1024*x+67840*x^2+417792*x^3) / (1-16*x)^5.
(End)
PROG
(PARI) a(n) = sum(k=-n, n, sum(l=-n, n, binomial(2*n, n+k)*binomial(2*n, n+l)*abs(k^2 - l^2)^4));
(PARI) concat(0, Vec(8*x*(1+1024*x+67840*x^2+417792*x^3)/(1-16*x)^5 + O(x^20))) \\ Colin Barker, Feb 11 2016
CROSSREFS
KEYWORD
easy,nonn
AUTHOR
Richard P. Brent, Jan 27 2016
STATUS
approved