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-10 of 10 results.

A001764 a(n) = binomial(3*n,n)/(2*n+1) (enumerates ternary trees and also noncrossing trees).

Original entry on oeis.org

1, 1, 3, 12, 55, 273, 1428, 7752, 43263, 246675, 1430715, 8414640, 50067108, 300830572, 1822766520, 11124755664, 68328754959, 422030545335, 2619631042665, 16332922290300, 102240109897695, 642312451217745, 4048514844039120, 25594403741131680, 162250238001816900
Offset: 0

Views

Author

Keywords

Comments

Smallest number of straight line crossing-free spanning trees on n points in the plane.
Number of dissections of some convex polygon by nonintersecting diagonals into polygons with an odd number of sides and having a total number of 2n+1 edges (sides and diagonals). - Emeric Deutsch, Mar 06 2002
Number of lattice paths of n East steps and 2n North steps from (0,0) to (n,2n) and lying weakly below the line y=2x. - David Callan, Mar 14 2004
With interpolated zeros, this has g.f. 2*sqrt(3)*sin(arcsin(3*sqrt(3)*x/2)/3)/(3*x) and a(n) = C(n+floor(n/2),floor(n/2))*C(floor(n/2),n-floor(n/2))/(n+1). This is the first column of the inverse of the Riordan array (1-x^2,x(1-x^2)) (essentially reversion of y-y^3). - Paul Barry, Feb 02 2005
Number of 12312-avoiding matchings on [2n].
Number of complete ternary trees with n internal nodes, or 3n edges.
Number of rooted plane trees with 2n edges, where every vertex has even outdegree ("even trees").
a(n) is the number of noncrossing partitions of [2n] with all blocks of even size. E.g.: a(2)=3 counts 12-34, 14-23, 1234. - David Callan, Mar 30 2007
Pfaff-Fuss-Catalan sequence C^{m}_n for m=3, see the Graham et al. reference, p. 347. eq. 7.66.
Also 3-Raney sequence, see the Graham et al. reference, p. 346-7.
The number of lattice paths from (0,0) to (2n,0) using an Up-step=(1,1) and a Down-step=(0,-2) and staying above the x-axis. E.g., a(2) = 3; UUUUDD, UUUDUD, UUDUUD. - Charles Moore (chamoore(AT)howard.edu), Jan 09 2008
a(n) is (conjecturally) the number of permutations of [n+1] that avoid the patterns 4-2-3-1 and 4-2-5-1-3 and end with an ascent. For example, a(4)=55 counts all 60 permutations of [5] that end with an ascent except 42315, 52314, 52413, 53412, all of which contain a 4-2-3-1 pattern and 42513. - David Callan, Jul 22 2008
Central terms of pendular triangle A167763. - Philippe Deléham, Nov 12 2009
With B(x,t)=x+t*x^3, the comp. inverse in x about 0 is A(x,t) = Sum_{j>=0} a(j) (-t)^j x^(2j+1). Let U(x,t)=(x-A(x,t))/t. Then DU(x,t)/Dt=dU/dt+U*dU/dx=0 and U(x,0)=x^3, i.e., U is a solution of the inviscid Burgers's, or Hopf, equation. Also U(x,t)=U(x-t*U(x,t),0) and dB(x,t)/dt = U(B(x,t),t) = x^3 = U(x,0). The characteristics for the Hopf equation are x(t) = x(0) + t*U(x(t),t) = x(0) + t*U(x(0),0) = x(0) + t*x(0)^3 = B(x(0),t). These results apply to all the Fuss-Catalan sequences with 3 replaced by n>0 and 2 by n-1 (e.g., A000108 with n=2 and A002293 with n=4), see also A086810, which can be generalized to A133437, for associahedra. - Tom Copeland, Feb 15 2014
Number of intervals (i.e., ordered pairs (x,y) such that x<=y) in the Kreweras lattice (noncrossing partitions ordered by refinement) of size n, see the Bernardi & Bonichon (2009) and Kreweras (1972) references. - Noam Zeilberger, Jun 01 2016
Number of sum-indecomposable (4231,42513)-avoiding permutations. Conjecturally, number of sum-indecomposable (2431,45231)-avoiding permutations. - Alexander Burstein, Oct 19 2017
a(n) is the number of topologically distinct endstates for the game Planted Brussels Sprouts on n vertices, see Ji and Propp link. - Caleb Ji, May 14 2018
Number of complete quadrillages of 2n+2-gons. See Baryshnikov p. 12. See also Nov 10 2014 comments in A134264. - Tom Copeland, Jun 04 2018
a(n) is the number of 2-regular words on the alphabet [n] that avoid the patterns 231 and 221. Equivalently, this is the number of 2-regular tortoise-sortable words on the alphabet [n] (see the Defant and Kravitz link). - Colin Defant, Sep 26 2018
a(n) is the number of Motzkin paths of length 3n with n steps of each type, with the condition that (1, 0) and (1, 1) steps alternate (starting with (1, 0)). - Helmut Prodinger, Apr 08 2019
a(n) is the number of uniquely sorted permutations of length 2n+1 that avoid the patterns 312 and 1342. - Colin Defant, Jun 08 2019
The compositional inverse o.g.f. pair in Copeland's comment above are related to a pair of quantum fields in Balduf's thesis by Theorem 4.2 on p. 92. - Tom Copeland, Dec 13 2019
The sequences of Fuss-Catalan numbers, of which this is the first after the Catalan numbers A000108 (the next is A002293), appear in articles on random matrices and quantum physics. See Banica et al., Collins et al., and Mlotkowski et al. Interpretations of these sequences in terms of the cardinality of specific sets of noncrossing partitions are provided by A134264. - Tom Copeland, Dec 21 2019
Call C(p, [alpha], g) the number of partitions of a cyclically ordered set with p elements, of cyclic type [alpha], and of genus g (the genus g Faa di Bruno coefficients of type [alpha]). This sequence counts the genus 0 partitions (non-crossing, or planar, partitions) of p = 3n into n parts of length 3: a(n) = C(3n, [3^n], 0). For genus 1 see A371250, for genus 2 see A371251. - Robert Coquereaux, Mar 16 2024
a(n) is the total number of down steps before the first up step in all 2_1-Dyck paths of length 3*n for n > 0. A 2_1-Dyck path is a lattice path with steps (1,2), (1,-1) that starts and ends at y = 0 and does not go below the line y = -1. - Sarah Selkirk, May 10 2020
a(n) is the number of pairs (A<=B) of noncrossing partitions of [n]. - Francesca Aicardi, May 28 2022
a(n) is the number of parking functions of size n avoiding the patterns 231 and 321. - Lara Pudwell, Apr 10 2023
Number of rooted polyominoes composed of n square cells of the hyperbolic regular tiling with Schläfli symbol {4,oo}. A rooted polyomino has one external edge identified, and chiral pairs are counted as two. A stereographic projection of the {4,oo} tiling on the Poincaré disk can be obtained via the Christensson link. - Robert A. Russell, Jan 27 2024
This is instance k = 3 of the family {C(k, n)}A130564.%20-%20_Wolfdieter%20Lang">{n>=0} given in a comment in A130564. - _Wolfdieter Lang, Feb 05 2024
The number of Apollonian networks (planar 3-trees) with n+3 vertices with a given base triangle. - Allan Bickle, Feb 20 2024
Number of rooted polyominoes composed of n tetrahedral cells of the hyperbolic regular tiling with Schläfli symbol {3,3,oo}. A rooted polyomino has one external face identified, and chiral pairs are counted as two. a(n) = T(n) in the second Beineke and Pippert link. - Robert A. Russell, Mar 20 2024

