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
a(n) = 2*
A001700(2*n-1) = (2*n+1)*C(2*n), n >= 1, C(n) :=
A000108(n) (Catalan).
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
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) 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
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
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