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.

A000139 a(n) = 2*(3*n)! / ((2*n+1)!*(n+1)!).

Original entry on oeis.org

2, 1, 2, 6, 22, 91, 408, 1938, 9614, 49335, 260130, 1402440, 7702632, 42975796, 243035536, 1390594458, 8038677054, 46892282815, 275750636070, 1633292229030, 9737153323590, 58392041019795, 352044769046880, 2132866978427640, 12980019040145352, 79319075627675556
Offset: 0

Views

Author

N. J. A. Sloane, entry revised Apr 24 2012

Keywords

Comments

This sequence arises in many different contexts, and over the years it has had several different definitions. I have now changed the definition back to one of the earlier ones, a self-contained formula. - N. J. A. Sloane, Apr 24 2012
The number of 2-stack sortable permutations on n letters (n >= 1).
The number of rooted non-separable planar maps with n+1 edges. - Valery A. Liskovets, Mar 17 2005
The shifted sequence starting with a(1): Number of quadrangular dissections of a square, counted by the number of vertices. Rooted, non-separable planar maps with no multiple edges, in which each non-root face has degree 4.
Number of left ternary trees having n nodes (n>=1). - Emeric Deutsch, Jul 23 2006
A combinatorial interpretation for this sequence in terms of a family of plane trees is given in [Schaeffer, Corollary 2 with k = 3]. - Peter Bala, Oct 12 2011
Number of canopy intervals in the Tamari lattices, see [Préville-Ratelle and Viennot, section 6]. - F. Chapoton, Apr 19 2015
The number of fighting fish (branching polyominoes). - David Bevan, Jan 10 2018
The number of 1324-avoiding dominoes (gridded permutations). - David Bevan, Jan 10 2018
For n > 0, a(n) is the number of simple strong triangulations of a fixed quadrilateral with n interior nodes. See A210664. - Andrew Howroyd, Feb 24 2021
Conjecture: a(n) is odd iff n is a term of A022341. - Peter Bala, Jul 24 2025

Examples

			G.f. = 2 + x + 2*x^2 + 6*x^3 + 22*x^4 + 91*x^5 + 408*x^6 + 1938*x^7 + ...

References

  • Miklos Bona, editor, Handbook of Enumerative Combinatorics, CRC Press, 2015, page 365.
  • Eric S. Egge, Defying God: The Stanley-Wilf Conjecture, Stanley-Wilf Limits, and a Two-Generation Explosion of Combinatorics, pp. 65-82 of "A Century of Advancing Mathematics", ed. S. F. Kennedy et al., MAA Press 2015.
  • J. L. Gross and J. Yellen, eds., Handbook of Graph Theory, CRC Press, 2004; p. 714.
  • S. Kitaev, Patterns in Permutations and Words, Springer-Verlag, 2011. See p. 399 Table A.7
  • W. F. Lunnon, Counting polyominoes, pp. 347-372 of A. O. L. Atkin and B. J. Birch, editors, Computers in Number Theory. Academic Press, NY, 1971.
  • 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).
  • R. P. Stanley, Enumerative Combinatorics, Cambridge, Vol. 2, 1999; see Problem 6.41.

Links

Crossrefs

Cf. A000142, A000309, A001764, A005802, A006335, A004677, A121545, A218473, A233389.
Row m=1 of A210664(n>0).

Programs

  • Haskell
    a000139 0 = 2
a000139 n = ((3 * n) `a007318` (2 * n + 1)) `div` a000217 n
-- Reinhard Zumkeller, Feb 17 2013
  • Magma
    [2*Factorial(3*n)/(Factorial(2*n+1)*Factorial(n+1)): n in [0..25]]; // Vincenzo Librandi, Apr 20 2015
  • Maple
    A000139 := n->2*(3*n)!/((2*n+1)!*((n+1)!)): seq(A000139(n), n=0..23);
  • Mathematica
    Table[(2(3n)!)/((2n+1)!(n+1)!),{n,0,30}] (* Harvey P. Dale, Mar 31 2013 *)
  • PARI
    a(n)=binomial(3*n,n)*2/((n+1)*(2*n+1)); \\ Joerg Arndt, Jul 21 2014
  • Python
    from sympy import binomial
def A000139(n): return (binomial(3*n, n)*2)//((n+1)*(2*n+1))
  • Python
    A000139_list = [2]
