|
|
A002896
|
|
Number of 2n-step polygons on cubic lattice.
(Formerly M4285 N1791)
|
|
32
|
|
|
1, 6, 90, 1860, 44730, 1172556, 32496156, 936369720, 27770358330, 842090474940, 25989269017140, 813689707488840, 25780447171287900, 825043888527957000, 26630804377937061000, 865978374333905289360, 28342398385058078078010, 932905175625150142902300
(list;
graph;
refs;
listen;
history;
text;
internal format)
|
|
|
OFFSET
|
0,2
|
|
COMMENTS
|
Number of walks with 2n steps on the cubic lattice Z^3 beginning and ending at (0,0,0).
If A is a random matrix in USp(6) (6 X 6 complex matrices that are unitary and symplectic) then a(n) is the 2n-th moment of tr(A^k) for all k >= 7. - Andrew V. Sutherland, Mar 24 2008
Diagonal of the rational function R(x,y,z,w) = 1/(1 - (w*x*y + w*x*z + w*y + x*z + y + z)). - Gheorghe Coserea, Jul 14 2016
Constant term in the expansion of (x + 1/x + y + 1/y + z + 1/z)^(2n). - Harry Richman, Apr 29 2020
|
|
REFERENCES
|
N. J. A. Sloane, A Handbook of Integer Sequences, Academic Press, 1973 (includes this sequence).
N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).
|
|
LINKS
|
|
|
FORMULA
|
a(n) = C(2*n, n)*Sum_{k=0..n} C(n, k)^2*C(2*k, k).
a(n) = (4^n*p(1/2, n)/n!)*hypergeom([-n, -n, 1/2], [1, 1], 4), where p(a, k) = Product_{i=0..k-1} (a+i).
E.g.f.: Sum_{n>=0} a(n)*x^(2*n)/(2*n)! = BesselI(0, 2*x)^3. - Corrected by Christopher J. Smyth, Oct 29 2012
D-finite with recurrence: n^3*a(n) = 2*(2*n-1)*(10*n^2-10*n+3)*a(n-1) - 36*(n-1)*(2*n-1)*(2*n-3)*a(n-2). - Vladeta Jovovic, Jul 16 2004
An asymptotic formula follows immediately from an observation of Bruce Richmond and me in SIAM Review - 31 (1989, 122-125). We use Hayman's method to find the asymptotic behavior of the sum of squares of the multinomial coefficients multi(n, k_1, k_2, ..., k_m) with m fixed. From this one gets a_n ~ (3/4)*sqrt(3)*6^(2*n)/(Pi*n)^(3/2). - Cecil C Rousseau (ccrousse(AT)memphis.edu), Mar 14 2006
G.f.: (1/sqrt(1+12*z)) * hypergeom([1/8,3/8],[1],64/81*z*(1+sqrt(1-36*z))^2*(2+sqrt(1-36*z))^4/(1+12*z)^4) * hypergeom([1/8, 3/8],[1],64/81*z*(1-sqrt(1-36*z))^2*(2-sqrt(1-36*z))^4/(1+12*z)^4). - Sergey Perepechko, Jan 26 2011
G.f.: (1/2)*(10-72*x-6*(144*x^2-40*x+1)^(1/2))^(1/2)*hypergeom([1/6, 1/3],[1],54*x*(108*x^2-27*x+1+(9*x-1)*(144*x^2-40*x+1)^(1/2)))^2. - Mark van Hoeij, Nov 12 2011
0 = (-x^2+40*x^3-144*x^4)*y''' + (-3*x+180*x^2-864*x^3)*y'' + (-1+132*x-972*x^2)*y' + (6-108*x)*y, where y is the g.f. - Gheorghe Coserea, Jul 14 2016
a(n) = (1/Pi)^3*Integral_{0 <= x, y, z <= Pi} (2*cos(x) + 2*cos(y) + 2*cos(z))^(2*n) dx dy dz. - Peter Bala, Feb 10 2022
a(n) = Sum_{i+j+k=n, 0<=i,j,k<=n} multinomial(2n [i,i,j,j,k,k]). - Shel Kaphan, Jan 16 2023
Sum_{k>=0} a(k)/36^k = A086231 = (sqrt(3)-1) * (Gamma(1/24) * Gamma(11/24))^2 / (32*Pi^3). - Vaclav Kotesovec, Apr 23 2023
|
|
EXAMPLE
|
1 + 6*x + 90*x^2 + 1860*x^3 + 44730*x^4 + 1172556*x^5 + 32496156*x^6 + ...
|
|
MAPLE
|
a := proc(n) local k; binomial(2*n, n)*add(binomial(n, k)^2 *binomial(2*k, k), k=0..n); end;
# second Maple program
a:= proc(n) option remember; `if`(n<2, 5*n+1,
(2*(2*n-1)*(10*n^2-10*n+3) *a(n-1)
-36*(n-1)*(2*n-1)*(2*n-3) *a(n-2)) /n^3)
end:
A002896 := n -> binomial(2*n, n)*hypergeom([1/2, -n, -n], [1, 1], 4):
|
|
MATHEMATICA
|
f[n_] := 4^n*Gamma[n + 1/2]*Sum[Binomial[n, k]^2 Binomial[2 k, k], {k, 0, n}]/(Sqrt[Pi]*n!); Array[f, 17, 0] (* Robert G. Wilson v, Oct 29 2011 *)
Table[Binomial[2n, n]Sum[Binomial[n, k]^2 Binomial[2k, k], {k, 0, n}], {n, 0, 20}] (* Harvey P. Dale, Jan 24 2012 *)
a[ n_] := If[ n < 0, 0, HypergeometricPFQ[ {-n, -n, 1/2}, {1, 1}, 4] Binomial[ 2 n, n]] (* Michael Somos, May 21 2013 *)
|
|
PROG
|
(PARI) a(n)=binomial(2*n, n)*sum(k=0, n, binomial(n, k)^2*binomial(2*k, k)) \\ Charles R Greathouse IV, Oct 31 2011
(Sage)
x, y, n = 1, 6, 1
while True:
yield x
n += 1
x, y = y, ((4*n-2)*((10*(n-1)*n+3)*y-18*(n-1)*(2*n-3)*x))//n^3
|
|
CROSSREFS
|
|
|
KEYWORD
|
nonn,easy,walk,nice
|
|
AUTHOR
|
|
|
STATUS
|
approved
|
|
|
|