A185650 -id:A185650 - cp's OEIS frontend

A000891 a(n) = (2n)!(2n+1)! / (n! (n+1)!)^2.

1, 3, 20, 175, 1764, 19404, 226512, 2760615, 34763300, 449141836, 5924217936, 79483257308, 1081724803600, 14901311070000, 207426250094400, 2913690606794775, 41255439318353700, 588272005095043500

Offset: 0

N. J. A. Sloane

nonn

Number of parallelogram polyominoes having n+1 columns and n+1 rows. - Emeric Deutsch, May 21 2003
Number of tilings of an  hexagon.
a(n) is the number of non-crossing partitions of [2n+1] into n+1 blocks. For example, a[1] counts 13-2, 1-23, 12-3. - David Callan, Jul 25 2005
The number of returning walks of length 2n on the upper half of a square lattice, since a(n) = Sum_{k=0..2n} binomial(2n,k)*A126120(k)*A126869(n-k). - Andrew V. Sutherland, Mar 24 2008
For sequences counting walks in the upper half-plane starting from the origin and finishing at the lattice points (0,m) see A145600 (m = 1), A145601 (m = 2), A145602 (m = 3) and A145603 (m = 4). - Peter Bala, Oct 14 2008
The number of proper mergings of two n-chains. - Henri Mühle, Aug 17 2012
a(n) is number of pairs of non-intersecting lattice paths from (0,0) to (n+1,n+1) using (1,0) and (0,1) as steps. Here, non-intersecting means two paths do not share a vertex except the origin and the destination. For example, a(1) = 3 because we have three such pairs from (0,0) to (2,2): {NNEE,EENN}, {NNEE,ENEN}, {NENE,EENN}. - Ran Pan, Oct 01 2015
Also the number of ordered rooted trees with 2(n+1) nodes and n+1 leaves, i.e., half of the nodes are leaves. These trees are ranked by A358579. The unordered version is A185650. - Gus Wiseman, Nov 27 2022
The number of secondary GL(2) invariants constructed from n+1 two component vectors. This number was evaluated by using the Molien-Weyl formula to compute the Hilbert series of the ring of invariants. - Jaco van Zyl, Jun 30 2025

			G.f. = 1 + 3*x + 20*x^2 + 175*x^3 + 1764*x^4 + 19404*x^5 + ...
From _Gus Wiseman_, Nov 27 2022: (Start)
The a(2) = 20 ordered rooted trees with 6 nodes and 3 leaves:
  (((o)oo))  (((o)o)o)  (((o))oo)
  (((oo)o))  (((oo))o)  ((o)(o)o)
  (((ooo)))  ((o)(oo))  ((o)o(o))
  ((o(o)o))  ((o(o))o)  (o((o))o)
  ((o(oo)))  ((oo)(o))  (o(o)(o))
  ((oo(o)))  (o((o)o))  (oo((o)))
             (o((oo)))
             (o(o(o)))
(End)

J. M. Borwein and P. B. Borwein, Pi and the AGM, Wiley, 1987, p. 8.
E. R. Hansen, A Table of Series and Products, Prentice-Hall, Englewood Cliffs, NJ, 1975, p. 94.

Vincenzo Librandi, Table of n, a(n) for n = 0..100
Abderrahim Arabi, Hacène Belbachir, and Jean-Philippe Dubernard, Enumeration of parallelogram polycubes, arXiv:2105.00971 [cs.DM], 2021.
Elena Barcucci, Andrea Frosini and Simone Rinaldi, On directed-convex polyominoes in a rectangle, Discr. Math., 298 (2005). 62-78.
Paul Barry, On Integer-Sequence-Based Constructions of Generalized Pascal Triangles, J. Integer Sequ., Vol. 9 (2006), Article 06.2.4.
Paul Barry, On a Generalization of the Narayana Triangle, J. Int. Seq. 14 (2011), Article 11.4.5.
William Y. C. Chen, Sabrina X. M. Pang, Ellen X. Y. Qu, and Richard P. Stanley, Pairs of Noncrossing Free Dyck Paths and Noncrossing Partitions, arXiv:0804.2930 [math.CO], 2008.
William Y. C. Chen, Sabrina X. M. Pang, Ellen X. Y. Qu, and Richard P. Stanley, Pairs of Noncrossing Free Dyck Paths and Noncrossing Partitions, Discrete Math., 309 (2009), 2834-2838.
Robert de Mello Koch, Animik Ghosh, and Hendrik J. R. Van Zyl, Bosonic Fortuity in Vector Models, arXiv:2504.14181 [hep-th], 2025. See p. 9; Journal of High Energy Physics 06 (2025) 246.
Ivan Marin and Emmanuel Wagner, A cubic defining algebra for the Links-Gould polynomial. arXiv preprint arXiv:1203.5981 [math.GT], 2012. - From _N. J. A. Sloane_, Sep 21 2012
Henri Mühle, Counting Proper Mergings of Chains and Antichains, arXiv:1206.3922 [math.CO], 2012.
Guoce Xin, Determinant formulas relating to tableaux of bounded height, Adv. Appl. Math. 45 (2010) 197-211.

