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

A130561 Numbers associated to partitions, used for combinatoric interpretation of Lah triangle numbers A105278; elementary Schur polynomials / functions.

Original entry on oeis.org

1, 2, 1, 6, 6, 1, 24, 24, 12, 12, 1, 120, 120, 120, 60, 60, 20, 1, 720, 720, 720, 360, 360, 720, 120, 120, 180, 30, 1, 5040, 5040, 5040, 5040, 2520, 5040, 2520, 2520, 840, 2520, 840, 210, 420, 42, 1, 40320, 40320, 40320, 40320, 20160, 20160, 40320, 40320, 20160
Offset: 1

Views

Author

Wolfdieter Lang, Jul 13 2007

Keywords

Comments

The order of this array is according to the Abramowitz-Stegun (A-St) ordering of partitions (see A036036).
The row lengths sequence is A000041 (partition numbers) [1, 2, 3, 5, 7, 11, 15, 22, 30, 42, ...].
These numbers are similar to M_0, M_1, M_2, M_3, M_4 given in A111786, A036038, A036039, A036040, A117506, respectively.
Combinatorial interpretation: a(n,k) counts the sets of lists (ordered subsets) obtained from partitioning the set {1..n}, with the lengths of the lists given by the k-th partition of n in A-St order. E.g., a(5,5) is computed from the number of sets of lists of lengths [1^1,2^2] (5th partition of 5 in A-St order). Hence a(5,5) = binomial(5,2)*binomial(3,2) = 5!/(1!*2!) = 60 from partitioning the numbers 1,2,...,5 into sets of lists of the type {[.],[..],[..]}.
This array, called M_3(2), is the k=2 member of a family of partition arrays generalizing A036040 which appears as M_3 = M_3(k=1). S2(2) = A105278 (unsigned Lah number triangle) is related to M_3(2) in the same way as S2(1), the Stirling2 number triangle, is related to M_3(1). - Wolfdieter Lang, Oct 19 2007
Another combinatorial interpretation: a(n,k) enumerates unordered forests of increasing binary trees which are described by the k-th partition of n in the Abramowitz-Stegun order. - Wolfdieter Lang, Oct 19 2007
A relation between partition polynomials formed from these "refined Lah numbers" and Lagrange inversion for an o.g.f. is presented in the link "Lagrange a la Lah" along with an e.g.f. and an umbral binary operator tree representation. - Tom Copeland, Apr 12 2011
With the indeterminates (x_1,x_2,x_3,...) = (t,-c_2*t,-c_3*t,...) with c_n >0, umbrally P(n,a.) = P(n,t)|{t^n = a_n} = 0 and P(j,a.)P(k,a.) = P(j,t)P(k,t)|{t^n =a_n} = d_{j,k} >= 0 is the coefficient of x^j/j!*y^k/k! in the Taylor series expansion of the formal group law FGL(x,y) = f[f^{-1}(x)+f^{-1}(y)], where a_n are the inversion partition polynomials for calculating f(x) from the coefficients of the series expansion of f^{-1}(x) given in A133437. - Tom Copeland, Feb 09 2018
Divided by n!, the row partition polynomials are the elementary homogeneous Schur polynomials presented on p. 44 of the Bracci et al. paper. - Tom Copeland, Jun 04 2018
Also presented (renormalized) as the Schur polynomials on p. 19 of the Konopelchenko and Schief paper with associations to differential operators related to the KP hierarchy. - Tom Copeland, Nov 19 2018
Through equation 4.8 on p. 26 of the Arbarello reference, these polynomials appear in the Hirota bilinear equations 4.7 related to tau-function solutions of the KP hierarchy. - Tom Copeland, Jan 21 2019
These partition polynomials appear as Feynman amplitudes in their Bell polynomial guise (put x_n = n!c_n in A036040 for the indeterminates of the Bell polynomials) in Kreimer and Yeats and Balduf (e.g., p. 27). - Tom Copeland, Dec 17 2019
From Tom Copeland, Oct 15 2020: (Start)
With a_n = n! * b_n = (n-1)! * c_n for n > 0, represent a function with f(0) = a_0 = b_0 = 1 as an
A) exponential generating function (e.g.f), or formal Taylor series: f(x) = e^{a.x} = 1 + Sum_{n > 0} a_n * x^n/n!
B) ordinary generating function (o.g.f.), or formal power series: f(x) = 1/(1-b.x) = 1 + Sum_{n > 0} b_n * x^n
C) logarithmic generating function (l.g.f): f(x) = 1 - log(1 - c.x) = 1 + Sum_{n > 0} c_n * x^n /n.
Expansions of log(f(x)) are given in
I) A127671 and A263634 for the e.g.f: log[ e^{a.*x} ] = e^{L.(a_1,a_2,...)x} = Sum_{n > 0} L_n(a_1,...,a_n) * x^n/n!, the logarithmic polynomials, cumulant expansion polynomials
II) A263916 for the o.g.f.: log[ 1/(1-b.x) ] = log[ 1 - F.(b_1,b_2,...)x ] = -Sum_{n > 0} F_n(b_1,...,b_n) * x^n/n, the Faber polynomials.
Expansions of exp(f(x)-1) are given in
III) A036040 for an e.g.f: exp[ e^{a.x} - 1 ] = e^{BELL.(a_1,...)x}, the Bell/Touchard/exponential partition polynomials, a.k.a. the Stirling partition polynomials of the second kind
IV) A130561 for an o.g.f.: exp[ b.x/(1-b.x) ] = e^{LAH.(b.,...)x}, the Lah partition polynomials
V) A036039 for an l.g.f.: exp[ -log(1-c.x) ] = e^{CIP.(c_1,...)x}, the cycle index polynomials of the symmetric groups S_n, a.k.a. the Stirling partition polynomials of the first kind.
Since exp and log are a compositional inverse pair, one can extract the indeterminates of the log set of partition polynomials from the exp set and vice versa. For a discussion of the relations among these polynomials and the combinatorics of connected and disconnected graphs/maps, see Novak and LaCroix on classical moments and cumulants and the two books on statistical mechanics referenced in A036040. (End)
These partition polynomials are referred to as Schur functions by Segal and Wilson, who present associations with Plucker coordinates, Grassmannians, and the tau functions of the KdV hierarchy. See pages 51 and 61. - Tom Copeland, Jan 08 2022

