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 41 results. Next

A125558 Central column of triangle A090181.

Original entry on oeis.org

1, 1, 6, 50, 490, 5292, 60984, 736164, 9202050, 118195220, 1551580888, 20734762776, 281248448936, 3863302870000, 53644719852000, 751920156592200, 10626401036545650, 151269944167296900, 2167317913508055000
Offset: 0

Views

Author

Philippe Deléham, Jan 01 2007, Oct 11 2007

Keywords

Comments

[1,6,50,490,5292,...] is a column in triangle of Narayana numbers A001263.
Number of Dyck 2n-paths with exactly n peaks. - Peter Luschny, May 10 2014
For n > 0, number of pairs of non-intersecting lattice paths with steps (1,0), (0,1), where one path goes from (0,0) to (n,n) and the other from (1,0) to (n+1,n). The proof is by switching intersecting path pairs after their first intersection, giving a(n) = binomial(2*n,n)^2 - binomial(2*n+1,n) * binomial(2*n-1,n). - Jeremy Tan, Apr 12 2021

Crossrefs

Equals A000888(n)/2 for n>0.
Cf. A090181.

Programs

  • Maple
    seq(ceil(1/2*(n+1)*((binomial(2*n,n)/(1+n))^2)), n=0..18); # Zerinvary Lajos, Jun 18 2007
  • Mathematica
    CoefficientList[
    Series[1 + (HypergeometricPFQ[{1/2, 1/2}, {2}, 16 x] - 1)/(2), {x, 0,
        20}], x]
    Join[{1},Table[CatalanNumber[n]^2 (n+1)/2,{n,20}]] (* Harvey P. Dale, Oct 19 2011 *)

Formula

a(0)=1, a(n) = Catalan(n)^2*(n+1)/2 = A000108(n)^2*(n+1)/2 for n>0.
a(n) = A090181(2*n, n).
G.f.: 1 + x*3F2( 1, 3/2, 3/2; 2, 3;16 x) = 1 + ( 2F1( 1/2, 1/2; 2;16*x) - 1)/2. - Olivier Gérard, Feb 16 2011
D-finite with recurrence n*(n+1)*a(n) -4*(2*n-1)^2*a(n-1)=0. - R. J. Mathar, Feb 08 2021
a(n) = binomial(2*n,n)^2 - binomial(2*n+1,n) * binomial(2*n-1,n). - Jeremy Tan, Apr 12 2021

A130749 Triangle A007318*A090181 (as infinite lower triangular matrices) .

Original entry on oeis.org

1, 1, 1, 1, 3, 1, 1, 7, 6, 1, 1, 15, 24, 10, 1, 1, 31, 80, 60, 15, 1, 1, 63, 240, 280, 125, 21, 1, 1, 127, 672, 1120, 770, 231, 28, 1, 1, 255, 1792, 4032, 3920, 1806, 392, 36, 1, 1, 511, 4608, 13440, 17472, 11340, 3780, 624, 45, 1
Offset: 0

Views

Author

Philippe Deléham, Jul 13 2007

Keywords

Examples

			Triangle begins:
  1;
  1,   1;
  1,   3,    1;
  1,   7,    6,     1;
  1,  15,   24,    10,     1;
  1,  31,   80,    60,    15,     1;
  1,  63,  240,   280,   125,    21,    1;
  1, 127,  672,  1120,   770,   231,   28,   1;
  1, 255, 1792,  4032,  3920,  1806,  392,  36,  1;
  1, 511, 4608, 13440, 17472, 11340, 3780, 624, 45,  1;
  ...
		

Crossrefs

Programs

  • Mathematica
    nmax = 9;
    T1[n_, k_] := Binomial[n, k];
    T2[n_, k_] := Sum[(-1)^(j-k) Binomial[2n-j, j] Binomial[j, k] CatalanNumber[n-j], {j, 0, n}];
    T[n_, k_] := Sum[T1[n, m] T2[m, k], {m, 0, n}];
    Table[T[n, k], {n, 0, nmax}, {k, 0, n}] // Flatten (* Jean-François Alcover, Nov 10 2018 *)
  • Maxima
    N(n, k):=(binomial(n, k-1)*binomial(n, k))/n;
    T(n, k):=if k=0 then 1 else sum(binomial(n, i)*N(i, k), i, 1, n); /* Vladimir Kruchinin, Jan 08 2022 */

Formula