Cf. A000356, A010370, A038535.
Cf. A145600, A145601, A145602, A145603. - Peter Bala, Oct 14 2008
Cf. A000142, A000894, A248045.
Equals half of A267981.
Counts the trees ranked by A358579.
A001263 counts ordered rooted trees by nodes and leaves.
A090181 counts ordered rooted trees by nodes and internals.
Cf. A000108, A065097, A080936, A185650, A358371, A358586, A358590.

a000891 n = a001263 (2 * n - 1) n  -- Reinhard Zumkeller, Oct 10 2013

[Factorial(2*n)*Factorial(2*n+1) / (Factorial(n) * Factorial(n+1))^2: n in [0..20]]; // Vincenzo Librandi, Aug 15 2011

with(combstruct): bin := {B=Union(Z,Prod(B,B))} :seq(1/2*binomial(2*i,i)*(count([B,bin,unlabeled],size=i)), i=1..18) ; # Zerinvary Lajos, Jun 06 2007

a[ n_] := If[ n == -1, 0, Binomial[2 n + 1, n]^2 / (2 n + 1)]; (* Michael Somos, May 28 2014 *)
a[ n_] := SeriesCoefficient[ (1 - Hypergeometric2F1[ -1/2, 1/2, 1, 16 x]) / (4 x), {x, 0, n}]; (* Michael Somos, May 28 2014 *)
a[ n_] := If[ n < 0, 0, (2 n)! SeriesCoefficient[ BesselI[0, 2 x] BesselI[1, 2 x] / x, {x, 0, 2 n}]]; (* Michael Somos, May 28 2014 *)
a[ n_] := SeriesCoefficient[ (1 - EllipticE[ 16 x] / (Pi/2)) / (4 x), {x, 0, n}]; (* Michael Somos, Sep 18 2016 *)
a[n_] := (2 n + 1) CatalanNumber[n]^2;
Array[a, 20, 0] (* Peter Luschny, Mar 03 2020 *)

{a(n) = binomial(2*n+1, n)^2 / (2*n + 1)}; /* Michael Somos, Jun 22 2005 */

a(n) = matdet(matrix(n, n, i, j, binomial(n+j+1,i+1))) \\ Hugo Pfoertner, Oct 22 2022

-4*a(n) = A010370(n+1).
G.f.: (1 - E(16*x)/(Pi/2))/(4*x) where E() is the elliptic integral of the second kind. [edited by Olivier Gérard, Feb 16 2011]
G.f.: 3F2(1, 1/2, 3/2; 2,2; 16*x)= (1 - 2F1(-1/2, 1/2; 1; 16*x)) / (4*x). - Olivier Gérard, Feb 16 2011
E.g.f.: Sum_{n>=0} a(n)*x^(2*n)/(2*n)! = BesselI(0, 2*x) * BesselI(1, 2*x) / x. - Michael Somos, Jun 22 2005
a(n) = A001700(n)*A000108(n) = (1/2)*A000984(n+1)*A000108(n). - Zerinvary Lajos, Jun 06 2007
For n > 0, a(n) = (n+2)*A000356(n) starting (1, 5, 35, 294, ...). - Gary W. Adamson, Apr 08 2011
a(n) = A001263(2*n+1,n+1) = binomial(2*n+1,n+1)*binomial(2*n+1,n)/(2*n+1) (central members of odd numbered rows of Narayana triangle).
G.f.: If G_N(x) = 1 + Sum_{k=1..N} ((2*k)!*(2*k+1)!*x^k)/(k!*(k+1)!)^2, G_N(x) = 1 + 12*x/(G(0) - 12*x); G(k) = 16*x*k^2 + 32*x*k + k^2 + 4*k + 12*x + 4 - 4*x*(2*k+3)*(2*k+5)*(k+2)^2/G(k+1); (continued fraction). - Sergei N. Gladkovskii, Nov 24 2011
D-finite with recurrence (n+1)^2*a(n) - 4*(2*n-1)*(2*n+1)*a(n-1) = 0. - R. J. Mathar, Dec 03 2012
a(n) = A005558(2n). - Mark van Hoeij, Aug 20 2014
a(n) = A000894(n) / (n+1) = A248045(n+1) / A000142(n+1). - Reinhard Zumkeller, Sep 30 2014
From Ilya Gutkovskiy, Feb 01 2017: (Start)
E.g.f.: 2F2(1/2,3/2; 2,2; 16*x).
a(n) ~ 2^(4*n+1)/(Pi*n^2). (End)
a(n) = A005408(n)*(A000108(n))^2. - Ivan N. Ianakiev, Nov 13 2019
a(n) = det(M(n)) where M(n) is the n X n matrix with m(i,j) = binomial(n+j+1,i+1). - Benoit Cloitre, Oct 22 2022
a(n) = Integral_{x=0..16} x^n*W(x) dx, where W(x) = (16*EllipticE(1 - x/16) - x*EllipticK(1 - x/16))/(8*Pi^2*sqrt(x)), n=>0. W(x) diverges at x=0, monotonically decreases for x>0, and vanishes at x=16. EllipticE and EllipticK are elliptic functions. This integral representation as n-th moment of a positive function W(x) on the interval [0, 16] is unique. - Karol A. Penson, Dec 20 2023