Examples

			a(2) = 3 because the only dissections with 5 edges are given by a square dissected by any of the two diagonals and the pentagon with no dissecting diagonal.
G.f. = 1 + x + 3*x^2 + 12*x^3 + 55*x^4 + 273*x^5 + 1428*x^6 + 7752*x^7 + 43263*x^8 + ...
		

References

  • Miklos Bona, editor, Handbook of Enumerative Combinatorics, CRC Press, 2015, page 23.
  • I. M. H. Etherington, On non-associative combinations, Proc. Royal Soc. Edinburgh, 59 (Part 2, 1938-39), 153-162.
  • I. M. H. Etherington, Some problems of non-associative combinations (I), Edinburgh Math. Notes, 32 (1940), pp. i-vi. Part II is by A. Erdelyi and I. M. H. Etherington, and is on pages vii-xiv of the same issue.
  • R. L. Graham, D. E. Knuth and O. Patashnik, Concrete Mathematics. Addison-Wesley, Reading, MA, 1990, pp. 200, 347. See also the Pólya-Szegő reference.
  • W. Kuich, Languages and the enumeration of planted plane trees. Nederl. Akad. Wetensch. Proc. Ser. A 73 = Indag. Math. 32, (1970), 268-280.
  • T. V. Narayana, Lattice Path Combinatorics with Statistical Applications. Univ. Toronto Press, 1979, p. 98.
  • G. Pólya and G. Szegő, Problems and Theorems in Analysis, Springer-Verlag, New York, Heidelberg, Berlin, 2 vols., 1972, Vol. 1, problem 211, p. 146 with solution on p. 348.
  • 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).

Crossrefs

Cf. A001762, A001763, A002294 - A002296, A006013, A025174, A063548, A064017, A072247, A072248, A134264, A143603, A258708, A256311, A188687 (binomial transform), A346628 (inverse binomial transform).
A column of triangle A102537.
Bisection of A047749 and A047761.
Row sums of triangles A108410 and A108767.
Second column of triangle A062993.
Mod 3 = A113047.
2D Polyominoes: A005034 (oriented), A005036 (unoriented), A369315 (chiral), A047749 (achiral), A000108 {3,oo}, A002293 {5,oo}.
3D Polyominoes: A007173 (oriented), A027610 (unoriented), A371350 (chiral), A371351 (achiral).
Cf. A130564 (for C(k, n) cases).

Programs

  • GAP
    List([0..25],n->Binomial(3*n,n)/(2*n+1)); # Muniru A Asiru, Oct 31 2018
    
  • Haskell
    a001764 n = a001764_list !! n
    a001764_list = 1 : [a258708 (2 * n) n | n <- [1..]]
    -- Reinhard Zumkeller, Jun 23 2015
    
  • Magma
    [Binomial(3*n,n)/(2*n+1): n in [0..30]]; // Vincenzo Librandi, Sep 04 2014
    
  • Maple
    A001764 := n->binomial(3*n,n)/(2*n+1): seq(A001764(n), n=0..25);
    with(combstruct): BB:=[T,{T=Prod(Z,F),F=Sequence(B),B=Prod(F,Z,F)}, unlabeled]:seq(count(BB,size=i),i=0..22); # Zerinvary Lajos, Apr 22 2007
    with(combstruct):BB:=[S, {B = Prod(S,S,Z), S = Sequence(B)}, labelled]: seq(count(BB, size=n)/n!, n=0..21); # Zerinvary Lajos, Apr 25 2008
    n:=30:G:=series(RootOf(g = 1+x*g^3, g),x=0,n+1):seq(coeff(G,x,k),k=0..n); # Robert FERREOL, Apr 03 2015
    alias(PS=ListTools:-PartialSums): A001764List := proc(m) local A, P, n;
    A := [1,1]; P := [1]; for n from 1 to m - 2 do P := PS(PS([op(P), P[-1]]));
    A := [op(A), P[-1]] od; A end: A001764List(25); # Peter Luschny, Mar 26 2022
  • Mathematica
    InverseSeries[Series[y-y^3, {y, 0, 24}], x] (* then a(n)=y(2n+1)=ways to place non-crossing diagonals in convex (2n+4)-gon so as to create only quadrilateral tiles *) (* Len Smiley, Apr 08 2000 *)
    Table[Binomial[3n,n]/(2n+1),{n,0,25}] (* Harvey P. Dale, Jul 24 2011 *)
  • PARI
    {a(n) = if( n<0, 0, (3*n)! / n! / (2*n + 1)!)};
    
  • PARI
    {a(n) = if( n<0, 0, polcoeff( serreverse( x - x^3 + O(x^(2*n + 2))), 2*n + 1))};
    
  • PARI
    {a(n) = my(A); if( n<0, 0, A = 1 + O(x); for( m=1, n, A = 1 + x * A^3); polcoeff(A, n))};
    
  • PARI
    b=vector(22);b[1]=1;for(n=2,22,for(i=1,n-1,for(j=1,n-1,for(k=1,n-1,if((i-1)+(j-1)+(k-1)-(n-2),NULL,b[n]=b[n]+b[i]*b[j]*b[k])))));a(n)=b[n+1]; print1(a(0));for(n=1,21,print1(", ",a(n))) \\ Gerald McGarvey, Oct 08 2008
    
  • PARI
    Vec(1 + serreverse(x / (1+x)^3 + O(x^30))) \\ Gheorghe Coserea, Aug 05 2015
    
  • Python
    from math import comb
    def A001764(n): return comb(3*n,n)//(2*n+1) # Chai Wah Wu, Nov 10 2022
  • Sage
    def A001764_list(n) :
        D = [0]*(n+1); D[1] = 1
        R = []; b = false; h = 1
        for i in range(2*n) :
            for k in (1..h) : D[k] += D[k-1]
            if not b : R.append(D[h])
            else : h += 1
            b = not b
        return R
    A001764_list(22) # Peter Luschny, May 03 2012
    