Sum_{k=0..n} T(n,k) = A007317(n+1).
G.f.: 1/(1-x-xy/(1-x/(1-x-xy/(1-x/(1-x-xy/(1-x.... (continued fraction); [Paul Barry, Jan 12 2009]
T(n,k) = Sum_{i=1..n} binomial(n, i)*N(i,k), T(n,0)=1, where N(n,k) is the triangle of Narayana numbers A001263. - Vladimir Kruchinin, Jan 08 2022

A001263 Triangle of Narayana numbers T(n,k) = C(n-1,k-1)*C(n,k-1)/k with 1 <= k <= n, read by rows. Also called the Catalan triangle.

Original entry on oeis.org

1, 1, 1, 1, 3, 1, 1, 6, 6, 1, 1, 10, 20, 10, 1, 1, 15, 50, 50, 15, 1, 1, 21, 105, 175, 105, 21, 1, 1, 28, 196, 490, 490, 196, 28, 1, 1, 36, 336, 1176, 1764, 1176, 336, 36, 1, 1, 45, 540, 2520, 5292, 5292, 2520, 540, 45, 1, 1, 55, 825, 4950, 13860, 19404, 13860, 4950, 825
Offset: 1

Views

Author

Keywords

Comments

Number of antichains (or order ideals) in the poset 2*(k-1)*(n-k) or plane partitions with rows <= k-1, columns <= n-k and entries <= 2. - Mitch Harris, Jul 15 2000
T(n,k) is the number of Dyck n-paths with exactly k peaks. a(n,k) = number of pairs (P,Q) of lattice paths from (0,0) to (k,n+1-k), each consisting of unit steps East or North, such that P lies strictly above Q except at the endpoints. - David Callan, Mar 23 2004
Number of permutations of [n] which avoid-132 and have k-1 descents. - Mike Zabrocki, Aug 26 2004
T(n,k) is the number of paths through n panes of glass, entering and leaving from one side, of length 2n with k reflections (where traversing one pane of glass is the unit length). - Mitch Harris, Jul 06 2006
Antidiagonal sums given by A004148 (without first term).
T(n,k) is the number of full binary trees with n internal nodes and k-1 jumps. In the preorder traversal of a full binary tree, any transition from a node at a deeper level to a node on a strictly higher level is called a jump. - Emeric Deutsch, Jan 18 2007
From Gary W. Adamson, Oct 22 2007: (Start)
The n-th row can be generated by the following operation using an ascending row of (n-1) triangular terms, (A) and a descending row, (B); e.g., row 6:
A: 1....3....6....10....15
B: 15...10....6.....3.....1
C: 1...15...50....50....15....1 = row 6.
Leftmost column of A,B -> first two terms of C; then followed by the operation B*C/A of current column = next term of row C, (e.g., 10*15/3 = 50). Continuing with the operation, we get row 6: (1, 15, 50, 50, 15, 1). (End)
The previous comment can be upgraded to: The ConvOffsStoT transform of the triangular series; and by rows, row 6 is the ConvOffs transform of (1, 3, 6, 10, 15). Refer to triangle A117401 as another example of the ConvOffsStoT transform, and OEIS under Maple Transforms. - Gary W. Adamson, Jul 09 2012
For a connection to Lagrange inversion, see A134264. - Tom Copeland, Aug 15 2008
T(n,k) is also the number of order-decreasing and order-preserving mappings (of an n-element set) of height k (height of a mapping is the cardinal of its image set). - Abdullahi Umar, Aug 21 2008
Row n of this triangle is the h-vector of the simplicial complex dual to an associahedron of type A_n [Fomin & Reading, p.60]. See A033282 for the corresponding array of f-vectors for associahedra of type A_n. See A008459 and A145903 for the h-vectors for associahedra of type B and type D respectively. The Hilbert transform of this triangle (see A145905 for the definition of this transform) is A145904. - Peter Bala, Oct 27 2008
T(n,k) is also the number of noncrossing set partitions of [n] into k blocks. Given a partition P of the set {1,2,...,n}, a crossing in P are four integers [a, b, c, d] with 1 <= a < b < c < d <= n for which a, c are together in a block, and b, d are together in a different block. A noncrossing partition is a partition with no crossings. - Peter Luschny, Apr 29 2011
Noncrossing set partitions are also called genus 0 partitions. In terms of genus-dependent Stirling numbers of the second kind S2(n,k,g) that count partitions of genus g of an n-set into k nonempty subsets, one has T(n,k) = S2(n,k,0). - Robert Coquereaux, Feb 15 2024
Diagonals of A089732 are rows of A001263. - Tom Copeland, May 14 2012
From Peter Bala, Aug 07 2013: (Start)
Let E(y) = Sum_{n >= 0} y^n/(n!*(n+1)!) = 1/sqrt(y)*BesselI(1,2*sqrt(y)). Then this triangle is the generalized Riordan array (E(y), y) with respect to the sequence n!*(n+1)! as defined in Wang and Wang.
Generating function E(y)*E(x*y) = 1 + (1 + x)*y/(1!*2!) + (1 + 3*x + x^2)*y^2/(2!*3!) + (1 + 6*x + 6*x^2 + x^3)*y^3/(3!*4!) + .... Cf. A105278 with a generating function exp(y)*E(x*y).
The n-th power of this array has a generating function E(y)^n*E(x*y). In particular, the matrix inverse A103364 has a generating function E(x*y)/E(y). (End)
T(n,k) is the number of nonintersecting n arches above the x axis, starting and ending on vertices 1 to 2n, with k being the number of arches starting on an odd vertice and ending on a higher even vertice. Example: T(3,2)=3 [16,25,34] [14,23,56] [12,36,45]. - Roger Ford, Jun 14 2014
Fomin and Reading on p. 31 state that the rows of the Narayana matrix are the h-vectors of the associahedra as well as its dual. - Tom Copeland, Jun 27 2017
The row polynomials P(n, x) = Sum_{k=1..n} T(n, k)*x^(k-1), together with P(0, x) = 1, multiplied by (n+1) are the numerator polynomials of the o.g.f.s of the diagonal sequences of the triangle A103371: G(n, x) = (n+1)*P(n, x)/(1 - x)^{2*n+1}, for n >= 0. This is proved with Lagrange's theorem applied to the Riordan triangle A135278 = (1/(1 - x)^2, x/(1 - x)). See an example below. - Wolfdieter Lang, Jul 31 2017
T(n,k) is the number of Dyck paths of semilength n with k-1 uu-blocks (pairs of consecutive up-steps). - Alexander Burstein, Jun 22 2020
In case you were searching for Narayama numbers, the correct spelling is Narayana. - N. J. A. Sloane, Nov 11 2020
Named after the Canadian mathematician Tadepalli Venkata Narayana (1930-1987). They were also called "Runyon numbers" after John P. Runyon (1922-2013) of Bell Telephone Laboratories, who used them in a study of a telephone traffic system. - Amiram Eldar, Apr 15 2021 The Narayana numbers were first studied by Percy Alexander MacMahon (see reference, Article 495) as pointed out by Bóna and Sagan (see link). - Peter Luschny, Apr 28 2022
From Andrea Arlette España, Nov 14 2022: (Start)
T(n,k) is the degree distribution of the paths towards synchronization in the transition diagram associated with the Laplacian system over the complete graph K_n, corresponding to ordered initial conditions x_1 < x_2 < ... < x_n.
T(n,k) for n=2N+1 and k=N+1 is the number of states in the transition diagram associated with the Laplacian system over the complete bipartite graph K_{N,N}, corresponding to ordered (x_1 < x_2 < ... < x_N and x_{N+1} < x_{N+2} < ... < x_{2N}) and balanced (Sum_{i=1..N} x_i/N = Sum_{i=N+1..2N} x_i/N) initial conditions. (End)
From Gus Wiseman, Jan 23 2023: (Start)
Also the number of unlabeled ordered rooted trees with n nodes and k leaves. See the link by Marko Riedel. For example, row n = 5 counts the following trees:
((((o)))) (((o))o) ((o)oo) (oooo)
(((o)o)) ((oo)o)
(((oo))) ((ooo))
((o)(o)) (o(o)o)
((o(o))) (o(oo))
(o((o))) (oo(o))
The unordered version is A055277. Leaves in standard ordered trees are counted by A358371. (End)

Examples

			The initial rows of the triangle are:
  [1] 1
  [2] 1,  1
  [3] 1,  3,   1
  [4] 1,  6,   6,    1
  [5] 1, 10,  20,   10,    1
  [6] 1, 15,  50,   50,   15,    1
  [7] 1, 21, 105,  175,  105,   21,   1
  [8] 1, 28, 196,  490,  490,  196,  28,  1
  [9] 1, 36, 336, 1176, 1764, 1176, 336, 36, 1;
  ...
For all n, 12...n (1 block) and 1|2|3|...|n (n blocks) are noncrossing set partitions.
Example of umbral representation:
  A007318(5,k)=[1,5/1,5*4/(2*1),...,1]=(1,5,10,10,5,1),
  so A001263(5,k)={1,b(5)/b(1),b(5)*b(4)/[b(2)*b(1)],...,1}
  = [1,30/2,30*20/(6*2),...,1]=(1,15,50,50,15,1).
  First = last term = b.(5!)/[b.(0!)*b.(5!)]= 1. - _Tom Copeland_, Sep 21 2011
Row polynomials and diagonal sequences of A103371: n = 4,  P(4, x) = 1 + 6*x + 6*x^2 + x^3, and the o.g.f. of fifth diagonal is G(4, x) = 5* P(4, x)/(1 - x)^9, namely [5, 75, 525, ...]. See a comment above. - _Wolfdieter Lang_, Jul 31 2017
		

References

  • Berman and Koehler, Cardinalities of finite distributive lattices, Mitteilungen aus dem Mathematischen Seminar Giessen, 121 (1976), pp. 103-124.
  • Miklos Bona, editor, Handbook of Enumerative Combinatorics, CRC Press, 2015, page 196.
  • P. A. MacMahon, Combinatory Analysis, Vols. 1 and 2, Cambridge University Press, 1915, 1916; reprinted by Chelsea, 1960, Sect. 495.
  • T. V. Narayana, Lattice Path Combinatorics with Statistical Applications. Univ. Toronto Press, 1979, pp. 100-101.
  • A. Nkwanta, Lattice paths and RNA secondary structures, in African Americans in Mathematics, ed. N. Dean, Amer. Math. Soc., 1997, pp. 137-147.
  • T. K. Petersen, Eulerian Numbers, Birkhäuser, 2015, Chapter 2.
  • J. Riordan, Combinatorial Identities, Wiley, 1968, p. 17.
  • R. P. Stanley, Enumerative Combinatorics, Cambridge, Vol. 2, 1999; see Problem 6.36(a) and (b).

Crossrefs

Other versions are in A090181 and A131198. - Philippe Deléham, Nov 18 2007
Cf. variants: A181143, A181144. - Paul D. Hanna, Oct 13 2010
Row sums give A000108 (Catalan numbers), n>0.
A008459 (h-vectors type B associahedra), A033282 (f-vectors type A associahedra), A145903 (h-vectors type D associahedra), A145904 (Hilbert transform). - Peter Bala, Oct 27 2008
Cf. A016098 and A189232 for numbers of crossing set partitions.
Cf. A243752.
Triangles of generalized binomial coefficients (n,k)_m (or generalized Pascal triangles) for m = 1,...,12: A007318 (Pascal), A001263, A056939, A056940, A056941, A142465, A142467, A142468, A174109, A342889, A342890, A342891.

Programs

  • GAP
    Flat(List([1..11],n->List([1..n],k->Binomial(n-1,k-1)*Binomial(n,k-1)/k))); # Muniru A Asiru, Jul 12 2018
  • Haskell
    a001263 n k = a001263_tabl !! (n-1) !! (k-1)
    a001263_row n = a001263_tabl !! (n-1)
    a001263_tabl = zipWith dt a007318_tabl (tail a007318_tabl) where
       dt us vs = zipWith (-) (zipWith (*) us (tail vs))
                              (zipWith (*) (tail us ++ [0]) (init vs))
    -- Reinhard Zumkeller, Oct 10 2013
    
  • Magma
    /* triangle */ [[Binomial(n-1,k-1)*Binomial(n,k-1)/k : k in [1..n]]: n in [1.. 15]]; // Vincenzo Librandi, Oct 19 2014
    
  • Maple
    A001263 := (n,k)->binomial(n-1,k-1)*binomial(n,k-1)/k;
    a:=proc(n,k) option remember; local i; if k=1 or k=n then 1 else add(binomial(n+i-1, 2*k-2)*a(k-1,i),i=1..k-1); fi; end:
    # Alternatively, as a (0,0)-based triangle:
    R := n -> simplify(hypergeom([-n, -n-1], [2], x)): Trow := n -> seq(coeff(R(n,x),x,j), j=0..n): seq(Trow(n), n=0..9); # Peter Luschny, Mar 19 2018
  • Mathematica
    T[n_, k_] := If[k==0, 0, Binomial[n-1, k-1] Binomial[n, k-1] / k];
    Flatten[Table[Binomial[n-1,k-1] Binomial[n,k-1]/k,{n,15},{k,n}]] (* Harvey P. Dale, Feb 29 2012 *)
    TRow[n_] := CoefficientList[Hypergeometric2F1[1 - n, -n, 2, x], x];
    Table[TRow[n], {n, 1, 11}] // Flatten (* Peter Luschny, Mar 19 2018 *)
    aot[n_]:=If[n==1,{{}},Join@@Table[Tuples[aot/@c],{c,Join@@Permutations/@IntegerPartitions[n-1]}]];
    Table[Length[Select[aot[n],Length[Position[#,{}]]==k&]],{n,2,9},{k,1,n-1}] (* Gus Wiseman, Jan 23 2023 *)
    T[1, 1] := 1; T[n_, k_]/;1<=k<=n := T[n, k] = (2n/k-1) T[n-1,k-1] + T[n-1, k]; T[n_, k_] := 0; Flatten@Table[T[n, k], {n, 1, 11}, {k, 1, n}] (* Oliver Seipel, Dec 31 2024 *)
  • PARI
    {a(n, k) = if(k==0, 0, binomial(n-1, k-1) * binomial(n, k-1) / k)};
    
  • PARI
    {T(n,k)=polcoeff(polcoeff(exp(sum(m=1,n,sum(j=0,m,binomial(m,j)^2*y^j)*x^m/m) +O(x^(n+1))),n,x),k,y)} \\ Paul D. Hanna, Oct 13 2010
    
  • Sage
    @CachedFunction
    def T(n, k):
        if k == n or k == 1: return 1
        if k <= 0 or k > n: return 0
        return binomial(n, 2) * (T(n-1, k)/((n-k)*(n-k+1)) + T(n-1, k-1)/(k*(k-1)))
    for n in (1..9): print([T(n, k) for k in (1..n)])  # Peter Luschny, Oct 28 2014
    

Formula

a(n, k) = C(n-1, k-1)*C(n, k-1)/k for k!=0; a(n, 0)=0.
Triangle equals [0, 1, 0, 1, 0, 1, ...] DELTA [1, 0, 1, 0, 1, 0, 1, 0, ...] where DELTA is Deléham's operator defined in A084938.
0Mike Zabrocki, Aug 26 2004
T(n, k) = C(n, k)*C(n-1, k-1) - C(n, k-1)*C(n-1, k) (determinant of a 2 X 2 subarray of Pascal's triangle A007318). - Gerald McGarvey, Feb 24 2005
T(n, k) = binomial(n-1, k-1)^2 - binomial(n-1, k)*binomial(n-1, k-2). - David Callan, Nov 02 2005
a(n,k) = C(n,2) (a(n-1,k)/((n-k)*(n-k+1)) + a(n-1,k-1)/(k*(k-1))) a(n,k) = C(n,k)*C(n,k-1)/n. - Mitch Harris, Jul 06 2006
Central column = A000891, (2n)!*(2n+1)! / (n!*(n+1)!)^2. - Zerinvary Lajos, Oct 29 2006
G.f.: (1-x*(1+y)-sqrt((1-x*(1+y))^2-4*y*x^2))/(2*x) = Sum_{n>0, k>0} a(n, k)*x^n*y^k.
From Peter Bala, Oct 22 2008: (Start)
Relation with Jacobi polynomials of parameter (1,1):
Row n+1 generating polynomial equals 1/(n+1)*x*(1-x)^n*Jacobi_P(n,1,1,(1+x)/(1-x)). It follows that the zeros of the Narayana polynomials are all real and nonpositive, as noted above. O.g.f for column k+2: 1/(k+1) * y^(k+2)/(1-y)^(k+3) * Jacobi_P(k,1,1,(1+y)/(1-y)). Cf. A008459.
T(n+1,k) is the number of walks of n unit steps on the square lattice (i.e., each step in the direction either up (U), down (D), right (R) or left (L)) starting from the origin and finishing at lattice points on the x axis and which remain in the upper half-plane y >= 0 [Guy]. For example, T(4,3) = 6 counts the six walks RRL, LRR, RLR, UDL, URD and RUD, from the origin to the lattice point (1,0), each of 3 steps. Compare with tables A145596 - A145599.
Define a functional I on formal power series of the form f(x) = 1 + ax + bx^2 + ... by the following iterative process. Define inductively f^(1)(x) = f(x) and f^(n+1)(x) = f(x*f^(n)(x)) for n >= 1. Then set I(f(x)) = lim_{n -> infinity} f^(n)(x) in the x-adic topology on the ring of formal power series; the operator I may also be defined by I(f(x)) := 1/x*series reversion of x/f(x).
The o.g.f. for this array is I(1 + t*x + t*x^2 + t*x^3 + ...) = 1 + t*x + (t + t^2)*x^2 + (t + 3*t^2 + t^3)*x^3 + ... = 1/(1 - x*t/(1 - x/(1 - x*t/(1 - x/(1 - ...))))) (as a continued fraction). Cf. A108767, A132081 and A141618. (End)
G.f.: 1/(1-x-xy-x^2y/(1-x-xy-x^2y/(1-... (continued fraction). - Paul Barry, Sep 28 2010
E.g.f.: exp((1+y)x)*Bessel_I(1,2*sqrt(y)x)/(sqrt(y)*x). - Paul Barry, Sep 28 2010
G.f.: A(x,y) = exp( Sum_{n>=1} [Sum_{k=0..n} C(n,k)^2*y^k] * x^n/n ). - Paul D. Hanna, Oct 13 2010
With F(x,t) = (1-(1+t)*x-sqrt(1-2*(1+t)*x+((t-1)*x)^2))/(2*x) an o.g.f. in x for the Narayana polynomials in t, G(x,t) = x/(t+(1+t)*x+x^2) is the compositional inverse in x. Consequently, with H(x,t) = 1/ (dG(x,t)/dx) = (t+(1+t)*x+x^2)^2 / (t-x^2), the n-th Narayana polynomial in t is given by (1/n!)*((H(x,t)*D_x)^n)x evaluated at x=0, i.e., F(x,t) = exp(x*H(u,t)*D_u)u, evaluated at u = 0. Also, dF(x,t)/dx = H(F(x,t),t). - Tom Copeland, Sep 04 2011
With offset 0, A001263 = Sum_{j>=0} A132710^j / A010790(j), a normalized Bessel fct. May be represented as the Pascal matrix A007318, n!/[(n-k)!*k!], umbralized with b(n)=A002378(n) for n>0 and b(0)=1: A001263(n,k)= b.(n!)/{b.[(n-k)!]*b.(k!)} where b.(n!) = b(n)*b(n-1)...*b(0), a generalized factorial (see example). - Tom Copeland, Sep 21 2011
With F(x,t) = {1-(1-t)*x-sqrt[1-2*(1+t)*x+[(t-1)*x]^2]}/2 a shifted o.g.f. in x for the Narayana polynomials in t, G(x,t)= x/[t-1+1/(1-x)] is the compositional inverse in x. Therefore, with H(x,t)=1/(dG(x,t)/dx)=[t-1+1/(1-x)]^2/{t-[x/(1-x)]^2}, (see A119900), the (n-1)-th Narayana polynomial in t is given by (1/n!)*((H(x,t)*d/dx)^n)x evaluated at x=0, i.e., F(x,t) = exp(x*H(u,t)*d/du) u, evaluated at u = 0. Also, dF(x,t)/dx = H(F(x,t),t). - Tom Copeland, Sep 30 2011
T(n,k) = binomial(n-1,k-1)*binomial(n+1,k)-binomial(n,k-1)*binomial(n,k). - Philippe Deléham, Nov 05 2011
A166360(n-k) = T(n,k) mod 2. - Reinhard Zumkeller, Oct 10 2013
Damped sum of a column, in leading order: lim_{d->0} d^(2k-1) Sum_{N>=k} T(N,k)(1-d)^N=Catalan(n). - Joachim Wuttke, Sep 11 2014
Multiplying the n-th column by n! generates the revert of the unsigned Lah numbers, A089231. - Tom Copeland, Jan 07 2016
Row polynomials: (x - 1)^(n+1)*(P(n+1,(1 + x)/(x - 1)) - P(n-1,(1 + x)/(x - 1)))/((4*n + 2)), n = 1,2,... and where P(n,x) denotes the n-th Legendre polynomial. - Peter Bala, Mar 03 2017
The coefficients of the row polynomials R(n, x) = hypergeom([-n,-n-1], [2], x) generate the triangle based in (0,0). - Peter Luschny, Mar 19 2018
Multiplying the n-th diagonal by n!, with the main diagonal n=1, generates the Lah matrix A105278. With G equal to the infinitesimal generator of A132710, the Narayana triangle equals Sum_{n >= 0} G^n/((n+1)!*n!) = (sqrt(G))^(-1) * I_1(2*sqrt(G)), where G^0 is the identity matrix and I_1(x) is the modified Bessel function of the first kind of order 1. (cf. Sep 21 2011 formula also.) - Tom Copeland, Sep 23 2020
T(n,k) = T(n,k-1)*C(n-k+2,2)/C(k,2). - Yuchun Ji, Dec 21 2020
From Sergii Voloshyn, Nov 25 2024: (Start)
G.f.: F(x,y) = (1-x*(1+y)-sqrt((1-x*(1+y))^2-4*y*x^2))/(2*x) is the solution of the differential equation x^3 * d^2(x*F(x,y))/dx^2 = y * d^2(x*F(x,y))/dy^2.
Let E be the operator x*D*D, where D denotes the derivative operator d/dx. Then (1/(n! (1 + n)!)) * E^n(x/(1 - x)) = (row n generating polynomial)/(1 - x)^(2*n+1) = Sum_{k >= 0} C(n-1, k-1)*C(n, k-1)/k*x^k. For example, when n = 4 we have (1/4!/5!)*E^3(x/(1 - x)) = x (1 + 6 x + 6 x^2 + x^3)/(1 - x)^9. (End)

Extensions

Deleted certain dangerous or potentially dangerous links. - N. J. A. Sloane, Jan 30 2021

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

Original entry on oeis.org

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

Views

Author

Keywords

Comments

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

Examples

			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)
		

References

  • 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.

Crossrefs

Cf. A145600, A145601, A145602, A145603. - Peter Bala, Oct 14 2008
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.

Programs

  • Haskell
    a000891 n = a001263 (2 * n - 1) n  -- Reinhard Zumkeller, Oct 10 2013
  • Magma
    [Factorial(2*n)*Factorial(2*n+1) / (Factorial(n) * Factorial(n+1))^2: n in [0..20]]; // Vincenzo Librandi, Aug 15 2011
    
  • Maple
    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
  • Mathematica
    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 *)
  • PARI
    {a(n) = binomial(2*n+1, n)^2 / (2*n + 1)}; /* Michael Somos, Jun 22 2005 */
    
  • PARI
    a(n) = matdet(matrix(n, n, i, j, binomial(n+j+1,i+1))) \\ Hugo Pfoertner, Oct 22 2022
    

Formula

-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

Extensions

More terms from Andrew V. Sutherland, Mar 24 2008

A243752 Number T(n,k) of Dyck paths of semilength n having exactly k (possibly overlapping) occurrences of the consecutive step pattern given by the binary expansion of n, where 1=U=(1,1) and 0=D=(1,-1); triangle T(n,k), n>=0, read by rows.

Original entry on oeis.org

1, 0, 1, 0, 1, 1, 1, 3, 1, 1, 11, 2, 9, 16, 12, 4, 1, 1, 57, 69, 5, 127, 161, 98, 35, 7, 1, 323, 927, 180, 1515, 1997, 1056, 280, 14, 4191, 5539, 3967, 1991, 781, 244, 64, 17, 1, 1, 10455, 25638, 18357, 4115, 220, 1, 20705, 68850, 77685, 34840, 5685, 246, 1
Offset: 0

Views

Author

Alois P. Heinz, Jun 09 2014

Keywords

Examples

			Triangle T(n,k) begins:
: n\k :    0     1     2     3    4    5  ...
+-----+----------------------------------------------------------
:  0  :    1;                                 [row  0 of A131427]
:  1  :    0,    1;                           [row  1 of A131427]
:  2  :    0,    1,    1;                     [row  2 of A090181]
:  3  :    1,    3,    1;                     [row  3 of A001263]
:  4  :    1,   11,    2;                     [row  4 of A091156]
:  5  :    9,   16,   12,    4,   1;          [row  5 of A091869]
:  6  :    1,   57,   69,    5;               [row  6 of A091156]
:  7  :  127,  161,   98,   35,   7,   1;     [row  7 of A092107]
:  8  :  323,  927,  180;                     [row  8 of A091958]
:  9  : 1515, 1997, 1056,  280,  14;          [row  9 of A135306]
: 10  : 4191, 5539, 3967, 1991, 781, 244, ... [row 10 of A094507]
		

Crossrefs

A185650 a(n) is the number of rooted trees with 2n vertices n of whom are leaves.

Original entry on oeis.org

1, 2, 8, 39, 214, 1268, 7949, 51901, 349703, 2415348, 17020341, 121939535, 885841162, 6511874216, 48359860685, 362343773669, 2736184763500, 20805175635077, 159174733727167, 1224557214545788, 9467861087020239, 73534456468877012, 573484090227222260
Offset: 1

Views

Author

Stepan Orevkov, Aug 29 2013

Keywords

Examples

			From _Gus Wiseman_, Nov 27 2022: (Start)
The a(1) = 1 through a(3) = 8 rooted trees:
  (o)  ((oo))  (((ooo)))
       (o(o))  ((o)(oo))
               ((o(oo)))
               ((oo(o)))
               (o((oo)))
               (o(o)(o))
               (o(o(o)))
               (oo((o)))
(End)
		

Crossrefs

The ordered version is A000891, ranked by A358579.
This is the central column of A055277.
These trees are ranked by A358578.
For height = internals we have A358587.
Square trees are counted by A358589.
A000081 counts rooted trees, ordered A000108.
A055277 counts rooted trees by nodes and leaves, ordered A001263.
A358575 counts rooted trees by nodes and internals, ordered A090181.

Programs

  • Mathematica
    terms = 23;
    m = 2 terms;
    T[, ] = 0;
    Do[T[x_, z_] = z x - x + x Exp[Sum[Series[1/k T[x^k, z^k], {x, 0, j}, {z, 0, j}], {k, 1, j}]] // Normal, {j, 1, m}];
    cc = CoefficientList[#, z]& /@ CoefficientList[T[x, z] , x];
    Table[cc[[2n+1, n+1]], {n, 1, terms}] (* Jean-François Alcover, Sep 14 2018 *)
    art[n_]:=If[n==1,{{}},Join@@Table[Select[Tuples[art/@c],OrderedQ],{c,Join@@Permutations/@IntegerPartitions[n-1]}]];
    Table[Length[Select[art[n],Count[#,{},{-2}]==n/2&]],{n,2,10,2}] (* Gus Wiseman, Nov 27 2022 *)
  • PARI
    \\ here R is A055277 as vector of polynomials
    R(n) = {my(A = O(x)); for(j=1, n, A = x*(y - 1  + exp( sum(i=1, j, 1/i * subst( subst( A + x * O(x^(j\i)), x, x^i), y, y^i) ) ))); Vec(A)};
    {my(A=R(2*30)); vector(#A\2, k, polcoeff(A[2*k],k))} \\ Andrew Howroyd, May 21 2018

Extensions

Terms a(20) and beyond from Andrew Howroyd, May 21 2018

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

Views

Author

Len Smiley, Nov 11 2001

Keywords

Comments

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

Examples

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

Crossrefs

Cf. A003150 (for analog with consecutive Fibonacci numbers).

Programs

  • Magma
    [Binomial(4*n-1, 2*n-1)/(2*n+1): n in [1..20]]; // Vincenzo Librandi, Dec 09 2015
  • Maple
    seq(binomial(4*n-1,2*n-1)/(2*n+1), n=0..30); # Robert Israel, Dec 08 2015
  • Mathematica
    a[ n_] := If[ n < 1, 0, Binomial[ 4 n - 1, 2 n - 1] / (2 n + 1)]; (* Michael Somos, Oct 25 2014 *)
  • MuPAD
    combinat::dyckWords::count(2*n)/2 $ n = 1..26 // Zerinvary Lajos, Apr 25 2007
    
  • PARI
    a(n) = { if(n==0, 1, (4*n - 1)!/((2*n + 1)!*(2*n - 1)!)) } \\ Harry J. Smith, Oct 07 2009
    
  • PARI
    vector(20, n, binomial(4*n-1, 2*n-1)/(2*n+1)) \\ Altug Alkan, Dec 08 2015
    
  • Sage
    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
    

Formula

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)

Extensions

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

A060693 Triangle (0 <= k <= n) read by rows: T(n, k) is the number of Schröder paths from (0,0) to (2n,0) having k peaks.

Original entry on oeis.org

1, 1, 1, 2, 3, 1, 5, 10, 6, 1, 14, 35, 30, 10, 1, 42, 126, 140, 70, 15, 1, 132, 462, 630, 420, 140, 21, 1, 429, 1716, 2772, 2310, 1050, 252, 28, 1, 1430, 6435, 12012, 12012, 6930, 2310, 420, 36, 1, 4862, 24310, 51480, 60060, 42042, 18018, 4620, 660, 45, 1, 16796
Offset: 0

Views

Author

F. Chapoton, Apr 20 2001

Keywords

Comments

The rows sum to A006318 (Schroeder numbers), the left column is A000108 (Catalan numbers); the next-to-left column is A001700, the alternating sum in each row but the first is 0.
T(n,k) is the number of Schroeder paths (i.e., consisting of steps U=(1,1), D=(1,-1), H=(2,0) and never going below the x-axis) from (0,0) to (2n,0), having k peaks. Example: T(2,1)=3 because we have UU*DD, U*DH and HU*D, the peaks being shown by *. E.g., T(n,k) = binomial(n,k)*binomial(2n-k,n-1)/n for n>0. - Emeric Deutsch, Dec 06 2003
A090181*A007318 as infinite lower triangular matrices. - Philippe Deléham, Oct 14 2008
T(n,k) is also the number of rooted plane trees with maximal degree 3 and k vertices of degree 2 (a node may have at most 2 children, and there are exactly k nodes with 1 child). Equivalently, T(n,k) is the number of syntactically different expressions that can be formed that use a unary operation k times, a binary operation n-k times, and nothing else (sequence of operands is fixed). - Lars Hellstrom (Lars.Hellstrom(AT)residenset.net), Dec 08 2009

Examples

			Triangle begins:
00: [    1]
01: [    1,     1]
02: [    2,     3,      1]
03: [    5,    10,      6,      1]
04: [   14,    35,     30,     10,      1]
05: [   42,   126,    140,     70,     15,      1]
06: [  132,   462,    630,    420,    140,     21,     1]
07: [  429,  1716,   2772,   2310,   1050,    252,    28,    1]
08: [ 1430,  6435,  12012,  12012,   6930,   2310,   420,   36,   1]
09: [ 4862, 24310,  51480,  60060,  42042,  18018,  4620,  660,  45,  1]
10: [16796, 92378, 218790, 291720, 240240, 126126, 42042, 8580, 990, 55, 1]
...
		

Crossrefs

Triangle in A088617 transposed.
T(2n,n) gives A007004.

Programs

  • Maple
    A060693 := (n,k) -> binomial(n,k)*binomial(2*n-k,n)/(n-k+1); # Peter Luschny, May 17 2011
  • Mathematica
    t[n_, k_] := Binomial[n, k]*Binomial[2 n - k, n]/(n - k + 1); Flatten[Table[t[n, k], {n, 0, 9}, {k, 0, n}]] (* Robert G. Wilson v, May 30 2011 *)
  • PARI
    T(n, k) = binomial(n, k)*binomial(2*n - k, n)/(n - k + 1);
    for(n=0, 10, for(k=0, n, print1(T(n, k),", ")); print); \\ Indranil Ghosh, Jul 28 2017
    
  • Python
    from sympy import binomial
    def T(n, k): return binomial(n, k) * binomial(2 * n - k, n) / (n - k + 1)
    for n in range(11): print([T(n, k) for k in range(n + 1)])  # Indranil Ghosh, Jul 28 2017

Formula

Triangle T(n, k) (0 <= k <= n) read by rows; given by [1, 1, 1, 1, 1, ...] DELTA [1, 0, 1, 0, 1, 0, ...] where DELTA is the operator defined in A084938. - Philippe Deléham, Aug 12 2003
If C_n(x) is the g.f. of row n of the Narayana numbers (A001263), C_n(x) = Sum_{k=1..n} binomial(n,k-1)*(binomial(n-1,k-1)/k) * x^k and T_n(x) is the g.f. of row n of T(n,k), then T_n(x) = C_n(x+1), or T(n,k) = [x^n]Sum_{k=1..n}(A001263(n,k)*(x+1)^k). - Mitch Harris, Jan 16 2007, Jan 31 2007
G.f.: (1 - t*y - sqrt((1-y*t)^2 - 4*y)) / 2.
T(n, k) = binomial(2n-k, n)*binomial(n, k)/(n-k+1). - Philippe Deléham, Dec 07 2003
A060693(n, k) = binomial(2*n-k, k)*A000108(n-k); A000108: Catalan numbers. - Philippe Deléham, Dec 30 2003
Sum_{k=0..n} T(n,k)*x^k = A000007(n), A000108(n), A006318(n), A047891(n+1), A082298(n), A082301(n), A082302(n), A082305(n), A082366(n), A082367(n), for x = -1, 0, 1, 2, 3, 4, 5, 6, 7, 8, respectively. - Philippe Deléham, Apr 01 2007
T(n,k) = Sum_{j>=0} A090181(n,j)*binomial(j,k). - Philippe Deléham, May 04 2007
Sum_{k=0..n} T(n,k)*x^(n-k) = (-1)^n*A107841(n), A080243(n), A000007(n), A000012(n), A006318(n), A103210(n), A103211(n), A133305(n), A133306(n), A133307(n), A133308(n), A133309(n) for x = -3, -2, -1, 0, 1, 2, 3, 4, 5, 6, 7, 8, respectively. - Philippe Deléham, Oct 18 2007
From Paul Barry, Jan 29 2009: (Start)
G.f.: 1/(1-xy-x/(1-xy-x/(1-xy-x/(1-xy-x/(1-xy-x/(1-.... (continued fraction);
G.f.: 1/(1-(x+xy)/(1-x/(1-(x+xy)/(1-x/(1-(x+xy)/(1-.... (continued fraction). (End)
T(n,k) = [k<=n]*(Sum_{j=0..n} binomial(n,j)^2*binomial(j,k))/(n-k+1). - Paul Barry, May 28 2009
T(n,k) = A104684(n,k)/(n-k+1). - Peter Luschny, May 17 2011
From Tom Copeland, Sep 21 2011: (Start)
With F(x,t) = (1-(2+t)*x-sqrt(1-2*(2+t)*x+(t*x)^2))/(2*x) an o.g.f. (nulling the n=0 term) in x for the A060693 polynomials in t,
G(x,t) = x/(1+t+(2+t)*x+x^2) is the compositional inverse in x.
Consequently, with H(x,t) = 1/(dG(x,t)/dx) = (1+t+(2+t)*x+x^2)^2 / (1+t-x^2), the n-th A060693 polynomial in t is given by (1/n!)*((H(x,t)*d/dx)^n) x evaluated at x=0, i.e., F(x,t) = exp(x*H(u,t)*d/d) u, evaluated at u = 0.
Also, dF(x,t)/dx = H(F(x,t),t). (End)
See my 2008 formulas in A033282 to relate this entry to A088617, A001263, A086810, and other matrices. - Tom Copeland, Jan 22 2016
Rows of this entry are non-vanishing antidiagonals of A097610. See p. 14 of Agapito et al. for a bivariate generating function and its inverse. - Tom Copeland, Feb 03 2016
From Werner Schulte, Jan 09 2017: (Start)
T(n,k) = A126216(n,k-1) + A126216(n,k) for 0 < k < n;
Sum_{k=0..n} (-1)^k*(1+x*(n-k))*T(n,k) = x + (1-x)*A000007(n).
(End)
Conjecture: Sum_{k=0..n} (-1)^k*T(n,k)*(n+1-k)^2 = 1+n+n^2. - Werner Schulte, Jan 11 2017

Extensions

More terms from Vladeta Jovovic, Apr 21 2001
New description from Philippe Deléham, Aug 12 2003
New name using a comment by Emeric Deutsch from Peter Luschny, Jul 26 2017

A034867 Triangle of odd-numbered terms in rows of Pascal's triangle.

Original entry on oeis.org

1, 2, 3, 1, 4, 4, 5, 10, 1, 6, 20, 6, 7, 35, 21, 1, 8, 56, 56, 8, 9, 84, 126, 36, 1, 10, 120, 252, 120, 10, 11, 165, 462, 330, 55, 1, 12, 220, 792, 792, 220, 12, 13, 286, 1287, 1716, 715, 78, 1, 14, 364, 2002, 3432, 2002, 364, 14, 15, 455, 3003, 6435, 5005, 1365, 105, 1
Offset: 0

Views

Author

Keywords

Comments

Also triangle of numbers of n-sequences of 0,1 with k subsequences of consecutive 01 because this number is C(n+1,2*k+1). - Roger Cuculiere (cuculier(AT)imaginet.fr), Nov 16 2002
From Gary W. Adamson, Oct 17 2008: (Start)
Received from Herb Conn:
Let T = tan x, then
tan x = T
tan 2x = 2T / (1 - T^2)
tan 3x = (3T - T^3) / (1 - 3T^2)
tan 4x = (4T - 4T^3) / (1 - 6T^2 + T^4)
tan 5x = (5T - 10T^3 + T^5) / (1 - 10T^2 + 5T^4)
tan 6x = (6T - 20T^3 + 6T^5) / (1 - 15T^2 + 15T^4 - T^6)
tan 7x = (7T - 35T^3 + 21T^5 - T^7) / (1 - 21T^2 + 35T^4 - 7T^6)
tan 8x = (8T - 56T^3 + 56T^5 - 8T^7) / (1 - 28T^2 + 70T^4 - 28T^6 + T^8)
tan 9x = (9T - 84T^3 + 126T^5 - 36T^7 + T^9) / (1 - 36 T^2 + 126T^4 - 84T^6 + 9T^8)
... To get the next one in the series, (tan 10x), for the numerator add:
9....84....126....36....1 previous numerator +
1....36....126....84....9 previous denominator =
10..120....252...120...10 = new numerator
For the denominator add:
......9.....84...126...36...1 = previous numerator +
1....36....126....84....9.... = previous denominator =
1....45....210...210...45...1 = new denominator
...where numerators = A034867, denominators = A034839
(End)
Column k is the sum of columns 2k and 2k+1 of A007318. - Philippe Deléham, Nov 12 2008
Triangle, with zeros omitted, given by (2, -1/2, 1/2, 0, 0, 0, 0, 0, 0, 0, ...) DELTA (0, 1/2, -1/2, 0, 0, 0, 0, 0, 0, 0, ...) where DELTA is the operator defined in A084938. - Philippe Deléham, Dec 12 2011
The row polynomials N(n,x) = Sum_{k=0..floor((n-1)/2)} T(n-1,k)*x^k, and D(n,x) = Sum_{k=0..floor(n/2)} A034839(n,k)*x^k, n >= 1, satisfy the recurrences N(n,x) = D(n-1,x) + N(n-1,x), D(n,x) = D(n-1,x) + x*N(n-1,x), with inputs N(1,x) = 1 = D(1,x). This is due to the Pascal triangle A007318 recurrence. Q(n,x) := tan(n*x)/tan(x) satisfies the recurrence Q(n,x) = (1 + Q(n-1,x))/(1 - v(x)*Q(n-1,x)) with input Q(1,x) = 1 and v = v(x) := (tan(x))^2. This recurrence is obtained from the addition theorem for tan(n*x) using n = 1 + (n-1). Therefore Q(n,x) = N(n,-v(x))/D(n,-v(x)). This proves the Gary W. Adamson contribution from above. See also A220673. This calculation was motivated by an e-mail of Thomas Olsen. The Oliver/Prodinger and Ma references resort to HAKEM Al Memo 239, Item 16, for the tan(n*x) formula in terms of tan(x). - Wolfdieter Lang, Jan 17 2013
The infinitesimal generator (infinigen) for the Narayana polynomials A090181/A001263 can be formed from the row polynomials P(n,y) of this entry. The resulting matrix is an instance of a matrix representation of the analytic infinigens presented in A145271 for general sets of binomial Sheffer polynomials and in A001263 and A119900 specifically for the Narayana polynomials. Given the column vector of row polynomials V = (1, P(1,x) = 2x, P(2,y) = 3x + x^2, P(3,y) = 4x + 4x^2, ...), form the lower triangular matrix M(n,k) = V(n-k,n-k), i.e., diagonally multiply the matrix with all ones on the diagonal and below by the components of V. Form the matrix MD by multiplying A132440^Transpose = A218272 = D (representing derivation of o.g.f.s) by M, i.e., MD = M*D. The non-vanishing component of the first row of (MD)^n * V / (n+1)! is the n-th Narayana polynomial. - Tom Copeland, Dec 09 2015
The diagonals of this entry are A078812 (also shifted A128908 and unsigned A053122, which are embedded in A030528, A102426, A098925, A109466, A092865). Equivalently, the antidiagonals of A078812 are the rows of A034867. - Tom Copeland, Dec 12 2015
Binomial(n,2k+1) is also the number of permutations avoiding both 132 and 213 with k peaks, i.e., positions with w[i]w[i+2]. - Lara Pudwell, Dec 19 2018
Binomial(n,2k+1) is also the number of permutations avoiding both 123 and 132 with k peaks, i.e., positions with w[i]w[i+2]. - Lara Pudwell, Dec 19 2018
The row polynomial P(n, x) = Sum_{0..floor(n/2)} T(n, k)*x^k appears as numerator polynomial of the diagonal sequence m of triangle A104698 as follows. G(m, x) = P(m, x^2)/(1 - x)^(m+1), for m >= 0. - Wolfdieter Lang, May 14 2025
Number of acyclic orientations of the path graph on n+1 vertices, with k-1 sinks. - Per W. Alexandersson, Aug 15 2025

Examples

			Triangle T starts:
  n\k   0   1   2   3   4  5 ...   ----------------------------------------
0:    1
1:    2
2:    3   1
3:    4   4
4:    5  10   1
5:    6  20   6
6:    7  35  21   1
7:    8  56  56   8
8:    9  84 126  36   1
9:   10 120 252 120  10
 10:   11 165 462 330  55  1
 11:   12 220 792 792 220 12
... ... reformatted and extended by - _Wolfdieter Lang_, May 14 2025
		

References

  • A. T. Benjamin and J. J. Quinn, Proofs that really count: the art of combinatorial proof, M.A.A. 2003, id. 136.

Crossrefs

From Wolfdieter Lang, May 14 2025:(Start)
Row length A008619. Row sums A000079. Alternating row sums A009545(n+1).
Column sequences (with certain offsets): A000027, A000292, A000389, A000580, A000582, A001288, ... (End)

Programs

  • Magma
    /* as a triangle */ [[Binomial(n+1,2*k+1): k in [0..Floor(n/2)]]: n in [0..20]]; // G. C. Greubel, Mar 06 2018
  • Maple
    seq(seq(binomial(n+1,2*k+1), k=0..floor(n/2)), n=0..14); # Emeric Deutsch, Apr 01 2005
  • Mathematica
    u[1, x_] := 1; v[1, x_] := 1; z = 12;
    u[n_, x_] := u[n - 1, x] + x*v[n - 1, x]
    v[n_, x_] := u[n - 1, x] + v[n - 1, x]
    cu = Table[CoefficientList[u[n, x], x], {n, 1, z}];
    TableForm[cu]  (* A034839 as a triangle *)
    cv = Table[CoefficientList[v[n, x], x], {n, 1, z}];
    TableForm[cv]  (* A034867 as a triangle *)
    (* Clark Kimberling, Feb 18 2012 *)
    Table[Binomial[n+1, 2*k+1], {n,0,20}, {k,0,Floor[n/2]}]//Flatten (* G. C. Greubel, Mar 06 2018 *)
  • PARI
    for(n=0,20, for(k=0,floor(n/2), print1(binomial(n+1,2*k+1), ", "))) \\ G. C. Greubel, Mar 06 2018
    

Formula

T(n,k) = C(n+1,2k+1) = Sum_{i=k..n-k} C(i,k) * C(n-i,k).
E.g.f.: 1+(exp(x)*sinh(x*sqrt(y)))/sqrt(y). - Vladeta Jovovic, Mar 20 2005
G.f.: 1/((1-z)^2-t*z^2). - Emeric Deutsch, Apr 01 2005
T(n,k) = Sum_{j = 0..n} A034839(j,k). - Philippe Deléham, May 18 2005
Pell(n+1) = A000129(n+1) = Sum_{k=0..n} T(n,k) * 2^k = (1/n!) Sum_{k=0..n} A131980(n,k) * 2^k. - Tom Copeland, Nov 30 2007
T(n,k) = A007318(n,2k) + A007318(n,2k+1). - Philippe Deléham, Nov 12 2008
O.g.f for column k, k>=0: (1/(1-x)^2)*(x/(1-x))^(2*k). See the G.f. of this array given above by Emeric Deutsch. - Wolfdieter Lang, Jan 18 2013
T(n,k) = (x^(2*k+1))*((1+x)^n-(1-x)^n)/2. - L. Edson Jeffery, Jan 15 2014

Extensions

More terms from Emeric Deutsch, Apr 01 2005

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

Views

Author

Gus Wiseman, Nov 23 2022

Keywords

Comments

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

Examples

			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))
		

Crossrefs

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.

Programs

  • Mathematica
    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}]
  • PARI
    \\ 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

Extensions

Terms a(19) and beyond from Andrew Howroyd, Jan 01 2023
Showing 1-10 of 41 results. Next