login

Reminder: The OEIS is hiring a new managing editor, and the application deadline is January 26.

The binomial sum a(n) = Sum_{k=0..n}(binomial(n,k)*binomial(n+1,k+1)*binomial(n+2,k+2)).
1

%I #13 Oct 06 2016 15:50:09

%S 1,7,43,281,1896,13112,92359,659941,4769758,34797170,255838760,

%T 1893389720,14091400480,105385445856,791504226943,5966958725021,

%U 45133376297922,342400478465678,2604549070175770,19860078537996958,151769147958738016

%N The binomial sum a(n) = Sum_{k=0..n}(binomial(n,k)*binomial(n+1,k+1)*binomial(n+2,k+2)).

%F a(n) = (1/2)*(n+1)^2*(n+2)*hypergeometric({-n,-n,-n},{2,3},-1).

%F From _Vaclav Kotesovec_, Oct 04 2016: (Start)

%F Recurrence: (n+2)^2*(6*n^3 + 3*n^2 - 3*n - 2)*a(n) = (42*n^5 + 147*n^4 + 147*n^3 - 8*n^2 - 60*n - 16)*a(n-1) + 8*(n-1)*n*(6*n^3 + 21*n^2 + 21*n + 4)*a(n-2).

%F a(n) ~ 2^(3*n+4)/(sqrt(3)*Pi*n).

%F (End)

%t Table[Sum[Binomial[n,k]Binomial[n+1,k+1]Binomial[n+2,k+2],{k,0,n}],{n,0,100}]

%o (Maxima) makelist((n+1)^2*(n+2)/2*hypergeometric([-n,-n,-n],[2,3],-1),n,0,12);

%o (PARI) a(n) = sum(k=0, n, (binomial(n,k)*binomial(n+1,k+1)*binomial(n+2,k+2))); \\ _Michel Marcus_, Oct 04 2016

%K nonn

%O 0,2

%A _Emanuele Munarini_, Oct 04 2016