Formula

From Karol A. Penson, Nov 08 2001: (Start)
G.f.: (2/sqrt(3*x))*sin((1/3)*arcsin(sqrt(27*x/4))).
E.g.f.: hypergeom([1/3, 2/3], [1, 3/2], 27/4*x).
Integral representation as n-th moment of a positive function on [0, 27/4]: a(n) = Integral_{x=0..27/4} (x^n*((1/12) * 3^(1/2) * 2^(1/3) * (2^(1/3)*(27 + 3 * sqrt(81 - 12*x))^(2/3) - 6 * x^(1/3))/(Pi * x^(2/3)*(27 + 3 * sqrt(81 - 12*x))^(1/3)))), n >= 0. This representation is unique. (End)
G.f. A(x) satisfies A(x) = 1+x*A(x)^3 = 1/(1-x*A(x)^2) [Cyvin (1998)]. - Ralf Stephan, Jun 30 2003
a(n) = n-th coefficient in expansion of power series P(n), where P(0) = 1, P(k+1) = 1/(1 - x*P(k)^2).
G.f. Rev(x/c(x))/x, where c(x) is the g.f. of A000108 (Rev=reversion of). - Paul Barry, Mar 26 2010
From Gary W. Adamson, Jul 07 2011: (Start)
Let M = the production matrix:
1, 1
2, 2, 1
3, 3, 2, 1
4, 4, 3, 2, 1
5, 5, 4, 3, 2, 1
...
a(n) = upper left term in M^n. Top row terms of M^n = (n+1)-th row of triangle A143603, with top row sums generating A006013: (1, 2, 7, 30, 143, 728, ...). (End)
Recurrence: a(0)=1; a(n) = Sum_{i=0..n-1, j=0..n-1-i} a(i)a(j)a(n-1-i-j) for n >= 1 (counts ternary trees by subtrees of the root). - David Callan, Nov 21 2011
G.f.: 1 + 6*x/(Q(0) - 6*x); Q(k) = 3*x*(3*k + 1)*(3*k + 2) + 2*(2*(k^2) + 5*k +3) - 6*x*(2*(k^2) + 5*k + 3)*(3*k + 4)*(3*k + 5)/Q(k+1); (continued fraction). - Sergei N. Gladkovskii, Nov 27 2011
D-finite with recurrence: 2*n*(2n+1)*a(n) - 3*(3n-1)*(3n-2)*a(n-1) = 0. - R. J. Mathar, Dec 14 2011
REVERT transform of A115140. BINOMIAL transform is A188687. SUMADJ transform of A188678. HANKEL transform is A051255. INVERT transform of A023053. INVERT transform is A098746. - Michael Somos, Apr 07 2012
(n + 1) * a(n) = A174687(n).
G.f.: F([2/3,4/3], [3/2], 27/4*x) / F([2/3,1/3], [1/2], (27/4)*x) where F() is the hypergeometric function. - Joerg Arndt, Sep 01 2012
a(n) = binomial(3*n+1, n)/(3*n+1) = A062993(n+1,1). - Robert FERREOL, Apr 03 2015
a(n) = A258708(2*n,n) for n > 0. - Reinhard Zumkeller, Jun 23 2015
0 = a(n)*(-3188646*a(n+2) + 20312856*a(n+3) - 11379609*a(n+4) + 1437501*a(n+5)) + a(n+1)*(177147*a(n+2) - 2247831*a(n+3) + 1638648*a(n+4) - 238604*a(n+5)) + a(n+2)*(243*a(n+2) + 31497*a(n+3) - 43732*a(n+4) + 8288*a(n+5)) for all integer n. - Michael Somos, Jun 03 2016
a(n) ~ 3^(3*n + 1/2)/(sqrt(Pi)*4^(n+1)*n^(3/2)). - Ilya Gutkovskiy, Nov 21 2016
Given g.f. A(x), then A(1/8) = -1 + sqrt(5), A(2/27) = (-1 + sqrt(3))*3/2, A(4/27) = 3/2, A(3/64) = -2 + 2*sqrt(7/3), A(5/64) = (-1 + sqrt(5))*2/sqrt(5), etc. A(n^2/(n+1)^3) = (n+1)/n if n > 1. - Michael Somos, Jul 17 2018
From Peter Bala, Sep 14 2021: (Start)
A(x) = exp( Sum_{n >= 1} (1/3)*binomial(3*n,n)*x^n/n ).
The sequence defined by b(n) := [x^n] A(x)^n = A224274(n) for n >= 1 and satisfies the congruence b(p) == b(1) (mod p^3) for prime p >= 3. Cf. A060941. (End)
G.f.: 1/sqrt(B(x)+(1-6*x)/(9*B(x))+1/3), with B(x):=((27*x^2-18*x+2)/54-(x*sqrt((-(4-27*x))*x))/(2*3^(3/2)))^(1/3). - Vladimir Kruchinin, Sep 28 2021
x*A'(x)/A(x) = (A(x) - 1)/(- 2*A(x) + 3) = x + 5*x^2 + 28*x^3 + 165*x^4 + ... is the o.g.f. of A025174. Cf. A002293 - A002296. - Peter Bala, Feb 04 2022
a(n) = hypergeom([1 - n, -2*n], [2], 1). Row sums of A108767. - Peter Bala, Aug 30 2023
G.f.: z*exp(3*z*hypergeom([1, 1, 4/3, 5/3], [3/2, 2, 2], (27*z)/4)) + 1.
- Karol A. Penson, Dec 19 2023
G.f.: hypergeometric([1/3, 2/3], [3/2], (3^3/2^2)*x). See the e.g.f. above. - Wolfdieter Lang, Feb 04 2024
a(n) = (3*n)! / (n!*(2*n+1)!). - Allan Bickle, Feb 20 2024
Sum_{n >= 0} a(n)*x^n/(1 + x)^(3*n+1) = 1. See A316371 and A346627. - Peter Bala, Jun 02 2024
G.f. A(x) satisfies A(x) = 1/A(-x*A(x)^5). - Seiichi Manyama, Jun 16 2025