for n in range(1,30):
    A000139_list.append(3*(3*n-2)*(3*n-1)*A000139_list[-1]//(2*n+2)//(2*n+1)) # Chai Wah Wu, Apr 02 2021
  • Sage
    def A000139(n): return (binomial(3*n, n)*2)//((n+1)*(2*n+1))
[A000139(n) for n in (0..23)]  # Peter Luschny, Jun 17 2013

Formula

a(n) = 2*binomial(3*n, 2*n+1)/(n*(n+1)), or 2*(3*n)!/((2*n+1)!*((n+1)!)).
Using Stirling's formula in A000142 it is easy to get the asymptotic expression a(n) ~ (27/4)^n / (sqrt(Pi*n / 3) * (2*n + 1) * (n + 1)). - Dan Fux (dan.fux(AT)OpenGaia.com or danfux(AT)OpenGaia.com), Apr 13 2001
G.f.: A(z) = 2 + z*B(z), where B(z) = 1 - 8*z + 2*z*(5-6*z)*B - 2*z^2*(1+3*z)*B^2 - z^4*B^3.
G.f.: (2/(3*x)) * (hypergeom([-2/3, -1/3],[1/2],(27/4)*x)-1). - Mark van Hoeij, Nov 02 2009
G.f.: (2-3*R)/(R-1)^2 where R := RootOf(x-t*(t-1)^2,t) is an algebraic function in Maple notation. - Mark van Hoeij, Nov 08 2011
G.f.: 2*Q(0), where Q(k) = 1 + 3*x*(3*k+1)*(6*k+1)/(2*(k+1)*(4*k+3) - 6*x*(k+1)*(3*k+2)*(4*k+3)*(6*k+5)/(3*x*(3*k+2)*(6*k+5) + (2*k+3)*(4*k+5)/Q(k+1))); (continued fraction). - Sergei N. Gladkovskii, Apr 30 2013
E.g.f.: 2*Q(0), where Q(k) = 1 + 3*x*(3*k+1)*(6*k+1)/(2*(k+1)*(2*k+1)*(4*k+3) - 6*x*(k+1)*(2*k+1)*(3*k+2)*(4*k+3)*(6*k+5)/(3*x*(3*k+2)*(6*k+5) + (2*k+2)*(2*k+3)*(4*k+5)/Q(k+1))); (continued fraction). - Sergei N. Gladkovskii, Apr 30 2013
a(n) = A007318(3*n, 2*n+1)/A000217(n) for n > 0. - Reinhard Zumkeller, Feb 17 2013
a(n) is the n-th Hausdorff moment of the positive function w(x) defined on (0,27) which is equal to w(x) = 3*sqrt(3)*2^(2/3)*(3-sqrt(81-12*x)/9)*(1+sqrt(81-12*x)/9)^(1/3)/(8*Pi*x^(2/3))-sqrt(3)*2^(1/3)*(3+sqrt(81-12*x)/9)*(1+sqrt(81-12*x)/9)^(-1/3)/(4*Pi*x^(1/3)), that is, a(n) is the integral Integral_{x=0..27/4} x^n*w(x) dx, n >= 0. The function w(x) is unique. - Karol A. Penson, Jun 17 2013
D-finite with recurrence 2*(n+1)*(2*n+1)*a(n) -3*(3*n-1)*(3*n-2)*a(n-1)=0. - R. J. Mathar, Aug 21 2014
G.f. A(z) is related to the g.f. M(z) of A000168 by M(z) = 1 + A(z*M(z)^2) (see Tutte 1963, equation 6.3). - Noam Zeilberger, Nov 02 2016
From Ilya Gutkovskiy, Jan 17 2017: (Start)
E.g.f.: 2*2F2(1/3,2/3; 3/2,2; 27*x/4).
Sum_{n>=0} 1/a(n) = (1/2)*3F2(1,3/2,2; 1/3,2/3; 4/27) = 2.226206199291261... (End)
G.f. A(z) is the solution to the initial value problem 4*A + 2*z*A' = 8 + 3*z*A + 9*z^2*A' + 2*z^2*A*A', A(0) = 2. - Bjarki Ágúst Guðmundsson, Jul 03 2017
a(n+1) = a(n)*3*(3*n+1)*(3*n+2)/(2*(n+2)*(2*n+3)). - Chai Wah Wu, Apr 02 2021
a(n) = 4*(3*n)!/(n!*(2*n+2)!). - Chai Wah Wu, Dec 15 2021
From Peter Bala, Feb 05 2022: (Start)
O.g.f.: A(x) = T(x)*(3 - T(x)), where T(x) = 1 + x*T(x)^3 is the o.g.f. of A001764.
(1/x)*(A(x) - 2)/(A(x) - 1) = 1 + x + 3*x^2 + 11*x^3 + 46*x^4 + 209*x^5 + ... is the o.g.f. of A233389.
1 + 2*x*A'(2*x)/A(2*x) = 1 + x + 7*x^2 + 61*x^3 + 591*x^4 + 6101*x^6 + ... is the o.g.f. of A218473.
Let B(x) = 1 + x*(A(x) - 1). Then x*B'(x)/B(x) = x + x^2 + 4*x^3 + 17*x^4 + 81*x^5 + ... is the o.g.f. of A121545. (End)

Started in 1964 by Neil J. A. Sloane | Maintained by The OEIS Foundation Inc.

Content is available under The OEIS End-User License Agreement.