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

A003480 a(0) = 1, a(1) = 2, for n > 1, a(n) = 4*a(n-1) - 2*a(n-2).

Original entry on oeis.org

1, 2, 7, 24, 82, 280, 956, 3264, 11144, 38048, 129904, 443520, 1514272, 5170048, 17651648, 60266496, 205762688, 702517760, 2398545664, 8189147136, 27959497216, 95459694592, 325919783936, 1112759746560, 3799199418368, 12971278180352, 44286713884672, 151204299177984
Offset: 0

Views

Author

Keywords

Comments

Gives the number of L-convex polyominoes with n cells, that is convex polyominoes where any two cells can be connected by a path internal to the polyomino and which has at most 1 change of direction (i.e., one of the four orientation of the L). - Simone Rinaldi (rinaldi(AT)unisi.it), Feb 19 2007
Joe Keane (jgk(AT)jgk.org) observes that this sequence (beginning at 2) is "size of raises in pot-limit poker, one blind, maximum raising".
Dimensions of the graded components of the Hopf algebra of noncommutative multi-symmetric functions of level 2. For level r, the sequence would be the INVERT transform of binomial(n+r-1,n). - Jean-Yves Thibon (jyt(AT)univ-mlv.fr), Jun 26 2008
The sum of the numbers in the n-th row of the summatory Pascal triangle (A059576). - Ron R. King, Jan 22 2009
(1 + 2x + 7x^2 + 24x^3 + ...) = 1 / (1 - 2x - 3x^2 - 4x^3 - ...). - Gary W. Adamson, Jul 27 2009
Let M be a triangle with the odd-indexed Fibonacci numbers (1, 2, 5, 13, ...) in every column, with the leftmost column shifted upwards one row. A003480 = lim_{n->oo} M^n, the left-shifted vector considered as a sequence. The analogous operation using the even-indexed Fibonacci numbers generates A001835 starting with offset 1. - Gary W. Adamson, Jul 27 2010
a(n) is the number of generalized compositions of n when there are i+1 different types of the part i, (i=1,2,...). - Milan Janjic, Sep 24 2010
Let h(t) = (1-t)^2/(2*(1-t)^2-1) = 1/(1-(2*t + 3*t^2 + 4*t^3 + ...)),
an o.g.f. for A003480, then
A001003(n) = (1/n!)*((h(t)*d/dt)^n) t, evaluated at t=0, with initial n=1. - Tom Copeland, Sep 06 2011
Excluding the initial 1, a(n) is the 2nd subdiagonal of A228405. - Richard R. Forberg, Sep 02 2013

References

  • G. Castiglione and A. Restivo, L-convex polyominoes: a survey, Chapter 2 of K. G. Subranian et al., eds., Formal Models, Languages and Applications, World Scientific, 2015.
  • N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).

Crossrefs

Row sums of A059576 and of A181289. Second differences of A007070.
Column k=2 of A261780.

Programs

  • Haskell
    a003480 n = a003480_list !! n
    a003480_list = 1 : 2 : 7 : (tail $ zipWith (-)
       (tail $ map (* 4) a003480_list) (map (* 2) a003480_list))
    -- Reinhard Zumkeller, Jan 16 2012, Oct 03 2011
  • Maple
    INVERT([seq(n+1,n=1..20)]); # Jean-Yves Thibon (jyt(AT)univ-mlv.fr), Jun 26 2008
  • Mathematica
    a[0]=1; a[1]=2; a[2]=7; a[n_]:=a[n]=4*a[n-1] - 2*a[n-2]; Table[a[n],{n,0,24}] (* Jean-François Alcover, Mar 22 2011 *)
    Join[{1},LinearRecurrence[{4,-2},{2,7},40]] (* Harvey P. Dale, Oct 23 2011 *)
  • PARI
    a(n)=polcoeff((1-x)^2/(1-4*x+2*x^2)+x*O(x^n),n)
    
  • PARI
    a(n)=local(x); if(n<1,n==0,x=(2+quadgen(8))^n; imag(x)+real(x)/2)
    

Formula

a(n) = (n+1)*a(0) + n*a(1) + ... + 3*a(n-2) + 2*a(n-1). - Amarnath Murthy, Aug 17 2002
G.f.: (1-x)^2/(1-4*x+2*x^2). - Simon Plouffe in his 1992 dissertation
a(n) = A007070(n)/2, n > 0.
G.f.: 1/( 1 - Sum_{k>=1} (k+1)*x^k ).
a(n+1)*a(n+1) - a(n+2)*a(n) = 2^n, n > 0. - D. G. Rogers, Jul 12 2004
For n > 0, a(n) = ((2+sqrt(2))^(n+1) - (2-sqrt(2))^(n+1))/(4*sqrt(2)). - Rolf Pleisch, Aug 03 2009
If the leading 1 is removed, 2, 7, 24, ... is the binomial transform of 2, 5, 12, 29, ..., which is A000129 without its first 2 terms, and the second binomial transform of 2, 3, 4, 6, ..., which is A029744, again without its leading 1. - Al Hakanson (hawkuu(AT)gmail.com), Aug 08 2009
a(n) = Sum((1+p_1)*(1+p_2)*...*(1+p_m)), summation being over all compositions (p_1, p_2, ..., p_m) of n. Example: a(3)=24; indeed, the compositions of 3 are (1,1,1), (1,2), (2,1), (3) and we have 2*2*2 + 2*3 + 3*2 + 4 = 24. - Emeric Deutsch, Oct 17 2010
a(n) = Sum_{k>=0} binomial(n+2*k-1,n) / 2^(k+1). - Vaclav Kotesovec, Dec 31 2013
E.g.f.: (1 + exp(2*x)*(cosh(sqrt(2)*x) + sqrt(2)*sinh(sqrt(2)*x)))/2. - Stefano Spezia, May 20 2024

A145271 Coefficients for expansion of (g(x)d/dx)^n g(x); refined Eulerian numbers for calculating compositional inverse of h(x) = (d/dx)^(-1) 1/g(x); iterated derivatives as infinitesimal generators of flows.

Original entry on oeis.org

1, 1, 1, 1, 1, 4, 1, 1, 11, 4, 7, 1, 1, 26, 34, 32, 15, 11, 1, 1, 57, 180, 122, 34, 192, 76, 15, 26, 16, 1, 1, 120, 768, 423, 496, 1494, 426, 294, 267, 474, 156, 56, 42, 22, 1, 1, 247, 2904, 1389, 4288, 9204, 2127, 496, 5946, 2829, 5142, 1206, 855, 768, 1344, 1038, 288, 56, 98, 64, 29, 1
Offset: 0

Views

Author

Tom Copeland, Oct 06 2008

Keywords

Comments