A006013 a(n) = binomial(3*n+1,n)/(n+1).

Original entry on oeis.org

1, 2, 7, 30, 143, 728, 3876, 21318, 120175, 690690, 4032015, 23841480, 142498692, 859515920, 5225264024, 31983672534, 196947587823, 1219199353190, 7583142491925, 47365474641870, 296983176369495, 1868545312633440, 11793499763070480
Offset: 0

Views

Author

Keywords

Comments

Enumerates pairs of ternary trees [Knuth, 2014]. - N. J. A. Sloane, Dec 09 2014
G.f. (offset 1) is series reversion of x - 2x^2 + x^3.
Hankel transform is A005156(n+1). - Paul Barry, Jan 20 2007
a(n) = number of ways to connect 2*n - 2 points labeled 1, 2, ..., 2*n-2 in a line with 0 or more noncrossing arcs above the line such that each maximal contiguous sequence of isolated points has even length. For example, with arcs separated by dashes, a(3) = 7 counts {} (no arcs), 12, 14, 23, 34, 12-34, 14-23. It does not count 13 because 2 is an isolated point. - David Callan, Sep 18 2007
In my 2003 paper I introduced L-algebras. These are K-vector spaces equipped with two binary operations > and < satisfying (x > y) < z = x > (y < z). In my arXiv paper math-ph/0709.3453 I show that the free L-algebra on one generator is related to symmetric ternary trees with odd degrees, so the dimensions of the homogeneous components are 1, 2, 7, 30, 143, .... These L-algebras are closely related to the so-called triplicial-algebras, 3 associative operations and 3 relations whose free object is related to even trees. - Philippe Leroux (ph_ler_math(AT)yahoo.com), Oct 05 2007
a(n-1) is also the number of projective dependency trees with n nodes. - Marco Kuhlmann (marco.kuhlmann(AT)lingfil.uu.se), Apr 06 2010
Number of subpartitions of [1^2, 2^2, ..., n^2]. - Franklin T. Adams-Watters, Apr 13 2011
a(n) = sum of (n+1)-th row terms of triangle A143603. - Gary W. Adamson, Jul 07 2011
Also the number of Dyck n-paths with up steps colored in two ways (N or A) and avoiding NA. The 7 Dyck 2-paths are NDND, ADND, NDAD, ADAD, NNDD, ANDD, and AADD. - David Scambler, Jun 24 2013
a(n) is also the number of permutations avoiding 213 in the classical sense which can be realized as labels on an increasing strict binary tree with 2n-1 nodes. See A245904 for more information on increasing strict binary trees. - Manda Riehl Aug 07 2014
With offset 1, a(n) is the number of ordered trees (A000108) with n non-leaf vertices and n leaf vertices such that every non-leaf vertex has a leaf child (and hence exactly one leaf child). - David Callan, Aug 20 2014
a(n) is the number of paths in the plane with unit east and north steps, never going above the line x=2y, from (0,0) to (2n+1,n). - Ira M. Gessel, Jan 04 2018
a(n) is the number of words on the alphabet [n+1] that avoid the patterns 231 and 221 and contain exactly one 1 and exactly two occurrences of every other letter. - Colin Defant, Sep 26 2018
a(n) is the number of Motzkin paths of length 3n with n of each type of step, such that (1, 1) and (1, 0) alternate (ignoring (-1, 1) steps). All paths start with a (1, 1) step. - Helmut Prodinger, Apr 08 2019
Hankel transform omitting a(0) is A051255(n+1). - Michael Somos, May 15 2022
If f(x) is the generating function for (-1)^n*a(n), a real solution of the equation y^3 - y^2 - x = 0 is given by y = 1 + x*f(x). In particular 1 + f(1) is Narayana's cow constant, A092526, aka the "supergolden" ratio. - R. James Evans, Aug 09 2023
This is instance k = 2 of the family {c(k, n+1)}A130564.%20_Wolfdieter%20Lang">{n>=0} described in A130564. _Wolfdieter Lang, Feb 04 2024
Also the number of quadrangulations of a (2n+4)-gon which do not contain any diagonals incident to a fixed vertex. - Esther Banaian, Mar 12 2025

Examples

			a(3) = 30 since the top row of Q^3 = (12, 12, 5, 1).
G.f. = 1 + 2*x + 7*x^2 + 30*x^3 + 143*x^4 + 728*x^5 + 3876*x^6 + 21318*x^7 + ... - _Michael Somos_, May 15 2022
		

References

  • N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).

Crossrefs

These are the odd indices of A047749.
Cf. A305574 (the same with offset 1 and the initial 1 replaced with 5).
Cf. A130564 (comment on c(k, n+1)).