Examples

			Triangle starts:
  [  1];
  [  2,   1];
  [  6,   6,   1];
  [ 24,  24,  12, 12,  1];
  [120, 120, 120, 60, 60, 20, 1];
  ...
a(5,6) = 20 = 5!/(3!*1!) because the 6th partition of 5 in A-St order is [1^3,2^1].
a(5,5) = 60 enumerates the unordered [1^1,2^2]-forest with 5 vertices (including the three roots) composed of three such increasing binary trees: 5*((binomial(4,2)*2)*(1*2))/2! = 5*12 = 60.
		

References

  • E. Arbarello, "Sketches of KdV", Contemp. Math. 312 (2002), p. 9-69.

Crossrefs

Cf. A105278 (unsigned Lah triangle |L(n, m)|) obtained by summing the numbers for given part number m.
Cf. A000262 (row sums), identical with row sums of unsigned Lah triangle A105278.
A134133(n, k) = A130561(n, k)/A036040(n, k) (division by the M_3 numbers). - Wolfdieter Lang, Oct 12 2007
Cf. A096162.
Cf. A133437.
Cf. A127671.

Formula

a(n,k) = n!/(Product_{j=1..n} e(n,k,j)!) with the exponent e(n,k,j) of j in the k-th partition of n in the A-St ordering of the partitions of n. Exponents 0 can be omitted due to 0!=1.
From Tom Copeland, Sep 18 2011: (Start)
Raising and lowering operators are given for the partition polynomials formed from A130561 in the Copeland link in "Lagrange a la Lah Part I" on pp. 22-23.
An e.g.f. for the partition polynomials is on page 3:
exp[t*:c.*x/(1-c.*x):] = exp[t*(c_1*x + c_2*x^2 + c_3*x^3 + ...)] where :(...): denotes umbral evaluation of the enclosed expression and c. is an umbral coefficient. (End)
From Tom Copeland, Sep 07 2016: (Start)
The row partition polynomials of this array P(n,x_1,x_2,...,x_n), given in the Lang link, are n! * S(n,x_1,x_2,...,x_n), where S(n,x_1,...,x_n) are the elementary Schur polynomials, for which d/d(x_m) S(n,x_1,...,x_n) = S(n-m,x_1,...,x_(n-m)) with S(k,...) = 0 for k < 0, so d/d(x_m) P(n,x_1,...,x_n) = (n!/(n-m)!) P(n-m,x_1,...,x_(n-m)), confirming that the row polynomials form an Appell sequence in the indeterminate x_1 with P(0,...) = 1. See p. 127 of the Ernst paper for more on these Schur polynomials.
With the e.g.f. exp[t * P(.,x_1,x_2,..)] = exp(t*x_1) * exp(x_2 t^2 + x_3 t^3 + ...), the e.g.f. for the partition polynomials that form the umbral compositional inverse sequence U(n,x_1,...,x_n) in the indeterminate x_1 is exp[t * U(.,x_1,x_2,...)] = exp(t*x_1) exp[-(x_2 t^2 + x_3 t^3 + ...)]; therefore, U(n,x_1,x_2,...,x_n) = P(n,x_1,-x_2,.,-x_n), so umbrally P[n,P(.,x_1,-x_2,-x_3,...),x_2,x_3,...,x_n] = (x_1)^n = P[n,P(.,x_1,x_2,...),-x_2,-x_3,...,-x_n]. For example, P(1,x_1) = x_1, P2(x_1,x_2) = 2 x_2 + x_1^2, and P(3,x_1,x_2,x_3) = 6 x_3 + 6 x_2 x_1 + x_1^3, then P[3,P(.,x_1,-x_2,...),x_2,x_3] = 6 x_3 + 6 x_2 P(1,x_1) + P(3,x_1,-x_2,-x_3) = 6 x_3 + 6 x_2 x_1 + 6 (-x_3) + 6 (-x_2) x_1 + x_1^3 = x_1^3.
From the Appell formalism, umbrally [P(.,0,x_2,x_3,...) + y]^n = P(n,y,x_2,x_3,...,x_n).
The indeterminates of the partition polynomials can also be extracted using the Faber polynomials of A263916 with -n * x_n = F(n,S(1,x_1),...,S(n,x_1,...,x_n)) = F(n,P(1,x_1),...,P(n,x_1,...,x_n)/n!). Compare with A263634.
Also P(n,x_1,...,x_n) = ST1(n,x_1,2*x_2,...,n*x_n), where ST1(n,...) are the row partition polynomials of A036039.
(End)