More terms from Andrew V. Sutherland, Mar 24 2008

A065097 a(n) = ((2n+1) + (2n-1) - 1)!/((2n+1)!*(2n-1)!).

Original entry on oeis.org

1, 1, 7, 66, 715, 8398, 104006, 1337220, 17678835, 238819350, 3282060210, 45741281820, 644952073662, 9183676536076, 131873975875180, 1907493251046152, 27767032438524099, 406472021074865382, 5979899192930226746, 88366931393503350700, 1311063521138246054410

Offset: 0

Len Smiley, Nov 11 2001

A Catalan-like formula using consecutive odd numbers. Recall that Catalan numbers (A000108) are given by ((n+1)+(n)-1)!/((n+1)!(n)!).
From David Callan, Jun 01 2006: (Start)
a(n) = number of Dyck (2n)-paths (i.e., semilength = 2n) all of whose interior returns to ground level (if any) occur at or before the (2n-2)-nd step, that is, they occur strictly before the midpoint of the path.
For example, a(2)=7 counts UUUUDDDD, UUUDUDDD, UUDUUDDD, UUDUDUDD, UUUDDUDD, UD.UUUDDD, UD.UUDUDD ("." denotes an interior return to ground level).
This result follows immediately from an involution on Dyck paths, due to Emeric Deutsch, defined by E->E, UPDQ -> UQDP (where E is the empty Dyck path; U=upstep, D=downstep and P,Q are arbitrary Dyck paths), because the involution is fixed-point-free on Dyck (2n)-paths and contains one path of the type being counted in each orbit.
a(n) = Sum_{k=0..n-1} C(2n-1-2k)*C(2k). This identity has the following combinatorial interpretation:
a(n) is the number of odd-GL-marked Dyck (2n-1)-paths. An odd-GL vertex is a vertex at location (2i,0) for some odd i >= 1 (path starts at origin). An odd-GL-marked Dyck path is a Dyck path with one of its odd-GL vertices marked. For example, a(2)=7 counts UUUDDD*, UUDUDD*, UD*UUDD, UDUUDD*, UD*UDUD, UDUDUD*, UUDDUD* (the * denotes the marked odd-GL vertex). (End)
a(n+1) = Sum_{k=0..n} C(k)*C(2*n+1-k), n >= 0, with C(n) = A000108(n), also gives the odd part of the bisection of the half-convolution of the Catalan sequence A000108 with itself. For the definition of the half-convolution of a sequence with itself see a comment on A201204. There one also finds the rule for the o.g.f. given below in the formula section. The even part of this bisection is found under A201205. - Wolfdieter Lang, Jan 05 2012
From Peter Bala, Dec 01 2015: (Start)
Let x = p/q be a positive rational in reduced form with p,q > 0. Define Cat(x) = (1/(2*p + q))*binomial(2*p + q, p). Then Cat(n) = Catalan(n). This sequence is Cat(n + 1/2) = (1/(4*n + 4))*binomial(4*n + 4, 2*n + 1). Cf. A265101 (Cat(n + 1/3)), A265102 (Cat(n + 1/4)) and A265103 (Cat(n + 1/5)).
Number of maximal faces of the rational associahedron Ass(2*n + 1, 2*n + 3). Number of lattice paths from (0, 0) to (2*n + 3, 2*n + 1) using steps of the form (1, 0) and (0, 1) and staying above the line y = (2*n + 1)/(2*n + 3)*x. See Armstrong et al. (End)
Also the number of ordered rooted trees with 2n nodes, most of which are leaves, i.e., the odd bisection of A358585. This follows from Callan's formula below. - Gus Wiseman, Nov 27 2022

			G.f.: 1 + x + 7*x^2 + 66*x^3 + 715*x^4 + 8398*x^5 + 104006*x^6 + ...