Programs

  • Haskell
    a006013 n = a007318 (3 * n + 1) n `div` (n + 1)
    a006013' n = a258708 (2 * n + 1) n
    -- Reinhard Zumkeller, Jun 22 2015
    
  • Magma
    [Binomial(3*n+1,n)/(n+1): n in [0..30]]; // Vincenzo Librandi, Mar 29 2015
    
  • Mathematica
    Binomial[3#+1,#]/(#+1)&/@Range[0,30]  (* Harvey P. Dale, Mar 16 2011 *)
  • PARI
    A006013(n) = binomial(3*n+1,n)/(n+1) \\ M. F. Hasler, Jan 08 2024
    
  • Python
    from math import comb
    def A006013(n): return comb(3*n+1,n)//(n+1) # Chai Wah Wu, Jul 30 2022
  • Sage
    def A006013_list(n) :
        D = [0]*(n+1); D[1] = 1
        R = []; b = false; h = 1
        for i in range(2*n) :
            for k in (1..h) : D[k] += D[k-1]
            if b : R.append(D[h]); h += 1
            b = not b
        return R
    A006013_list(23) # Peter Luschny, May 03 2012
    

Formula

G.f. is square of g.f. for ternary trees, A001764 [Knuth, 2014]. - N. J. A. Sloane, Dec 09 2014
Convolution of A001764 with itself: 2*C(3*n + 2, n)/(3*n + 2), or C(3*n + 2, n+1)/(3*n + 2).
G.f.: (4/(3*x)) * sin((1/3)*arcsin(sqrt(27*x/4)))^2.
G.f.: A(x)/x with A(x)=x/(1-A(x))^2. - Vladimir Kruchinin, Dec 26 2010
From Gary W. Adamson, Jul 14 2011: (Start)
a(n) is the top left term in M^n, where M is the infinite square production matrix:
2, 1, 0, 0, 0, ...
3, 2, 1, 0, 0, ...
4, 3, 2, 1, 0, ...
5, 4, 3, 2, 1, ...
... (End)
From Gary W. Adamson, Aug 11 2011: (Start)
a(n) is the sum of top row terms in Q^n, where Q is the infinite square production matrix as follows:
1, 1, 0, 0, 0, ...
2, 2, 1, 0, 0, ...
3, 3, 2, 1, 0, ...
4, 4, 3, 2, 1, ...
... (End)
D-finite with recurrence: 2*(n+1)*(2n+1)*a(n) = 3*(3n-1)*(3n+1)*a(n-1). - R. J. Mathar, Dec 17 2011
a(n) = 2*A236194(n)/n for n > 0. - Bruno Berselli, Jan 20 2014
a(n) = A258708(2*n+1, n). - Reinhard Zumkeller, Jun 22 2015
From Ilya Gutkovskiy, Dec 29 2016: (Start)
E.g.f.: 2F2([2/3, 4/3]; [3/2,2]; 27*x/4).
a(n) ~ 3^(3*n+3/2)/(sqrt(Pi)*4^(n+1)*n^(3/2)). (End)
a(n) = A110616(n+1, 1). - Ira M. Gessel, Jan 04 2018
0 = v0*(+98415*v2 -122472*v3 +32340*v4) +v1*(+444*v3 -2968*v4) +v2*(-60*v2 +56*v3 +64*v4) where v0=a(n)^2, v1=a(n)*a(n+1), v2=a(n+1)^2, v3=a(n+1)*a(n+2), v4=a(n+2)^2 for all n in Z if a(-1)=-2/3 and a(n)=0 for n<-1. - Michael Somos, May 15 2022
a(n) = (1/4^n) * Product_{1 <= i <= j <= 2*n} (2*i + j + 2)/(2*i + j - 1). Cf. A000260. - Peter Bala, Feb 21 2023
From Karol A. Penson, Jun 02 2023: (Start)
a(n) = Integral_{x=0..27/4} x^n*W(x) dx, where
W(x) = (((9 + sqrt(81 - 12*x))^(2/3) - (9 - sqrt(81 - 12*x))^(2/3)) * 2^(1/3) * 3^(1/6)) / (12 * Pi * x^(1/3)), for x in (0, 27/4).
This integral representation is unique as W(x) is the solution of the Hausdorff power moment problem. Using only the definition of a(n), W(x) can be proven to be positive. W(x) is singular at x = 0, with the singularity x^(-1/3), and for x > 0 is monotonically decreasing to zero at x = 27/4. At x = 27/4 the first derivative of W(x) is infinite. (End)
G.f.: hypergeometric([2/3,1,4/3], [3/2,2], (3^3/2^2)*x). See the e.g.f. above. - Wolfdieter Lang, Feb 04 2024
a(n) = A024485(n+1)/3. - Michael Somos, Oct 14 2024
G.f.: (Sum_{n >= 0} binomial(3*n+2, n)*x^n) / (Sum_{n >= 0} binomial(3*n, n)*x^n) = (B(x) - 1)/(x*B(x)), where B(x) = Sum_{n >= 0} binomial(3*n, n)/(2*n+1) * x^n is the g.f. of A001764. - Peter Bala, Dec 13 2024
The g.f. A(x) is uniquely determined by the conditions A(0) = 1 and [x^n] A(x)^(-n) = -2 for all n >= 1. Cf. A006632. - Peter Bala, Jul 24 2025

Extensions

Edited by N. J. A. Sloane, Feb 21 2008

A005809 a(n) = binomial(3n,n).

Original entry on oeis.org

1, 3, 15, 84, 495, 3003, 18564, 116280, 735471, 4686825, 30045015, 193536720, 1251677700, 8122425444, 52860229080, 344867425584, 2254848913647, 14771069086725, 96926348578605, 636983969321700, 4191844505805495, 27619435402363035
Offset: 0

Views

Author

Keywords

Comments

Number of paths in Z X Z starting at (0,0) and ending at (3n,0) using steps in {(1,1),(1,-2)}.
Number of even trees with 2n edges and one distinguished vertex. Even trees are rooted plane trees where every vertex (including root) has even degree.
Hankel transform is 3^n*A051255(n), where A051255 is the Hankel transform of C(3n,n)/(2n+1). - Paul Barry, Jan 21 2007
a(n) is the number of stack polyominoes inscribed in an (n+1) X (n+1) box. Equivalently, a(n) is the number of unimodal compositions with n+1 parts in which the maximum value of the parts is n+1. For instance, for n = 2, we have the following compositions: (3,3,3), (2,3,3), (1,3,3), (3,3,1), (3,3,2), (2,2,3), (1,2,3), (2,3,1), (1,1,3), (1,3,1), (3,1,1), (2,3,2), (1,3,2), (3,2,1), (3,2,2). - Emanuele Munarini, Apr 07 2011
Conjecture: a(n)==3 (mod n^3) iff n is an odd prime. - Gary Detlefs, Mar 23 2013. The congruence a(p) = binomial(3*p,p) = 3 (mod p^3) for odd prime p is a known generalization of Wolstenholme's theorem. See Mestrovic, Section 6, equation 35. - Peter Bala, Dec 28 2014
In general, C(k*n,n) = C(k*n-1,n-1)*C((k*n)^2,2)/(3*n*C(k*n+1,3)), n>0. - Gary Detlefs, Jan 02 2014

Examples

			G.f. = 1 + 3*x + 15*x^2 + 84*x^3 + 495*x^4 + 3003*x^5 + 18564*x^6 + ... - _Michael Somos_, Jan 30 2019
		

References

  • M. Abramowitz and I. A. Stegun, eds., Handbook of Mathematical Functions, National Bureau of Standards Applied Math. Series 55, 1964 (and various reprintings), p. 828.
  • N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).