Extensions

Name augmented by Tom Copeland, Dec 08 2022

A187535 Central Lah numbers: a(n) = A105278(2*n,n) = A008297(2*n,n).

Original entry on oeis.org

1, 2, 36, 1200, 58800, 3810240, 307359360, 29682132480, 3339239904000, 428906814336000, 61934143990118400, 9931984545324441600, 1751339941492209868800, 336796142594655744000000, 70149825129001153536000000, 15732267448930658699673600000
Offset: 0

Views

Author

Emanuele Munarini, Mar 11 2011

Keywords

Comments

a(n) is the number of Lah partitions of a set of size 2n with n blocks.

Crossrefs

Programs

  • Maple
    A187535:= n -> if n=0 then 1 else binomial(2*n-1,n-1)*(2*n)!/n! fi;
    seq(A187535(n),n=0..12);
  • Mathematica
    a[n_]:=If[n==0,1,Binomial[2n-1,n-1](2n)!/n!]
    Table[a[n],{n,0,12}]
    (* Alternative: *)
    a[n_] := Binomial[2*n, n] FactorialPower[2*n - 1, n];
    Table[a[n], {n, 0, 15}] (* Peter Luschny, Jun 15 2022 *)
  • Maxima
    a(n) := if n=0 then 1 else binomial(2*n-1,n-1)*(2*n)!/n!;
    makelist(a(n),n,0,12);
    
  • Sage
    [catalan_number(n)*binomial(2*n-1,n)*factorial(n+1) for n in range(15)] # Peter Luschny, Oct 07 2014

Formula

a(n) = binomial(2n-1,n-1)*(2n)!/n! (for n>0).
D-finite with recurrence (n+1)*a(n+1) = 4*(2n+1)^2*a(n) - delta(n,0).
a(n) ~ 2^(4*n)*n^n*exp(-n)/sqrt(2*n*Pi).
a(n)*a(n+2) - a(n+1)^2 is >= 0 and is a multiple of 2^(n+3) for all nonnegative n.
a(n) == 0 (mod 10) for n>3.
E.g.f.: 1/2 + K(16x)/Pi, where K(z) is the complete elliptic integral of the first kind, which can also be written as a Legendre function of the second kind.
a(n) = Catalan(n)*C(2*n-1,n)*(n+1)!. - Peter Luschny, Oct 07 2014
a(n) = A125558(n)*(n+1)! = A090181(2*n,n)*(n+1)!. - Peter Luschny, Oct 07 2014
a(n) = (2/n)*(Gamma(2*n)^2/Gamma(n)^3) for n>0. - Peter Luschny, Oct 17 2014

A132710 Infinitesimal generator for a diagonally-shifted Lah matrix, unsigned A105278, related to n! Laguerre(n,-x,1).

Original entry on oeis.org

0, 2, 0, 0, 6, 0, 0, 0, 12, 0, 0, 0, 0, 20, 0, 0, 0, 0, 0, 30, 0, 0, 0, 0, 0, 0, 42, 0, 0, 0, 0, 0, 0, 0, 56, 0, 0, 0, 0, 0, 0, 0, 0, 72, 0, 0, 0, 0, 0, 0, 0, 0, 0, 90, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 110, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 132, 0
Offset: 0

Views

Author

