login
a(n) = binomial(2n+1, n+1)*binomial(n+2, 2).
7

%I #24 May 14 2022 04:33:22

%S 1,9,60,350,1890,9702,48048,231660,1093950,5080790,23279256,105462084,

%T 473227300,2106121500,9307051200,40873466520,178520875830,

%U 775924068150,3357800061000,14473885526100,62168784497820,266168518910580

%N a(n) = binomial(2n+1, n+1)*binomial(n+2, 2).

%H G. C. Greubel, <a href="/A085373/b085373.txt">Table of n, a(n) for n = 0..1000</a>

%F From _David Callan_, Nov 20 2003: (Start)

%F a(n-1) = Sum_{1<=i1<=i2<=...<=in<=n} (i1 + i2 + ... + in).

%F G.f.: (1 - x)/(1 - 4*x)^(5/2). (End)

%F a(n) = A119578(n+1)/2. - _Zerinvary Lajos_, Jun 19 2008

%F From _Amiram Eldar_, May 14 2022: (Start)

%F Sum_{n>=0} 1/a(n) = 4*Pi^2/9 - 4*Pi/sqrt(3) + 4.

%F Sum_{n>=0} (-1)^n/a(n) = 8*sqrt(5)*log(phi) - 16*log(phi)^2 - 4, where phi = A001622 is the golden ratio. (End)

%t Table[Binomial[2n+1, n+1]Binomial[n+2, 2], {n, 0, 30}]

%o (PARI) a(n)=binomial(2*n+1,n+1)*binomial(n+2,2)

%o (Magma) [Binomial(2*n+1, n+1)*Binomial(n+2, 2): n in [0..30]]; // _G. C. Greubel_, Feb 12 2019

%o (Sage) [binomial(2*n+1, n+1)*binomial(n+2, 2) for n in (0..30)] # _G. C. Greubel_, Feb 12 2019

%Y Cf. A000984, A002544, A119578.

%K easy,nonn

%O 0,2

%A Mario Catalani (mario.catalani(AT)unito.it), Jun 26 2003