Crossrefs

binomial(k*n,n): A000984 (k = 2), A005810 (k = 4), A001449 (k = 5), A004355 (k = 6), A004368 (k = 7), A004381 (k = 8), A169958 - A169961 (k = 9 thru 12).

Programs

  • Haskell
    a005809 n = a007318 (3*n) n  -- Reinhard Zumkeller, May 06 2012
    
  • Magma
    [ Binomial(3*n,n): n in [0..150] ]; // Vincenzo Librandi, Apr 21 2011
    
  • Maple
    A005809:=n->binomial(3*n,n); seq(A005809(n), n=0..40); # Wesley Ivan Hurt, Mar 21 2014
  • Mathematica
    R[ z_ ] := ((2-18*z + 27*z^2 + 3^(3/2)*z^(3/2)*(27*z-4)^(1/2))/2)^(1/3); f[ z_ ] := ( (R[ z ])^3 + (1-3*z)*(R[ z ])^2 + (1-6*z)*R[ z ] )/( (R[ z ])^4 + (1-6*z)*(R[ z ])^2 + (6*z-1)^2 )
    Table[Binomial[3*n,n],{n,0,40}] (* Vladimir Joseph Stephan Orlovsky, Mar 03 2011 *)
  • Maxima
    makelist(binomial(3*n,n),n,0,100); /* Emanuele Munarini, Apr 07 2011 */
    
  • Maxima
    B(x):=(2/sqrt(3*x))*sin((1/3)*asin(sqrt(27*x/4)))-1;
    taylor(x*diff(B(x),x)/B(x),x,0,10); /* Vladimir Kruchinin, Oct 02 2015 */
  • PARI
    a(n)=binomial(3*n,n) \\ Charles R Greathouse IV, Nov 20 2012
    
  • Sage
    [binomial(3*n,n) for n in range(0, 22)] # Zerinvary Lajos, Dec 16 2009
    

Formula

The g.f. R[ z_ ] below (in the Mathematica field) was found by Kurt Persson (kurt(AT)math.chalmers.se) and communicated by Einar Steingrimsson (einar(AT)math.chalmers.se).
Using Stirling's formula in A000142, it is easy to get the asymptotic expression a(n) ~ (1/2) * (27/4)^n / sqrt(Pi*n / 3). - Dan Fux (dan.fux(AT)OpenGaia.com or danfux(AT)OpenGaia.com), Apr 07 2001
a(n) = Sum_{k=0..n} C(n, k)*C(2n, k). - Paul Barry, May 15 2003
G.f.: 1/(1-3zg^2), where g=g(z) is given by g=1+zg^3, g(0)=1, i.e., (in Maple notation) g := 2*sin(arcsin(3*sqrt(3*z)/2)/3)/sqrt(3*z). - Emeric Deutsch, May 22 2003
G.f.: x*B'(x)/B(x), where B(x)+1 is the g.f. for A001764. - Vladimir Kruchinin, Oct 02 2015
a(n) ~ (1/2)*3^(1/2)*Pi^(-1/2)*n^(-1/2)*2^(-2*n)*3^(3*n)*(1 - 7/72*n^-1 + 49/10368*n^-2 + 6425/2239488*n^-3 - ...). - Joe Keane (jgk(AT)jgk.org), Nov 07 2003
a(n) = A006480(n)/A000984(n). - Lior Manor, May 04 2004
a(n) = Sum_{i_1=0..n, i_2=0..n} binomial(n, i_1)*binomial(n, i_2)*binomial(n, i_1+i_2). - Benoit Cloitre, Oct 14 2004
a(n) = Sum_{k=0..n} A109971(k)*3^k; a(0)=1, a(n) = Sum_{k=0..n} 3^k*C(3n-k,n-k)2k/(3n-k), n>0. - Paul Barry, Jan 21 2007
a(n) = A085478(2n,n). - Philippe Deléham, Sep 17 2009
E.g.f.: 2F2(1/3,2/3;1/2,1;27*x/4), where F(a1,a2;b1,b2;z) is a hypergeometric series. - Emanuele Munarini, Apr 12 2011
a(n) = Sum_{k=0..n} binomial(2*n+k-1,k). - Arkadiusz Wesolowski, Apr 02 2012
G.f.: cos((1/3)*asin(sqrt(27x/4)))/sqrt(1-27x/4). - Tom Copeland, May 24 2012
G.f.: A(x) = 1 + 6*x/(G(0)-6*x) where G(k) = (2*k+2)*(2*k+1) + 3*x*(3*k+1)*(3*k+2) - 6*x*(k+1)*(2*k+1)*(3*k+4)*(3*k+5)/G(k+1); (continued fraction). - Sergei N. Gladkovskii, Jun 30 2012
D-finite with recurrence: 2*n*(2*n-1)*a(n) - 3*(3*n-1)*(3*n-2)*a(n-1) = 0. - R. J. Mathar, Feb 05 2013
a(n) = (2n+1)*A001764(n). - Johannes W. Meijer, Aug 22 2013
a(n) = C(3*n-1,n-1)*C(9*n^2,2)/(3*n*C(3*n+1,3)), n>0. - Gary Detlefs, Jan 02 2014
a(n) = [x^n] 1/(1 - x)^(2*n+1). - Ilya Gutkovskiy, Oct 03 2017
a(n) = hypergeom([-2*n, -n], [1], 1). - Peter Luschny, Mar 19 2018
a(n) = Sum_{k=0..n} binomial(n, k) * binomial(2*n, n-k) = row sums of A110608. - Michael Somos, Jan 30 2019
0 = a(n)*(-3188646*a(n+2) +7322076*a(n+3) -2805111*a(n+4) +273585*a(n+5)) +a(n+1)*(+413343*a(n+2) -1252017*a(n+3) +538344*a(n+4) -55940*a(n+5)) +a(n+2)*(-4131*a(n+2) +38733*a(n+3) -21628*a(n+4) +2528*a(n+5)) for all n in Z. - Michael Somos, Jan 30 2019
Sum_{n>=1} 1/a(n) = A229705. - Amiram Eldar, Nov 14 2020
From Peter Bala, Feb 20 2022: (Start)
The o.g.f. A(x) satisfies the differential equation (4*x - 27*x^2)*A''(x) + (2 - 54*x)*A'(x) - 6*A(x) = 0, with A(0) = 1 and A'(0) = 3.
Algebraic equation: (1 - A(x))*(1 + 2*A(x))^2 + 27*x*A(x)^3 = 0.
Sum_{n >= 1} a(n)*( x*(2*x + 3)^2/(27*(1 + x)^3) )^n = x. (End)
From Vaclav Kotesovec, May 13 2022: (Start)
Sum_{n>=0} a(n) / 3^(2*n) = 2*cos(Pi/9).
Sum_{n>=0} a(n) / (27/2)^n = (1 + sqrt(3))/2.
Sum_{n>=0} a(n) / 3^(3*n) = 2*cos(Pi/18) / sqrt(3).
In general, for k > 27/4, Sum_{n>=0} a(n)/k^n = 2*cos(arccos(1 - 27/(2*k))/6) / sqrt(4 - 27/k). (End)
G.f.: hypergeom([1/3, 2/3], [1/2], 27*z/4), the Gauss hypergeometric function 2F1. - Karol A. Penson, Dec 12 2023
a(n) = 1/4^n * Sum_{k = n..3*n} binomial(k, n)*binomial(3*n, k). - Peter Bala, Jun 29 2025
a(n) = Sum_{k = 0..n} (-1)^(n+k)*binomial(4*n+1, k)*binomial(2*n-k, n-k). - Peter Bala, Sep 04 2025