Tom Copeland, Nov 15 2007, Nov 16 2007, Nov 27 2007

Keywords

Comments

Analogous to the infinitesimal generators of A132681 and A132792.
The matrix T begins
0;
2, 0;
0, 6, 0;
0, 0, 12 0;
0, 0, 0, 20, 0;
Along the nonvanishing diagonal the n-th term is (n+2)*(n+1).
Let LM(t) = exp(t*T) = lim_{n->infinity} (1 + t*T/n)^n.
Shifted Lah matrix = [bin(n+1,k+1)*(n)!/(k)! ] = LM(1) = exp(T). Truncating the series gives the n X n submatrices. In fact, the submatrices of T are nilpotent with [Tsub_n]^(n+1) = 0 for n=0,1,2,....
Inverse shifted Lah matrix = LM(-1) = exp(-T)
Umbrally shifted Lah[b(.)] = exp(b(.)*T) = [ binomial(n+1,k+1)*(n)!/(k)! * b(n-k) ]
A(j) = T^j / j! equals the matrix [binomial(n+1,k+1)*(n)!/(k)! * delta(n-k-j)] where delta(n) = 1 if n=0 and vanishes otherwise (Kronecker delta); i.e. A(j) is a matrix with all the terms 0 except for the j-th lower (or main for j=0) diagonal which equals that of the Lah matrix. Hence the A(j)'s form a linearly independent basis for all matrices of the form [binomial(n+1,k+1) * (n)! / (k)! * d(n-k)].
For sequences with b(0) = 1, umbrally,
LM[b(.)] = exp(b(.)*T) = [ bin(n+1,k+1)*(n)!/(k)! * b(n-k) ] .
[LM[b(.)]]^(-1) = exp(c(.)*T) = [ bin(n+1,k+1)*(n)!/(k)! * c(n-k) ] where c = LPT(b) with LPT the list partition transform of A133314. Or,
[LM[b(.)]]^(-1) = exp[LPT(b(.))*T] = LPT[LM(b(.))] = LM[LPT(b(.))] = LM[c(.)] .
The matrix operation b = T*a can be characterized in several ways in terms of the coefficients a(n) and b(n), their o.g.f.'s A(x) and B(x), or e.g.f.'s EA(x) and EB(x).
1) b(0) = 0, b(n) = (n+1)*(n) * a(n-1),
2) B(x) = x * D^2 * x^2 A(x)
3) B(x) = x * 2 *Lag(2,-:xD:,0) A(x)
4) EB(x) = D * x^2 EA(x)
where D is the derivative w.r.t. x, (:xD:)^j = x^j*D^j and Lag(n,x,m) is the associated Laguerre polynomial of order m.
The exponentiated operator can be characterized (with loose notation) as
5) exp(t*T) * a = LM(t) * a = [sum(k=0,...,n) bin(n+1,k+1) * n!/k! t^(n-k) * a(k)] = [ t^n * n! * Lag(n,-a(.)/t,1) ], a vector array.
With t=1 and a(k) = (-x)^k, then LM(1) * a = [ n! * Laguerre(n,x,1) ], a vector array with index n .
6) exp(t*T) EA(x) = EB(x) = EA[ x / (1-x*t) ] / (1-x*t)^2

Programs

  • Mathematica
    Table[PadLeft[{n*(n-1), 0}, n], {n, 0, 12}] // Flatten (* Jean-François Alcover, Apr 30 2014 *)

Formula

Given a polynomial sequence p_n(x) with p_0(x)=1 and the lowering and raising operators L and R defined by L P_n(x) = n * P_(n-1)(x) and R P_n(x) = P_(n+1)(x), the matrix T represents the action of R*L^2*R^2 in the p_n(x) basis. For p_n(x) = x^n, L = D = d/dx and R = x. For p_n(x) = x^n/n!, L = DxD and R = D^(-1). - Tom Copeland, Oct 25 2012

A079621 Matrix square of unsigned Lah triangle abs(A008297(n,k)) or A105278(n,k).

Original entry on oeis.org

1, 4, 1, 24, 12, 1, 192, 144, 24, 1, 1920, 1920, 480, 40, 1, 23040, 28800, 9600, 1200, 60, 1, 322560, 483840, 201600, 33600, 2520, 84, 1, 5160960, 9031680, 4515840, 940800, 94080, 4704, 112, 1, 92897280, 185794560, 108380160, 27095040, 3386880, 225792
Offset: 1

Views

Author

Vladeta Jovovic, Jan 29 2003

Keywords

Comments

Also the Bell transform of A002866(n+1). For the definition of the Bell transform see A264428. - Peter Luschny, Jan 26 2016
Also the number of k-dimensional flats of the extended Shi arrangement of dimension n consisting of hyperplanes x_i - x_j = d (1 <= i < j <= n, -1 <= d <= 2). - Shuhei Tsujie, Apr 26 2019

