A001448 a(n) = binomial(4n,2n) or (4*n)!/((2*n)!*(2*n)!).
1, 6, 70, 924, 12870, 184756, 2704156, 40116600, 601080390, 9075135300, 137846528820, 2104098963720, 32247603683100, 495918532948104, 7648690600760440, 118264581564861424, 1832624140942590534, 28453041475240576740, 442512540276836779204, 6892620648693261354600
Offset: 0
Examples
a(n) = (1/Pi)*Integral_{x=0..4} x^(2n)/sqrt(4-(x-2)^2) dx. - _Paul Barry_, Sep 17 2010
G.f. = 1 + 6*x + 70*x^2 + 924*x^3 + 12870*x^4 + 184756*x^5 + 2704156*x^6 + ...
Links
- T. D. Noe, Table of n, a(n) for n = 0..100
- R. J. Chapman, T. Y. Chowa, A. Khetana, D. P. Moulton and R. J. Waters, Simple formulas for lattice paths avoiding certain periodic staircase boundaries, Journal of Combinatorial Theory, Series A 116 (2009) 205-214.
- Maciej Dziemianczuk, On Directed Lattice Paths With Additional Vertical Steps, arXiv preprint arXiv:1410.5747 [math.CO], 2014.
- Quy Nhan, A finite sum of binomial coefficients, MathStackExchange, 2025-06-21.
- Tito Piezas III, Fermat primes and Binomial sums, tpiezas, A Collection of Algebraic IdentitiesJune 3, 2012.
- K. H. Pilehrood and T. H. Pilehrood, Jacobi Polynomials and Congruences Involving Some Higher-Order Catalan Numbers and Binomial Coefficients, J. Int. Seq. 18 (2015) 15.11.7.
- Ricardo A. Podestá, New identities for binary Krawtchouk polynomials, binomial coefficients and Catalan numbers, arXiv:1603.09156 [math.CO], 2016.
- Renzo Sprugnoli, Sums of reciprocals of the central binomial coefficients, INTEGERS 6 (2006) #A27.
Crossrefs
Programs
-
Magma
[Factorial(4*n)/(Factorial(2*n)*Factorial(2*n)): n in [0..20]]; // Vincenzo Librandi, Sep 13 2011
-
Maple
A001448 := n-> binomial(4*n,2*n) ;
-
Mathematica
Table[Binomial[4n,2n],{n,0,20}] (* Harvey P. Dale, Apr 26 2014 *) a[ n_] := If[ n < 0, 0, HypergeometricPFQ[ {-2 n, -2 n}, {1}, 1]]; (* Michael Somos, Oct 22 2014 *) -
PARI
a(n)=binomial(4*n,2*n) \\ Charles R Greathouse IV, Sep 13 2011
-
Python
from math import comb def A001448(n): return comb(n<<2,n<<1) # Chai Wah Wu, Aug 10 2023
Formula
a(n) = A000984(2*n).
Using Stirling's formula in sequence A000142 it is easy to get the asymptotic expression a(n) ~ 16^n / sqrt(2 * Pi * n). - Dan Fux (dan.fux(AT)OpenGaia.com or danfux(AT)OpenGaia.com), Apr 07 2001
From Wolfdieter Lang, Dec 13 2001: (Start)
G.f.: (1-y*((1+4*y)*c(y)-(1-4*y)*c(-y)))/(1-(4*y)^2) with y^2=x, c(y) = g.f. for A000108 (Catalan). (End)
a(n) ~ 2^(-1/2)*Pi^(-1/2)*n^(-1/2)*2^(4*n)*{1 - (1/16)*n^-1 + ...}. - Joe Keane (jgk(AT)jgk.org), Jun 11 2002
a(n) = (1/Pi)*Integral_{x=-2..2} (2+x)^(2*n)/sqrt((2-x)*(2+x)) dx. Peter Luschny, Sep 12 2011
G.f.: (1/2) * (1/sqrt(1+4*sqrt(x)) + 1/sqrt(1-4*sqrt(x))). - Mark van Hoeij, Oct 25 2011
Sum_{n>=1} 1/a(n) = 16/15 + Pi*sqrt(3)/27 - 2*sqrt(5)*log(phi)/25, [T. Trif, Fib Quart 38 (2000) 79] with phi=A001622. - R. J. Mathar, Jul 18 2012
D-finite with recurrence n*(2*n-1)*a(n) -2*(4*n-1)*(4*n-3)*a(n-1)=0. - R. J. Mathar, Dec 02 2012
G.f.: sqrt((1 + sqrt(1-16*x))/(2*(1-16*x))) = 1 + 6*x/(G(0)-6*x), where G(k) = 2*x*(4*k+3)*(4*k+1) + (2*k+1)*(k+1) - 2*x*(k+1)*(2*k+1)*(4*k+5)*(4*k+7)/G(k+1); (continued fraction). - Sergei N. Gladkovskii, Jun 30 2013
a(n) = hypergeom([1-2*n,-2*n],[2],1)*(2*n+1). - Peter Luschny, Sep 22 2014
From Michael Somos, Oct 22 2014: (Start)
0 = a(n)*(+65536*a(n+2) - 16896*a(n+3) + 858*a(n+4)) + a(n+1)*(-3584*a(n+2) + 1176*a(n+3) - 66*a(n+4)) + a(n+2)*(+14*a(n+2) - 14*a(n+3) + a(n+4)) for all n in Z.
0 = a(n)^2*(+196608*a(n+1)^2 - 40960*a(n+1)*a(n+2) + 2100*a(n+2)^2) + a(n)*a(n+1)*(-12288*a(n+1)^2 + 2840*a(n+1)*a(n+2) - 160*a(n+2)^2) + a(n+1)^2*(+180*a(n+1)^2 - 48*a(n+1)*a(n+2) + 3*a(n+2)^2) for all n in Z. (End)
a(n) = [x^n] ( (1 + x)^4/(1 - x)^2 )^n; exp( Sum_{n >= 1} a(n)*x^n/n ) = 1 + 6*x + 53*x^2 + 554*x^3 + ... = Sum_{n >= 0} A066357(n+1)*x^n. - Peter Bala, Jun 23 2015
a(n) = Sum_{i = 0..n} binomial(4*n,i) * binomial(3*n-i-1,n-i). - Peter Bala, Sep 29 2015
a(n) = A000984(n)*Product_{j=0..n} (2^j/(j!*(2*j-1)!!))*A068424(n, j)^2, with A068424 the falling factorial. See (5.4) in Podestá link. - Michel Marcus, Mar 31 2016
a(n) = GegenbauerC(2*n, -2*n, -1). - Peter Luschny, May 07 2016
a(n) = [x^n] 1/sqrt(1 - 4*x)^(2*n+1). - Ilya Gutkovskiy, Oct 10 2017
a(n) is the n-th moment of the positive weight function w(x) on (0,16), i.e. a(n) = Integral_{x=0..16} x^n*w(x) dx, n = 0,1,..., where w(x) = (1/(2*Pi))/(sqrt(4 - sqrt(x))*x^(3/4)). The function w(x) is the solution of the Hausdorff moment problem and is unique. - Karol A. Penson, Mar 06 2018
a(n) = (16^n*(Beta(2*n - 1/2, 1/2) - Beta(2*n - 1/2, 3/2)))/Pi. - Peter Luschny, Mar 06 2018
E.g.f.: hypergeom([1/4,3/4],[1/2,1],16*x). - Karol A. Penson, Mar 08 2018
From Peter Bala, Feb 16 2020: (Start)
a(m*p^k) == a(m*p^(k-1)) ( mod p^(3*k) ) for prime p >= 5 and positive integers m and k.
a(n) = [(x*y)^(2*n)] (1 + x + y)^(4*n). (End)
a(n) = (2^n/n!)*Product_{k = n..2*n-1} (2*k + 1). - Peter Bala, Feb 26 2023
a(n) = Sum_{k = 0..2*n} binomial(2*n+k-1, k). - Peter Bala, Nov 02 2024
Sum_{n>=0} (-1)^n/a(n) = 16/17 + 4*sqrt(34)*(sqrt(17)-2)*arctan(sqrt(2/(sqrt(17)-1)))/(289*sqrt(sqrt(17)-1)) + 2*sqrt(34)*(sqrt(17)+2)*log((sqrt(sqrt(17)+1)-sqrt(2))/(sqrt(sqrt(17)+1)+sqrt(2)))/(289*sqrt(sqrt(17)+1)) (Sprugnoli, 2006, Theorem 3.8, p. 11; Piezas, 2012). - Amiram Eldar, Nov 03 2024
For n >= 1, a(n) = Sum_{k=1}^n a(n-k) * A337350(n) = Sum_{k=1}^n a(n-k) * a(k) * (8k + 1) / (8k^2 + 2k - 1). For proof, see the Quy Nhan link. - Lucas A. Brown, Jun 26 2025
From Seiichi Manyama, Aug 09 2025: (Start)
a(n) = [x^n] 1/((1-x)^(n+1) * (1-2*x)^(2*n)).
a(n) = Sum_{k=0..n} 2^k * (-1)^(n-k) * binomial(4*n,k) * binomial(2*n-k,n-k).
a(n) = Sum_{k=0..n} 2^k * binomial(2*n+k-1,k) * binomial(2*n-k,n-k).
a(n) = 4^n * binomial((4*n-1)/2,n).
a(n) = [x^n] (1+4*x)^((4*n-1)/2). (End)
Comments