A004747 -id:A004747 - cp's OEIS frontend

A015735 Row sums of triangle A004747.

1, 3, 17, 145, 1661, 23931, 415773, 8460257, 197360985, 5192853011, 152137882601, 4911873672113, 173268075672277, 6630323916472075, 273555262963272501, 12105084133976359361, 571897644855277242673, 28731255563712689630627, 1529450942687399074134465

Offset: 1

Wolfdieter Lang

Vincenzo Librandi, Table of n, a(n) for n = 1..200
Wolfdieter Lang, On generalizations of Stirling number triangles, J. Integer Seqs., Vol. 3 (2000), #00.2.4.

Cf. A001515, A004747, A016036, A028575.

R:=PowerSeriesRing(Rationals(), 30); Coefficients(R!(Laplace( Exp(1-(1-3*x)^(1/3)) - 1 ))); // G. C. Greubel, Oct 02 2023

a[1]=1; a[n_]:= 1 +(n-1)!*Sum[Binomial[k, n-m-k]*Binomial[k+n-1,n-1]*(-1/3)^(n-m-k)/(m-1)!, {m,n}, {k,n-m}]; Table[a[n], {n,20}] (* Jean-François Alcover, Jul 05 2013, after Vladimir Kruchinin *)
Rest@With[{m=30}, CoefficientList[Series[Exp[1-Surd[1-3*x,3]] -1, {x, 0,m}], x]*Range[0,m]!] (* G. C. Greubel, Oct 02 2023 *)

a(n):=if n=1 then 1 else (n-1)!*sum(sum(binomial(k,n-m-k)* (-1/3)^(n-m-k)*binomial(k+n-1,n-1),k,1,n-m)/(m-1)!,m,1,n)+1; /* Vladimir Kruchinin, Aug 08 2010 */

def A015735_list(prec):
    P. = PowerSeriesRing(QQ, prec)
    return P( exp(1-(1-3*x)^(1/3)) -1 ).egf_to_ogf().list()
a=A015735_list(40); a[1:] # G. C. Greubel, Oct 02 2023

E.g.f.: exp(1-(1-3*x)^(1/3)) - 1, if one takes a(0)=0.
a(n) = 6*(n-2)*a(n-1) - (3*n-8)*(3*n-7)*a(n-2) + a(n-3), a(0)=1, a(1)=1, a(2)=3.
a(n) = 1 + (n-1)!*Sum_{m=1..n} ( Sum_{k=1..n-m} C(k, n-m-k)*C(k+n-1, n-1)*(-1/3)^(n-m-k) ) / (m-1)!, n > 1. - Vladimir Kruchinin, Aug 08 2010
a(n) = D^n(exp(x)) evaluated at x = 0, where D is the operator 1/(1-x)^2*d/dx. Cf. A001515, A016036 and A028575. - Peter Bala, Nov 25 2011
E.g.f. with offset 0: exp(1-(1-3*x)^(1/3))/(1-3*x)^(2/3). - Sergei N. Gladkovskii, Jul 07 2012.
a(n) ~ sqrt(2*Pi)*3^(n-1)*exp(1-n)*n^(n-5/6)/Gamma(2/3) * (1-sqrt(3)*Gamma(2/3)^2/(2*Pi*n^(1/3))). - Vaclav Kotesovec, Aug 10 2013
From Seiichi Manyama, Jan 20 2025: (Start)
a(n) = Sum_{k=0..n} (-1)^k * 3^(n-k) * |Stirling1(n,k)| * A000587(k).
a(n) = e * (-3)^n * n! * Sum_{k>=0} (-1)^k * binomial(k/3,n)/k!. (End)

A143172 Partition number array, called M32(-2), related to A004747(n,m) = |S2(-2;n,m)| (generalized Stirling triangle).

Original entry on oeis.org

1, 2, 1, 10, 6, 1, 80, 40, 12, 12, 1, 880, 400, 200, 100, 60, 20, 1, 12320, 5280, 2400, 1000, 1200, 1200, 120, 200, 180, 30, 1, 209440, 86240, 36960, 28000, 18480, 16800, 7000, 4200, 2800, 4200, 840, 350, 420, 42, 1, 4188800, 1675520, 689920, 492800, 224000, 344960

Offset: 1

Wolfdieter Lang, Oct 09 2008

Each partition of n, ordered as in Abramowitz-Stegun (A-St order; for the reference see A134278), is mapped to a nonnegative integer a(n,k)=:M32(-2;n,k) with the k-th partition of n in A-St order.
The sequence of row lengths is A000041 (partition numbers) [1, 2, 3, 5, 7, 11, 15, 22, 30, 42, ...].
a(n,k) enumerates special unordered forests related to the k-th partition of n in the A-St order. The k-th partition of n is given by the exponents enk =(e(n,k,1),...,e(n,k,n)) of 1,2,...n. The number of parts is m = sum(e(n,k,j),j=1..n). The special (enk)-forest is composed of m rooted increasing (r+1)-ary trees if the outdegree is r>=0.
If M32(-2;n,k) is summed over those k with fixed number of parts m one obtains triangle A004747(n,m)= |S2(-2;n,m)|, a generalization of Stirling numbers of the second kind. For S2(K;n,m), K from the integers, see the reference under A035342.

			a(4,3)=12. The relevant partition of 4 is (2^2). The 12 unordered (0,2,0,0)-forests are composed of the following 2 rooted increasing trees 1--2,3--4; 1--3,2--4 and 1--4,2--3. The trees are binary because r=1 vertices are binary (2-ary) and for the leaves (r=0) the arity does not matter. Each of the three differently labeled forests comes therefore in 4 versions due to the two binary root vertices.

W. Lang, First 10 rows of the array and more.
W. Lang, Combinatorial Interpretation of Generalized Stirling Numbers, J. Int. Seqs. Vol. 12 (2009) 09.3.3.

Cf. A143171 (M32(-1) array), A143173 (M32(-3) array).