Examples

			Triangle begins:
     1;
     4,    1;
    24,   12,   1;
   192,  144,  24,  1;
  1920, 1920, 480, 40, 1;
  ...
		

Crossrefs

Cf. A002866 (first column), A025168 (row sums).

Programs

  • Maple
    # The function BellMatrix is defined in A264428.
    # Adds (1, 0, 0, 0, ..) as column 0.
    BellMatrix(n -> 2^n*(n+1)!, 9); # Peter Luschny, Jan 26 2016
  • Mathematica
    BellMatrix[f_, len_] := With[{t = Array[f, len, 0]}, Table[BellY[n, k, t], {n, 0, len - 1}, {k, 0, len - 1}]];
    B = BellMatrix[2^#*(#+1)!&, rows = 12];
    Table[B[[n, k]], {n, 2, rows}, {k, 2, n}] // Flatten (* Jean-François Alcover, Jun 28 2018, after Peter Luschny *)

Formula

E.g.f.: exp(x*y/(1-2*x)).
T(n, k) = n!/k!*binomial(n-1, k-1)*2^(n-k). - Vladeta Jovovic, Sep 24 2003
The n-th row polynomial equals x o (x + 2) o (x + 4) o ... o (x + 2*n), where o is the deformed Hadamard product of power series defined in Bala, section 3.1. - Peter Bala, Jan 18 2018

A276965 Square row sums of the triangle of Lah numbers (A105278).

Original entry on oeis.org

1, 1, 5, 73, 2017, 86801, 5289301, 430814665, 45052534913, 5868875082817, 930114039075301, 175964489469769001, 39125942325820605025, 10092849114680961297553, 2987365449592984040715317, 1005030253302269078318250601
Offset: 0

Views

Author

Emanuele Munarini, Sep 27 2016

Keywords

Crossrefs

Programs

  • Mathematica
    Table[HypergeometricPFQ[{1-n,1-n,-n,-n},{1},1],{n,0,100}]
  • Maxima
    makelist(hypergeometric([-n+1,-n+1,-n,-n],[1],1),n,0,12);
    
  • PARI
    concat([1], for(n=1,25, print1(sum(k=0,n, binomial(n,k)^2*binomial(n-1,k-1)^2*((n-k)!)^2), ", "))) \\ G. C. Greubel, Jun 05 2017
  • Perl
    use ntheory ":all"; for my $n (0..20) { say "$n ",vecsum(map{my $l=stirling($n,$,3); vecprod($l,$l); } 0..$n) } # _Dana Jacobsen, Mar 16 2017
    

Formula

a(n) = Sum_{k=0..n} lah(n,k)^2.
a(n) = Sum_{k=0..n} binomial(n,k)^2*binomial(n-1,k-1)^2*((n-k)!)^2.
a(n) = hypergeometric([-n+1,-n+1,-n,-n],[1],1).
a(n) = (n!)^2 * hypergeometric([-n+1,-n+1],[1,2,2],1) for n > 0.
Recurrence: n*(16*n^3 - 96*n^2 + 185*n - 116)*a(n) = 2*(32*n^6 - 272*n^5 + 930*n^4 - 1668*n^3 + 1670*n^2 - 867*n + 164)*a(n-1) - (n-2)*(96*n^7 - 1056*n^6 + 4646*n^5 - 10500*n^4 + 12990*n^3 - 8644*n^2 + 2827*n - 364)*a(n-2) + 2*(n-3)*(n-2)^3*(32*n^6 - 336*n^5 + 1410*n^4 - 2978*n^3 + 3268*n^2 - 1731*n + 353)*a(n-3) - (n-4)^2*(n-3)^3*(n-2)^4*(16*n^3 - 48*n^2 + 41*n - 11)*a(n-4). - Vaclav Kotesovec, Sep 27 2016
a(n) ~ n^(2*n - 3/4) * exp(4*sqrt(n) - 2*n - 1) / (2^(3/2) * sqrt(Pi)) * (1 + 31/(96*sqrt(n)) + 937/(18432*n)). - Vaclav Kotesovec, Sep 27 2016

A308281 The third power of the unsigned Lah triangular matrix A105278.

Original entry on oeis.org

1, 6, 1, 54, 18, 1, 648, 324, 36, 1, 9720, 6480, 1080, 60, 1, 174960, 145800, 32400, 2700, 90, 1, 3674160, 3674160, 1020600, 113400, 5670, 126, 1, 88179840, 102876480, 34292160, 4762800, 317520, 10584, 168, 1, 2380855680, 3174474240, 1234517760, 205752960, 17146080, 762048, 18144, 216, 1
Offset: 1

Views

Author

Shuhei Tsujie, May 18 2019

Keywords

Comments

Also the number of k-dimensional flats of the extended Shi arrangement of dimension n consisting of hyperplanes x_i - x_j = d (1 <= i < j <= n, -2 <= d <= 3).

Examples

			Triangle begins:
     1;
     6,    1;
    54,   18,    1;
   648,  324,   36,  1;
  9720, 6480, 1080, 60, 1;
  ...
		

Crossrefs

Cf. A105278.

Programs

  • Mathematica
    Table[3^(n - k) * Binomial[n - 1, k - 1] * n! / k!, {n, 1, 10}, {k, 1, n}] // Flatten (* Amiram Eldar, Jul 13 2019 *)

Formula

E.g.f.: exp(x*y/(1-3*x)).
T(n,k) = 3^(n-k)*binomial(n-1, k-1)*n!/k! = 3^(n-k)*A105278.

A308282 The fifth power of the unsigned Lah triangular matrix A105278.

Original entry on oeis.org

1, 10, 1, 150, 30, 1, 3000, 900, 60, 1, 75000, 30000, 3000, 100, 1, 2250000, 1125000, 150000, 7500, 150, 1, 78750000, 47250000, 7875000, 525000, 15750, 210, 1, 3150000000, 2205000000, 441000000, 36750000, 1470000, 29400, 280, 1, 141750000000, 113400000000, 26460000000, 2646000000, 132300000, 3528000, 50400, 360, 1
Offset: 1

Views

Author

Shuhei Tsujie, May 18 2019

Keywords

Comments

Also the number of k-dimensional flats of the extended Shi arrangement of dimension n consisting of hyperplanes x_i - x_j = d (1 <= i < j <= n, -4 <= d <= 5).

Examples

			Triangle begins:
      1;
     10,     1;
    150,    30,    1;
   3000,   900,   60,   1;
  75000, 30000, 3000, 100, 1;
  ...
		

Crossrefs

Cf. A105278.

Programs

  • Mathematica
    Table[5^(n - k) * Binomial[n - 1, k - 1] * n! / k!, {n, 1, 10}, {k, 1, n}] // Flatten (* Amiram Eldar, Jul 13 2019 *)

Formula

E.g.f.: exp(x*y/(1-5*x)).
T(n,k) = 5^(n-k)*binomial(n-1, k-1)*n!/k! = 5^(n-k)*A105278.

A308440 Matrix product of triangle of Stirling numbers of second kind A008277 and square of unsigned Lah triangle A105278.

Original entry on oeis.org

1, 5, 1, 37, 15, 1, 365, 223, 30, 1, 4501, 3675, 745, 50, 1, 66605, 68071, 18450, 1865, 75, 1, 1149877, 1411515, 479101, 64750, 3920, 105, 1, 22687565, 32512663, 13260030, 2244501, 181650, 7322, 140, 1, 503589781, 825175275, 393017185, 79948050, 8103711, 436590, 12558, 180, 1
Offset: 1

Views

Author

Shuhei Tsujie, May 27 2019

Keywords

Comments

Also the number of k-dimensional flats of the extended Catalan arrangement of dimension n consisting of hyperplanes x_i - x_j = d (1 <= i < j <= n, -2 <= d <= 2).

Examples

			Triangle begins:
     1;
     5,    1;
    37,   15,   1;
   365,  223,  30,  1;
  4501, 3675, 745, 50, 1;
  ...
		

Crossrefs

Cf. A008277, A105278, A050351 (first column), A109092 (row sums).

Formula

E.g.f.: exp((exp(x)-1)*y/(3-2exp(x))).

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

A008297 Triangle of Lah numbers.

Original entry on oeis.org

-1, 2, 1, -6, -6, -1, 24, 36, 12, 1, -120, -240, -120, -20, -1, 720, 1800, 1200, 300, 30, 1, -5040, -15120, -12600, -4200, -630, -42, -1, 40320, 141120, 141120, 58800, 11760, 1176, 56, 1, -362880, -1451520, -1693440, -846720, -211680, -28224, -2016, -72, -1, 3628800, 16329600, 21772800, 12700800
Offset: 1

Views

Author

Keywords

Comments

|a(n,k)| = number of partitions of {1..n} into k lists, where a list means an ordered subset.
Let N be a Poisson random variable with parameter (mean) lambda, and Y_1,Y_2,... independent exponential(theta) variables, independent of N, so that their density is given by (1/theta)*exp(-x/theta), x > 0. Set S=Sum_{i=1..N} Y_i. Then E(S^n), i.e., the n-th moment of S, is given by (theta^n) * L_n(lambda), n >= 0, where L_n(y) is the Lah polynomial Sum_{k=0..n} |a(n,k)| * y^k. - Shai Covo (green355(AT)netvision.net.il), Feb 09 2010
For y = lambda > 0, formula 2) for the Lah polynomial L_n(y) dated Feb 02 2010 can be restated as follows: L_n(lambda) is the n-th ascending factorial moment of the Poisson distribution with parameter (mean) lambda. - Shai Covo (green355(AT)netvision.net.il), Feb 10 2010
See A111596 for an expression of the row polynomials in terms of an umbral composition of the Bell polynomials and relation to an inverse Mellin transform and a generalized Dobinski formula. - Tom Copeland, Nov 21 2011
Also the Bell transform of the sequence (-1)^(n+1)*(n+1)! without column 0. For the definition of the Bell transform see A264428. - Peter Luschny, Jan 28 2016
Named after the Slovenian mathematician and actuary Ivo Lah (1896-1979). - Amiram Eldar, Jun 13 2021