For more detail, including connections to Legendre transformations, rooted trees, A139605, A139002 and A074060, see Mathemagical Forests p. 9.
For connections to the h-polynomials associated to the refined f-polynomials of permutohedra see my comments in A008292 and A049019.
From Tom Copeland, Oct 14 2011: (Start)
Given analytic functions F(x) and FI(x) such that F(FI(x))=FI(F(x))=x about 0, i.e., they are compositional inverses of each other, then, with g(x) = 1/dFI(x)/dx, a flow function W(s,x) can be defined with the following relations:
W(s,x) = exp(s g(x)d/dx)x = F(s+FI(x)) ,
W(s,0) = F(s) ,
W(0,x) = x ,
dW(0,x)/ds = g(x) = F'[FI(x)] , implying
dW(0,F(x))/ds = g(F(x)) = F'(x) , and
W(s,W(r,x)) = F(s+FI(F(r+FI(x)))) = F(s+r+FI(x)) = W(s+r,x) . (See MF link below.) (End)
dW(s,x)/ds - g(x)dW(s,x)/dx = 0, so (1,-g(x)) are the components of a vector orthogonal to the gradient of W and, therefore, tangent to the contour of W, at (s,x) . - Tom Copeland, Oct 26 2011
Though A139605 contains A145271, the op. of A145271 contains that of A139605 in the sense that exp(s g(x)d/dx) w(x) = w(F(s+FI(x))) = exp((exp(s g(x)d/dx)x)d/du)w(u) evaluated at u=0. This is reflected in the fact that the forest of rooted trees assoc. to (g(x)d/dx)^n, FOR_n, can be generated by removing the single trunk of the planted rooted trees of FOR_(n+1). - Tom Copeland, Nov 29 2011
Related to formal group laws for elliptic curves (see Hoffman). - Tom Copeland, Feb 24 2012
The functional equation W(s,x) = F(s+FI(x)), or a restriction of it, is sometimes called the Abel equation or Abel's functional equation (see Houzel and Wikipedia) and is related to Schröder's functional equation and Koenigs functions for compositional iterates (Alexander, Goryainov and Kudryavtseva). - Tom Copeland, Apr 04 2012
g(W(s,x)) = F'(s + FI(x)) = dW(s,x)/ds = g(x) dW(s,x)/dx, connecting the operators here to presentations of the Koenigs / Königs function and Loewner / Löwner evolution equations of the Contreras et al. papers. - Tom Copeland, Jun 03 2018
The autonomous differential equation above also appears with a change in variable of the form x = log(u) in the renormalization group equation, or Beta function. See Wikipedia, Zinn-Justin equations 2.10 and 3.11, and Krajewski and Martinetti equation 21. - Tom Copeland, Jul 23 2020
A variant of these partition polynomials appears on p. 83 of Petreolle et al. with the indeterminates e_n there related to those given in the examples below by e_n = n!*(n'). The coefficients are interpreted as enumerating certain types of trees. See also A190015. - Tom Copeland, Oct 03 2022

Examples

			From _Tom Copeland_, Sep 19 2014: (Start)
Let h(x) = log((1+a*x)/(1+b*x))/(a-b); then, g(x) = 1/(dh(x)/dx) = (1+ax)(1+bx), so (0')=1, (1')=a+b, (2')=2ab, evaluated at x=0, and higher order derivatives of g(x) vanish. Therefore, evaluated at x=0,
R^0 g(x) =  1
R^1 g(x) =  a+b
R^2 g(x) = (a+b)^2 + 2ab = a^2 + 4 ab + b^2
R^3 g(x) = (a+b)^3 + 4*(a+b)*2ab = a^3 + 11 a^2*b + 11 ab^2 + b^3
R^4 g(x) = (a+b)^4 + 11*(a+b)^2*2ab + 4*(2ab)^2
         =  a^4 + 26 a^3*b + 66 a^2*b^2 + 26 ab^3 + b^4,
etc., and these bivariate Eulerian polynomials (A008292) are the first few coefficients of h^(-1)(x) = (e^(ax) - e^(bx))/(a*e^(bx) - b*e^(ax)), the inverse of h(x). (End)
Triangle starts:
  1;
  1;
  1,   1;
  1,   4,    1;
  1,  11,    4,    7,    1;
  1,  26,   34,   32,   15,   11,    1;
  1,  57,  180,  122,   34,  192,   76,  15,   26,   16,    1;
  1, 120,  768,  423,  496, 1494,  426, 294,  267,  474,  156,   56,  42,  22,    1;
  1, 247, 2904, 1389, 4288, 9204, 2127, 496, 5946, 2829, 5142, 1206, 855, 768, 1344, 1038, 288, 56, 98, 64, 29, 1;
		

References

  • D. S. Alexander, A History of Complex Dynamics: From Schröder to Fatou to Julia, Friedrich Vieweg & Sohn, 1994.
  • T. Mansour and M. Schork, Commutation Relations, Normal Ordering, and Stirling Numbers, Chapman and Hall/CRC, 2015.

Crossrefs

Cf. (A133437, A086810, A181289) = (LIF, reduced LIF, associated g(x)), where LIF is a Lagrange inversion formula. Similarly for (A134264, A001263, A119900), (A134685, A134991, A019538), (A133932, A111999, A007318).
Second column is A000295, subdiagonal is A000124, row sums are A000142, row lengths are A000041. - Peter Luschny, Jul 21 2016

Programs

  • Maple
    with(LinearAlgebra): with(ListTools):
    A145271_row := proc(n) local b, M, V, U, G, R, T;
    if n < 2 then return 1 fi;
    b := (n,k) -> `if`(k=1 or k>n+1,0,binomial(n-1,k-2)*g[n-k+1]);
    M := n -> Matrix(n, b):
    V := n -> Vector[row]([1, seq(0,i=2..n)]):
    U := n -> VectorMatrixMultiply(V(n), M(n)^(n-1)):
    G := n -> Vector([seq(g[i], i=0..n-1)]);
    R := n -> VectorMatrixMultiply(U(n), G(n)):
    T := Reverse([op(sort(expand(R(n+1))))]);
    seq(subs({seq(g[i]=1, i=0..n)},T[j]),j=1..nops(T)) end:
    for n from 0 to 9 do A145271_row(n) od; # Peter Luschny, Jul 20 2016

Formula

Let R = g(x)d/dx; then
R^0 g(x) = 1 (0')^1
R^1 g(x) = 1 (0')^1 (1')^1
R^2 g(x) = 1 (0')^1 (1')^2 + 1 (0')^2 (2')^1
R^3 g(x) = 1 (0')^1 (1')^3 + 4 (0')^2 (1')^1 (2')^1 + 1 (0')^3 (3')^1
R^4 g(x) = 1 (0')^1 (1')^4 + 11 (0')^2 (1')^2 (2')^1 + 4 (0')^3 (2')^2 + 7 (0')^3 (1')^1 (3')^1 + 1 (0')^4 (4')^1
R^5 g(x) = 1 (0') (1')^5 + 26 (0')^2 (1')^3 (2') + (0')^3 [34 (1') (2')^2 + 32 (1')^2 (3')] + (0')^4 [ 15 (2') (3') + 11 (1') (4')] + (0')^5 (5')
R^6 g(x) = 1 (0') (1')^6 + 57 (0')^2 (1')^4 (2') + (0')^3 [180 (1')^2 (2')^2 + 122 (1')^3 (3')] + (0')^4 [ 34 (2')^3 + 192 (1') (2') (3') + 76 (1')^2 (4')] + (0')^5 [15 (3')^2 + 26 (2') (4') + 16 (1') (5')] + (0')^6 (6')
where (j')^k = ((d/dx)^j g(x))^k. And R^(n-1) g(x) evaluated at x=0 is the n-th Taylor series coefficient of the compositional inverse of h(x) = (d/dx)^(-1) 1/g(x), with the integral from 0 to x.
The partitions are in reverse order to those in Abramowitz and Stegun p. 831. Summing over coefficients with like powers of (0') gives A008292.
Confer A190015 for another way to compute numbers for the array for each partition. - Tom Copeland, Oct 17 2014
Equivalent matrix computation: Multiply the n-th diagonal (with n=0 the main diagonal) of the lower triangular Pascal matrix by g_n = (d/dx)^n g(x) to obtain the matrix VP with VP(n,k) = binomial(n,k) g_(n-k). Then R^n g(x) = (1, 0, 0, 0, ...) [VP * S]^n (g_0, g_1, g_2, ...)^T, where S is the shift matrix A129185, representing differentiation in the divided powers basis x^n/n!. - Tom Copeland, Feb 10 2016 (An evaluation removed by author on Jul 19 2016. Cf. A139605 and A134685.)
Also, R^n g(x) = (1, 0, 0, 0, ...) [VP * S]^(n+1) (0, 1, 0, ...)^T in agreement with A139605. - Tom Copeland, Jul 21 2016
A recursion relation for computing each partition polynomial of this entry from the lower order polynomials and the coefficients of the cycle index polynomials of A036039 is presented in the blog entry "Formal group laws and binomial Sheffer sequences". - Tom Copeland, Feb 06 2018
A formula for computing the polynomials of each row of this matrix is presented as T_{n,1} on p. 196 of the Ihara reference in A139605. - Tom Copeland, Mar 25 2020
Indeterminate substitutions as illustrated in A356145 lead to [E] = [L][P] = [P][E]^(-1)[P] = [P][RT] and [E]^(-1) = [P][L] = [P][E][P] = [RT][P], where [E] contains the refined Eulerian partition polynomials of this entry; [E]^(-1), A356145, the inverse set to [E]; [P], the permutahedra polynomials of A133314; [L], the classic Lagrange inversion polynomials of A134685; and [RT], the reciprocal tangent polynomials of A356144. Since [L]^2 = [P]^2 = [RT]^2 = [I], the substitutional identity, [L] = [E][P] = [P][E]^(-1) = [RT][P], [RT] = [E]^(-1)[P] = [P][L][P] = [P][E], and [P] = [L][E] = [E][RT] = [E]^(-1)[L] = [RT][E]^(-1). - Tom Copeland, Oct 05 2022

Extensions

Title amplified by Tom Copeland, Mar 17 2014
R^5 and R^6 formulas and terms a(19)-a(29) added by Tom Copeland, Jul 11 2016
More terms from Peter Luschny, Jul 20 2016

A033282 Triangle read by rows: T(n, k) is the number of diagonal dissections of a convex n-gon into k+1 regions.

Original entry on oeis.org

1, 1, 2, 1, 5, 5, 1, 9, 21, 14, 1, 14, 56, 84, 42, 1, 20, 120, 300, 330, 132, 1, 27, 225, 825, 1485, 1287, 429, 1, 35, 385, 1925, 5005, 7007, 5005, 1430, 1, 44, 616, 4004, 14014, 28028, 32032, 19448, 4862, 1, 54, 936, 7644, 34398, 91728, 148512, 143208, 75582, 16796
Offset: 3

Views

Author

Keywords

Comments

T(n+3, k) is also the number of compatible k-sets of cluster variables in Fomin and Zelevinsky's cluster algebra of finite type A_n. Take a row of this triangle regarded as a polynomial in x and rewrite as a polynomial in y := x+1. The coefficients of the polynomial in y give a row of the triangle of Narayana numbers A001263. For example, x^2 + 5*x + 5 = y^2 + 3*y + 1. - Paul Boddington, Mar 07 2003
Number of standard Young tableaux of shape (k+1,k+1,1^(n-k-3)), where 1^(n-k-3) denotes a sequence of n-k-3 1's (see the Stanley reference).
Number of k-dimensional 'faces' of the n-dimensional associahedron (see Simion, p. 168). - Mitch Harris, Jan 16 2007
Mirror image of triangle A126216. - Philippe Deléham, Oct 19 2007
For relation to Lagrange inversion or series reversion and the geometry of associahedra or Stasheff polytopes (and other combinatorial objects) see A133437. - Tom Copeland, Sep 29 2008
Row generating polynomials 1/(n+1)*Jacobi_P(n,1,1,2*x+1). Row n of this triangle is the f-vector of the simplicial complex dual to an associahedron of type A_n [Fomin & Reading, p. 60]. See A001263 for the corresponding array of h-vectors for associahedra of type A_n. See A063007 and A080721 for the f-vectors for associahedra of type B and type D respectively. - Peter Bala, Oct 28 2008
f-vectors of secondary polytopes for Grobner bases for optimization and integer programming (see De Loera et al. and Thomas). - Tom Copeland, Oct 11 2011
From Devadoss and O'Rourke's book: The Fulton-MacPherson compactification of the configuration space of n free particles on a line segment with a fixed particle at each end is the n-Dim Stasheff associahedron whose refined f-vector is given in A133437 which reduces to A033282. - Tom Copeland, Nov 29 2011
Diagonals of A132081 are rows of A033282. - Tom Copeland, May 08 2012
The general results on the convolution of the refined partition polynomials of A133437, with u_1 = 1 and u_n = -t otherwise, can be applied here to obtain results of convolutions of these polynomials. - Tom Copeland, Sep 20 2016
The signed triangle t(n, k) =(-1)^k* T(n+2, k-1), n >= 1, k = 1..n, seems to be obtainable from the partition array A111785 (in Abramowitz-Stegun order) by adding the entries corresponding to the partitions of n with the number of parts k. E.g., triangle t, row n=4: -1, (6+3) = 9, -21, 14. - Wolfdieter Lang, Mar 17 2017
The preceding conjecture by Lang is true. It is implicit in Copeland's 2011 comments in A086810 on the relations among a gf and its compositional inverse for that entry and inversion through A133437 (a differently normalized version of A111785), whose integer partitions are the same as those for A134685. (An inversion pair in Copeland's 2008 formulas below can also be used to prove the conjecture.) In addition, it follows from the relation between the inversion formula of A111785/A133437 and the enumeration of distinct faces of associahedra. See the MathOverflow link concernimg Loday and the Aguiar and Ardila reference in A133437 for proofs of the relations between the partition polynomials for inversion and enumeration of the distinct faces of the A_n associahedra, or Stasheff polytopes. - Tom Copeland, Dec 21 2017
The rows seem to give (up to sign) the coefficients in the expansion of the integer-valued polynomial (x+1)*(x+2)^2*(x+3)^2*...*(x+n)^2*(x+n+1)/(n!*(n+1)!) in the basis made of the binomial(x+i,i). - F. Chapoton, Oct 07 2022
Chapoton's observation above is correct: the precise expansion is (x+1)*(x+2)^2*(x+3)^2*...*(x+n)^2*(x+n+1)/ (n!*(n+1)!) = Sum_{k = 0..n-1} (-1)^k*T(n+2,n-k-1)*binomial(x+2*n-k,2*n-k), as can be verified using the WZ algorithm. For example, n = 4 gives (x+1)*(x+2)^2*(x+3)^2*(x+4)^2*(x+5)/(4!*5!) = 14*binomial(x+8,8) - 21*binomial(x+7,7) + 9*binomial(x+6,6) - binomial(x+5,5). - Peter Bala, Jun 24 2023

Examples

			The triangle T(n, k) begins:
n\k  0  1   2    3     4     5      6      7     8     9
3:   1
4:   1  2
5:   1  5   5
6:   1  9  21   14
7:   1 14  56   84    42
8:   1 20 120  300   330   132
9:   1 27 225  825  1485  1287    429
10:  1 35 385 1925  5005  7007   5005   1430
11:  1 44 616 4004 14014 28028  32032  19448  4862
12:  1 54 936 7644 34398 91728 148512 143208 75582 16796
... reformatted. - _Wolfdieter Lang_, Mar 17 2017
		

References

  • S. Devadoss and J. O'Rourke, Discrete and Computational Geometry, Princeton Univ. Press, 2011 (See p. 241.)
  • Ronald L. Graham, Donald E. Knuth, Oren Patashnik, Concrete Mathematics, 2nd ed., Addison-Wesley, 1994. Exercise 7.50, pages 379, 573.
  • T. K. Petersen, Eulerian Numbers, Birkhauser, 2015, Section 5.8.

Crossrefs

Cf. diagonals: A000012, A000096, A033275, A033276, A033277, A033278, A033279; A000108, A002054, A002055, A002056, A007160, A033280, A033281; row sums: A001003 (Schroeder numbers, first term omitted). See A086810 for another version.
A007160 is a diagonal. Cf. A001263.
With leading zero: A086810.
Cf. A019538 'faces' of the permutohedron.
Cf. A063007 (f-vectors type B associahedra), A080721 (f-vectors type D associahedra), A126216 (mirror image).
Cf. A248727 for a relation to f-polynomials of simplices.
Cf. A111785 (contracted partition array, unsigned; see a comment above).
Antidiagonal sums give A005043. - Jordan Tirrell, Jun 01 2017

Programs

  • Magma
    [[Binomial(n-3, k)*Binomial(n+k-1, k)/(k+1): k in [0..(n-3)]]: n in [3..12]];  // G. C. Greubel, Nov 19 2018
    
  • Maple
    T:=(n,k)->binomial(n-3,k)*binomial(n+k-1,k)/(k+1): seq(seq(T(n,k),k=0..n-3),n=3..12); # Muniru A Asiru, Nov 24 2018
  • Mathematica
    t[n_, k_] = Binomial[n-3, k]*Binomial[n+k-1, k]/(k+1);
    Flatten[Table[t[n, k], {n, 3, 12}, {k, 0, n-3}]][[1 ;; 52]] (* Jean-François Alcover, Jun 16 2011 *)
  • PARI
    Q=(1+z-(1-(4*w+2+O(w^20))*z+z^2+O(z^20))^(1/2))/(2*(1+w)*z);for(n=3,12,for(m=1,n-2,print1(polcoef(polcoef(Q,n-2,z),m,w),", "))) \\ Hugo Pfoertner, Nov 19 2018
    
  • PARI
    for(n=3,12, for(k=0,n-3, print1(binomial(n-3,k)*binomial(n+k-1,k)/(k+1), ", "))) \\ G. C. Greubel, Nov 19 2018
    
  • Sage
    [[ binomial(n-3,k)*binomial(n+k-1,k)/(k+1) for k in (0..(n-3))] for n in (3..12)] # G. C. Greubel, Nov 19 2018

Formula

G.f. G = G(t, z) satisfies (1+t)*G^2 - z*(1-z-2*t*z)*G + t*z^4 = 0.
T(n, k) = binomial(n-3, k)*binomial(n+k-1, k)/(k+1) for n >= 3, 0 <= k <= n-3.
From Tom Copeland, Nov 03 2008: (Start)
Two g.f.s (f1 and f2) for A033282 and their inverses (x1 and x2) can be derived from the Drake and Barry references.
1. a: f1(x,t) = y = {1 - (2t+1) x - sqrt[1 - (2t+1) 2x + x^2]}/[2x (t+1)] = t x + (t + 2 t^2) x^2 + (t + 5 t^2 + 5 t^3) x^3 + ...
b: x1 = y/[t + (2t+1)y + (t+1)y^2] = y {1/[t/(t+1) + y] - 1/(1+y)} = (y/t) - (1+2t)(y/t)^2 + (1+ 3t + 3t^2)(y/t)^3 +...
2. a: f2(x,t) = y = {1 - x - sqrt[(1-x)^2 - 4xt]}/[2(t+1)] = (t/(t+1)) x + t x^2 + (t + 2 t^2) x^3 + (t + 5 t^2 + 5 t^3) x^4 + ...
b: x2 = y(t+1) [1- y(t+1)]/[t + y(t+1)] = (t+1) (y/t) - (t+1)^3 (y/t)^2 + (t+1)^4 (y/t)^3 + ...
c: y/x2(y,t) = [t/(t+1) + y] / [1- y(t+1)] = t/(t+1) + (1+t) y + (1+t)^2 y^2 + (1+t)^3 y^3 + ...
x2(y,t) can be used along with the Lagrange inversion for an o.g.f. (A133437) to generate A033282 and show that A133437 is a refinement of A033282, i.e., a refinement of the f-polynomials of the associahedra, the Stasheff polytopes.
y/x2(y,t) can be used along with the indirect Lagrange inversion (A134264) to generate A033282 and show that A134264 is a refinement of A001263, i.e., a refinement of the h-polynomials of the associahedra.
f1[x,t](t+1) gives a generator for A088617.
f1[xt,1/t](t+1) gives a generator for A060693, with inverse y/[1 + t + (2+t) y + y^2].
f1[x(t-1),1/(t-1)]t gives a generator for A001263, with inverse y/[t + (1+t) y + y^2].
The unsigned coefficients of x1(y t,t) are A074909, reverse rows of A135278. (End)
G.f.: 1/(1-x*y-(x+x*y)/(1-x*y/(1-(x+x*y)/(1-x*y/(1-(x+x*y)/(1-x*y/(1-.... (continued fraction). - Paul Barry, Feb 06 2009
Let h(t) = (1-t)^2/(1+(u-1)*(1-t)^2) = 1/(u + 2*t + 3*t^2 + 4*t^3 + ...), then a signed (n-1)-th row polynomial of A033282 is given by u^(2n-1)*(1/n!)*((h(t)*d/dt)^n) t, evaluated at t=0, with initial n=2. The power series expansion of h(t) is related to A181289 (cf. A086810). - Tom Copeland, Sep 06 2011
With a different offset, the row polynomials equal 1/(1 + x)*Integral_{0..x} R(n,t) dt, where R(n,t) = Sum_{k = 0..n} binomial(n,k)*binomial(n+k,k)*t^k are the row polynomials of A063007. - Peter Bala, Jun 23 2016
n-th row polynomial = ( LegendreP(n-1,2*x + 1) - LegendreP(n-3,2*x + 1) )/((4*n - 6)*x*(x + 1)), n >= 3. - Peter Bala, Feb 22 2017
n*T(n+1, k) = (4n-6)*T(n, k-1) + (2n-3)*T(n, k) - (n-3)*T(n-1, k) for n >= 4. - Fang Lixing, May 07 2019

Extensions

Missing factor of 2 for expansions of f1 and f2 added by Tom Copeland, Apr 12 2009

A086810 Triangle obtained by adding a leading diagonal 1,0,0,0,... to A033282.

Original entry on oeis.org

1, 0, 1, 0, 1, 2, 0, 1, 5, 5, 0, 1, 9, 21, 14, 0, 1, 14, 56, 84, 42, 0, 1, 20, 120, 300, 330, 132, 0, 1, 27, 225, 825, 1485, 1287, 429, 0, 1, 35, 385, 1925, 5005, 7007, 5005, 1430, 0, 1, 44, 616, 4004, 14014, 28028, 32032, 19448, 4862, 0, 1, 54, 936, 7644, 34398, 91728
Offset: 0

Views

Author

Philippe Deléham, Aug 05 2003

Keywords

Comments

Mirror image of triangle A133336. - Philippe Deléham, Dec 10 2008
From Tom Copeland, Oct 09 2011: (Start)
With polynomials
P(0,t) = 0
P(1,t) = 1
P(2,t) = t
P(3,t) = t + 2 t^2
P(4,t) = t + 5 t^2 + 5 t^3
P(5,t) = t + 9 t^2 + 21 t^3 + 14 t^4
The o.g.f. A(x,t) = {1+x-sqrt[(1-x)^2-4xt]}/[2(1+t)] (see Drake et al.).
B(x,t)= x-t x^2/(1-x)= x-t(x^2+x^3+x^4+...) is the comp. inverse in x.
Let h(x,t) = 1/(dB/dx) = (1-x)^2/(1+(1+t)*x*(x-2)) = 1/(1-t(2x+3x^2+4x^3+...)), an o.g.f. in x for row polynomials in t of A181289. Then P(n,t) is given by (1/n!)(h(x,t)*d/dx)^n x, evaluated at x=0, A = exp(x*h(y,t)*d/dy) y, eval. at y=0, and dA/dx = h(A(x,t),t). These results are a special case of A133437 with u(x,t) = B(x,t), i.e., u_1=1 and (u_n)=-t for n > 1. See A001003 for t = 1. (End)
Let U(x,t) = [A(x,t)-x]/t, then U(x,0) = -dB(x,t)/dt and U satisfies dU/dt = UdU/dx, the inviscid Burgers' equation (Wikipedia), also called the Hopf equation (see Buchstaber et al.). Also U(x,t) = U(A(x,t),0) = U(x+tU,0) since U(x,0) = [x-B(x,t)]/t. - Tom Copeland, Mar 12 2012
Diagonals of A132081 are essentially rows of this sequence. - Tom Copeland, May 08 2012
T(r, s) is the number of [0,r]-covering hierarchies with s segments (see Kreweras). - Michel Marcus, Nov 22 2014
From Yu Hin Au, Dec 07 2019: (Start)
T(n,k) is the number of small Schröder n-paths (lattice paths from (0,0) to (2n,0) using steps U=(1,1), F=(2,0), D=(1,-1) with no F step on the x-axis) that has exactly k U steps.
T(n,k) is the number of Schröder trees (plane rooted tree where each internal node has at least two children) with exactly n+1 leaves and k internal nodes. (End)

Examples

			Triangle starts:
  1;
  0,  1;
  0,  1,  2;
  0,  1,  5,  5;
  0,  1,  9, 21, 14;
  ...
		

Crossrefs

The diagonals (except for A000007) are also the diagonals of A033282.
Row sums: A001003 (Schroeder numbers).

Programs

  • Mathematica
    Table[Boole[n == 2] + If[# == -1, 0, Binomial[n - 3, #] Binomial[n + # - 1, #]/(# + 1)] &[k - 1], {n, 2, 12}, {k, 0, n - 2}] // Flatten (* after Jean-François Alcover at A033282, or *)
    Table[If[n == 0, 1, Binomial[n, k] Binomial[n + k, k - 1]/n], {n, 0, 10}, {k, 0, n}] // Flatten (* Michael De Vlieger, Aug 22 2016 *)
  • PARI
    t(n, k) = if (n==0, 1, binomial(n, k)*binomial(n+k, k-1)/n);
    tabl(nn) = {for (n=0, nn, for (k=0, n, print1(t(n,k), ", ");); print(););} \\ Michel Marcus, Nov 22 2014

Formula

Triangle T(n, k) read by rows; given by [0, 1, 0, 1, 0, 1, ...] DELTA [1, 1, 1, 1, 1, 1, 1, 1, 1, ...] where DELTA is Deléham's operator defined in A084938.
For k>0, T(n, k) = binomial(n-1, k-1)*binomial(n+k, k)/(n+1); T(0, 0) = 1 and T(n, 0) = 0 if n > 0. [corrected by Marko Riedel, May 04 2023]
Sum_{k>=0} T(n, k)*2^k = A107841(n). - Philippe Deléham, May 26 2005
Sum_{k>=0} T(n-k, k) = A005043(n). - Philippe Deléham, May 30 2005
T(n, k) = A108263(n+k, k). - Philippe Deléham, May 30 2005
Sum_{k=0..n} T(n,k)*x^k = A000007(n), A001003(n), A107841(n), A131763(n), A131765(n), A131846(n), A131926(n), A131869(n), A131927(n) for x = 0, 1, 2, 3, 4, 5, 6, 7, 8, respectively. - Philippe Deléham, Nov 05 2007
Sum_{k=0..n} T(n,k)*5^k*(-2)^(n-k) = A152601(n). - Philippe Deléham, Dec 10 2008
Sum_{k=0..n} T(n,k)*(-1)^k*3^(n-k) = A154825(n). - Philippe Deléham, Jan 17 2009
Umbrally, P(n,t) = Lah[n-1,-t*a.]/n! = (1/n)*Sum_{k=1..n-1} binomial(n-2,k-1)a_k t^k/k!, where (a.)^k = a_k = (n-1+k)!/(n-1)!, the rising factorial, and Lah(n,t) = n!*Laguerre(n,-1,t) are the Lah polynomials A008297 related to the Laguerre polynomials of order -1. - Tom Copeland, Oct 04 2014
T(n, k) = binomial(n, k)*binomial(n+k, k-1)/n, for k >= 0; T(0, 0) = 1 (see Kreweras, p. 21). - Michel Marcus, Nov 22 2014
P(n,t) = Lah[n-1,-:Dt:]/n! t^(n-1) with (:Dt:)^k = (d/dt)^k t^k = k! Laguerre(k,0,-:tD:) with (:tD:)^j = t^j D^j. The normalized Laguerre polynomials of 0 order are given in A021009. - Tom Copeland, Aug 22 2016

Extensions

Typo in a(60) corrected by Michael De Vlieger, Nov 21 2019

A133437 Irregular triangle of coefficients of a partition transform for direct Lagrange inversion of an o.g.f., complementary to A134685 for an e.g.f.; normalized by the factorials, these are signed, refined face polynomials of the associahedra.

Original entry on oeis.org

1, -2, 12, -6, -120, 120, -24, 1680, -2520, 360, 720, -120, -30240, 60480, -20160, -20160, 5040, 5040, -720, 665280, -1663200, 907200, 604800, -60480, -362880, -181440, 20160, 40320, 40320, -5040, -17297280, 51891840, -39916800, -19958400, 6652800, 19958400, 6652800, -1814400, -1814400, -3628800, -1814400, 362880, 362880, 362880, -40320
Offset: 1

Views

Author

Tom Copeland, Jan 27 2008

Keywords

Comments

Let f(t) = u(t) - u(0) = Sum_{n>=1} u_n * t^n.
If u_1 is not equal to 0, then the compositional inverse for f(t) is given by g(t) = Sum_{j>=1} P(n,t) where, with u_n denoted by (n'),
P(1,t) = (1')^(-1) * [ 1 ] * t
P(2,t) = (1')^(-3) * [ -2 (2') ] * t^2 / 2!
P(3,t) = (1')^(-5) * [ 12 (2')^2 - 6 (1')(3') ] * t^3 / 3!
P(4,t) = (1')^(-7) * [ -120 (2')^3 + 120 (1')(2')(3') - 24 (1')^2 (4') ] * t^4 / 4!
P(5,t) = (1')^(-9) * [ 1680 (2')^4 - 2520 (1') (2')^2 (3') + 360 (1')^2 (3')^2 + 720 (1')^2 (2') (4') - 120 (1')^3 (5') ] * t^5 / 5!
P(6,t) = (1')^(-11) * [ -30240 (2')^5 + 60480 (1') (2')^3 (3') - 20160 (1')^2 (2') (3')^2 - 20160 (1')^2 (2')^2 (4') + 5040 (1')^3 (3')(4') + 5040 (1')^3 (2')(5') - 720 (1')^4 (6') ] * t^6 / 6!
P(7,t) = (1')^(-13) * [ 665280 (2')^6 - 1663200 (1')(2')^4(3') + (1')^2 [907200 (2')^2(3')^2 + 604800 (2')^3(4')] - (1')^3 [60480 (3')^3 + 362880 (2')(3')(4') + 181440 (2')^2(5')] + (1')^4 [20160 (4')^2 + 40320 (3')(5') + 40320 (2')(6')] - 5040 (1')^5(7')] * t^7 / 7!
P(8,t) = (1')^(-15) * [ -17297280 (2')^7 + 51891840 (1')(2')^5(3') - (1')^2 [39916800 (2')^3(3')^2 + 19958400 (2')^4(4')] + (1')^3 [6652800 (2')(3')^3 + 19958400 (2')^2(3')(4') + 6652800 (2')^3(5')] - (1')^4 [1814400 (3')^2(4') + 1814400 (2')(4')^2 + 3628800 (2')(3')(5') + 1814400 (2')^2(6')] + (1')^5 [362880 (4')(5') + 362880 (3')(6') + 362880 (2')(7')] - 40320 (1')^6(8')] * t^8 / 8!
...
See A134685 for more information.
A111785 is obtained from A133437 by dividing through the bracketed terms of the P(n,t) by n! and unsigned A111785 is a refinement of A033282 and A126216. - Tom Copeland, Sep 28 2008
For relation to the geometry of associahedra or Stasheff polytopes (and other combinatorial objects) see the Loday and McCammond links. E.g., P(5,t) = (1')^(-9) * [ 14 (2')^4 - 21 (1') (2')^2 (3') + 6 (1')^2 (2') (4')+ 3 (1')^2 (3')^2 - 1 (1')^3 (5') ] * t^5 is related to the 3-D associahedron with 14 vertices (0-D faces), 21 edges (1-D faces), 6 pentagons (2-D faces), 3 rectangles (2-D faces), 1 3-D polytope (3-D faces). Summing over faces of the same dimension gives A033282 or A126216. - Tom Copeland, Sep 29 2008
A relation between this Lagrange inversion for an o.g.f. and partition polynomials formed from the "refined Lah numbers" A130561 is presented in the link "Lagrange a la Lah" along with umbral binary tree representations.
With f(x,t) = x + t*Sum_{n>=2} u_n*x^n, the compositional inverse in x is related to the velocity profile of particles governed by the inviscid Burgers's, or Hopf, eqn. See A001764 and A086810. - Tom Copeland, Feb 15 2014
Newton was aware of this power series expansion for series reversion. See the Ferraro reference p. 75 eqn. 52. - Tom Copeland, Jan 22 2017
The coefficients of the partition polynomials divided by the associated factorial enumerate the faces of the convex, bounded polytopes called the associahedra, and the absolute value of the sum of the renormalized coefficients gives the Euler characteristic of unity for each polytope; i.e., the absolute value of the sum of each row of the array is either n! (unnormalized) or unity (normalized). In addition, the signs of the faces alternate with dimension, and the coefficients of faces with the same dimension for each polytope have the same sign. - Tom Copeland, Nov 13 2019
With u_1 = 1 and the other u_n replaced by suitably signed partition polynomials of A263633, the partition polynomials enumerating the refined noncrossing partitions of A134264 with a shift in indices are obtained (cf. In the Realm of Shadows). - Tom Copeland, Nov 16 2019
Relations between associahedra and oriented n-simplices are presented by Halvorson and by Street. - Tom Copeland, Dec 08 2019
Let f(x;t,n) = x - t * x^(n+1), giving u_1 = (1') = 1 and u_(n+1) = (n+1) = -t. Then inverting in x with t a parameter gives finv(x;t,n) = Sum_{j>=0} {binomial((n+1)*j,j) / (n*j + 1)} * t^j * x^(n*j + 1), which gives the Catalan numbers for n=1, and the Fuss-Catalan sequences for n>1 (see A001764, n=2). Comparing this with the same result in A134264 gives relations between the faces of associahedra and noncrossing partitions (and other combinatorial constructs related to these inversion formulas and those listed in A145271). - Tom Copeland, Jan 27 2020
From Tom Copeland, Mar 24 2020: (Start)
There is a mapping between the faces of K_n, the associahedron of dimension (n-1), and polygon dissections. The dissecting noncrossing diagonals (i.e., nonintersecting in the interior) form subpolygons. Assign the indeterminate x_k to a subpolygon where k = (number of vertices of the subpolygon) - 1. Multiply the x_k together to form the monomials for the inversion formula.
For the 3-dimensional associahedron K_4, the fundamental polygon is the hexagon, which can be dissected into pentagons, associated to x_4; tetragons, to x_3; and triangles, to x_2; for example, there are six distinguished partitions of the hexagon into one triangle and one pentagon, sharing two vertices, associated to the monomial 6 * x_2 * x_4 since the unshared vertex of the triangle can be moved consecutively from one vertex of the hexagon to the next. This term corresponds to 720 (1')^2 (2') (4') / 5! in P(5,t) above, denumerating the six pentagonal facets of K_4. (End)

References

  • G. Ferraro, The Rise and Development of the Theory of Series up to the Early 1820s, Springer Science and Business Media, 2007.
  • H. Halvorson (editor), Deep Beauty: Understanding the Quantum World Through Innovation, Cambridge Univ. Press, 2011.
  • H. Turnbull (editor), The Correspondence of Isaac Newton Vol. II 1676-1687, Cambridge Univ. Press, 1960, p. 147.

Crossrefs

Cf. A145271, (A086810, A181289) = (reduced array, associated g(x)).

Programs

  • Mathematica
    rows[nn_] := {{1}}~Join~With[{s = InverseSeries[t (1 + Sum[u[k] t^k, {k, nn}] + O[t]^(nn+1))]}, Table[(n+1)! Coefficient[s, t^(n+1) Product[u[w], {w, p}]], {n, nn}, {p, Reverse[Sort[Sort /@ IntegerPartitions[n]]]}]];
    rows[7] // Flatten (* Andrey Zabolotskiy, Mar 07 2024 *)

Formula

The bracketed partitions of P(n,t) are of the form (u_1)^e(1) (u_2)^e(2) ... (u_n)^e(n) with coefficients given by (-1)^(n-1+e(1)) * [2*(n-1)-e(1)]! / [ (e(2))! * (e(3))! * ... * (e(n))! ].
From Tom Copeland, Sep 06 2011: (Start)
Let h(t) = 1/(df(t)/dt)
= 1/Ev[u./(1-u.t)^2]
= 1/((u_1) + 2*(u_2)*t + 3*(u_3)*t^2 + 4*(u_4)*t^3 + ...),
where Ev denotes umbral evaluation.
Then for the partition polynomials of A133437,
n!*P(n,t) = ((t*h(y)*d/dy)^n) y evaluated at y=0,
and the compositional inverse of f(t) is
g(t) = exp(t*h(y)*d/dy) y evaluated at y=0.
Also, dg(t)/dt = h(g(t)). (End)
From Tom Copeland, Oct 20 2011: (Start)
With exp[x* PS(.,t)] = exp[t*g(x)] = exp[x*h(y)d/dy] exp(t*y) eval. at y=0, the raising/creation and lowering/annihilation operators defined by R PS(n,t)=PS(n+1,t) and L PS(n,t) = n*PS(n-1,t) are
R = t*h(d/dt) = t* 1/[(u_1) + 2*(u_2)*d/dt + 3*(u_3)*(d/dt)^2 + ...] and
L = f(d/dt) = (u_1)*d/dt + (u_2)*(d/dt)^2 + (u_3)*(d/dt)^3 + ....
Then P(n,t) = (t^n/n!) dPS(n,z)/dz eval. at z=0. (Cf. A139605, A145271, and link therein to Mathemagical Forests for relation to planted trees on p. 13.) (End)
The bracketed partition polynomials of P(n,t) are also given by (d/dx)^(n-1) 1/[u_1 + u_2 * x + u_3 * x^2 + ... + u_n * x^(n-1)]^n evaluated at x=0. - Tom Copeland, Jul 07 2015
From Tom Copeland, Sep 20 2016: (Start)
Let PS(n,u1,u2,...,un) = P(n,t) / t^n, i.e., the square-bracketed part of the partition polynomials in the expansion for the inverse in the comment section, with u_k = uk.
Also let PS(n,u1=1,u2,...,un) = PB(n,b1,b2,...,bK,...) where each bK represents the partitions of PS, with u1 = 1, that have K components or blocks, e.g., PS(5,1,u2,...,u5) = PB(5,b1,b2,b3,b4) = b1 + b2 + b3 + b4 with b1 = -u5, b2 = 6 u2 u4 + 3 u3^2, b3 = -21 u2^2 u3, and b4 = 14 u2^4.
The relation between solutions of the inviscid Burgers' equation and compositional inverse pairs (cf. A086810) implies that, for n > 2, PB(n, 0 * b1, 1 * b2, ..., (K-1) * bK, ...) = [(n+1)/2] * Sum_{k = 2..n-1} PS(n-k+1,u_1=1,u_2,...,u_(n-k+1)) * PS(k,u_1=1,u_2,...,u_k).
For example, PB(5,0 * b1, 1 * b2, 2 * b3, 3 * b4) = 3 * 14 u2^4 - 2 * 21 u2^2 u3 + 1 * 6 u2 u4 + 1 * 3 u3^2 - 0 * u5 = 42 u2^4 - 42 u2^2 u3 + 6 u2 u4 + 3 u3^2 = 3 * [2 * PS(2,1,u2) * PS(4,1,u2,...,u4) + PS(3,1,u2,u3)^2] = 3 * [ 2 * (-u2) (-5 u2^3 + 5 u2 u3 - u4) + (2 u2^2 - u3)^2].
Also, PB(n,0*b1,1*b2,...,(K-1)*bK,...) = d/dt t^(n-2)*PS(n,u1=1/t,u2,...,un)|{t=1} = d/dt (1/t)*PS(n,u1=1,t*u2,...,t*un)|{t=1}.
(End)
From Tom Copeland, Sep 22 2016: (Start)
Equivalent matrix computation: Multiply the m-th diagonal (with m=1 the index of the main diagonal) of the lower triangular Pascal matrix A007318 by f_m = m!*u_m = (d/dx)^m f(x) evaluated at x=0 to obtain the matrix UP with UP(n,k) = binomial(n,k) f_{n+1-k}, or equivalently multiply the diagonals of A132159 by u_m. Then P(n,t) = (1, 0, 0, 0, ...) [UP^(-1) * S]^(n-1) FC * t^n/n!, where S is the shift matrix A129185, representing differentiation in the basis x^n//n!, and FC is the first column of UP^(-1), the inverse matrix of UP. These results follow from A145271 and A133314.
Also, P(n,t) = (1, 0, 0, 0, ...) [UP^(-1) * S]^n (0, 1, 0, ...)^T * t^n/n! in agreement with A139605. (End)
A recursion relation for computing each partition polynomial of this entry from the lower order polynomials and the coefficients of the refined Lah polynomials of A130561 is presented in the blog entry "Formal group laws and binomial Sheffer sequences." - Tom Copeland, Feb 06 2018
The derivative of the partition polynomials of A350499 with respect to a distinguished indeterminate give polynomials proportional to those of this entry. The connection of this derivative relation to the inviscid Burgers-Hopf evolution equation is given in a reference for that entry. - Tom Copeland, Feb 19 2022

Extensions

Missing coefficient in P(6,t) replaced by Tom Copeland, Nov 06 2008
P(7,t) and P(8,t) data added by Tom Copeland, Jan 14 2016
Title modified by Tom Copeland, Jan 13 2020
Terms ordered according to the reversed Abramowitz-Stegun ordering of partitions (with every k' replaced by (k-1)') by Andrey Zabolotskiy, Mar 07 2024

A126216 Triangle read by rows: T(n,k) is the number of Schroeder paths of semilength n containing exactly k peaks but no peaks at level one (n >= 1; 0 <= k <= n-1).

Original entry on oeis.org

1, 2, 1, 5, 5, 1, 14, 21, 9, 1, 42, 84, 56, 14, 1, 132, 330, 300, 120, 20, 1, 429, 1287, 1485, 825, 225, 27, 1, 1430, 5005, 7007, 5005, 1925, 385, 35, 1, 4862, 19448, 32032, 28028, 14014, 4004, 616, 44, 1, 16796, 75582, 143208, 148512, 91728, 34398, 7644, 936, 54, 1
Offset: 1

Views

Author

Emeric Deutsch, Dec 20 2006

Keywords

Comments

A Schroeder path of semilength n is a lattice path in the first quadrant, from the origin to the point (2n,0) and consisting of steps U=(1,1), D=(1,-1) and H=(2,0).
Also number of Schroeder paths of semilength n containing exactly k doublerises but no (2,0) steps at level 0 (n >= 1; 0 <= k <= n-1). Also number of doublerise-bicolored Dyck paths (doublerises come in two colors; also called marked Dyck paths) of semilength n and having k doublerises of a given color (n >= 1; 0 <= k <= n-1). Also number of 12312- and 121323-avoiding matchings on [2n] with exactly k crossings.
Essentially the triangle given by [1,1,1,1,1,1,1,1,...] DELTA [0,1,0,1,0,1,0,1,0,1,...] where DELTA is the operator defined in A084938. - Philippe Deléham, Oct 20 2007
Mirror image of triangle A033282. - Philippe Deléham, Oct 20 2007
For relation to Lagrange inversion, or series reversion and the geometry of associahedra, or Stasheff polytopes (and other combinatorial objects), see A133437. - Tom Copeland, Sep 29 2008
First column (k=0) gives the Catalan numbers (A000108). - Alexander Karpov, Jun 10 2018
T(n,k) is the multiplicity of the k-hook representation of the symmetric group in the (n-1)st parking space representation (see Pak and Postnikov, 1995). - Joshua Mundinger, Jul 18 2025

Examples

			T(3,1)=5 because we have HUUDD, UUDDH, UUUDDD, UHUDD and UUDHD.
Triangle starts:
   n\k  0      1      2      3      4     5    6   7  8
   1    1;
   2    2,     1;
   3    5,     5;     1;
   4   14,    21,     9,     1;
   5   42,    84,    56,    14,     1;
   6  132,   330,   300,   120,    20,    1;
   7  429,  1287,  1485,   825,   225,   27,   1;
   8 1430,  5005,  7007,  5005,  1925,  385,  35,  1;
   9 4862, 19448, 32032, 28028, 14014, 4004, 616, 44, 1;
  10 ...
Triangle [1,1,1,1,1,1,1,...] DELTA [0,1,0,1,0,1,0,1,...] begins:
   1;
   1,  0;
   2,  1,  0;
   5,  5,  1,  0;
  14, 21,  9,  1,  0;
  42, 84, 56, 14,  1,  0;
  ...
		

Crossrefs

Programs

  • Maple
    T:=(n,k)->binomial(n,k)*binomial(2*n-k,n+1)/n: for n from 1 to 11 do seq(T(n,k),k=0..n-1) od; # yields sequence in triangular form
  • Mathematica
    Table[Binomial[n, k] Binomial[2 n - k, n + 1]/n, {n, 10}, {k, 0, n - 1}] // Flatten (* Michael De Vlieger, Jan 09 2016 *)
  • PARI
    tabl(nn) = {mP = matrix(nn, nn, n, k, binomial(n-1, k-1)); mN = matrix(nn, nn, n, k, binomial(n-1, k-1) * binomial(n, k-1) / k); mprod = mN*mP; for (n=1, nn, for (k=1, n, print1(mprod[n, k], ", ");); print(););} \\ Michel Marcus, Apr 16 2015
    
  • PARI
    t(n,k) = binomial(n,k)*binomial(2*n-k,n+1)/n;
    concat(vector(10, n, vector(n, k, t(n,k-1))))  \\ Gheorghe Coserea, Apr 24 2016

Formula

T(n,k) = C(n,k)*C(2*n-k,n+1)/n (0 <= k <= n-1).
G.f.: G(t,z) = (1-2*z-t*z-sqrt(1-4*z-2*t*z+t^2*z^2))/(2*(1+t)*z).
Equals N * P, where N = the Narayana triangle (A001263) and P = Pascal's triangle, as infinite lower triangular matrices. A126182 = P * N. - Gary W. Adamson, Nov 30 2007
G.f.: 1/(1-x-(x+xy)/(1-xy/(1-(x+xy)/(1-xy/(1-(x+xy)/(1-xy/(1-.... (continued fraction). - Paul Barry, Feb 06 2009
Let h(t) = (1-t)^2/(1+(u-1)*(1-t)^2) = 1/(u + 2*t + 3*t^2 + 4*t^3 + ...), then a signed (n-1)-th row polynomial of A126216 is given by u^(2n-1)*(1/n!)*((h(t)*d/dt)^n) t, evaluated at t=0, with initial n=2. The power series expansion of h(t) is related to A181289 (cf. A086810). - Tom Copeland, Oct 09 2011
From Tom Copeland, Oct 10 2011: (Start)
With polynomials
P(0,t) = 0
P(1,t) = 1
P(2,t) = 1
P(3,t) = 2 + t
P(4,t) = 5 + 5 t + t^2
P(5,t) = 14 + 21 t + 9 t^2 + t^3
The o.g.f. A(x,t) = (1+x*t-sqrt((1-x*t)^2-4x))/(2(1+t)), and
B(x,t) = x - x^2/(1-t*x) = x - x^2 - ((t*x)^3 + (t*x)^4 + ...)/t^2 is the compositional inverse in x. [series corrected by Tom Copeland, Dec 10 2019]
Let h(x,t) = 1/(dB/dx) = (1-tx)^2/(1-(t+1)(2x-tx^2)) = 1/(1 - 2x - 3tx^2 + 4t^2x^3 + ...). Then P(n,t) = (1/n!)(h(x,t)*d/dx)^n x, evaluated at x=0, A = exp(x*h(u,t)*d/du) u, evaluated at u=0, and dA/dx = h(A(x,t),t). (End)
From Tom Copeland, Dec 09 2019: (Start)
The polynomials in my 2011 formula entry above evaluate to an aerated, alternating sign sequence of the Catalan numbers A000108 with t = -2. The first few are P(2,-2) = 1, P(3,-2) = 0, P(4,t) = -1, P(5,-2) = 0, P(6,-2) = 2, P(7,-2) = 0, P(8,-2) = -5, P(9,-2) = 0, P(10,-2) = 14.
Generalizing the relations between w = theta and u = phi in Mizera on pp. 32-34, modify the inverse pair above to w = i * B(-i*u,t) = u + i * u^2/(1+i*t*u), where i is the imaginary number, and u = i*A(-i*w,t) = i*(1 - i*w*t - sqrt((1 + i*w*t)^2 + i*4*w))/(2(1+t)). Then the expression for V'(w) in Mizera generalizes to V'(w) = -i*(w - u) and reduces to V'(w) = (1 - sqrt(1-4 w^2))/2 when evaluated at t = -2, which is an o.g.f. for A126120. Cf. also A086810. (End)
Sum_{k = 0..n-1} (-1)^k*T(n,k)*binomial(x + 2*n - k, 2*n - k) = ( (x + 1) * ( Product_{k = 2..n} (x + k)^2 ) * (x + n + 1) )/(n!*(n + 1)!) for n >= 1. Cf. A243660 and A243661. - Peter Bala, Oct 08 2022

A181290 The sum of the lengths of the 2-compositions of n. A 2-composition of n is a nonnegative matrix with two rows, such that each column has at least one nonzero entry and whose entries sum up to n. The length of the 2-composition is the number of columns.

Original entry on oeis.org

0, 2, 11, 52, 227, 944, 3800, 14944, 57748, 220128, 829968, 3101376, 11502704, 42393088, 155392768, 566918144, 2059768384, 7456496128, 26905720576, 96804463616, 347386161920, 1243665567744, 4442849839104, 15840448094208, 56375692407808, 200307512532992, 710622022258688, 2517475213557760
Offset: 0

Views

Author

Emeric Deutsch, Oct 12 2010

Keywords

Crossrefs

Cf. A181289.

Programs

  • Maple
    g := z*(1-z)^2*(2-z)/(1-4*z+2*z^2)^2: gser := series(g, z = 0, 28): seq(coeff(gser, z, n), n = 0 .. 25);

Formula

a(n) = Sum_{k=0..n} k * A181289(n,k).
G.f.: z*(2-z)*(1-z)^2/(1-4*z+2*z^2)^2.

A382820 Number of columns in all n-compositions of n.

Original entry on oeis.org

1, 11, 163, 3019, 67251, 1753877, 52468711, 1772042699, 66708748963, 2770212058261, 125812351808551, 6203908746628501, 330108021642012407, 18853083403505443593, 1150352428059538611663, 74685045367715777653195, 5140745255774277374241411, 373950591013899715795929605
Offset: 1

Views

Author

John Tyler Rascoe, Apr 05 2025

Keywords

Comments

A k-composition of n is a rectangular array of nonnegative integers with k rows, at least one nonzero entry in each column, and having the sum of all entries equal to n.

Crossrefs

C.f. A001792, A145839, A181289, A181290, (main diagonal of A382818).

Programs

  • PARI
    A382818_Column(k,N) = {my(x='x+O('x^N)); Vec(-(((1 - x)^k - 1)*(1 - x)^k)/( ((1 - x)^k - 1) + (1 - x)^k)^2)}
    A382820(n)= {A382818_Column(n,n+1)[n]}

Formula

a(n) = [x^n] -((1 - x)^n - 1)*(1 - x)^n/(((1 - x)^n - 1) + (1 - x)^n)^2

A382818 Square array A(n,k), n > 0, k > 0, read by downward antidiagonals: A(n,k) is the number of columns in all k-compositions of n.

Original entry on oeis.org

1, 2, 3, 3, 11, 8, 4, 24, 52, 20, 5, 42, 163, 227, 48, 6, 65, 372, 1017, 944, 112, 7, 93, 710, 3019, 6030, 3800, 256, 8, 126, 1208, 7095, 23256, 34563, 14944, 576, 9, 164, 1897, 14340, 67251, 173076, 193392, 57748, 1280, 10, 207, 2808, 26082, 161394, 615630, 1256936, 1062756, 220128, 2816
Offset: 1

Views

Author

John Tyler Rascoe, Apr 05 2025

Keywords

Comments

A k-composition of n is a rectangular array of nonnegative integers with k rows, at least one nonzero entry in each column, and having the sum of all entries equal to n.

Examples

			Square array begins:
   1,   2,    3,     4,     5,      6, ...
   3,  11,   24,    42,    65,     93, ...
   8,  52,  163,   372,   710,   1208, ...
  20, 227, 1017,  3019,  7095,  14340, ...
  48, 944, 6030, 23256, 67251, 161394, ...
  ...
A(2,2) = 11 counts the columns in the 2-compositions of 2:
 [2]   [0]   [1]   [1,0]   [0,1]   [0,0]   [1,1]
 [0],  [2],  [1],  [0,1],  [1,0],  [1,1],  [0,0].
		

Crossrefs

C.f. A001792 (column k=1), A005475 (row n=2), A145839, A181289, A181290 (column k=2), A382820 (main diagonal).

Programs

  • PARI
    A382818_Column(k,N) = {my(x='x+O('x^N)); Vec(-(((1 - x)^k - 1)*(1 - x)^k)/( ((1 - x)^k - 1) + (1 - x)^k)^2)}
    A382818_array(max_row) = {my(m=matrix(0)); for(n=1,max_row, m=matconcat([m,A382818_Column(n,max_row)~])); m}
    A382818_array(10)

Formula

Column k has g.f.: -((1 - x)^k - 1)*(1 - x)^k/(((1 - x)^k - 1) + (1 - x)^k)^2.

A006636 a(n) = (n + 1)*(n + 2)*(n + 4)*(n + 8)*(n + 15)/120.

Original entry on oeis.org

8, 36, 102, 231, 456, 819, 1372, 2178, 3312, 4862, 6930, 9633, 13104, 17493, 22968, 29716, 37944, 47880, 59774, 73899, 90552, 110055, 132756, 159030, 189280, 223938, 263466, 308357, 359136, 416361, 480624, 552552, 632808, 722092, 821142, 930735
Offset: 0

Views

Author

Keywords

Comments

Former name: From generalized Catalan numbers. - G. C. Greubel, Sep 03 2025

References

  • H. M. Finucan, Some decompositions of generalized Catalan numbers, pp. 275-293 of Combinatorial Mathematics IX. Proc. Ninth Australian Conference (Brisbane, August 1981). Ed. E. J. Billington, S. Oates-Williams and A. P. Street. Lecture Notes Math., 952. Springer-Verlag, 1982.
  • N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).

Crossrefs

Cf. A181289.

Programs

  • Magma
    A006636:= func< n | (n+1)*(n+2)*(n+4)*(n+8)*(n+15)/120 >;
    [A006636(n): n in [0..40]]; // G. C. Greubel, Sep 03 2025
    
  • Mathematica
    Table[(n+1)*(n+2)*(n+4)*(n+8)*(n+15)/120, {n,0,40}] (* G. C. Greubel, Sep 03 2025 *)
  • SageMath
    def A006636(n): return (n+1)*(n+2)*(n+4)*(n+8)*(n+15)//120
    print([A006636(n) for n in range(41)]) # G. C. Greubel, Sep 03 2025

Formula

From Sean A. Irvine, May 30 2017: (Start)
a(n) = (n + 1)*(n + 2)*(n + 4)*(n + 8)*(n + 15)/120.
G.f.: (2 - x)^3/(1 - x)^6. (End)
E.g.f.: exp(x)*(960 + 3360*x + 2280*x^2 + 500*x^3 + 40*x^4 + x^5)/120. - Stefano Spezia, Oct 15 2022

Extensions

a(6) and a(8) corrected and more terms from Sean A. Irvine, May 30 2017
New name by G. C. Greubel, Sep 03 2025
Showing 1-10 of 12 results. Next