Alois P. Heinz, Table of n, a(n) for n = 0..834 (terms n = 1..100 from Harry J. Smith)
D. Armstrong, B. Rhoades, and N. Williams, Rational associahedra and noncrossing partitions arxiv:1305.7286 [math.CO], 2008.
Emeric Deutsch, An involution on Dyck paths and its consequences, Discrete Math., 204 (1999), no. 1-3, 163-166.

Cf. A003150 (for analog with consecutive Fibonacci numbers).
Equals (1/2) A048990. Cf. A000108, A001795, A201204, A201205, A265101, A265102, A265103.
Cf. A000891, A001263, A001622, A090181, A185650, A358585.

[Binomial(4*n-1, 2*n-1)/(2*n+1): n in [1..20]]; // Vincenzo Librandi, Dec 09 2015

seq(binomial(4*n-1,2*n-1)/(2*n+1), n=0..30); # Robert Israel, Dec 08 2015

a[ n_] := If[ n < 1, 0, Binomial[ 4 n - 1, 2 n - 1] / (2 n + 1)]; (* Michael Somos, Oct 25 2014 *)

combinat::dyckWords::count(2*n)/2 $ n = 1..26 // Zerinvary Lajos, Apr 25 2007

a(n) = { if(n==0, 1, (4*n - 1)!/((2*n + 1)!*(2*n - 1)!)) } \\ Harry J. Smith, Oct 07 2009

vector(20, n, binomial(4*n-1, 2*n-1)/(2*n+1)) \\ Altug Alkan, Dec 08 2015

A065097 = lambda n: hypergeometric([1-2*n,-2*n],[2],1)/2
[Integer(A065097(n).n(500)) for n in (1..20)] # Peter Luschny, Sep 22 2014

a(n) = binomial(4*n-1, 2*n-1)/(2*n+1).
a(n) = C(2n)/2 where C(n) is the Catalan number A000108. - David Callan, Jun 01 2006
G.f.: 1/2 + (sqrt(2)/2)/sqrt(1+sqrt(1-16*x)). - Vladeta Jovovic, Sep 26 2003
G.f.: 1 + 3F2([1, 5/4, 7/4], [2, 5/2], 16*x). - Olivier Gérard, Feb 16 2011
O.g.f.: (1 + (cata(sqrt(x)) + cata(-sqrt(x)))/2)/2, with the o.g.f. cata(x) of the Catalan numbers. See the W. Lang comment above. - Wolfdieter Lang, Jan 05 2012
a(n) = hypergeometric([1-2*n,-2*n],[2],1)/2. - Peter Luschny, Sep 22 2014
a(n) = A001448(n) / (4*n + 2) if n>0. - Michael Somos, Oct 25 2014
n*(2*n+1)*a(n) - 2*(4*n-1)*(4*n-3)*a(n-1) = 0. - R. J. Mathar, Oct 31 2015
O.g.f. is 1 + Revert( x*(1 + x)/(1 + 2*x)^4 ). - Peter Bala, Dec 01 2015
Sum_{n>=0} 1/a(n) = 39/25 + 4*Pi/(9*sqrt(3)) - 24*log(phi)/(25*sqrt(5)), where phi is the golden ratio (A001622). - Amiram Eldar, Mar 02 2023
From Peter Bala, Apr 29 2024: (Start)
For n >= 1, a(n) = (1/8)*Sum_{k = 0..2*n-1} (-1)^k * 4^(2*n-k)*binomial(2*n-1, k)*Catalan(k+1).
For n >= 1, a(n) = (1/8)*(16^n)*hypergeom([1 - 2*n, 3/2], [3], 1). (End)

a(0)=1 prepended by Alois P. Heinz, Nov 28 2021

A358577 Matula-Goebel numbers of "square" rooted trees, i.e., whose height equals their number of leaves.

Original entry on oeis.org

1, 4, 12, 14, 18, 19, 21, 27, 40, 52, 60, 68, 70, 74, 78, 86, 89, 90, 91, 92, 95, 100, 102, 105, 107, 111, 117, 119, 122, 129, 130, 134, 135, 138, 146, 150, 151, 153, 161, 163, 169, 170, 175, 176, 181, 183, 185, 195, 201, 206, 207, 215, 219, 221, 225, 227, 230

Offset: 1

Gus Wiseman, Nov 25 2022

nonn

The Matula-Goebel number of a rooted tree is the product of primes indexed by the Matula-Goebel numbers of the branches of its root, which gives a bijective correspondence between positive integers and unlabeled rooted trees.

			The terms together with their corresponding rooted trees begin:
   1: o
   4: (oo)
  12: (oo(o))
  14: (o(oo))
  18: (o(o)(o))
  19: ((ooo))
  21: ((o)(oo))
  27: ((o)(o)(o))
  40: (ooo((o)))
  52: (oo(o(o)))
  60: (oo(o)((o)))
  68: (oo((oo)))
  70: (o((o))(oo))
  74: (o(oo(o)))
  78: (o(o)(o(o)))
  86: (o(o(oo)))
  89: ((ooo(o)))
  90: (o(o)(o)((o)))