Examples

			|a(2,1)| = 2: (12), (21); |a(2,2)| = 1: (1)(2). |a(4,1)| = 24: (1234) (24 ways); |a(4,2)| = 36: (123)(4) (6*4 ways), (12)(34) (3*4 ways); |a(4,3)| = 12: (12)(3)(4) (6*2 ways); |a(4,4)| = 1: (1)(2)(3)(4) (1 way).
Triangle:
    -1;
     2,    1;
    -6,   -6,   -1;
    24,   36,   12,   1;
  -120, -240, -120, -20, -1; ...
		

References

  • Louis Comtet, Advanced Combinatorics, Reidel, 1974, p. 156.
  • Shai Covo, The moments of a compound Poisson process with exponential or centered normal jumps, J. Probab. Stat. Sci., Vol. 7, No. 1 (2009), pp. 91-100.
  • Theodore S. Motzkin, Sorting numbers for cylinders and other classification numbers, in Combinatorics, Proc. Symp. Pure Math. 19, AMS, 1971, pp. 167-176; the sequence called {!}^{n+}. For a link to this paper see A000262.
  • John Riordan, An Introduction to Combinatorial Analysis, Wiley, 1958, p. 44.
  • S. Gill Williamson, Combinatorics for Computer Science, Computer Science Press, 1985; see p. 176.