A047098 a(n) = 2*binomial(3*n, n) - Sum_{k=0..n} binomial(3*n, k).

Original entry on oeis.org

1, 2, 8, 38, 196, 1062, 5948, 34120, 199316, 1181126, 7080928, 42860534, 261542752, 1607076200, 9934255472, 61732449648, 385393229460, 2415935640198, 15200964233864, 95962904716402, 607640599286276, 3858198001960438, 24559243585545644, 156692889782067712
Offset: 0

Views

Author

Clark Kimberling, Aug 15 1998

Keywords

Comments

T(2n,n), array T as in A047089. [Corrected Dec 08 2006]
Let B_3^+ denote the semigroup with presentation . Let D=aba be the 'fundamental word'. Then this sequence is also equal to the number of words in B_3^+ equal in B_3^+ to D^n, n >= 0. - Stephen P. Humphries, Jan 20 2004
In the language of Riordan arrays, row sums of (1/(1+x), x/(1+x)^3)^-1, where (1/(1+x), x/(1+x)^3) has general term (-1)^(n-k)*binomial(n+2k, 3k). - Paul Barry, May 09 2005
Hankel transform is 2^n*A051255(n) where A051255 is the Hankel transform of C(3n,n)/(2n+1). - Paul Barry, Jan 21 2007

Crossrefs

Column k=2 of A213028.

Programs

  • Maple
    A047098 := n -> 2*binomial(3*n, n)-add(binomial(3*n, k), k=0..n);
  • Mathematica
    Table[2Binomial[3n,n]-Sum[Binomial[3n,k],{k,0,n}],{n,0,35}] (* Harvey P. Dale, Jul 27 2011 *)
  • PARI
    a(n)=if(n<0,0,polcoeff((((1+10*x-2*x^2)+(1-4*x)*sqrt(1-4*x+x*O(x^n)))/2)^n,n))
    
  • PARI
    a(n)=if(n<0,0, 2*binomial(3*n,n)-sum(k=0,n,binomial(3*n,k)))

Formula

G.f. A(x)=y satisfies (8x-1)y^3-y^2+y+1=0. - Michael Somos, Jan 28 2004
Coefficient of x^n in ((1+10x-2x^2+(1-4x)^(3/2))/2)^n. - Michael Somos, Sep 25 2003
a(n) = Sum_{k = 0..n} A109971(k)*2^k; a(0) = 1, a(n) = Sum_{k = 0..n} 2^k*C(3n-k,n-k)*2*k/(3*n-k), n > 0. - Paul Barry, Jan 21 2007
Conjecture: 2*n*(2*n-1)*a(n) +(-71*n^2+112*n-48)*a(n-1) +3*(131*n^2-391*n+296)*a(n-2) -72*(3*n-7)*(3*n-8)*a(n-3)=0. - R. J. Mathar, Nov 30 2012
a(n) = A321957(n) + 2*binomial(3*n, n) - 8^n. - Peter Luschny, Nov 22 2018
The Gauss congruences a(n*p^k) == a(n*p^(k-1)) (mod p^k) hold for all primes p and positive integers n and k. - Peter Bala, Mar 05 2022

Extensions

Clark Kimberling, Dec 08 2006, changed "T(3n,2n)" to "T(2n,n)" in the comment line, but observes that some of the other comments seem to apply to the sequence T(3n,2n) rather than to the sequence T(2n,n).
Edited by N. J. A. Sloane, Dec 21 2006, replacing the old definition in terms of A047089 by an explicit formula supplied by Benoit Cloitre, Oct 25 2003.

A005161 Number of alternating sign 2n+1 X 2n+1 matrices symmetric with respect to both horizontal and vertical axes (VHSASM's).

Original entry on oeis.org

1, 1, 1, 2, 6, 33, 286, 4420, 109820, 4799134, 340879665, 42235307100, 8564558139000, 3012862604463000, 1742901718473961200, 1742218029490675762080, 2873822682985675809192288, 8167157387273280570395662320, 38402596062535617548517706584760, 310388509293255836481583597538626504
Offset: 0

Views

Author

Keywords

References

  • N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).
  • R. P. Stanley, A baker's dozen of conjectures concerning plane partitions, pp. 285-293 of "Combinatoire Enumerative (Montreal 1985)", Lect. Notes Math. 1234, 1986.

