%I #16 Aug 07 2024 09:19:28
%S 0,4,16,436,3200,78004,675808,15919320,150266880,3450748180,
%T 34461586016,774842070600,8061900244736,178065876017176,
%U 1912172640160960,41596867935469936,458156035085377536
%N a(n) = (1/4)*(binomial(4*n, 2*n) - (-1)^n*binomial(2*n, n) + (1-(-1)^n)*binomial(2*n, n)^2).
%D The right-hand side of a binomial coefficient identity in H. W. Gould, Combinatorial Identities, Morgantown, 1972. (See (3.75) on page 31.)
%H G. C. Greubel, <a href="/A037976/b037976.txt">Table of n, a(n) for n = 0..825</a>
%F From _G. C. Greubel_, Jun 22 2022: (Start)
%F a(n) = Sum_{k=0..floor((n-1)/2)} binomial(2*n, 2*k+1)^2.
%F a(n) = (1/4)*( (2*n+1)*A000108(2*n) - (-1)^n*(n+1)*A000108(n) + (1-(-1)^n)*(n+1)^2*A000108(n)^2 ).
%F G.f.: (1/4)*(sqrt(1 + sqrt(1-16*x))/(sqrt(2)*sqrt(1-16*x)) - 1/sqrt(1+4*x)) + (1/(2*Pi))*( EllipticK(16*x) - EllipticK(-16*x)). (End)
%p A037976 := proc(n)
%p binomial(4*n,2*n)/4-(-1)^n*binomial(2*n,n)/4+(1-(-1)^n)*binomial(2*n,n)^2/4 ;
%p end proc:
%p seq(A037976(n),n=0..30) ; # _R. J. Mathar_, Jul 26 2015
%t With[{B=Binomial}, Table[(1/4)*(B[4*n,2*n] +B[2*n,n]^2 -2*(-1)^n*B[B[2*n,n] +1, 2]), {n,0,30}]] (* _G. C. Greubel_, Jun 22 2022 *)
%o (Magma) [(1/4)*((2*n+1)*Catalan(2*n) -(-1)^n*(n+1)*Catalan(n) +(1-(-1)^n)*(n+1)^2*Catalan(n)^2): n in [0..30]]; // _G. C. Greubel_, Jun 22 2022
%o (SageMath) b=binomial; [(1/4)*(b(4*n, 2*n) -(-1)^n*b(2*n, n) +(1-(-1)^n)*b(2*n, n)^2) for n in (0..30)] # _G. C. Greubel_, Jun 22 2022
%Y Cf. A000108, A037964, A037972.
%K nonn,easy
%O 0,2
%A _N. J. A. Sloane_