Crossrefs

Same as A066667 and A105278 except for signs. Other variants: A111596 (differently signed triangle and (0,0)-based), A271703 (unsigned and (0,0)-based), A089231.
A293125 (row sums) and A000262 (row sums of unsigned triangle).
Columns 1-6 (unsigned): A000142, A001286, A001754, A001755, A001777, A001778.
A002868 gives maximal element (in magnitude) in each row.
A248045 (central terms, negated). A130561 is a natural refinement.

Programs

  • Haskell
    a008297 n k = a008297_tabl !! (n-1) !! (k-1)
    a008297_row n = a008297_tabl !! (n-1)
    a008297_tabl = [-1] : f [-1] 2 where
       f xs i = ys : f ys (i + 1) where
         ys = map negate $
              zipWith (+) ([0] ++ xs) (zipWith (*) [i, i + 1 ..] (xs ++ [0]))
    -- Reinhard Zumkeller, Sep 30 2014
    
  • Maple
    A008297 := (n,m) -> (-1)^n*n!*binomial(n-1,m-1)/m!;
  • Mathematica
    a[n_, m_] := (-1)^n*n!*Binomial[n-1, m-1]/m!; Table[a[n, m], {n, 1, 10}, {m, 1, n}] // Flatten (* Jean-François Alcover, Dec 12 2012, after Maple *)
    T[n_, n_] := (-1)^n; T[n_, k_]/;0Oliver Seipel, Dec 06 2024 *)
  • PARI
    T(n, m) = (-1)^n*n!*binomial(n-1, m-1)/m!
    for(n=1,9, for(m=1,n, print1(T(n,m)", "))) \\ Charles R Greathouse IV, Mar 09 2016
    
  • Perl
    use bigint; use ntheory ":all"; my @L; for my $n (1..9) { push @L, map { stirling($n,$,3)*(-1)**$n } 1..$n; } say join(", ",@L); # _Dana Jacobsen, Mar 16 2017
  • Sage
    def A008297_triangle(dim): # computes unsigned T(n, k).
        M = matrix(ZZ,dim,dim)
        for n in (0..dim-1): M[n,n] = 1
        for n in (1..dim-1):
            for k in (0..n-1):
                M[n,k] = M[n-1,k-1]+(2+2*k)*M[n-1,k]+((k+1)*(k+2))*M[n-1,k+1]
        return M
    A008297_triangle(9) # Peter Luschny, Sep 19 2012
    