Gus Wiseman, The first 64 square rooted trees.

Internals instead of leaves: A358576, counted by A358587, ordered A358588.
Internals instead of height: A358578, counted by A185650, ordered A358579.
These trees are counted by A358589, ordered A358590.
A000081 counts rooted trees, ordered A000108.
A034781 counts trees by nodes and height.
A055277 counts trees by nodes and leaves, ordered A001263.
MG statistics: A061775, A109082, A109129, A196050, A342507, A358552.
Cf. A000040, A000720, A001222, A007097, A056239, A112798.
Cf. A206487, A358379, A358580, A358581-A358586, A358592.

MGTree[n_]:=If[n==1,{},MGTree/@Flatten[Cases[FactorInteger[n],{p_,k_}:>Table[PrimePi[p],{k}]]]];
Select[Range[100],Count[MGTree[#],{},{0,Infinity}]==Depth[MGTree[#]]-1&]

A358552(a(n)) = A109129(a(n)).

A358578 Matula-Goebel numbers of rooted trees whose number of leaves equals their number of internal (non-leaf) nodes.

Original entry on oeis.org

2, 6, 7, 18, 20, 21, 26, 34, 37, 43, 54, 60, 63, 67, 70, 78, 88, 91, 92, 95, 102, 111, 116, 119, 122, 129, 142, 146, 151, 162, 164, 173, 180, 181, 189, 200, 201, 202, 210, 227, 234, 236, 239, 245, 260, 264, 269, 273, 276, 278, 285, 306, 308, 314, 322, 333, 337

Offset: 1

Gus Wiseman, Nov 25 2022

nonn

			The terms together with their corresponding rooted trees begin:
   2: (o)
   6: (o(o))
   7: ((oo))
  18: (o(o)(o))
  20: (oo((o)))
  21: ((o)(oo))
  26: (o(o(o)))
  34: (o((oo)))
  37: ((oo(o)))
  43: ((o(oo)))
  54: (o(o)(o)(o))
  60: (oo(o)((o)))
  63: ((o)(o)(oo))
  67: (((ooo)))
  70: (o((o))(oo))
  78: (o(o)(o(o)))
  88: (ooo(((o))))
  91: ((oo)(o(o)))

Gus Wiseman, The first 64 rooted trees whose number of leaves equals their number of internal nodes.

These trees are counted by A185650, ordered A358579.
Height instead of leaves: A358576, counted by A358587, ordered A358588.
Height instead of internals: A358577, counted by A358589, ordered A358590.
Positions of 0's in A358580.
A000081 counts rooted trees, ordered A000108.
A034781 counts trees by nodes and height.
A055277 counts trees by nodes and leaves, ordered A001263.
MG statistics: A061775, A109082, A109129, A196050, A342507, A358552.
Cf. A000040, A000720, A001222, A007097, A056239, A112798.
Cf. A206487, A358371, A358581-A358586, A358592.

MGTree[n_]:=If[n==1,{},MGTree/@Flatten[Cases[FactorInteger[n],{p_,k_}:>Table[PrimePi[p],{k}]]]];
Select[Range[100],Count[MGTree[#],{},{0,Infinity}]==Count[MGTree[#],[_],{0,Infinity}]&]

A342507(a(n)) = A109129(a(n)).

A055290 Triangle of trees with n nodes and k leaves, 2 <= k <= n.

Original entry on oeis.org

1, 1, 0, 1, 1, 0, 1, 1, 1, 0, 1, 2, 2, 1, 0, 1, 3, 4, 2, 1, 0, 1, 4, 8, 6, 3, 1, 0, 1, 5, 14, 14, 9, 3, 1, 0, 1, 7, 23, 32, 26, 12, 4, 1, 0, 1, 8, 36, 64, 66, 39, 16, 4, 1, 0, 1, 10, 54, 123, 158, 119, 60, 20, 5, 1, 0, 1, 12, 78, 219, 350, 325, 202, 83, 25, 5, 1, 0

Offset: 2

Christian G. Bower, May 09 2000

			Triangle begins:
  n=2:  1
  n=3:  1   0
  n=4:  1   1   0
  n=5:  1   1   1   0
  n=6:  1   2   2   1   0
  n=7:  1   3   4   2   1   0
  n=8:  1   4   8   6   3   1   0
  n=9:  1   5  14  14   9   3   1   0
  n=10: 1   7  23  32  26  12   4   1   0
  n=11: 1   8  36  64  66  39  16   4   1   0
  n=12: 1  10  54 123 158 119  60  20   5   1   0
  n=13: 1  12  78 219 350 325 202  83  25   5   1   0

F. Harary and E. M. Palmer, Graphical Enumeration, Academic Press, NY, 1973, p. 80, Problem 3.9.

Andrew Howroyd, Table of n, a(n) for n = 2..1226 (rows 2..50)
Index entries for sequences related to trees

Row sums give A000055, row sums with weight k give A003228.
Columns 3 through 12: A001399(n-4), A055291, A055292, A055293, A055294, A055295, A055296, A055297, A055298, A055299.
Cf. A055300, A055301, A238416, A304222.
The labeled version is A055314.
Central column is A358107.
Left of central column is A359398.
Cf. A000088, A000272, A001349, A001433, A002494, A185650.

EulerMT(u)={my(n=#u, p=x*Ser(u), vars=variables(p)); Vec(exp( sum(i=1, n, substvec(p + O(x*x^(n\i)), vars, apply(v->v^i,vars))/i ))-1)}
T(n)={my(u=[y]); for(n=2, n, u=concat([y], EulerMT(u))); my(r=x*Ser(u), v=Vec(r*(1-x+x*y) + (substvec(r,[x,y],[x^2,y^2]) - r^2)/2)); vector(n-1, k, Vecrev(v[1+k]/y^2, k))}
{ my(A=T(10)); for(n=1, #A, print(A[n])) }

G.f.: A(x, y)=(1-x+x*y)*B(x, y)+(1/2)*(B(x^2, y^2)-B(x, y)^2), where B(x, y) is g.f. of A055277.

A358589 Number of square rooted trees with n nodes.

Original entry on oeis.org

1, 0, 1, 0, 3, 2, 11, 17, 55, 107, 317, 720, 1938, 4803, 12707, 32311, 85168, 220879, 581112, 1522095, 4014186, 10568936, 27934075, 73826753, 195497427, 517927859, 1373858931, 3646158317, 9684878325, 25737819213, 68439951884, 182070121870, 484583900955, 1290213371950

Offset: 1

Gus Wiseman, Nov 23 2022

nonn

We say that a tree is square if it has the same height as number of leaves.

			The a(1) = 1 through a(7) = 11 trees:
  o  .  (oo)  .  ((ooo))  ((o)(oo))  (((oooo)))
                 (o(oo))  (o(o)(o))  ((o(ooo)))
                 (oo(o))             ((oo(oo)))
                                     ((ooo(o)))
                                     (o((ooo)))
                                     (o(o(oo)))
                                     (o(oo(o)))
                                     (oo((oo)))
                                     (oo(o(o)))
                                     (ooo((o)))
                                     ((o)(o)(o))

Andrew Howroyd, Table of n, a(n) for n = 1..200

For internals instead of height we have A185650 aerated, ranked by A358578.
These trees are ranked by A358577.
For internals instead of leaves we have A358587, ranked by A358576.
The ordered version is A358590.
A000081 counts rooted trees, ordered A000108.
A034781 counts rooted trees by nodes and height, ordered A080936.
A055277 counts rooted trees by nodes and leaves, ordered A001263.
A358575 counts rooted trees by nodes and internal nodes, ordered A090181.
Cf. A000891, A065097, A109129, A358552, A358588, A358591.

art[n_]:=If[n==1,{{}},Join@@Table[Select[Tuples[art/@c],OrderedQ],{c,Join@@Permutations/@IntegerPartitions[n-1]}]];
Table[Length[Select[art[n],Count[#,{},{0,Infinity}]==Depth[#]-1&]],{n,1,10}]

\\ R(n,f) enumerates trees by height(h), nodes(x) and leaves(y).
R(n,f) = {my(A=O(x*x^n), Z=0); for(h=1, n, my(p = A); A = x*(y - 1  + exp( sum(i=1, n-1, 1/i * subst( subst( A + O(x*x^((n-1)\i)), x, x^i), y, y^i) ) )); Z += f(h, A-p)); Z}
seq(n) = {Vec(R(n, (h,p)->polcoef(p,h,y)), -n)} \\ Andrew Howroyd, Jan 01 2023

Terms a(19) and beyond from Andrew Howroyd, Jan 01 2023

A358576 Matula-Goebel numbers of rooted trees whose node-height equals their number of internal (non-leaf) nodes.

Original entry on oeis.org

9, 15, 18, 21, 23, 30, 33, 35, 36, 39, 42, 46, 47, 49, 51, 57, 60, 61, 66, 70, 72, 73, 77, 78, 83, 84, 87, 91, 92, 93, 94, 95, 98, 102, 111, 113, 114, 119, 120, 122, 123, 129, 132, 133, 137, 140, 144, 146, 149, 151, 154, 156, 159, 166, 167, 168, 174, 177, 181

Offset: 1

Gus Wiseman, Nov 25 2022

nonn

Node-height is the number of nodes in the longest path from root to leaf.

			The terms together with their corresponding rooted trees begin:
   9: ((o)(o))
  15: ((o)((o)))
  18: (o(o)(o))
  21: ((o)(oo))
  23: (((o)(o)))
  30: (o(o)((o)))
  33: ((o)(((o))))
  35: (((o))(oo))
  36: (oo(o)(o))
  39: ((o)(o(o)))
  42: (o(o)(oo))
  46: (o((o)(o)))
  47: (((o)((o))))
  49: ((oo)(oo))
  51: ((o)((oo)))
  57: ((o)(ooo))
  60: (oo(o)((o)))
  61: ((o(o)(o)))

Gus Wiseman, The first 64 rooted trees whose node-height equals their number of internal nodes.

The version for edge-height is A209638.
Square trees are A358577, counted by A358589, ordered A358590.
The version for leaves instead of height is A358578, counted by A185650.
These trees are counted by A358587, ordered A358588.
A000081 counts rooted trees, ordered A000108.
A034781 counts rooted trees by nodes and height.
A055277 counts rooted trees by leaves, ordered A001263.
MG statistics: A061775, A109082, A109129, A196050, A342507, A358552.
Cf. A000040, A000720, A001222, A007097, A056239, A112798.
Cf. A206487, A358580, A358581-A358586, A358592, A358729.

MGTree[n_]:=If[n==1,{},MGTree/@Flatten[Cases[FactorInteger[n],{p_,k_}:>Table[PrimePi[p],{k}]]]];
Select[Range[100],Count[MGTree[#],[_],{0,Infinity}]==Depth[MGTree[#]]-1&]

A358552(a(n)) = A342507(a(n)).

A358580 Difference between the number of leaves and the number of internal (non-leaf) nodes in the rooted tree with Matula-Goebel number n.

Original entry on oeis.org

1, 0, -1, 1, -2, 0, 0, 2, -1, -1, -3, 1, -1, 1, -2, 3, -1, 0, 1, 0, 0, -2, -2, 2, -3, 0, -1, 2, -2, -1, -4, 4, -3, 0, -1, 1, 0, 2, -1, 1, -2, 1, 0, -1, -2, -1, -3, 3, 1, -2, -1, 1, 2, 0, -4, 3, 1, -1, -2, 0, -1, -3, 0, 5, -2, -2, 0, 1, -2, 0, -1, 2, -1, 1, -3

Offset: 1

Gus Wiseman, Nov 25 2022

sign

			The Matula-Goebel number of ((ooo(o))) is 89, and it has 4 leaves and 3 internal nodes, so a(89) = 1.

Zeros are A358578, counted by A185650 (ordered A358579).
Positions of positive terms are counted by A358581, negative A358582.
Positions of nonnegative terms are counted by A358583, nonpositive A358584.
A000081 counts rooted trees, ordered A000108.
A034781 counts trees by nodes and height.
A055277 counts trees by nodes and leaves, ordered A001263.
MG statistics: A061775, A109082, A109129, A196050, A342507, A358552.
Cf. A000040, A000720, A001222, A007097, A056239, A112798.
Cf. A206487, A358371, A358576, A358577, A358592.

MGTree[n_]:=If[n==1,{},MGTree/@Flatten[Cases[FactorInteger[n],{p_,k_}:>Table[PrimePi[p],{k}]]]];
Table[Count[MGTree[n],{},{0,Infinity}]-Count[MGTree[n],[_],{0,Infinity}],{n,100}]

a(n) = A109129(n) - A342507(n).

A358590 Number of square ordered rooted trees with n nodes.

Original entry on oeis.org

1, 0, 1, 0, 6, 5, 36, 84, 309, 890, 3163, 9835, 32979, 108252, 360696, 1192410, 3984552, 13276769, 44371368, 148402665, 497072593, 1665557619, 5586863093, 18750662066, 62968243731, 211565969511, 711187790166, 2391640404772, 8045964959333, 27077856222546

Offset: 1

Gus Wiseman, Nov 25 2022

nonn

We say that a tree is square if it has the same height as number of leaves.

			The a(1) = 1 through a(6) = 5 ordered trees:
  o  .  (oo)  .  ((o)oo)  ((o)(o)o)
                 ((oo)o)  ((o)(oo))
                 ((ooo))  ((o)o(o))
                 (o(o)o)  ((oo)(o))
                 (o(oo))  (o(o)(o))
                 (oo(o))

Andrew Howroyd, Table of n, a(n) for n = 1..200

For internals instead of height we have A000891, unordered A185650 aerated.
For internals instead of leaves we have A358588, unordered A358587.
The unordered version is A358589, ranked by A358577.
A000108 counts ordered rooted trees, unordered A000081.
A001263 counts ordered rooted trees by nodes and leaves, unordered A055277.
A080936 counts ordered rooted trees by nodes and height, unordered A034781.
A090181 counts ordered rooted trees by nodes and internals, unord. A358575.
Cf. A065097, A109129, A358371, A358552, A358579, A358586.

aot[n_]:=If[n==1,{{}},Join@@Table[Tuples[aot/@c],{c,Join@@Permutations/@IntegerPartitions[n-1]}]];
Table[Length[Select[aot[n],Count[#,{},{0,Infinity}]==Depth[#]-1&]],{n,1,10}]

\\ R(n,f) enumerates trees by height(h), nodes(x) and leaves(y).
R(n,f) = {my(A=O(x*x^n), Z=0); for(h=1, n, my(p = A); A = x*(y - 1  + 1/(1 - A + O(x^n))); Z += f(h, A-p)); Z}
seq(n) = {Vec(R(n, (h,p)->polcoef(p,h,y)), -n)} \\ Andrew Howroyd, Jan 01 2023

Terms a(16) and beyond from Andrew Howroyd, Jan 01 2023

A358586 Number of ordered rooted trees with n nodes, at least half of which are leaves.

Original entry on oeis.org

1, 1, 1, 4, 7, 31, 66, 302, 715, 3313, 8398, 39095, 104006, 484706, 1337220, 6227730, 17678835, 82204045, 238819350, 1108202513, 3282060210, 15195242478, 45741281820, 211271435479, 644952073662, 2971835602526, 9183676536076, 42217430993002, 131873975875180, 604834233372884

Offset: 1

Gus Wiseman, Nov 24 2022

nonn

			The a(1) = 1 through a(5) = 7 ordered trees:
  o  (o)  (oo)  (ooo)   (oooo)
                ((o)o)  ((o)oo)
                ((oo))  ((oo)o)
                (o(o))  ((ooo))
                        (o(o)o)
                        (o(oo))
                        (oo(o))

Andrew Howroyd, Table of n, a(n) for n = 1..1000

For equality we have A000891, unordered A185650.
Odd-indexed terms appear to be A065097.
The unordered version is A358583.
The opposite is the same, unordered A358584.
The strict case is A358585, unordered A358581.
A000108 counts ordered rooted trees, unordered A000081.
A001263 counts ordered rooted trees by nodes and leaves, unordered A055277.
A080936 counts ordered rooted trees by nodes and height, unordered A034781.
A090181 counts ordered rooted trees by nodes and internals, unord. A358575.
A358590 counts square ordered trees, unordered A358589 (ranked by A358577).
Cf. A109129, A342507, A358371, A358579, A358580, A358582, A358584, A358588.

aot[n_]:=If[n==1,{{}},Join@@Table[Tuples[aot/@c],{c,Join@@Permutations/@IntegerPartitions[n-1]}]];
Table[Length[Select[aot[n],Count[#,{},{0,Infinity}]>=Count[#,[_],{0,Infinity}]&]],{n,1,10}]

a(n) = if(n==1, 1, n--; (binomial(2*n,n)/(n+1) + if(n%2, binomial(n, (n-1)/2)^2 / n))/2) \\ Andrew Howroyd, Jan 13 2024

From Andrew Howroyd, Jan 13 2024: (Start)
a(n) = Sum_{k=1..floor(n/2)} A001263(n-1, k) for n >= 2.
a(2*n) = (A000108(2*n-1) + A000891(n-1))/2 for n >= 1;
a(2*n+1) = A000108(2*n)/2 for n >= 1. (End)

a(16) onwards from Andrew Howroyd, Jan 13 2024

cp's OEIS Frontend

A000891 a(n) = (2*n)!*(2*n+1)! / (n! * (n+1)!)^2.

Views

Author

Keywords

Comments

Examples

References

Links

Crossrefs

Programs

Haskell

Magma

Maple

Mathematica

PARI

PARI

Formula

Extensions

A065097 a(n) = ((2n+1) + (2n-1) - 1)!/((2n+1)!*(2n-1)!).

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Magma

Maple

Mathematica

MuPAD

PARI

PARI

Sage

Formula

Extensions

A358577 Matula-Goebel numbers of "square" rooted trees, i.e., whose height equals their number of leaves.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Mathematica

Formula

A358578 Matula-Goebel numbers of rooted trees whose number of leaves equals their number of internal (non-leaf) nodes.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Mathematica

Formula

A055290 Triangle of trees with n nodes and k leaves, 2 <= k <= n.

Views

Author

Keywords

Examples

References

Links

Crossrefs

Programs

PARI

Formula

A358589 Number of square rooted trees with n nodes.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Mathematica

A000891 a(n) = (2n)!(2n+1)! / (n! (n+1)!)^2.