Using Stirling's formula in
A000142 it is easy to get the asymptotic expression a(n) ~ 1/2 * sqrt(3) * 27^n / (Pi*n) - Dan Fux (dan.fux(AT)OpenGaia.com or danfux(AT)OpenGaia.com), Apr 07 2001
O.g.f.: hypergeom([1/3, 2/3], [1], 27*x).
E.g.f.: hypergeom([1/3, 2/3], [1, 1], 27*x).
Integral representation as n-th moment of a positive function on [0, 27]:
a(n) = int( x^n*(-1/24*(3*sqrt(3)*hypergeom([2/3, 2/3], [4/3], 1/27*x)* Gamma(2/3)^6*x^(1/3) - 8*hypergeom([1/3, 1/3], [2/3], 1/27*x)*Pi^3)/Pi^3 /x^(2/3)/Gamma(2/3)^3), x=0..27). This representation is unique. (End)
a(n) = Sum_{k=-n..n} (-1)^k*binomial(2*n, n+k)^3. -
Benoit Cloitre, Mar 02 2005
G.f. satisfies: A(x^3) = A( x*(1+3*x+9*x^2)/(1+6*x)^3 )/(1+6*x). -
Paul D. Hanna, Oct 29 2010
D-finite with recurrence: n^2*a(n) - 3*(3*n-1)*(3*n-2)*a(n-1) = 0. -
R. J. Mathar, Dec 04 2012
0 = a(n)^2*(472392*a(n+1)^2 - 83106*a(n+1)*a(n+2) + 3600*a(n+2)^2) + a(n)*a(n+1)*(-8748*a(n+1)^2 + 1953*a(n+1)*a(n+2) - 120*a(n+2)^2) + a(n+1)^2*(36*a(n+1)^2 - 12*a(n+1)*a(n+2) + a(n+2)^2) for all n in Z. -
Michael Somos, Oct 22 2014
0 = x*(27*x-1)*y'' + (54*x-1)*y' + 6*y, where y is g.f. -
Gheorghe Coserea, Jul 06 2016
a(n) = 3*binomial(2*n - 1,n)*binomial(3*n - 1,n) = 3*[x^n] 1/(1 - x)^n * [x^n] 1/(1 - x)^(2*n) for n >= 1.
a(n) = binomial(2*n,n)*binomial(3*n,n) = ([x^n](1 + x)^(2*n)) *([x^n](1 + x)^(3*n)) = [x^n](F(x)^(6*n)), where F(x) = 1 + x + 2*x^2 + 14*x^3 + 127*x^4 + 1364*x^5 + 16219*x^6 + ... appears to have integer coefficients. Cf.
A002894.
This sequence occurs as the right-hand side of several binomial sums:
Sum_{k = 0..2*n} (-1)^(n+k)*binomial(2*n,k)^3 = a(n) (Dixon's identity).
Sum_{k = 0..n} binomial(n,k)*binomial(2*n,n - k)*binomial(3*n + k,k) = a(n) (Gould, Vol. 4, 6.86)
Sum_{k = 0..n} (-1)^(n+k)*binomial(n,k)*binomial(2*n + k,n)*binomial(3*n + k,n) = a(n).
Sum_{k = 0..n} binomial(n,k)*binomial(2*n + k,k)*binomial(3*n,n - k) = a(n).
Sum_{k = 0..n} (-1)^(k)*binomial(n,k)*binomial(3*n - k,n)*binomial(4*n - k,n) = a(n).
Sum_{k = 0..2*n} (-1)^(n+k)*binomial(2*n + k,2*n - k)*binomial(2*k,k)*binomial(4*n - k,2*n) = a(n) (see Gould, Vol.5, 9.23).
Sum_{k = 0..2*n} (-1)^k*binomial(3*n,k)*binomial(3*n - k,n)^3 = a(n) (Sprugnoli, Section 2.9, Table 10, p. 123). (End)
G.f.: F(x) = 1/(2*Pi) Integral_{z=0..2*Pi} 2F1(1/3,2/3; 1/2; 27*x*sin^2(z)) dz.
With G(x) = x*2F1(1/3,2/3; 2; 27*x): F(x) = d/dx G(x). (Cf.
A007004) (End)
F(x)*G(1/27-x) + F(1/27-x)*G(x) = 1/(4*Pi*sqrt(3)). -
Bradley Klee, Sep 29 2018
a(n) = Sum_{k = n..2*n} binomial(2*n,k)^2 * binomial(k,n). Cf.
A001459.
a(n*p^k) == a(n*p^(k-1)) ( mod p^(3*k) ) for any prime p >= 5 and any positive integers n and k (write a(n) as C(3*n,2*n)*C(2*n,n) and apply Mestrovic, equation 39, p. 12). (End)
Occurs on the right-hand side of the binomial sum identities Sum_{k = -n..n} (-1)^k * (n + x - k) * binomial(2*n, n+k)^3 = (x + n)*a(n) and Sum_{k = -n..n} (-1)^k * (n + x - k)^3 * binomial(2*n, n+k)^3 = x*(x + n)*(x + 2*n)*a(n) (x arbitrary). Compare with Dixon's identity: Sum_{k = -n..n} (-1)^k * binomial(2*n, n+k)^3 = a(n). -
Peter Bala, Jul 31 2023
a(n) = (-1)^n * [x^(2*n)] ( (1 - x)^(4*n) * Legendre_P(2*n, (1 + x)/(1 - x)) ).
The following hold for n >= 1:
a(n) = Sum_{k = 0.. 2*n} (-1)^(n+k) * k/n * binomial(2*n, k)^3 = 3/2 * Sum_{k = 0.. 2*n} (-1)^(n+k) * (k/n)^2 * binomial(2*n, k)^3.
a(n) = 3/2 * Sum_{0..2*n-1} (-1)^(n+k) * k/n * binomial(2*n, k)^2*binomial(2*n-1, k).
a(n) = 3 * Sum_{0..2*n-1} (-1)^(n+k) * k/n * binomial(2*n, k)*binomial(2*n-1, k)^2. (End)
a(n) = Sum_{k = 0..n} (-1)^(n+k) * binomial(n, k) *
A108625(2*n, k) (verified using the MulZeil procedure in Doron Zeilberger's MultiZeilberger package). -
Peter Bala, Oct 15 2024
Comments