a(n,k)= (n!/product(e(n,k,j)!*j!^(e(n,k,j),j=1..n))*product(|S2(-2,j,1)|^e(n,k,j),j=1..n) = M3(n,k)*product(|S2(-2,j,1)|^e(n,k,j),j=1..n), with |S2(-2,n,1)|= A008544(n-1) = (3*n-4)(!^3) (3-factorials) for n>=2 and 1 if n=1 and 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. M3(n,k):= A036040(n,k), k=1..p(n), p(n):= A000041(n).

A144274 Partition number array, called M32hat(-2)= 'M32(-2)/M3'= 'A143172/A036040', related to A004747(n,m)= |S2(-2;n,m)| (generalized Stirling triangle).

Original entry on oeis.org

1, 2, 1, 10, 2, 1, 80, 10, 4, 2, 1, 880, 80, 20, 10, 4, 2, 1, 12320, 880, 160, 100, 80, 20, 8, 10, 4, 2, 1, 209440, 12320, 1760, 800, 880, 160, 100, 40, 80, 20, 8, 10, 4, 2, 1, 4188800, 209440, 24640, 8800, 6400, 12320, 1760, 800, 320, 200, 880, 160, 100, 40, 16, 80, 20

Offset: 1

Wolfdieter Lang, Oct 09 2008

Each partition of n, ordered as in Abramowitz-Stegun (A-St order; for the reference see A134278), is mapped to a nonnegative integer a(n,k) =: M32hat(-2;n,k) with the k-th partition of n in A-St order.
The sequence of row lengths is A000041 (partition numbers) [1, 2, 3, 5, 7, 11, 15, 22, 30, 42,...].
If M32hat(-2;n,k) is summed over those k with fixed number of parts m one obtains triangle S2hat(-2):= A144275(n,m).

			a(4,3) = 4 = |S2(-2,2,1)|^2. The relevant partition of 4 is (2^2).

Wolfdieter Lang, First 10 rows of the array and more.
Wolfdieter Lang, Combinatorial Interpretation of Generalized Stirling Numbers, J. Int. Seqs. Vol. 12 (2009) 09.3.3.

Cf. A144269 (M32hat(-1) array). A144279 (M32hat(-3) array).

a(n,k) = Product_{j=1..n} |S2(-2,j,1)|^e(n,k,j) with |S2(-2,n,1)|= A008544(n-1) = (3*n-4)(!^3) (3-factorials) for n>=2 and 1 if n=1 and 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.

Formally a(n,k)= 'M32(-2)/M3' = 'A143172/A036040' (elementwise division of arrays).

A144345 Second column (m=2) of triangle S2p(-2) = A004747.

Original entry on oeis.org

1, 6, 52, 600, 8680, 151200, 3082240, 71998080, 1896294400, 55601145600, 1796277683200, 63397990656000, 2427084884224000, 100175046107136000, 4434284662872064000, 209554432423784448000, 10530302071553904640000, 560682451860226375680000

Offset: 0

Wolfdieter Lang, Oct 09 2008

Cf. A004747, A008544(n-1) (m=1 column), A144346 (m=3 column).

a(n) = A004747(n+2,2), n>=0.

A144346 Third column (m=3) of triangle S2p(-2) = A004747.

Original entry on oeis.org

1, 12, 160, 2520, 46480, 987840, 23826880, 643843200, 19280060800, 634002969600, 22718375680000, 881259515136000, 36796205974528000, 1645615697037312000, 78486991029551104000, 3976930001842237440000, 213353732165508014080000, 12081783988797659136000000

Offset: 0

Wolfdieter Lang, Oct 09 2008

Cf. A004747, A144345 (m=2 column).

a(n) = A004747(n+3,3), n>=0.

A008544 Triple factorial numbers: Product_{k=0..n-1} (3*k+2).

Original entry on oeis.org

1, 2, 10, 80, 880, 12320, 209440, 4188800, 96342400, 2504902400, 72642169600, 2324549427200, 81359229952000, 3091650738176000, 126757680265216000, 5577337931669504000, 262134882788466688000

Offset: 0

Joe Keane (jgk(AT)jgk.org)

a(n-1), n >= 1, enumerates increasing plane (aka ordered) trees with n vertices (one of them a root labeled 1) where each vertex with outdegree r >= 0 comes in r+1 types (like an (r+1)-ary vertex). See the increasing tree comments under A004747. - Wolfdieter Lang, Oct 12 2007
An example for the case of 3 vertices is shown below. For the enumeration of non-plane trees of this type see A029768. - Peter Bala, Aug 30 2011
a(n) is the product of the positive integers k <= 3*n that have k modulo 3 = 2. - Peter Luschny, Jun 23 2011
See A094638 for connections to differential operators. - Tom Copeland, Sep 20 2011
Partial products of A016789. - Reinhard Zumkeller, Sep 20 2013
The Mathar conjecture is true. Generally from the factorial form, the last term is the "extra" product beyond the prior term, from k=n-1 and 3k+2 evaluates to 3*(n-1)+2 = 3n-1, yielding a(n) = a(n-1)*(3n-1) (eqn1). Similarly, a(n) = a(n-2)*(3n-1)*(3(n-2)+2) = a(n-2)*(3n-1)*(3n-4) (eqn2) and a(n) = a(n-3)*(3n-1)*(3n-4)*(3*(n-2)+2) = a(n-3)*(3n-1)*(3n-4)*(3n-7) (eqn3). We equate (eqn2) and (eqn3) to get a(n-2)*(3n-1)*(3n-4) = a(n-3)*(3n-1)*(3n-4)*(3n-7) or a(n-2)+(7-3n)*a(n-3) = 0 (eqn4). From (eqn1) we have a(n)+(1-3n)*a(n-1) = 0 (eqn5). Combining (eqn4) and (eqn5) yields a(n)+(1-3n)*a(n-1)+a(n-2)+(7-3n)*a(n-3) = 0. - Bill McEachen, Jan 01 2016
a(n-1), n>=1, is the dimension of the n-th component of the operad encoding the multilinearization of the following identity in nonassociative algebras: s*(a,a,b)-(s+t)*(a,b,a)+t*(b,a,a)=0, for any given pair of scalars (s,t). Here (a,b,c) is the associator (ab)c-a(bc). This is proved in the referenced article on associator dependent algebras by Bremner and me. - Vladimir Dotsenko, Mar 22 2022

			a(2) = 10 from the described trees with 3 vertices: there are three trees with a root vertex (label 1) with outdegree r=2 (like the three 3-stars each with one different ray missing) and the four trees with a root (r=1 and label 1) a vertex with (r=1) and a leaf (r=0). Assigning labels 2 and 3 yields 2*3+4=10 such trees.
a(2) = 10. The 10 possible plane increasing trees on 3 vertices, where vertices of outdegree 1 come in 2 colors (denoted a or b) and vertices of outdegree 2 come in 3 colors (a, b or c), are:
.
   1a    1b    1a    1b        1a       1b       1c
   |     |     |     |        / \      / \      / \
   2a    2b    2b    2a      2   3    2   3    2   3
   |     |     |     |
   3     3     3     3         1a       1b       1c
                              / \      / \      / \
                             3   2    3   2    3   2

T. D. Noe, Table of n, a(n) for n = 0..100
Murray Bremner and Vladimir Dotsenko, Associator dependent algebras and Koszul duality, arXiv:2203.11142 [math.RA], 2022.
Wolfdieter Lang, On generalizations of Stirling number triangles, J. Integer Seq., Vol. 3 (2000), Article 00.2.4.
Keiichi Shigechi, On the lattice of weighted partitions, arXiv:2212.14666 [math.CO], 2022. See p. 27.
Garrett Southwood and Hua Wang, The Statistics and Combinatorics of Increasing Trees and Colored Increasing Trees, J. Combin. Math. Combin. Comput. (2024) Vol. 123, 475-487. See p. 485.

a(n) = A004747(n+1, 1) (first column of triangle). Cf. A051141.
Cf. A000165, A001813, A047055, A047657, A084947, A084948, A084949.
Cf. A029768, A034724, A049308.
Cf. A024396, A193534.
Cf. A225470, A290596 (first columns).
Subsequence of A007661.

a008544 n = a008544_list !! n
a008544_list = scanl (*) 1 a016789_list
-- Reinhard Zumkeller, Sep 20 2013

[Round((Gamma(2*n-5/3)/Gamma(n-5/6)*Gamma(2/3)/Gamma(5/6) )/ Sqrt(3)*3^n/4^(n-1)): n in [1..20]]; // Vincenzo Librandi, Feb 21 2015

[Round(3^n*Gamma(n+2/3)/Gamma(2/3)): n in [0..20]]; // G. C. Greubel, Mar 31 2019

a := n -> mul(3*k-1, k = 1..n);
A008544 := n -> mul(k, k = select(k-> k mod 3 = 2, [$1 .. 3*n])): seq(A008544(n), n = 0 .. 16); # Peter Luschny, Jun 23 2011

k = 3; b[1]=2; b[n_]:= b[n] = b[n-1]+k; a[0]=1; a[1]=2; a[n_]:= a[n] = a[n-1]*b[n]; Table[a[n], {n,0,20}] (* Roger L. Bagula, Sep 17 2008 *)
Product[3 k + 2, {k, 0, # - 1}] & /@ Range[0, 16] (* Michael De Vlieger, Jan 02 2016 *)
Table[3^n*Pochhammer[2/3, n], {n,0,20}] (* G. C. Greubel, Mar 31 2019 *)

a(n):=((n)!*sum(binomial(k,n-k)*binomial(n+k,k)*3^(-n+k)*(-1)^(n-k),k,floor(n/2),n)); /* Vladimir Kruchinin, Sep 28 2013 */

PARI
```
a(n) = prod(k=0,n-1, 3*k+2 );
```

vector(20, n, n--; round(3^n*gamma(n+2/3)/gamma(2/3))) \\ G. C. Greubel, Mar 31 2019

@CachedFunction
def A008544(n): return 1 if n == 0 else (3*n-1)*A008544(n-1)
[A008544(n) for n in (0..16)]  # Peter Luschny, May 20 2013

[3^n*rising_factorial(2/3, n) for n in (0..20)] # G. C. Greubel, Mar 31 2019

a(n) = Product_{k=0..n-1} (3*k+2) = A007661(3*n-1) (with A007661(-1) = 1).
E.g.f.: (1-3*x)^(-2/3).
a(n) = 2*A034000(n), n >= 1, a(0) = 1.
a(n) ~ 2^(1/2)*Pi^(1/2)*Gamma(2/3)^-1*n^(1/6)*3^n*e^-n*n^n*{1 - 1/36*n^-1 + ...}. - Joe Keane (jgk(AT)jgk.org), Nov 22 2001
a(n) = (Gamma(2*n-5/3)/Gamma(n-5/6)*Gamma(2/3)/Gamma(5/6))/sqrt(3)*3^n/4^(n-1). - Jeremy L. Martin, Mar 31 2002 (typo fixed by Vincenzo Librandi, Feb 21 2015)
From Daniel Dockery (peritus(AT)gmail.com), Jun 13 2003: (Start)
a(n) = A084939(n)/A000142(n)*A000079(n).
a(n) = 3^n*Pochhammer(2/3, n) = 3^n*Gamma(n+2/3)/Gamma(2/3). (End)
Let T = A094638 and c(t) = column vector(1, t, t^2, t^3, t^4, t^5,...), then A008544 = unsigned [ T * c(-3) ] and the list partition transform A133314 of [1,T * c(-3)] gives [1,T * c(3)] with all odd terms negated, which equals a signed version of A007559; i.e., LPT[(1,signed A008544)] = signed A007559. Also LPT[A007559] = (1,-A008544) and e.g.f. [1,T * c(t)] = (1-x*t)^(-1/t) for t = 3 or -3. Analogous results hold for the double factorial, quadruple factorial and so on. - Tom Copeland, Dec 22 2007
G.f.: 1/(1-2x/(1-3x/(1-5x/(1-6x/(1-8x/(1-9x/(1-11x/(1-12x/(1-...))))))))) (continued fraction). - Philippe Deléham, Jan 08 2012
a(n) = (-1)^n*Sum_{k=0..n} 3^k*s(n+1,n+1-k), where s(n,k) are the Stirling numbers of the first kind, A048994. - Mircea Merca, May 03 2012
G.f.: 1/Q(0) where Q(k) = 1 - x*(3*k+2)/(1 - x*(3*k+3)/Q(k+1) ); (continued fraction). - Sergei N. Gladkovskii, Mar 20 2013
G.f.: G(0)/2, where G(k) = 1 + 1/(1 - x*(3*k+2)/(x*(3*k+2) + 1/G(k+1))); (continued fraction). - Sergei N. Gladkovskii, May 25 2013
D-finite with recurrence: a(n) = (9*(n-2)*(n-1)+2)*a(n-2) + 4*a(n-1), n>=2. - Ivan N. Ianakiev, Aug 09 2013
a(n) = n!*Sum_{k=floor(n/2)..n} binomial(k,n-k)*binomial(n+k,k)*3^(-n+k)*(-1)^(n-k). - Vladimir Kruchinin, Sep 28 2013
Recurrence equation: a(n) = 3*a(n-1) + (3*n - 4)^2*a(n-2) with a(0) = 1 and a(1) = 2. A024396 satisfies the same recurrence (but with different initial conditions). This observation leads to a continued fraction expansion for the constant A193534 due to Euler. - Peter Bala, Feb 20 2015
a(n) = A225470(n, 0), n >= 0. - Wolfdieter Lang, May 29 2017
G.f.: Hypergeometric2F0(1, 2/3; -; 3*x). - G. C. Greubel, Mar 31 2019
D-finite with recurrence: a(n) + (-3*n+1)*a(n-1)=0. - R. J. Mathar, Jan 17 2020
G.f.: 1/(1-2*x-6*x^2/(1-8*x-30*x^2/(1-14*x-72*x^2/(1-20*x-132*x^2/(1-...))))) (Jacobi continued fraction). - Nikolaos Pantelidis, Feb 28 2020
G.f.: 1/G(0), where G(k) = 1 - (6*k+2)*x - 3*(k+1)*(3*k+2)*x^2/G(k+1). - Nikolaos Pantelidis, Feb 28 2020
Sum_{n>=0} 1/a(n) = 1 + (e/3)^(1/3) * (Gamma(2/3) - Gamma(2/3, 1/3)). - Amiram Eldar, Mar 01 2022

A105278 Triangle read by rows: T(n,k) = binomial(n,k)*(n-1)!/(k-1)!.

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

Offset: 1

Miklos Kristof, Apr 25 2005

T(n,k) is the number of partially ordered sets (posets) on n elements that consist entirely of k chains. For example, T(4, 3)=12 since there are exactly 12 posets on {a,b,c,d} that consist entirely of 3 chains. Letting ab denote a<=b and using a slash "/" to separate chains, the 12 posets can be given by a/b/cd, a/b/dc, a/c/bd, a/c/db, a/d/bc, a/d/cb, b/c/ad, b/c/da, b/d/ac, b/d/ca, c/d/ab and c/d/ba, where the listing of the chains is arbitrary (e.g., a/b/cd = a/cd/b =...cd/b/a). - Dennis P. Walsh, Feb 22 2007
Also the matrix product |S1|.S2 of Stirling numbers of both kinds.
This Lah triangle is a lower triangular matrix of the Jabotinsky type. See the column e.g.f. and the D. E. Knuth reference given in A008297. - Wolfdieter Lang, Jun 29 2007
The infinitesimal matrix generator of this matrix is given in A132710. See A111596 for an interpretation in terms of circular binary words and generalized factorials. - Tom Copeland, Nov 22 2007
Three combinatorial interpretations: T(n,k) is (1) the number of ways to split [n] = {1,...,n} into a collection of k nonempty lists ("partitions into sets of lists"), (2) the number of ways to split [n] into an ordered collection of n+1-k nonempty sets that are noncrossing ("partitions into lists of noncrossing sets"), (3) the number of Dyck n-paths with n+1-k peaks labeled 1,2,...,n+1-k in some order. - David Callan, Jul 25 2008
Given matrices A and B with A(n,k) = T(n,k)*a(n-k) and B(n,k) = T(n,k)*b(n-k), then A*B = D where D(n,k) = T(n,k)*[a(.)+b(.)]^(n-k), umbrally. - Tom Copeland, Aug 21 2008
An e.g.f. for the row polynomials of A(n,k) = T(n,k)*a(n-k) is exp[a(.)* D_x * x^2] exp(x*t) = exp(x*t) exp[(.)!*Lag(.,-x*t,1)*a(.)*x], umbrally, where [(.)! Lag(.,x,1)]^n = n! Lag(n,x,1) is a normalized Laguerre polynomial of order 1. - Tom Copeland, Aug 29 2008
Triangle of coefficients from the Bell polynomial of the second kind for f = 1/(1-x). B(n,k){x1,x2,x3,...} = B(n,k){1/(1-x)^2,...,(j-1)!/(1-x)^j,...} = T(n,k)/(1-x)^(n+k). - Vladimir Kruchinin, Mar 04 2011
The triangle, with the row and column offset taken as 0, is the generalized Riordan array (exp(x), x) with respect to the sequence n!*(n+1)! as defined by Wang and Wang (the generalized Riordan array (exp(x), x) with respect to the sequence n! is Pascal's triangle A007318, and with respect to the sequence n!^2 is A021009 unsigned). - Peter Bala, Aug 15 2013
For a relation to loop integrals in QCD, see p. 33 of Gopakumar and Gross and Blaizot and Nowak. - Tom Copeland, Jan 18 2016
Also the Bell transform of (n+1)!. For the definition of the Bell transform see A264428. - Peter Luschny, Jan 27 2016
Also the number of k-dimensional flats of the n-dimensional Shi arrangement. - Shuhei Tsujie, Apr 26 2019
The numbers T(n,k) appear as coefficients when expanding the rising factorials (x)^k = x(x+1)...(x+k-1) in the basis of falling factorials (x)k = x(x-1)...(x-k+1). Specifically, (x)^n = Sum{k=1..n} T(n,k) (x)k. - _Jeremy L. Martin, Apr 21 2021

			T(1,1) = C(1,1)*0!/0! = 1,
T(2,1) = C(2,1)*1!/0! = 2,
T(2,2) = C(2,2)*1!/1! = 1,
T(3,1) = C(3,1)*2!/0! = 6,
T(3,2) = C(3,2)*2!/1! = 6,
T(3,3) = C(3,3)*2!/2! = 1,
Sheffer a-sequence recurrence: T(6,2)= 1800 = (6/3)*120 + 6*240.
B(n,k) =
   1/(1-x)^2;
   2/(1-x)^3,  1/(1-x)^4;
   6/(1-x)^4,  6/(1-x)^5,  1/(1-x)^6;
  24/(1-x)^5, 36/(1-x)^6, 12/(1-x)^7, 1/(1-x)^8;
The triangle T(n,k) begins:
  n\k      1       2       3      4      5     6    7  8  9 ...
  1:       1
  2:       2       1
  3:       6       6       1
  4:      24      36      12      1
  5:     120     240     120     20      1
  6:     720    1800    1200    300     30     1
  7:    5040   15120   12600   4200    630    42    1
  8:   40320  141120  141120  58800  11760  1176   56  1
  9:  362880 1451520 1693440 846720 211680 28224 2016 72  1
  ...
Row n=10: [3628800, 16329600, 21772800, 12700800, 3810240, 635040, 60480, 3240, 90, 1]. - _Wolfdieter Lang_, Feb 01 2013
From _Peter Bala_, Feb 24 2025: (Start)
The array factorizes as an infinite product (read from right to left):
  /  1                \        /1             \^m /1           \^m /1           \^m
  |  2    1            |      | 0   1          |  |0  1         |  |1  1         |
  |  6    6   1        | = ...| 0   0   1      |  |0  1  1      |  |0  2  1      |
  | 24   36  12   1    |      | 0   0   1  1   |  |0  0  2  1   |  |0  0  3  1   |
  |120  240 120  20   1|      | 0   0   0  2  1|  |0  0  0  3  1|  |0  0  0  4  1|
  |...                 |      |...             |  |...          |  |...          |
where m = 2. Cf. A008277 (m = 1), A035342 (m = 3), A035469 (m = 4), A049029 (m = 5) A049385 (m = 6), A092082 (m = 7), A132056 (m = 8), A223511 - A223522 (m = 9 through 20), A001497 (m = -1), A004747 (m = -2), A000369 (m = -3), A011801 (m = -4), A013988 (m = -5). (End)

Reinhard Zumkeller, Rows n = 1..100 of triangle, flattened
Peter Bala, Factorising (r,b)-Stirling arrays
J. Fernando Barbero G., Jesús Salas, and Eduardo J. S. Villaseñor, Bivariate Generating Functions for a Class of Linear Recurrences. I. General Structure, arXiv:1307.2010 [math.CO], 2013.
Paul Barry, Eulerian polynomials as moments, via exponential Riordan arrays, arXiv preprint arXiv:1105.3043 [math.CO], 2011, J. Int. Seq. 14 (2011) # 11.9.5
Jean-Paul Blaizot and Maciej A. Nowak, Large N_c confinement and turbulence, arXiv:0801.1859 [hep-th], 2008.
David Callan, Sets, Lists and Noncrossing Partitions, arXiv:0711.4841 [math.CO], 2007-2008.
Pietro Codara, Ottavio M. D'Antona, and Pavol Hell, A simple combinatorial interpretation of certain generalized Bell and Stirling numbers, arXiv preprint arXiv:1308.1700 [cs.DM], 2013.
Tom Copeland, Mathemagical Forests, Addendum to Mathemagical Forests, The Inverse Mellin Transform, Bell Polynomials, a Generalized Dobinski Relation, and the Confluent Hypergeometric Functions, A Class of Differential Operators and the Stirling Numbers
Siad Daboul, Jan Mangaldan, Michael Z. Spivey and Peter Taylor, The Lah Numbers and the n-th Derivative of exp(1/x), Math. Mag., 86 (2013), 39-47.
Bérénice Delcroix-Oger and Clément Dupont, Lie-operads and operadic modules from poset cohomology, arXiv:2505.06094 [math.CO], 2025. See p. 33.
G. H. E. Duchamp et al., Feynman graphs and related Hopf algebras, J. Phys. (Conf Ser) 30 (2006) 107-118.
Rajesh Gopakumar and David J. Gross, Mastering the master field, arXiv:hep-th/9411021, 1994.
Gábor Hetyei, Meixner polynomials of the second kind and quantum algebras representing su(1,1), arXiv preprint arXiv:0909.4352 [math.QA], 2009, p. 4. - From _Tom Copeland_, Oct 01 2015
Milan Janjic, Some classes of numbers and derivatives, JIS 12 (2009) 09.8.3.
Donald E. Knuth, Convolution polynomials, The Mathematica J., 2 (1992), 67-78.
Shi-Mei Ma, Some combinatorial sequences associated with context-free grammars, arXiv:1208.3104v2 [math.CO], 2012. - From _N. J. A. Sloane_, Aug 21 2012
MacTutor History of Mathematics archive: Biography of Ivo Lah.
Robert S. Maier, Boson Operator Ordering Identities from Generalized Stirling and Eulerian Numbers, arXiv:2308.10332 [math.CO], 2023. See p. 19.
Norihiro Nakashima and Shuhei Tsujie, Enumeration of Flats of the Extended Catalan and Shi Arrangements with Species, arXiv:1904.09748 [math.CO], 2019. See p. 18.
Michael Penn, Lah Numbers and an appearance of exponential generating functions, YouTube video, 2025.
Mathias Pétréolle and Alan D. Sokal, Lattice paths and branched continued fractions. II. Multivariate Lah polynomials and Lah symmetric functions, arXiv:1907.02645 [math.CO], 2019. See p. 4.
Tilman Piesk, Illustration of the first four rows
Kornelia Ufniarz and Grzegorz Siudem, Combinatorial origins of the canonical ensemble, arXiv:2008.00244 [math-ph], 2020. See p. 5.
Weiping Wang and Tianming Wang, Generalized Riordan arrays, Discrete Mathematics, Vol. 308, No. 24, 6466-6500.
Wikipedia, Lah number

Cf. A000262, A103194, A105220.
Triangle of Lah numbers (A008297) unsigned.
Cf. A111596 (signed triangle with extra n=0 row and m=0 column).
Cf. A130561 (for a natural refinement).
Cf. A094638 (for differential operator representation).
Cf. A248045 (central terms), A002868 (row maxima).
Cf, A059110.
Cf. A001263, A008277.
Cf. A089231 (triangle with mirrored rows).
Cf. A271703 (triangle with extra n=0 row and m=0 column).

Flat(List([1..10],n->List([1..n],k->Binomial(n,k)*Factorial(n-1)/Factorial(k-1)))); # Muniru A Asiru, Jul 25 2018

a105278 n k = a105278_tabl !! (n-1) !! (k-1)
a105278_row n = a105278_tabl !! (n-1)
a105278_tabl = [1] : f [1] 2 where
   f xs i = ys : f ys (i + 1) where
     ys = zipWith (+) ([0] ++ xs) (zipWith (*) [i, i + 1 ..] (xs ++ [0]))
-- Reinhard Zumkeller, Sep 30 2014, Mar 18 2013

/* As triangle */ [[Binomial(n,k)*Factorial(n-1)/Factorial(k-1): k in [1..n]]: n in [1.. 15]]; // Vincenzo Librandi, Oct 31 2014

The triangle: for n from 1 to 13 do seq(binomial(n,k)*(n-1)!/(k-1)!,k=1..n) od;
the sequence: seq(seq(binomial(n,k)*(n-1)!/(k-1)!,k=1..n),n=1..13);
# The function BellMatrix is defined in A264428.
# Adds (1, 0, 0, 0, ...) as column 0.
BellMatrix(n -> (n+1)!, 9); # Peter Luschny, Jan 27 2016

nn = 9; a = x/(1 - x); f[list_] := Select[list, # > 0 &]; Flatten[Map[f, Drop[Range[0, nn]! CoefficientList[Series[Exp[y a], {x, 0, nn}], {x, y}], 1]]] (* Geoffrey Critzer, Dec 11 2011 *)
nn = 9; Flatten[Table[(j - k)! Binomial[j, k] Binomial[j - 1, k - 1], {j, nn}, {k, j}]] (* Jan Mangaldan, Mar 15 2013 *)
rows = 10;
t = Range[rows]!;
T[n_, k_] := BellY[n, k, t];
Table[T[n, k], {n, 1, rows}, {k, 1, n}] // Flatten (* Jean-François Alcover, Jun 23 2018, after Peter Luschny *)
T[n_, n_] := 1; T[n_, k_] /;0Oliver Seipel, Dec 06 2024 *)

use ntheory ":all"; say join ", ", map { my $n=$; map { stirling($n,$,3) } 1..$n; } 1..9; # Dana Jacobsen, Mar 16 2017

T(n,k) = Sum_{m=n..k} |S1(n,m)|*S2(m,k), k>=n>=1, with Stirling triangles S2(n,m):=A048993 and S1(n,m):=A048994.
T(n,k) = C(n,k)*(n-1)!/(k-1)!.
Sum_{k=1..n} T(n,k) = A000262(n).
n*Sum_{k=1..n} T(n,k) = A103194(n) = Sum_{k=1..n} T(n,k)*k^2.
E.g.f. column k: (x^(k-1)/(1-x)^(k+1))/(k-1)!, k>=1.
Recurrence from Sheffer (here Jabotinsky) a-sequence [1,1,0,...] (see the W. Lang link under A006232): T(n,k)=(n/k)*T(n-1,m-1) + n*T(n-1,m). - Wolfdieter Lang, Jun 29 2007
The e.g.f. is, umbrally, exp[(.)!* L(.,-t,1)*x] = exp[t*x/(1-x)]/(1-x)^2 where L(n,t,1) = Sum_{k=0..n} T(n+1,k+1)*(-t)^k = Sum_{k=0..n} binomial(n+1,k+1)* (-t)^k / k! is the associated Laguerre polynomial of order 1. - Tom Copeland, Nov 17 2007
For this Lah triangle, the n-th row polynomial is given umbrally by
  n! C(B.(x)+1+n,n) = (-1)^n C(-B.(x)-2,n), where C(x,n)=x!/(n!(x-n)!),
  the binomial coefficient, and B_n(x)= exp(-x)(xd/dx)^n exp(x), the n-th Bell / Touchard / exponential polynomial (cf. A008277). E.g.,
  2! C(-B.(-x)-2,2) = (-B.(x)-2)(-B.(x)-3) = B_2(x) + 5*B_1(x) + 6 = 6 + 6x + x^2.
  n! C(B.(x)+1+n,n) = n! e^(-x) Sum_{j>=0} C(j+1+n,n)x^j/j! is a corresponding Dobinski relation. See the Copeland link for the relation to inverse Mellin transform. - Tom Copeland, Nov 21 2011
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 (D = (1+x)*d/dx), A035342 (D = (1+x)^3*d/dx), A035469 (D = (1+x)^4*d/dx) and A049029 (D = (1+x)^5*d/dx). - Peter Bala, Nov 25 2011
T(n,k) = Sum_{i=k..n} A130534(n-1,i-1)*A008277(i,k). - Reinhard Zumkeller, Mar 18 2013
Let E(x) = Sum_{n >= 0} x^n/(n!*(n+1)!). Then a generating function is exp(t)*E(x*t) = 1 + (2 + x)*t + (6 + 6*x + x^2)*t^2/(2!*3!) + (24 + 36*x + 12*x^2 + x^3)*t^3/(3!*4!) + ... . - Peter Bala, Aug 15 2013
P_n(x) = L_n(1+x) = n!*Lag_n(-(1+x);1), where P_n(x) are the row polynomials of A059110; L_n(x), the Lah polynomials of A105278; and Lag_n(x;1), the Laguerre polynomials of order 1. These relations follow from the relation between the iterated operator (x^2 D)^n and ((1+x)^2 D)^n with D = d/dx. - Tom Copeland, Jul 23 2018
Dividing each n-th diagonal by n!, where the main diagonal is n=1, generates the Narayana matrix A001263. - Tom Copeland, Sep 23 2020
T(n,k) = A089231(n,n-k). - Ron L.J. van den Burg, Dec 12 2021
T(n,k) = T(n-1,k-1) + (n+k-1)*T(n-1,k). - Bérénice Delcroix-Oger, Jun 25 2025

Stirling comments and e.g.f.s from Wolfdieter Lang, Apr 11 2007

A000369 Triangle of numbers related to triangle A049213; generalization of Stirling numbers of second kind A008277, Bessel triangle A001497.

Original entry on oeis.org

1, 3, 1, 21, 9, 1, 231, 111, 18, 1, 3465, 1785, 345, 30, 1, 65835, 35595, 7650, 825, 45, 1, 1514205, 848925, 196245, 24150, 1680, 63, 1, 40883535, 23586255, 5755050, 775845, 62790, 3066, 84, 1, 1267389585, 748471185, 190482705, 27478710

Offset: 1

Wolfdieter Lang

a(n,m) := S2p(-3; n,m), a member of a sequence of triangles including S2p(-1; n,m) := A001497(n-1,m-1) (Bessel triangle) and ((-1)^(n-m))*S2p(1; n,m) := A008277(n,m) (Stirling 2nd kind). a(n,1)= A008545(n-1).
a(n,m), n>=m>=1, enumerates unordered n-vertex m-forests composed of m increasing plane (aka ordered) trees, with one vertex of out-degree r=0 (leafs or a root) and each vertex with out-degree r>=1 comes in r+2 types (like for an (r+2)-ary vertex). Proof from the e.g.f. of the first column Y(z):=1-(1-4*x)^(1/4) and the F. Bergeron et al. reference given in A001498, eq. (8), Y'(z)= phi(Y(z)), Y(0)=0, with out-degree o.g.f. phi(w)=1/(1-w)^3. - Wolfdieter Lang, Oct 12 2007
Also the Bell transform of the quadruple factorial numbers Product_{k=0..n-1} (4*k+3) (A008545) adding 1,0,0,0,... as column 0. For the definition of the Bell transform see A264428 and for cross-references A265606. - Peter Luschny, Dec 31 2015

			Triangle begins:
  1;
  3, 1;
  21, 9, 1;
  231, 111, 18, 1;
  3465, 1785, 345, 30, 1; ...
Tree combinatorics for a(3,2)=9: there are three m=2 forests each with one tree a root (with out-degree r=0) and the other tree a root and a leaf coming in three versions (like for a 3-ary vertex). Each such forest can be labeled increasingly in three ways (like (1,(23)), (2,(13)) and (3,(12))) yielding 9 such forests. - _Wolfdieter Lang_, Oct 12 2007

Vincenzo Librandi, Rows n = 1..50, flattened
P. Blasiak, K. A. Penson and A. I. Solomon, The general boson normal ordering problem, arXiv:quant-ph/0212072, 2002.
Tom Copeland, A Class of Differential Operators and the Stirling Numbers
M. Janjic, Some classes of numbers and derivatives, JIS 12 (2009) 09.8.3
Wolfdieter Lang, First ten rows.
Wolfdieter Lang, On generalizations of Stirling number triangles, J. Integer Seqs., Vol. 3 (2000), #00.2.4.
Toufik Mansour, Matthias Schork and Mark Shattuck, The Generalized Stirling and Bell Numbers Revisited, Journal of Integer Sequences, Vol. 15 (2012), #12.8.3.
Mathias Pétréolle, Alan D. Sokal, Lattice paths and branched continued fractions. II. Multivariate Lah polynomials and Lah symmetric functions, arXiv:1907.02645 [math.CO], 2019.
Index entries for sequences related to Bessel functions or polynomials

Row sums give A016036. Cf. A004747.
Columns include A008545.
Alternating row sums A132163.

a[n_, m_] /; n >= m >= 1 := a[n, m] = (4(n-1) - m)*a[n-1, m] + a[n-1, m-1]; a[n_, m_] /; n < m = 0; a[, 0] = 0; a[1, 1] = 1; Flatten[Table[a[n, m], {n, 1, 9}, {m, 1, n}]] (* _Jean-François Alcover, Jul 22 2011 *)

# uses[bell_transform from A264428]
# Adds a column 1,0,0,0,... at the left side of the triangle.
def A000369_row(n):
    multifact_4_3 = lambda n: prod(4*k + 3 for k in (0..n-1))
    mfact = [multifact_4_3(k) for k in (0..n)]
    return bell_transform(n, mfact)
[A000369_row(n) for n in (0..9)] # Peter Luschny, Dec 31 2015

a(n, m) = n!*A049213(n, m)/(m!*4^(n-m)); a(n+1, m) = (4*n-m)*a(n, m) + a(n, m-1), n >= m >= 1; a(n, m) := 0, n

E.g.f. of m-th column: ((1-(1-4*x)^(1/4))^m)/m!.

From Peter Bala, Jun 08 2016: (Start)

With offset 0, the e.g.f. is 1/(1 - 4*x)^(3/4)*exp(t*(1 - (1 - 4*x)^(1/4))) = 1 + (3 + t)*x + (21 + 9*t + t^2)*x^2/2! + ....

Thus with row and column numbering starting at 0, this triangle is the exponential Riordan array [d/dx(F(x)), F(x)], belonging to the Derivative subgroup of the exponential Riordan group, where F(x) = 1 - (1 - 4*x)^(1/4).

Row polynomial recurrence: R(n+1,t) = t*Sum_{k = 0..n} binomial(n,k)*A008545(k)*R(n-k,t) with R(0,t) = 1. (End)

A011801 Triangle read by rows, the inverse Bell transform of n!*binomial(4,n) (without column 0).

Original entry on oeis.org

1, 4, 1, 36, 12, 1, 504, 192, 24, 1, 9576, 3960, 600, 40, 1, 229824, 100656, 17160, 1440, 60, 1, 6664896, 3048192, 563976, 54600, 2940, 84, 1, 226606464, 107255232, 21095424, 2256576, 142800, 5376, 112, 1, 8837652096, 4302305280, 887785920, 102332160, 7254576, 325584, 9072, 144, 1

Offset: 1

Wolfdieter Lang

Previous name was: Triangle of numbers related to triangle A049223; generalization of Stirling numbers of second kind A008277, Bessel triangle A001497.
T(n, m) = S2p(-4; n,m), a member of a sequence of triangles including S2p(-1; n,m) = A001497(n-1,m-1) (Bessel triangle) and ((-1)^(n-m))*S2p(1; n,m) = A008277(n,m) (Stirling 2nd kind). T(n, 1) = A008546(n-1).
For the definition of the Bell transform see A264428 and the link. - Peter Luschny, Jan 16 2016

			Triangle starts:
          1;
          4,         1;
         36,        12,        1;
        504,       192,       24,       1;
       9576,      3960,      600,      40,      1;
     229824,    100656,    17160,    1440,     60,     1;
    6664896,   3048192,   563976,   54600,   2940,    84,    1;
  226606464, 107255232, 21095424, 2256576, 142800,  5376,  112,   1;

G. C. Greubel, Rows n = 1..50 of the triangle, flattened
P. Blasiak, K. A. Penson and A. I. Solomon, The general boson normal ordering problem, arXiv:quant-ph/0402027, 2004.
Wolfdieter Lang, On generalizations of Stirling number triangles, J. Integer Seqs., Vol. 3 (2000), #00.2.4.
Peter Luschny, The Bell transform
Index entries for sequences related to Bessel functions or polynomials

Cf. A008277, A008546, A049223, A264428.
Cf. A028575 (row sums).
Triangles with the recurrence T(n,k) = (m*(n-1)-k)*T(n-1,k) + T(n-1,k-1): A010054 (m=1), A001497 (m=2), A004747 (m=3), A000369 (m=4), this sequence (m=5), A013988 (m=6).

function T(n,k) // T = A011801
  if k eq 0 then return 0;
  elif k eq n then return 1;
  else return (5*(n-1)-k)*T(n-1,k) + T(n-1,k-1);
  end if;
end function;
[T(n,k): k in [1..n], n in [1..12]]; // G. C. Greubel, Oct 03 2023

(* First program *)
T[n_, m_] /; n>=m>=1:= T[n, m]= (5*(n-1)-m)*T[n-1, m] + T[n-1, m-1]; T[n_, m_] /; nJean-François Alcover, Jun 20 2018 *)
(* Second program *)
rows = 10;
b[n_, m_]:= BellY[n, m, Table[k! Binomial[4, k], {k, 0, rows}]];
T= Table[b[n, m], {n,rows}, {m,rows}]//Inverse//Abs;
A011801= Table[T[[n, m]], {n,rows}, {m,n}]//Flatten (* Jean-François Alcover, Jun 22 2018 *)

# uses[inverse_bell_matrix from A264428]
# Adds 1,0,0,0, ... as column 0 at the left side of the triangle.
inverse_bell_matrix(lambda n: factorial(n)*binomial(4, n), 8) # Peter Luschny, Jan 16 2016

T(n, m) = n!*A049223(n, m)/(m!*5^(n-m)).
T(n+1, m) = (5*n-m)*T(n, m) + T(n, m-1), for n >= m >= 1, with T(n, m) = 0, for n < m, and T(n, 0) = 0, T(1, 1) = 1.
E.g.f. of n-th column: (1/n!)*( 1 - (1-5*x)^(1/5) )^n.
Sum_{k=1..n} T(n, k) = A028575(n).

New name from Peter Luschny, Jan 16 2016

A051141 Triangle read by rows: a(n, m) = S1(n, m)*3^(n-m), where S1 are the signed Stirling numbers of first kind A008275 (n >= 1, 1 <= m <= n).

Original entry on oeis.org

1, -3, 1, 18, -9, 1, -162, 99, -18, 1, 1944, -1350, 315, -30, 1, -29160, 22194, -6075, 765, -45, 1, 524880, -428652, 131544, -19845, 1575, -63, 1, -11022480, 9526572, -3191076, 548289, -52920, 2898, -84, 1, 264539520, -239660208

Offset: 1

Wolfdieter Lang

Previous name was: Generalized Stirling number triangle of first kind.
a(n,m) = R_n^m(a=0,b=3) in the notation of the given reference.
a(n,m) is a Jabotinsky matrix, i.e., the monic row polynomials E(n,x) := Sum_{m=1..n} a(n,m)*x^m = Product_{j=0..n-1} (x - 3*j), n >= 1 and E(0,x) := 1 are exponential convolution polynomials (see A039692 for the definition and a Knuth reference).
This is the signed Stirling1 triangle with diagonals d>=0 (main diagonal d=0) scaled with 3^d.
Exponential Riordan array [1/(1 + 3*x), log(1 + 3*x)/3]. The unsigned triangle is [1/(1 - 3*x), log(1/(1 - 3*x)^(1/3))]. - Paul Barry, Apr 29 2009
Also the Bell transform of the triple factorial numbers A032031 which adds a first column (1, 0, 0 ...) on the left side of the triangle and computes the unsigned values. For the definition of the Bell transform, see A264428. See A004747 for the triple factorial numbers A008544 and A203412 for the triple factorial numbers A007559 as well as A039683 and A132062 for the case of double factorial numbers. - Peter Luschny, Dec 21 2015

			Triangle starts:
       1;
      -3,       1;
      18,      -9,      1;
    -162,      99,    -18,      1;
    1944,   -1350,    315,    -30,    1;
  -29160,   22194,  -6075,    765,  -45,   1;
  524880, -428652, 131544, -19845, 1575, -63, 1;
---
Row polynomial E(3,x) = 18*x-9*x^2+x^3.
From _Paul Barry_, Apr 29 2009: (Start)
The unsigned array [1/(1 - 3*x), log(1/(1 - 3*x)^(1/3))] has production matrix
    3,    1;
    9,    6,    1;
   27,   27,    9,   1;
   81,  108,   54,  12,   1;
  243,  405,  270,  90,  15,  1;
  729, 1458, 1215, 540, 135, 18, 1;
  ...
which is A007318^{3} beheaded (by viewing A007318 as a lower triangular matrix). See the comment above. (End)

G. C. Greubel, Table of n, a(n) for the first 50 rows, flattened
Richell O. Celeste, Roberto B. Corcino, and Ken Joffaniel M. Gonzales, Two Approaches to Normal Order Coefficients, Journal of Integer Sequences, Vol. 20 (2017), Article 17.3.5.
Wolfdieter Lang, First 10 rows.
D. S. Mitrinovic and M. S. Mitrinovic, Tableaux d'une classe de nombres reliés aux nombres de Stirling, Univ. Beograd. Publ. Elektrotehn. Fak. Ser. Mat. Fiz. 77 (1962), 1-77.

First (m=1) column sequence is: A032031(n-1).
Row sums (signed triangle): A008544(n-1)*(-1)^(n-1).
Row sums (unsigned triangle): A007559(n).
Cf. A008275 (Stirling1 triangle, b=1), A039683 (b=2), A051142 (b=4).
Cf. A039683, A132062, A264428.

a[n_, m_] /; n >= m >= 1 := a[n, m] = a[n-1, m-1] - 3(n-1)*a[n-1, m]; a[n_, m_] /; n < m = 0; a[, 0] = 0; a[1, 1] = 1; Flatten[Table[a[n, m], {n, 1, 9}, {m, 1, n}]][[1 ;; 38]] (* _Jean-François Alcover, Jun 01 2011, after formula *)
Table[StirlingS1[n, m]*3^(n - m), {n, 1, 10}, {m, 1, n}]//Flatten (* G. C. Greubel, Oct 24 2017 *)

for(n=1,10, for(m=1,n, print1(stirling(n,m,1)*3^(n-m), ", "))) \\ G. C. Greubel, Oct 24 2017

# uses[bell_transform from A264428]
triplefactorial = lambda n: 3^n*factorial(n)
def A051141_row(n):
    trifact = [triplefactorial(k) for k in (0..n)]
    return bell_transform(n, trifact)
[A051141_row(n) for n in (0..8)] # Peter Luschny, Dec 21 2015

a(n, m) = a(n-1, m-1) - 3*(n-1)*a(n-1, m) for n >= m >= 1; a(n, m) = 0 for n < m; a(n, 0) = 0 for n >= 1; a(0, 0) = 1.
E.g.f. for the m-th column of the signed triangle: (log(1 + 3*x)/3)^m/m!.
|a(n,1)| = A032031(n-1). - Peter Luschny, Dec 23 2015

Name clarified using a formula of the author by Peter Luschny, Dec 23 2015

Showing 1-10 of 16 results. Next

cp's OEIS Frontend

A015735 Row sums of triangle A004747.

Views

Author

Keywords

Links

Crossrefs

Programs

Magma

Mathematica

Maxima

SageMath

Formula

A143172 Partition number array, called M32(-2), related to A004747(n,m) = |S2(-2;n,m)| (generalized Stirling triangle).

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Formula

A144274 Partition number array, called M32hat(-2)= 'M32(-2)/M3'= 'A143172/A036040', related to A004747(n,m)= |S2(-2;n,m)| (generalized Stirling triangle).

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Formula

A144345 Second column (m=2) of triangle S2p(-2) = A004747.

Views

Author

Keywords

Crossrefs

Formula

A144346 Third column (m=3) of triangle S2p(-2) = A004747.

Views

Author

Keywords

Crossrefs

Formula

A008544 Triple factorial numbers: Product_{k=0..n-1} (3*k+2).

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Haskell

Magma

Magma

Maple

Mathematica

Maxima

PARI

PARI

Sage

Sage

Formula

A105278 Triangle read by rows: T(n,k) = binomial(n,k)*(n-1)!/(k-1)!.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

GAP

Haskell

Magma

Maple

Mathematica

Perl

Formula

Extensions