cp's OEIS Frontend

This is a front-end for the Online Encyclopedia of Integer Sequences, made by Christian Perfect. The idea is to provide OEIS entries in non-ancient HTML, and then to think about how they're presented visually. The source code is on GitHub.

Showing 1-2 of 2 results.

A018224 a(n) = binomial(n, floor(n/2))^2 = A001405(n)^2.

Original entry on oeis.org

1, 1, 4, 9, 36, 100, 400, 1225, 4900, 15876, 63504, 213444, 853776, 2944656, 11778624, 41409225, 165636900, 590976100, 2363904400, 8533694884, 34134779536, 124408576656, 497634306624, 1828114918084, 7312459672336, 27043120090000, 108172480360000, 402335398890000
Offset: 0

Views

Author

Keywords

Comments

a(n) is also the number of rooted two-vertex (or, dually, two-face) regular planar maps of valency n+1. - Valery A. Liskovets, Oct 19 2005
If A is a random matrix in USp(4) (4 X 4 complex matrices that are unitary and symplectic), then a(n)=(-1)^n*E[(tr(A^4))^n]. - Andrew V. Sutherland, Apr 01 2008
Number of square lattice walks with unit steps in all four directions (NSWE), starting at the origin, ending on the y-axis, and never going below the x-axis. Row sums of A378061. - Peter Luschny, Dec 08 2024

Examples

			The 9 lattice walks defined in the comments: 'NNN', 'NNS', 'NSN', 'NWE', 'NEW', 'WNE', 'WEN', 'ENW', 'EWN'.
		

Crossrefs

Programs

  • Maple
    s := x -> (1+x)*EllipticK(x)/(x*Pi/2)-1/x:
    seq(4^n*coeff(series(s(x),x,n+2),x,n),n=0..23); # Peter Luschny, Oct 14 2015
  • Mathematica
    (* Note that Mathematica uses a different definition of the EllipticK function. *)
    CoefficientList[Series[(-Pi + (2 + 8 x) EllipticK[16 x^2])/(4 Pi x), {x,0,23}], x] (* Peter Luschny, Oct 14 2015 *)
    Table[Binomial[n,Floor[n/2]]^2,{n,0,30}] (* Harvey P. Dale, Dec 02 2022 *)
  • PARI
    vector(50, n, n--; binomial(n, n\2)^2) \\ Altug Alkan, Oct 14 2015

Formula

E.g.f.: BesselI(0, 2*x)*(BesselI(0, 2*x)+BesselI(1, 2*x)). - Vladeta Jovovic, Jun 12 2005
G.f. (1+1/(4*x))*hypergeom([1/2, 1/2],[1],16*x^2)-1/(4*x). - Mark van Hoeij, Oct 13 2009
a(n) = (n!/(floor(n/2)!*floor((n+1)/2)!))^2. - Peter Luschny, Apr 29 2014
a(n) = A056040(n) * A056040(n+1) / (n+1). - Peter Luschny, Apr 29 2014
a(n) = 4^n*[x^n]((1+x)*EllipticK(x)/(x*Pi/2)-1/x). - Peter Luschny, Oct 14 2015
a(n) ~ 4^n*((2*n+3)/(2*n+1))^((-1)^n/2)/((n+1)*Pi/2). - Peter Luschny, Oct 14 2015
a(n) = Sum_{k=0..n} (-1)^k*binomial(n,k)*C(k)*binomial(2*n-2*k,n-k) where C(k) are Catalan numbers (A000108), see Prodinger. - Michel Marcus, Nov 19 2019
From Peter Bala, Jul 03 2023: (Start)
Right hand side of the binomial sum identity (1/2)*Sum_{k = 0..n+1} (-1)^k*4^(n+1-k)*binomial(n+1,k)*binomial(n+k,k)*binomial(2*k,k) = a(n).
a(n) = (1/2)*4^(n+1) * hypergeom([n+1, -n-1, 1/2], [1, 1], 1).
P-recursive:
(2*n - 1)*(n + 1)^2*a(n) = 4*(2*n^2 - 1)*a(n-1) + 16*(2*n + 1)*(n - 1)^2*a(n-2) with a(0) = a(1) = 1. (End)

A138544 Moment sequence of tr(A^4) in USp(6).

Original entry on oeis.org

1, -1, 4, -9, 42, -130, 660, -2415, 12810, -51786, 281736, -1216446, 6727644, -30440124, 170316432, -798126615, 4504487130, -21692469370, 123255492360, -606672653730, 3465702008340, -17366224451940, 99645553785960, -506814533253210, 2918768920720380, -15034038412333500
Offset: 0

Views

Author

Andrew V. Sutherland, Mar 24 2008

Keywords

Comments

If A is a random matrix in the compact group USp(6) (6 X 6 complex matrices which are unitary and symplectic), then a(n) = E[(tr(A^4))^n] is the n-th moment of the trace of A^4. See A138545 for central moments.

Examples

			a(3) = -9 because E[(tr(A^4))^3] = -9 for a random matrix A in USp(6).
		

Crossrefs

Formula

mgf is A(z) = det[F_{i+j-2}(z)], 1<=i,j<=3, where F_m(z) = Sum_j binomial(m,j)(B_{(2j-m)/4}(z)-B_{(2j-m+2)/4}(z)) and B_v(z)=0 for non-integer v and otherwise B_v(z)=I_v(2z), with I_v(z) the hyperbolic Bessel function (of the first kind) of order v.
Showing 1-2 of 2 results.