Crossrefs

Programs

  • PARI
    \\ here b(n) and c(n) are A005156 and A051255.
    b(n) = prod(k=0, n-1, (3*k+2)*(6*k+3)!*(2*k+1)!/((4*k+2)!*(4*k+3)!));
    c(n) = prod(k=0, n-1, (3*k+1)*(6*k)!*(2*k)!/((4*k)!*(4*k+1)!));
    a(n) = b(n\2) * c((n+1)\2) \\ Andrew Howroyd, May 09 2023

Formula

Robbins gives a simple (conjectured) formula, which was proven by Okada.
a(2*n) = A005156(n)*A051255(n); a(2*n+1) = A005156(n)*A051255(n+1). - Paul Zinn-Justin, May 05 2023
a(n) = A005156(floor(n/2)) * A051255(ceiling(n/2)). - Andrew Howroyd, May 09 2023

Extensions

More terms (from the P. Pyatov paper) from Vladeta Jovovic, Aug 15 2008
Terms a(13) and beyond from Andrew Howroyd, May 09 2023

A181119 Number of transpose-complementary plane partitions of n.

Original entry on oeis.org

1, 2, 84, 81796, 1844536720, 962310111888300, 11608208114358751650000, 3236574482779383546336417240000, 20853456581643133066208521560263633137920, 3104385823530881109001458753652585998600603921849920, 10676554307318599842868990948461304923921623250562199975300214736
Offset: 0

Views

Author

Arvind Ayyer, Jan 21 2011

Keywords

Comments

The complement of a plane partition inside an m X m X m cube consists of the boxes which are within the cube, but not in the plane partition, rotated in an appropriate way.
a(n) is the number of plane partitions inside an 2n X 2n X 2n cube whose (matrix) transpose when written as an 2n X 2n array is the same as its complement.

Examples

			When n=2, there are two transpose-complementary plane partitions,
[1 1] and [2 1], both of whose transpose and complement is equal to themselves.
[1 1]     [1 0]
		

Crossrefs

Programs

  • Mathematica
    Table[Binomial[3n-1,n]Product[(2n+i+j+1)/(i+j+1),{i,1,2n-2}, {j,i,2n-2}], {n,0,10}] (* Harvey P. Dale, Jan 27 2012 *)
  • PARI
    a(n) = binomial(3*n-1,n)*prod(i=1,2*n-2,prod(j=i,2*n-2,(2*n+i+j+1)/(i+j+1))); \\ Michel Marcus, Jun 18 2015

Formula

a(n) = binomial(3n-1,n)*Product(i=1..2n-2,Product(j=i..2n-2,(2n+i+j+1)/(i+j+1))).
a(n) ~ exp(1/24) * 3^(9*n^2 - 3*n/2 - 1/24) / (sqrt(A) * n^(1/24) * 2^(12*n^2 - n - 1/3)), where A = A074962 = 1.2824271291... is the Glaisher-Kinkelin constant. - Vaclav Kotesovec, Feb 28 2015

A049504 a(n) = Product_{i = 0..n-1} ((3*i+1)!*(6*i)!*(2*i)!)/((4*i)!*(4*i+1)!).

Original entry on oeis.org

1, 1, 12, 47520, 266499072000, 5578457158440714240000, 903833169262981594760400076800000000, 2035652583056655211566004660439314466655436800000000000, 103962610930356904475854868257296244089884364267142052118842572800000000000000
Offset: 0

Views

Author

Keywords

Comments

Given in first printing of Bressoud book as number of cyclically symmetric transpose complement plane partitions. For correct version see A051255.

References

  • D. M. Bressoud, Proofs and Confirmations, Camb. Univ. Press, 1999; Eq. (6.15), p. 199.

Programs

  • Maple
    a := proc(n) local i; mul((3*i+1)!*(6*i)!*(2*i)!/((4*i)!*(4*i+1)!),i = 0..n-1); end;
  • Mathematica
    Table[Product[((3i+1)!(6i)!(2i)!)/((4i)!(4i+1)!),{i,0,n-1}],{n,0,10}] (* Harvey P. Dale, Apr 25 2016 *)

Formula

a(n) ~ A^(-1/2) * Gamma(1/3) * 2^(-1/9 + 3*n/2 - 4*n^2) * 3^(-1/24 - 5*n/2 + 9*n^2/2) * exp(1/24 + n - 9*n^2/4) * n^(1/8 - n + 3*n^2/2) * Pi^((n-1)/2), where A = A074962 is the Glaisher-Kinkelin constant. - Vaclav Kotesovec, Apr 25 2016

Extensions

Definition corrected by Harvey P. Dale, Apr 25 2016

A127636 Hankel transform of C(3n,n), A005809.

Original entry on oeis.org

1, 6, 99, 4590, 601749, 223671780, 236051526780, 707876792166456, 6035175245831686020, 146337688589606021286324, 10094030504239191283895044707, 1981045079056219549583047340549844
Offset: 0

Views

Author

Paul Barry, Jan 21 2007

Keywords

Formula

a(n)=3^n*A051255(n).

A127635 Hankel transform of A047098.

Original entry on oeis.org

1, 4, 44, 1360, 118864, 29454720, 20723316480, 41430374667264, 235483137163985920, 3806579106735674587136, 175045931960590896961989632, 22902901668710944230193460535296, 8527272133354589357030560193109508096
Offset: 0

Views

Author

Paul Barry, Jan 21 2007

Keywords

Crossrefs

Formula

a(n) = 2^n*A051255(n).

A133402 Hankel transform of A006588.

Original entry on oeis.org

1, 96, 405504, 77007421440, 661630022502580224, 257876005135691663309537280, 4565900567740737406606787243126292480, 3675506444195600567841408683430769715388692299776
Offset: 0

Views

Author

Paul Barry, Jan 08 2008

Keywords

Formula

a(n)=4^(n(n+1))*A127636(n); a(n)=A000244(n)*A051255(n)*A053765(n+1);
Showing 1-10 of 10 results.