Formula

a(n, m) = (-1)^n*n!*A007318(n-1, m-1)/m!, n >= m >= 1.
a(n+1, m) = (n+m)*a(n, m)+a(n, m-1), a(n, 0) := 0; a(n, m) := 0, n < m; a(1, 1)=1.
a(n, m) = ((-1)^(n-m+1))*L(1, n-1, m-1) where L(1, n, m) is the triangle of coefficients of the generalized Laguerre polynomials n!*L(n, a=1, x). These polynomials appear in the radial l=0 eigen-functions for discrete energy levels of the H-atom.
|a(n, m)| = Sum_{k=m..n} |A008275(n, k)|*A008277(k, m), where A008275 = Stirling numbers of first kind, A008277 = Stirling numbers of second kind. - Wolfdieter Lang
If L_n(y) = Sum_{k=0..n} |a(n, k)|*y^k (a Lah polynomial) then the e.g.f. for L_n(y) is exp(x*y/(1-x)). - Vladeta Jovovic, Jan 06 2001
E.g.f. for the k-th column (unsigned): x^k/(1-x)^k/k!. - Vladeta Jovovic, Dec 03 2002
a(n, k) = (n-k+1)!*N(n, k) where N(n, k) is the Narayana triangle A001263. - Philippe Deléham, Jul 20 2003
From Shai Covo (green355(AT)netvision.net.il), Feb 02 2010: (Start)
We have the following expressions for the Lah polynomial L_n(y) = Sum_{k=0..n} |a(n, k)|*y^k -- exact generalizations of results in A000262 for A000262(n) = L_n(1):
1) L_n(y) = y*exp(-y)*n!*M(n+1,2,y), n >= 1, where M (=1F1) is the confluent hypergeometric function of the first kind;
2) L_n(y) = exp(-y)* Sum_{m>=0} y^m*[m]^n/m!, n>=0, where [m]^n = m*(m+1)*...*(m+n-1) is the rising factorial;
3) L_n(y) = (2n-2+y)L_{n-1}(y)-(n-1)(n-2)L_{n-2}(y), n>=2;
4) L_n(y) = y*(n-1)!*Sum_{k=1..n} (L_{n-k}(y) k!)/((n-k)! (k-1)!), n>=1. (End)
The row polynomials are given by D^n(exp(-x*t)) evaluated at x = 0, where D is the operator (1-x)^2*d/dx. Cf. A008277 and A035342. - Peter Bala, Nov 25 2011
n!C(-xD,n) = Lah(n,:xD:) where C(m,n) is the binomial coefficient, xD= x d/dx, (:xD:)^k = x^k D^k, and Lah(n,x) are the row polynomials of this entry. E.g., 2!C(-xD,2)= 2 xD + x^2 D^2. - Tom Copeland, Nov 03 2012
From Tom Copeland, Sep 25 2016: (Start)
The Stirling polynomials of the second kind A048993 (A008277), i.e., the Bell-Touchard-exponential polynomials B_n[x], are umbral compositional inverses of the Stirling polynomials of the first kind signed A008275 (A130534), i.e., the falling factorials, (x)_n = n! binomial(x,n); that is, umbrally B_n[(x).] = x^n = (B.[x])_n.
An operational definition of the Bell polynomials is (xD_x)^n = B_n[:xD:], where, by definition, (:xD_x:)^n = x^n D_x^n, so (B.[:xD_x:])_n = (xD_x)_n = :xD_x:^n = x^n (D_x)^n.
Let y = 1/x, then D_x = -y^2 D_y; xD_x = -yD_y; and P_n(:yD_y:) = (-yD_y)_n = (-1)^n (1/y)^n (y^2 D_y)^n, the row polynomials of this entry in operational form, e.g., P_3(:yD_y:) = (-yD_y)_3 = (-yD_y) (yD_y-1) (yD_y-2) = (-1)^3 (1/y)^3 (y^2 D_y)^3 = -( 6 :yD_y: + 6 :yD_y:^2 + :yD_y:^3 ) = - ( 6 y D_y + 6 y^2 (D_y)^2 + y^3 (D_y)^3).
Therefore, P_n(y) = e^(-y) P_n(:yD_y:) e^y = e^(-y) (-1/y)^n (y^2 D_y)^n e^y = e^(-1/x) x^n (D_x)^n e^(1/x) = P_n(1/x) and P_n(x) = e^(-1/x) x^n (D_x)^n e^(1/x) = e^(-1/x) (:x D_x:)^n e^(1/x). (Cf. also A094638.) (End)
T(n,k) = Sum_{j=k..n} (-1)^j*A008296(n,j)*A360177(j,k). - Mélika Tebni, Feb 02 2023
Showing 1-10 of 43 results. Next