login
a(n) = Sum_{k=0..n} binomial(n,k)*binomial(n+k,k)^2.
5

%I #36 Jan 12 2019 20:11:57

%S 1,5,55,749,11251,178835,2949115,49906925,860905315,15071939255,

%T 266982872905,4774722189275,86070844191775,1561948324845095,

%U 28507384046515555,522867506128197869,9631571375362268515,178094411589895650815,3304192479145474141741,61487420580006795749999

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

%C Diagonal of rational function 1/(1 - x - y - z - x*y + x*z + x*y*z). - _Gheorghe Coserea_, Jul 01 2018

%H Vincenzo Librandi, <a href="/A112019/b112019.txt">Table of n, a(n) for n = 0..200</a>

%H C. Elsner, <a href="http://www.fq.math.ca/Papers1/43-1/paper43-1-5.pdf">On recurrence formulas for sums involving binomial coefficients</a>, Fib. Q., 43,1 (2005), 31-45.

%H V. Kotesovec, <a href="https://oeis.org/wiki/User:Vaclav_Kotesovec">Asymptotic of generalized Apery sequences with powers of binomial coefficients</a>, Nov 04 2012

%F a(n) = 3F2( {-n, 1 + n, 1 + n} ; {1, 1} )(-1). - _Olivier Gérard_, Apr 23 2009

%F a(n) ~ (1+r)^(4*n+5/2)/r^(3*n+5/2)/(2*Pi*n)*sqrt((1-r)/(3-r)), where r is positive real root of the equation (1-r)*(1+r)^2=r^3, r = 1/6*((44-3*sqrt(177))^(1/3)+(44+3*sqrt(177))^(1/3)-1) = 0.82948354095849... - _Vaclav Kotesovec_, Nov 04 2012

%F Recurrence: 2*n^2*(59*n - 83)*a(n) = (2301*n^3 - 5538*n^2 + 3797*n - 800)*a(n-1) + 5*(59*n^3 - 201*n^2 + 213*n - 64)*a(n-2) + (59*n - 24)*(n-2)^2*a(n-3). - _Vaclav Kotesovec_, Nov 04 2012

%F G.f. y=A(x) satisfies: 0 = x*(5*x + 8)*(x^3 + 5*x^2 + 39*x - 2)*y'' + (15*x^4 + 82*x^3 + 315*x^2 + 624*x - 16)*y' + (5*x^3 + 21*x^2 + 80)*y. - _Gheorghe Coserea_, Jul 01 2018

%p seq(add((multinomial(n+k,n-k,k,k))*binomial(n+k,k),k=0..n),n=0..19); # _Zerinvary Lajos_, Oct 18 2006

%p ogf := hypergeom([1/12,5/12],[1], -1728*(x^3+5*x^2+39*x-2)*x^4 / (x^4+4*x^3+30*x^2-20*x+1)^3 ) / (x^4+4*x^3+30*x^2-20*x+1)^(1/4);

%p series(ogf, x=0, 30); # _Mark van Hoeij_, Jan 22 2013

%t Table[HypergeometricPFQ[{-n, 1 + n, 1 + n}, {1, 1}, -1], {n, 0, 20}] (* _Olivier Gérard_, Apr 23 2009 *)

%t Table[Sum[Binomial[n,k]*Binomial[n+k,k]^2,{k,0,n}],{n,0,20}] (* _Vaclav Kotesovec_, Nov 04 2012 *)

%o (PARI) a(n) = sum(k=0, n, binomial(n,k)*binomial(n+k,k)^2); \\ _Michel Marcus_, Mar 09 2016

%Y Cf. A218693, A111968, A014178, A014180, A218689, A218692.

%K nonn

%O 0,2

%A _N. J. A. Sloane_, Nov 28 2005