A008277 -id:A008277 - cp's OEIS frontend

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

0, 1, 10, 45, 136, 325, 666, 1225, 2080, 3321, 5050, 7381, 10440, 14365, 19306, 25425, 32896, 41905, 52650, 65341, 80200, 97461, 117370, 140185, 166176, 195625, 228826, 266085, 307720, 354061, 405450, 462241, 524800, 593505, 668746, 750925, 840456, 937765

Offset: 0

Aaron Gulliver (gulliver(AT)elec.canterbury.ac.nz)

Sum of first n^2 positive integers.
Start from xanthene and attach amino acids according to the reaction scheme that describes the reaction between the active sites. See the hyperlink below on chemistry. - Robert G. Wilson v, Aug 02 2002; Amarnath Murthy, Aug 01 2002
Sum of the next n multiples of n. - Amarnath Murthy, Aug 01 2002
The sum of the terms in an n X n spiral. These are also triangular numbers. - William A. Tedeschi, Feb 27 2008
Hypotenuse of Pythagorean triangles with smallest side a cube: A000578(n)^2 + A083374(n)^2 = a(n)^2. - Martin Renner, Nov 12 2011
For n>1, triangular numbers that can be represented as a sum of a square and a triangular number. For example, a(2)=10=4+6=9+1. - Ivan N. Ianakiev, Apr 24 2012
A037270 can be constructed in the following manner: Take A000217 and for every n not in A000290 delete the corresponding A000217(n). - Ivan N. Ianakiev, Apr 26 2012
Starting at a(1)=1 simply take 1*1=1, a(2)= 2*(2+3)=10, a(3)= 3*(4+5+6)=45, a(4)=4*(7+8+9+10) and so on. - J. M. Bergot, May 01 2015
Observation: The digital roots of the terms repeat in the sequence 1, 1, 9; e.g., the digital roots of 1, 10, 45, 136, 325, and 666 are 1, 1, 9, 1, 1, and 9. Verified for the first 10000 terms. - Rob Barton, Mar 28 2018
The above observation is easily explained and proved given that the digital root of a positive number equals the number modulo 9, and a(n + 9k) == a(n) (mod 9). - M. F. Hasler, Apr 05 2018
Number of unoriented rows of length 4 using up to n colors. For a(0)=0, there are no rows using no colors. For a(1)=1, there is one row using that one color for all positions. For a(2)=10, there are 4 achiral (AAAA, ABBA, BAAB, BBBB) and 6 chiral pairs (AAAB-BAAA, AABA-ABAA, AABB-BBAA, ABAB-BABA, ABBB-BBBA, BABB-BBAB). - Robert A. Russell, Nov 14 2018
For n > 0, a(2n+1) is the number of non-isomorphic 6C_m-snakes, where m = 2n+1 or m = 2n (for n>=2). A kC_n-snake is a connected graph in which the k>=2 blocks are isomorphic to the cycle C_n and the block-cutpoint graph is a path. - Christian Barrientos, May 15 2019
Number of achiral colorings of the edges of a tetrahedron with n available colors. - Robert A. Russell, Sep 07 2019

C. Alsina and R. B. Nelson, Charming Proofs: A Journey into Elegant Mathematics, MAA, 2010. See p. 5.
C. Barrientos, Graceful labelings of cyclic snakes, Ars Combin., 60(2001), 85-96.
Albert H. Beiler, Recreations in the theory of numbers, New York: Dover, (2nd ed.) 1966, p. 106, table 55.
T. A. Gulliver, Sequences from Arrays of Integers, Int. Math. Journal, Vol. 1, No. 4, pp. 323-332, 2002.
T. A. Gulliver, Sequences from Cubes of Integers, Int. Math. Journal, 4 (2003), 439-445.
R. A. Wilson, Cosmic Trigger, epilogue of S.-P. Sirag.

T. D. Noe, Table of n, a(n) for n = 0..1000
J. D. Bell, A translation of Leonhard Euler's "De Quadratis Magicis", E795, arXiv:math/0408230 [math.CO], 2004-2005.
N. G de Bruijn, Some classes of integer-valued functions, Nederl. Akad. Wetensch. Proc, Ser. A, 58 (1955), 363-367. See page 363.
Th. Gruner, A. Kerber, R. Laue, and M. Meringer, Mathematics for Combinatorial Chemistry, In: F. Keil, W. Mackens, H. Voß and J. Wether, Scientific Computing in Chemical Engineering II, Springer, 1999, 74-81.
Index entries for linear recurrences with constant coefficients, signature (5,-10,10,-5,1).

Cf. A000217, A236770 (see crossrefs).
Row 4 of A277504.
Cf. A000583 (oriented), A083374 (chiral), A000290 (achiral).
Cf. A317617.
Row 3 of A327086 (achiral simplex edge colorings).

a:=List([0..30],n->n^2*(n^2+1)/2); # Muniru A Asiru, Mar 28 2018

[n^2*(n^2 + 1)/2: n in [0..30]] // Stefano Spezia, Jan 15 2019

seq(n^2*(n^2+1)/2,n=0..30); # Muniru A Asiru, Mar 28 2018

Table[ n^2*((n^2 + 1)/2), {n, 0, 30} ]
Table[(1/8) Round[N[Sinh[2 ArcSinh[n]]^2, 100]], {n, 0, 30}] (* Artur Jasinski, Feb 10 2010 *)
LinearRecurrence[{5,-10,10,-5,1},{0,1,10,45,136},30] (* Harvey P. Dale, Aug 03 2014 *)

a(n)=binomial(n^2+1,2) \\ Charles R Greathouse IV, Apr 25 2012

for n in range(0,30): print(n**2*(n**2+1)/2, end=', ') # Stefano Spezia, Jan 10 2019

a(n) = a(n-1) + n^3 + (n-1)^3.
a(n) = A000537(n)+A000537(n-1), i.e., square of sum of first n integers plus square of sum of first n-1 integers. - Henry Bottomley, Oct 15 2001
a(n) = Sum_{k=0..n^2} k. - William A. Tedeschi, Feb 27 2008
a(n) = (1/8)*sinh(2*arcsinh(n)). - Artur Jasinski, Feb 10 2010
G.f.: x*(1+x)*(1+4*x+x^2)/(1-x)^5. - Colin Barker, Mar 22 2012
a(n) = a(n-1) + A005898(n-1). - Ivan N. Ianakiev, May 13 2012
a(n) = 2 * A000217(n-1) * A000217(n) + A000290(n). - Ivan N. Ianakiev, May 26 2012
a(n) = A000217(n^2). - J. M. Bergot, Jun 07 2012
a(n) = 5*a(n-1) -10*a(n-2) +10*a(n-3) -5*a(n-4) +a(n-5) n>4, a(0)=0, a(1)=1, a(2)=10, a(3)=45, a(4)=136. - Yosu Yurramendi, Sep 02 2013
For n>0, a(n) = A000217(n)^2 + A000217(n-1)^2. - Richard R. Forberg, Dec 25 2013
a(n) = T(T(n)) + T(T(n-1)) + T(T(n)-1) + T(T(n-1)-1), where T(n) = A000217(n). - Charlie Marion, Sep 10 2016
a(n) = t(n-3)*t(n)+t(n-1)*t(n+2), with t(n)=A000217(n). - J. M. Bergot, Apr 07 2018
From Robert A. Russell, Nov 14 2018: (Start)
a(n) = (A000583(n) + A000290(n)) / 2 = (n^4 + n^2) / 2.
a(n) = A000583(n) - A083374(n) = A083374(n) + A000290(n).
G.f.: (Sum_{j=1..4} S2(4,j)*j!*x^j/(1-x)^(j+1) + Sum_{j=1..2} S2(2,j)*j!*x^j/(1-x)^(j+1)) / 2, where S2 is the Stirling subset number A008277.
G.f.: Sum_{k=1..4} A145882(4,k) * x^k / (1-x)^5.
E.g.f.: (Sum_{k=1..4} S2(4,k)*x^k + Sum_{k=1..2} S2(2,k)*x^k) * exp(x) / 2, where S2 is the Stirling subset number A008277.
For n>4, a(n) = Sum_{j=1..5} -binomial(j-6,j) * a(n-j). (End)
a(n) = n*A006003(n). - Kritsada Moomuang, Dec 16 2018
For n > 0, a(n) = Sum_{k=1..n} A317617(n,k). - Stefano Spezia, Jan 10 2019
Sum_{n>=1} 1/a(n) = 1 + Pi^2/3 - Pi*coth(Pi) = 1.13652003875929052467672874379... - Vaclav Kotesovec, Jan 21 2019
Sum_{n>=1} (-1)^(n+1)/a(n) = Pi*csch(Pi) + Pi^2/6 - 1. - Amiram Eldar, Nov 02 2021

A026898 a(n) = Sum_{k=0..n} (n-k+1)^k.

Original entry on oeis.org

1, 2, 4, 9, 23, 66, 210, 733, 2781, 11378, 49864, 232769, 1151915, 6018786, 33087206, 190780213, 1150653921, 7241710930, 47454745804, 323154696185, 2282779990495, 16700904488706, 126356632390298, 987303454928973, 7957133905608837, 66071772829247410

Offset: 0

N. J. A. Sloane

nonn

Row sums of A004248, A009998, A009999.
First differences are in A047970.
First differences of A103439.
Antidiagonal sums of array A003992.
a(n-1), for n>=1, is the number of length-n restricted growth strings (RGS) [s(0),s(1),...,s(n-1)] where s(0)=0 and s(k)<=1+max(prefix) for k>=1, that are simultaneously projections as maps f: [n] -> [n] where f(x)<=x and f(f(x))=f(x); see example and the two comments (Arndt, Apr 30 2011 Jan 04 2013) in A000110. - Joerg Arndt, Mar 07 2015
Number of finite sequences s of length n+1 whose discriminator sequence is s itself. Here the discriminator sequence of s is the one where the n-th term (n>=1) is the least positive integer k such that the first n terms are pairwise incongruent, modulo k. - Jeffrey Shallit, May 17 2016
From Gus Wiseman, Jan 08 2019: (Start)
Also the number of set partitions of {1,...,n+1} whose minima form an initial interval of positive integers. For example, the a(3) = 9 set partitions are:
  {{1},{2},{3},{4}}
  {{1},{2},{3,4}}
  {{1},{2,4},{3}}
  {{1,4},{2},{3}}
  {{1},{2,3,4}}
  {{1,3},{2,4}}
  {{1,4},{2,3}}
  {{1,3,4},{2}}
  {{1,2,3,4}}
Missing from this list are:
  {{1},{2,3},{4}}
  {{1,2},{3},{4}}
  {{1,3},{2},{4}}
  {{1,2},{3,4}}
  {{1,2,3},{4}}
  {{1,2,4},{3}}
(End)
a(n) is the number of m-tuples of nonnegative integers less than or equal to n-m (including the "0-tuple"). - Mathew Englander, Apr 11 2021

			G.f.: A(x) = 1 + 2*x + 4*x^2 + 9*x^3 + 23*x^4 + 66*x^5 + 210*x^6 + ...
where we have the identity:
A(x) = 1/(1-x) + x/(1-2*x) + x^2/(1-3*x) + x^3/(1-4*x) + x^4/(1-5*x) + ...
is equal to
A(x) = 1/(1-x) + x/((1-x)^2*(1+x)) + 2!*x^2/((1-x)^3*(1+x)*(1+2*x)) + 3!*x^3/((1-x)^4*(1+x)*(1+2*x)*(1+3*x)) + 4!*x^4/((1-x)^5*(1+x)*(1+2*x)*(1+3*x)*(1+4*x)) + ...
From _Joerg Arndt_, Mar 07 2015: (Start)
The a(5-1) = 23 RGS described in the comment are (dots denote zeros):
01:  [ . . . . . ]
02:  [ . 1 . . . ]
03:  [ . 1 . . 1 ]
04:  [ . 1 . 1 . ]
05:  [ . 1 . 1 1 ]
06:  [ . 1 1 . . ]
07:  [ . 1 1 . 1 ]
08:  [ . 1 1 1 . ]
09:  [ . 1 1 1 1 ]
10:  [ . 1 2 . . ]
11:  [ . 1 2 . 1 ]
12:  [ . 1 2 . 2 ]
13:  [ . 1 2 1 . ]
14:  [ . 1 2 1 1 ]
15:  [ . 1 2 1 2 ]
16:  [ . 1 2 2 . ]
17:  [ . 1 2 2 1 ]
18:  [ . 1 2 2 2 ]
19:  [ . 1 2 3 . ]
20:  [ . 1 2 3 1 ]
21:  [ . 1 2 3 2 ]
22:  [ . 1 2 3 3 ]
23:  [ . 1 2 3 4 ]
(End)

Reinhard Zumkeller, Table of n, a(n) for n = 0..500
Fufa Beyene, Jörgen Backelin, Roberto Mantaci, and Samuel A. Fufa, Set Partitions and Other Bell Number Enumerated Objects, J. Int. Seq., Vol. 26 (2023), Article 23.1.8.
Giulio Cerbai, Pattern-avoiding modified ascent sequences, arXiv:2401.10027 [math.CO], 2024. See p. 12.
Giulio Cerbai and Anders Claesson, Counting fixed-point-free Cayley permutations, arXiv:2507.09304 [math.CO], 2025. See p. 25.
Sajed Haque, Discriminators of Integer Sequences, 2017, See p. 33 Corollary 29.
Mathematics Stack Exchange, Asymptotics of ..., 2011.
Chunyan Yan and Zhicong Lin, Inversion sequences avoiding pairs of patterns, arXiv:1912.03674 [math.CO], 2019. See p. 3.
Robin D. P. Zhou, Pattern avoidance in revised ascent sequences, arXiv:2505.05171 [math.CO], 2025. See p. 4.

Cf. A000110, A000258, A000670, A003101, A008277, A038125, A062810.

Cf. A105795, A179928, A287215, A287216.

a026898 n = sum $ zipWith (^) [n + 1, n .. 1] [0 ..]
-- Reinhard Zumkeller, Sep 14 2014

[(&+[(n-k+1)^k: k in [0..n]]): n in [0..50]]; // Stefano Spezia, Jan 09 2019

a:= n-> add((n+1-j)^j, j=0..n): seq(a(n), n=0..23); # Zerinvary Lajos, Apr 18 2009

Table[Sum[(n-k+1)^k, {k,0,n}], {n, 0, 25}] (* Michael De Vlieger, Apr 01 2015 *)

{a(n)=polcoeff(sum(m=0,n,x^m/(1-(m+1)*x+x*O(x^n))),n)} /* Paul D. Hanna, Sep 13 2011 */

{INTEGRATE(n,F)=local(G=F);for(i=1,n,G=intformal(G));G}
{a(n)=local(A=1+x);A=sum(k=0,n,INTEGRATE(k,exp((k+1)*x+x*O(x^n))));n!*polcoeff(A,n)} \\ Paul D. Hanna, Dec 28 2013
for(n=0,30,print1(a(n),", "))

{a(n)=polcoeff( sum(m=0, n, m!*x^m/(1-x +x*O(x^n))^(m+1)/prod(k=1, m, 1+k*x +x*O(x^n))), n)}  /* From o.g.f. (Paul D. Hanna, Jul 20 2014) */
for(n=0, 25, print1(a(n), ", "))

[sum((n-j+1)^j for j in (0..n)) for n in (0..30)] # G. C. Greubel, Jun 15 2021

a(n) = A003101(n) + 1.
G.f.: Sum_{n>=0} x^n/(1 - (n+1)*x). - Paul D. Hanna, Sep 13 2011
G.f.: G(0) where G(k) = 1 + x*(2*k*x-1)/((2*k*x+x-1) - x*(2*k*x+x-1)^2/(x*(2*k*x+x-1) + (2*k*x+2*x-1)/G(k+1) )); (continued fraction). - Sergei N. Gladkovskii, Jan 26 2013
E.g.f.: Sum_{n>=0} Integral^n exp((n+1)*x) dx^n, where Integral^n F(x) dx^n is the n-th integration of F(x) with no constant of integration. - Paul D. Hanna, Dec 28 2013
O.g.f.: Sum_{n>=0} n! * x^n/(1-x)^(n+1) / Product_{k=1..n} (1 + k*x). - Paul D. Hanna, Jul 20 2014
a(n) = A101494(n+1,0). - Vladimir Kruchinin, Apr 01 2015
a(n-1) = Sum_{k = 1..n} k^(n-k). - Gus Wiseman, Jan 08 2019
log(a(n)) ~ (1 - 1/LambertW(exp(1)*n)) * n * log(1 + n/LambertW(exp(1)*n)). - Vaclav Kotesovec, Jun 15 2021
a(n) ~ sqrt(2*Pi/(n+1 + (n+1)/w(n))) * ((n+1)/w(n))^(n+2 - (n+1)/w(n)), where w(n) = LambertW(exp(1)*(n+1)). - Vaclav Kotesovec, Jun 25 2021, after user "leonbloy", see Mathematics Stack Exchange link.

a(23)-a(25) from Paul D. Hanna, Dec 28 2013

A319193 Irregular triangle where T(n,k) is the number of permutations of the integer partition with Heinz number A215366(n,k).

Original entry on oeis.org

1, 1, 1, 1, 1, 2, 1, 1, 1, 2, 3, 1, 1, 2, 2, 3, 3, 4, 1, 1, 2, 2, 1, 1, 3, 6, 6, 4, 5, 1, 1, 2, 2, 2, 6, 3, 3, 3, 4, 4, 12, 10, 5, 6, 1, 1, 2, 2, 1, 3, 2, 3, 6, 6, 3, 1, 12, 4, 12, 6, 10, 5, 20, 15, 6, 7, 1, 1, 2, 2, 2, 3, 2, 6, 3, 3, 4, 6, 6, 1, 12, 12, 4, 12

Offset: 0

Gus Wiseman, Sep 13 2018

A refinement of Pascal's triangle, these are the unsigned coefficients appearing in the expansion of homogeneous symmetric functions in terms of elementary symmetric functions.

			Triangle begins:
  1
  1
  1  1
  1  2  1
  1  1  2  3  1
  1  2  2  3  3  4  1
  1  2  2  1  1  3  6  6  4  5  1
The fourth row corresponds to the symmetric function identity: h(4) = -e(4) + e(22) + 2 e(31) - 3 e(211) + e(1111).

Alois P. Heinz, Rows n = 0..33, flattened

A different row ordering is A072811.

Cf. A000041, A005651, A008277, A008480, A056239, A124794, A215366, A319182, A319191, A319192.

b:= proc(n, i) option remember; `if`(n=0 or i<2, [2^n], [seq(
      map(p-> p*ithprime(i)^j, b(n-i*j, i-1))[], j=0..n/i)])
    end:
T:= n-> map(m-> (l-> add(i, i=l)!/mul(i!, i=l))(map(
        i-> i[2], ifactors(m)[2])), sort(b(n$2)))[]:
seq(T(n), n=0..10);  # Alois P. Heinz, Feb 14 2020

primeMS[n_]:=If[n==1,{},Flatten[Cases[FactorInteger[n],{p_,k_}:>Table[PrimePi[p],{k}]]]];
Table[Length[Permutations[primeMS[k]]],{n,6},{k,Sort[Times@@Prime/@#&/@IntegerPartitions[n]]}]
(* Second program: *)
b[n_, i_] := b[n, i] = If[n == 0 || i < 2, {2^n}, Flatten[Table[ #*Prime[i]^j& /@ b[n - i*j, i - 1], {j, 0, n/i}]]];
T[n_] := Map[Function[m, Function[l, Total[l]!/Times @@ (l!)][ FactorInteger[m][[All, 2]]]], Sort[b[n, n]]];
T /@ Range[0, 10] // Flatten (* Jean-François Alcover, May 10 2021, after Alois P. Heinz *)

T(n,k) = A008480(A215366(n,k)).

T(0,1)=1 prepended by Alois P. Heinz, Feb 14 2020

A001497 Triangle of coefficients of Bessel polynomials (exponents in decreasing order).

Original entry on oeis.org

1, 1, 1, 3, 3, 1, 15, 15, 6, 1, 105, 105, 45, 10, 1, 945, 945, 420, 105, 15, 1, 10395, 10395, 4725, 1260, 210, 21, 1, 135135, 135135, 62370, 17325, 3150, 378, 28, 1, 2027025, 2027025, 945945, 270270, 51975, 6930, 630, 36, 1, 34459425, 34459425, 16216200, 4729725, 945945, 135135, 13860, 990, 45, 1

Offset: 0

N. J. A. Sloane

The (reverse) Bessel polynomials P(n,x):=Sum_{m=0..n} a(n,m)*x^m, the row polynomials, called Theta_n(x) in the Grosswald reference, solve x*(d^2/dx^2)P(n,x) - 2*(x+n)*(d/dx)P(n,x) + 2*n*P(n,x) = 0.
With the related Sheffer associated polynomials defined by Carlitz as
B(0,x) = 1
B(1,x) = x
B(2,x) = x + x^2
B(3,x) = 3 x + 3 x^2 + x^3
B(4,x) = 15 x + 15 x^2 + 6 x^3 + x^4
... (see Mathworld reference), then P(n,x) = 2^n * B(n,x/2) are the Sheffer polynomials described in A119274. - Tom Copeland, Feb 10 2008
Exponential Riordan array [1/sqrt(1-2x), 1-sqrt(1-2x)]. - Paul Barry, Jul 27 2010
From Vladimir Kruchinin, Mar 18 2011: (Start)
For B(n,k){...} the Bell polynomial of the second kind we have
B(n,k){f', f'', f''', ...} = T(n-1,k-1)*(1-2*x)^(k/2-n), where f(x) = 1-sqrt(1-2*x).
The expansions of the first few rows are:
1/sqrt(1-2*x);
1/(1-2*x)^(3/2), 1/(1-2*x);
3/(1-2*x)^(5/2), 3/(1-2*x)^2, 1/(1-2*x)^(3/2);
15/(1-2*x)^(7/2), 15/(1-2*x)^3, 6/(1-2*x)^(5/2), 1/(1-2*x)^2. (End)
Also the Bell transform of A001147 (whithout column 0 which is 1,0,0,...). For the definition of the Bell transform see A264428. - Peter Luschny, Jan 19 2016
Antidiagonals of A099174 are rows of this entry. Dividing each diagonal by its first element generates A054142. - Tom Copeland, Oct 04 2016
The row polynomials p_n(x) of A107102 are (-1)^n B_n(1-x), where B_n(x) are the modified Carlitz-Bessel polynomials above, e.g., (-1)^2 B_2(1-x) = (1-x) + (1-x)^2 = 2 - 3 x + x^2 = p_2(x). - Tom Copeland, Oct 10 2016
a(n-1,m-1) counts rooted unordered binary forests with n labeled leaves and m roots. - David desJardins, Feb 23 2019
From Jianing Song, Nov 29 2021: (Start)
The polynomials P_n(x) = Sum_{k=0..n} T(n,k)*x^k satisfy: P_n(x) - (d/dx)P_n(x) = x*P_{n-1}(x) for n >= 1.
{P(n,x)} are related to the Fourier transform of 1/(1+x^2)^(n+1) and x/(1+x^2)^(n+2):
(i) For n >= 0, real number t, we have Integral_{x=-oo..oo} exp(-i*t*x)/(1+x^2)^(n+1) dx = Pi/(2^n*n!) * P_n(|t|) * exp(-|t|);
(ii) For n >= 0, real number t, we have Integral_{x=-oo..oo} x*exp(-i*t*x)/(1+x^2)^(n+2) dx = Pi/(2^(n+1)*(n+1)!) * ((-t)*P_n(-|t|)) * exp(-|t|). (End)
Suppose that f(x) is an n-times differentiable function defined on (a,b) for 0 <= a < b <= +oo, then for n >= 1, the n-th derivative of f(sqrt(x)) on (a^2,b^2) is Sum_{k=1..n} ((-1)^(n-k)*T(n-1,k-1)*f^(k)(sqrt(x))) / (2^n*x^(n-(k/2))), where f^(k) is the k-th derivative of f. - Jianing Song, Nov 30 2023

			Triangle begins
        1,
        1,       1,
        3,       3,      1,
       15,      15,      6,      1,
      105,     105,     45,     10,     1,
      945,     945,    420,    105,    15,    1,
    10395,   10395,   4725,   1260,   210,   21,   1,
   135135,  135135,  62370,  17325,  3150,  378,  28,  1,
  2027025, 2027025, 945945, 270270, 51975, 6930, 630, 36, 1
Production matrix begins
       1,      1,
       2,      2,      1,
       6,      6,      3,     1,
      24,     24,     12,     4,     1,
     120,    120,     60,    20,     5,    1,
     720,    720,    360,   120,    30,    6,   1,
    5040,   5040,   2520,   840,   210,   42,   7,  1,
   40320,  40320,  20160,  6720,  1680,  336,  56,  8, 1,
  362880, 362880, 181440, 60480, 15120, 3024, 504, 72, 9, 1
This is the exponential Riordan array A094587, or [1/(1-x),x], beheaded.
- _Paul Barry_, Mar 18 2011

J. Riordan, Combinatorial Identities, Wiley, 1968, p. 77.

T. D. Noe, Rows n = 0..50 of triangle, flattened
Peter Bala, Generalized Dobinski formulas
Paul Barry, Riordan Arrays, Orthogonal Polynomials as Moments, and Hankel Transforms, J. Int. Seq. 14 (2011) # 11.2.2, chapter 8.
Tom Copeland, A Class of Differential Operators and the Stirling Numbers
E. Deutsch, L. Ferrari and S. Rinaldi, Production Matrices, Advances in Applied Mathematics, 34 (2005) pp. 101-122.
O. Frink and H. L. Krall, A new class of orthogonal polynomials, Trans. Amer. Math. Soc. 65,100-115, 1945. [From _Roger L. Bagula_, Feb 15 2009]
E. Grosswald, Bessel Polynomials, Lecture Notes Math. vol. 698 1978 p. 18.
Milan Janjic, Some classes of numbers and derivatives, JIS 12 (2009) #09.8.3.
Wolfdieter Lang, On generalizations of Stirling number triangles, J. Integer Seqs., Vol. 3 (2000), #00.2.4.
B. Leclerc, Powers of staircase Schur functions and symmetric analogues of Bessel polynomials, Discrete Math., 153 (1996), 213-227.
Robert S. Maier, Boson Operator Ordering Identities from Generalized Stirling and Eulerian Numbers, arXiv:2308.10332 [math.CO], 2023. See p. 19.
Toufik Mansour, Matthias Schork and Mark Shattuck, The Generalized Stirling and Bell Numbers Revisited, Journal of Integer Sequences, Vol. 15 (2012), #12.8.3.
Toufik Mansour, Matthias Schork and Mark Shattuck, On the Stirling numbers associated with the meromorphic Weyl algebra, Applied Mathematics Letters, Volume 25, Issue 11, November 2012, Pages 1767-1771. - From _N. J. A. Sloane_, Sep 15 2012
W. Mlotkowski and A. Romanowicz, A family of sequences of binomial type, Probability and Mathematical Statistics, Vol. 33, Fasc. 2 (2013), pp. 401-408.
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.
Feng Qi and Bai-Ni Guo, "Some Properties and Generalizations of the Catalan, Fuss, and Fuss-Catalan Numbers", Mathematical Analysis and Applications : Selected Topics (2018), Wiley, Ch. 5, 101-133.
Feng Qi, X.-T. Shi and F.-F. Liu, Several formulas for special values of the Bell polynomials of the second kind and applications, Preprint 2015.
Alexander Stoimenow, On the number of chord diagrams, Discr. Math. 218 (2000), 209-233. Lemma 2.2.
Eric Weisstein's World of Mathematics, Bessel Polynomial
Index entries for sequences related to Bessel functions or polynomials

Reflected version of A001498 which is considered the main entry.
Other versions of this same triangle are given in A144299, A111924 and A100861.
Row sums give A001515. a(n, 0)= A001147(n) (double factorials).
Cf. A104556 (matrix inverse). A039683, A122850.
Cf. A245066 (central terms).
Cf. A054142, A099174, A107102.

a001497 n k = a001497_tabl !! n !! k
a001497_row n = a001497_tabl !! n
a001497_tabl = [1] : f [1] 1 where
   f xs z = ys : f ys (z + 2) where
     ys = zipWith (+) ([0] ++ xs) (zipWith (*) [z, z-1 ..] (xs ++ [0]))
-- Reinhard Zumkeller, Jul 11 2014

/* As triangle */ [[Factorial(2*n-k)/(Factorial(k)*Factorial(n-k)*2^(n-k)): k in [0..n]]: n in [0.. 15]]; // Vincenzo Librandi, Aug 12 2015

f := proc(n) option remember; if n <=1 then (1+x)^n else expand((2*n-1)*x*f(n-1)+f(n-2)); fi; end;
row := n -> seq(coeff(f(n), x, n - k), k = 0..n): seq(row(n), n = 0..9);

m = 9; Flatten[ Table[(n + k)!/(2^k*k!*(n - k)!), {n, 0, m}, {k, n, 0, -1}]] (* Jean-François Alcover, Sep 20 2011 *)
y[n_, x_] := Sqrt[2/(Pi*x)]*E^(1/x)*BesselK[-n-1/2, 1/x]; t[n_, k_] := Coefficient[y[n, x], x, k]; Table[t[n, k], {n, 0, 9}, {k, n, 0, -1}] // Flatten (* Jean-François Alcover, Mar 01 2013 *)

T(k, n) = if(n>k||k<0||n<0,0,(2*k-n)!/(n!*(k-n)!*2^(k-n))) /* Ralf Stephan */

{T(n, k) = if( k<0 || k>n, 0, binomial(n, k)*(2*n-k)!/2^(n-k)/n!)}; /* Michael Somos, Oct 03 2006 */

# uses[bell_matrix from A264428]
# Adds a column 1,0,0,0, ... at the left side of the triangle.
bell_matrix(lambda n: A001147(n), 9) # Peter Luschny, Jan 19 2016

a(n, m) = (2*n-m)!/(m!*(n-m)!*2^(n-m)) if n >= m >= 0 else 0 (from Grosswald, p. 7).
a(n, m)= 0, n= m >= 0 (from Grosswald p. 23, (19)).
E.g.f. for m-th column: ((1-sqrt(1-2*x))^m)/(m!*sqrt(1-2*x)).
G.f.: 1/(1-xy-x/(1-xy-2x/(1-xy-3x/(1-xy-4x/(1-.... (continued fraction). - Paul Barry, Jan 29 2009
T(n,k) = if(k<=n, C(2n-k,2(n-k))*(2(n-k)-1)!!,0) = if(k<=n, C(2n-k,2(n-k))*A001147(n-k),0). - Paul Barry, Mar 18 2011
Row polynomials for n>=1 are given by 1/t*D^n(exp(x*t)) evaluated at x = 0, where D is the operator 1/(1-x)*d/dx. - Peter Bala, Nov 25 2011
The matrix product A039683*A008277 gives a signed version of this triangle. Dobinski-type formula for the row polynomials: R(n,x) = (-1)^n*exp(x)*Sum_{k = 0..inf} k*(k-2)*(k-4)*...*(k-2*(n-1))*(-x)^k/k!. Cf. A122850. - Peter Bala, Jun 23 2014

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

A132440 Infinitesimal Pascal matrix: generator (lower triangular matrix representation) of the Pascal matrix, the classical operator xDx, iterated Laguerre transforms, associated matrices of the list partition transform and general Euler transformation for sequences.

Original entry on oeis.org

0, 1, 0, 0, 2, 0, 0, 0, 3, 0, 0, 0, 0, 4, 0, 0, 0, 0, 0, 5, 0, 0, 0, 0, 0, 0, 6, 0, 0, 0, 0, 0, 0, 0, 7, 0, 0, 0, 0, 0, 0, 0, 0, 8, 0, 0, 0, 0, 0, 0, 0, 0, 0, 9, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 10, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 11, 0

Offset: 0

Tom Copeland, Nov 13 2007, Nov 15 2007, Nov 22 2007, Dec 02 2007

Let M(t) = exp(t*T) = lim_{n->oo} (1 + t*T/n)^n.
Pascal matrix = [ binomial(n,k) ] = M(1) = exp(T), truncating the series gives the n X n submatrices.
Inverse Pascal matrix = M(-1) = exp(-T) = matrix for inverse binomial transform.
A(j) = T^j / j! equals the matrix [binomial(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 Pascal triangle. Hence the A(j)'s form a linearly independent basis for all matrices of the form [binomial(n,k) * d(n-k)] which include as a subset the invertible associated matrices of the list partition transform (LPT) of A133314.
For sequences with b(0) = 1, umbrally,
M[b(.)] = exp(b(.)*T) = [ binomial(n,k) * b(n-k) ] = matrices associated to b by LPT.
[M[b(.)]]^(-1) = exp(c(.)*T) = [ binomial(n,k) * c(n-k) ] = matrices associated to c, where c = LPT(b) . Or,
[M[b(.)]]^(-1) = exp[LPT(b(.))*T] = LPT[M(b(.))] = M[LPT(b(.))]= M[c(.)].
This is related to xDx, the iterated Laguerre transform and the general Euler transformation of a sequence through the comments in A132013 and A132014 and the relation [Sum_{k=0..n} binomial(n,k) * b(n-k) * d(k)] = M(b)*d, (n-th term). See also A132382.
If b(n,x) is a binomial type Sheffer sequence, then M[b(.,x)]*s(y) = s(x+y) when s(y) = (s(0,y),s(1,y),s(2,y),...) is an array for a Sheffer sequence with the same delta operator as b(n,x) and [M[b(.,x)]]^(-1) is given by the formulas above with b(n) replaced by b(n,x) as b(0,x)=1 for a binomial-type Sheffer sequence.
T = I - A132013 and conversely A132013 = I - T, which is the matrix representation for the iterated mixed order Laguerre transform characterized in A132013 (and A132014).
(I-T)^m generates the group [A132013]^m for m = 0,1,2,... discussed in A132014.
The inverse is 1/(I-T) = I + T + T^2 + T^3 + ... = [A132013]^(-1) = A094587 with the associated sequence (0!,1!,2!,3!,...) under the LPT.
And 1/(I-T)^2 = I + 2*T + 3*T^2 + 4*T^3 + ... = [A132013]^(-2) = A132159 with the associated sequence (1!,2!,3!,4!,...) under the LPT.
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 * a(n-1),
2) B(x) = xDx A(x)
3) B(x) = x * Lag(1,-:xD:) A(x)
4) EB(x) = x * EA(x) where D is the derivative w.r.t. x, (:xD:)^j = x^j*D^j and Lag(n,x) is the Laguerre polynomial.
So the exponentiated operator can be characterized as
5) exp(t*T) A(x) = exp(t*xDx) A(x) = [Sum_{n=0,1,...} (t*x)^n * Lag(n,-:xD:)] A(x) = [exp{[t*u/(1-t*u)]*:xD:} / (1-t*u) ] A(x) (eval. at u=x) = A[x/(1-t*x)]/(1-t*x), a generalized Euler transformation for an o.g.f.,
6) exp(t*T) EA(x) = exp(t*x)*EA(x) = exp[(t+a(.))*x], gen. Euler trf. for an e.g.f.
7) exp(t*T) * a = M(t) * a = [Sum_{k=0..n} binomial(n,k) * t^(n-k) * a(k)].
The umbral extension of formulas 5, 6 and 7 gives formally
8) exp[c(.)*T] A(x) = exp(c(.)*xDx) A(x) = [Sum_{n>=0} (c(.)*x)^n * Lag(n,-:xD:)] A(x) = [exp{[c(.)*u/(1-c(.)*u)]*:xD:} / (1-c(.)*u) ] A(x) (eval. at u=x) = A[x/(1-c(.)*x)]/(1-c(.)*x), where the umbral evaluation should be applied only after a power series in c is obtained,
9) exp[c(.)*T] EA(x) = exp(c(.)*x)*EA(x) = exp[(c(.)+a(.))*x]
10) exp[c(.)*T] * a = M[c(.)] * a = [Sum_{k=0..n} binomial(n,k) * c(n-k) * a(k)] .
The n X n principal submatrix of T is nilpotent, in particular, [Tsub_n]^(n+1) = 0, n=0,1,2,3,....
Note (xDx)^n = x^n D^n x^n = x^n n! (:Dx:)^n/n! = x^n n! Lag(n,-:xD:).
The operator xDx is an important, classical operator explored by among others Dattoli, Al-Salam, Carlitz and Stokes and even earlier investigators.
For a recent treatment of xDx, DxD and more general operators see the paper "Laguerre-type derivatives: Dobinski relations and combinatorial identities". - Karol A. Penson, Sep 15 2009
See Copeland's link for generalized Laguerre functions and connection to fractional differ-integrals in exercises through (:Dx:)^a/a!=(D^a x^a)/a!. - Tom Copeland, Nov 17 2011
From Tom Copeland, Apr 25 2014: (Start)
Conjugation or "similarity" transformations of [dP]=A132440 have an operator interpretation (cf. A074909 and A238363):
In general, select two operators A and B such that A^n = F1(n,B) and B^n = F2(n,A); then A^n = F1(n,F2(.,A)) and B^n = F2(n,F1(.,B)), evaluated umbrally, i.e., F1(n,F2(.,x))=F2(n,F1(.,x))=x^n, implying the polynomials F1 and F2 are an umbral compositional inverse pair.
One such pair are the Bell polynomials Bell(n,x) and falling factorials (x)_n with Bell(n,:xD:)=(xD)^n and (xD)_n=:xD:^n (cf. A074909). Another are the Laguerre polynomials LN(n,x)= n!*Lag(n,x) (A021009), which are umbrally self-inverse, with LN(n,-:xD:)=:Dx:^n and LN(n,:Dx:)= (-:xD:)^n with :Dx:^n=D^n*x^n.
Evaluating, for n>=0, the operator derivative d(B^n)/dA = d(F2(n,A))/dA in the basis B^n, i.e., with A^n finally replaced by F1(n,B), or A^n=F1(.,B)^n=F1(n,B), is equivalent to the matrix conjugation
A)   [F2]*[dP]*[F1]
B) = [F2]*[dP]*[F2]^(-1)
C) = [F1]^(-1)*[dP]*[F1],
where [F1] is the lower triangular matrix with the n-th row the coefficients of F1(n,x) and analogously for [F2].
So, given the row vector Rv=(c0 c1 c2 c3 ...) and the column vector Cv(x)=(1 x x^2 x^3 ...)^Transpose, form the power series V(x)=Rv*Cv(x).
D) dV(B)/dA = Rv * [F2]*[dP]*[F1] * Cv(B).
E) With A=D and B=D, F1(n,x)=F2(n,x)=x^n and [F1]=[F2]=I. Then d(B^n)/dA = d(D^n)/dD = n * D^(n-1); therefore, consistently [F2]*[dP]*[F1] = [dP] and dV(D)/dD = Rv * [dP] * Cv(D). (End)

			Matrix T begins
  0;
  1,0;
  0,2,0;
  0,0,3,0;
  0,0,0,4,0;
  ...

T. Mansour and M. Schork, Commutation Relations, Normal Ordering, and Stirling Numbers, Chapman and Hall/CRC, 2015, (x^n D^n x^n on p. 187).

Robert Israel, Table of n, a(n) for n = 0..10000
W. A. Al-Salam, Operational representations for the Laguerre and other polynomials, Duke Math. Jour., vol. 31 (1964), pp. 127-142.
Tom Copeland, Fractional Calculus, Gamma Classes, the Riemann Zeta Function, and an Appell Pair of Sequences.
Tom Copeland, Goin' with the Flow: Logarithm of the Derivative Operator.
Tom Copeland, Infinigens, the Pascal Pyramid, and the Witt and Virasoro Algebras.
Tom Copeland, The Inverse Mellin Transform, Bell Polynomials, a Generalized Dobinski Relation, and the Confluent Hypergeometric Functions (pdf).
Tom Copeland, Mathemagical Forests.
Tom Copeland, Mellin Interpolation of Differential Ops and Associated Infinigens and Appell Polynomials: The Ordered, Laguerre, and Scherk-Witt-Lie Diff Ops.
G. Hetyei, Meixner polynomials of the second kind and quantum algebras representing su(1,1), arXiv preprint arXiv:0909.4352 [math.QA], 2009. (Cf. Viennot's Laguerre histoires.)
K. A. Penson, P. Blasiak, A. Horzela, G.H.E. Duchamp and A. I. Solomon, Laguerre-type derivatives: Dobinski relations and combinatorial identities, arXiv:0904.0369 [math-ph], 2009.
K. A. Penson, P. Blasiak, A. Horzela, G.H.E. Duchamp and A. I. Solomon, Laguerre-type derivatives: Dobinski relations and combinatorial identities, Journal of Mathematical Physics vol. 50, (2009) 083512.
G. Stokes, Note on certain formulae in the calculus of operations, Proceedings of the Royal Society of Edinburgh, IX, pp. 101-102, 1876.

seq(op([0$i,i]),i=1..20); # Robert Israel, Oct 02 2015

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

T = log(P) with the Pascal matrix P:=A007318. This should be read as T_N = log(P_N) with P_N the N X N matrix P, N>=2. Because P_N is lower triangular with all diagonal elements 1, the series log(1_N-(1_N-P_N)) stops after N-1 terms because (1_N-P_N)^N is the 0_N-matrix. - Wolfdieter Lang, Oct 14 2010
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*R 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
From Tom Copeland, Apr 26 2014: (Start)
A) T = exp(A238385-I) - I
B)   = [St1]*P*[St2] - I
C)   = [St1]*P*[St1]^(-1) - I
D)   = [St2]^(-1)*P*[St2] - I
E)   = [St2]^(-1)*P*[St1]^(-1) - I
where P=A007318, [St1]=padded A008275 just as [St2]=A048993=padded A008277, and I=identity matrix. (End)
From Robert Israel, Oct 02 2015: (Start)
G.f. Sum_{k >= 1} k x^((k+3/2)^2/2 - 17/8) is related to Jacobi theta functions.
If 8*n+17 = y^2 is a square, then a(n) = (y-3)/2, otherwise a(n) = 0. (End)

Missing zero added in table by Tom Copeland, Feb 25 2014

A049029 Triangle read by rows, the Bell transform of the quartic factorial numbers A007696(n+1) without column 0.

Original entry on oeis.org

1, 5, 1, 45, 15, 1, 585, 255, 30, 1, 9945, 5175, 825, 50, 1, 208845, 123795, 24150, 2025, 75, 1, 5221125, 3427515, 775845, 80850, 4200, 105, 1, 151412625, 108046575, 27478710, 3363045, 219450, 7770, 140, 1, 4996616625, 3824996175, 1069801425

Offset: 1

Wolfdieter Lang

Previous name was: Triangle of numbers related to triangle A048882; generalization of Stirling numbers of second kind A008277, Lah-numbers A008297, ...
a(n,m) enumerates unordered n-vertex m-forests composed of m plane increasing quintic (5-ary) trees. Proof based on the a(n,m) recurrence. See also the F. Bergeron et al. reference, especially Table 1, first row and Example 1 for the e.g.f. for m=1. - Wolfdieter Lang, Sep 14 2007
Also the Bell transform of A007696(n+1). For the definition of the Bell transform see A264428. - Peter Luschny, Jan 28 2016

			Triangle starts:
{1};
{5,1};
{45,15,1};
{585,255,30,1};
{9945,5175,825,50,1};
...

F. Bergeron, Ph. Flajolet and B. Salvy, Varieties of Increasing Trees, in Lecture Notes in Computer Science vol. 581, ed. J.-C. Raoult, Springer 1992, pp. 24-48.
P. Blasiak, K. A. Penson and A. I. Solomon, The general boson normal ordering problem, arXiv:quant-ph/0402027, 2004.
T. Copeland, Mathemagical Forests
T. Copeland, Addendum to Mathemagical Forests
T. Copeland, A Class of Differential Operators and the Stirling Numbers
M. Janjic, Some classes of numbers and derivatives, JIS 12 (2009) #09.8.3.
W. Lang, On generalizations of Stirling number triangles, J. Integer Seqs., Vol. 3 (2000), #00.2.4.
W. Lang, First 10 rows.
Shi-Mei Ma, Some combinatorial sequences associated with context-free grammars, arXiv:1208.3104v2 [math.CO]. - From N. J. A. Sloane, Aug 21 2012
E. Neuwirth, Recursively defined combinatorial Functions: Extending Galton's board, Discr. Maths. 239 (2001) 33-51.
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.

a(n, m) := S2(5, n, m) is the fifth triangle of numbers in the sequence S2(1, n, m) := A008277(n, m) (Stirling 2nd kind), S2(2, n, m) := A008297(n, m) (Lah), S2(3, n, m) := A035342(n, m), S2(4, n, m) := A035469(n, m). a(n, 1)= A007696(n). A007559(n).
Cf. A048882, A007696. Row sums: A049120(n), n >= 1.
Cf. A094638

# The function BellMatrix is defined in A264428.
# Adds (1,0,0,0, ..) as column 0.
BellMatrix(n -> mul(4*k+1, k=0..n), 9); # Peter Luschny, Jan 28 2016

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 *)
rows = 9;
a[n_, m_] := BellY[n, m, Table[Product[4k+1, {k, 0, j}], {j, 0, rows}]];
Table[a[n, m], {n, 1, rows}, {m, 1, n}] // Flatten (* Jean-François Alcover, Jun 22 2018 *)

a(n, m) = n!*A048882(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

a(n, m) = sum(|A051142(n, j)|*S2(j, m), j=m..n) (matrix product), with S2(j, m) := A008277(j, m) (Stirling2 triangle). Priv. comm. to W. Lang by E. Neuwirth, Feb 15 2001; see also the 2001 Neuwirth reference. See the general comment on products of Jabotinsky matrices given under A035342.

From Peter Bala, Nov 25 2011: (Start)

E.g.f.: G(x,t) = exp(t*A(x)) = 1+t*x+(5*t+t^2)*x^2/2!+(45*t+15*t^2+t^3)*x^3/3!+..., where A(x) = -1+(1-4*x)^(-1/4) satisfies the autonomous differential equation A'(x) = (1+A(x))^5.

The generating function G(x,t) satisfies the partial differential equation t*(dG/dt+G) = (1-4*x)*dG/dx, from which follows the recurrence given above.

The row polynomials are given by D^n(exp(x*t)) evaluated at x = 0, where D is the operator (1+x)^5*d/dx. Cf. A008277 (D = (1+x)*d/dx), A105278 (D = (1+x)^2*d/dx), A035342 (D = (1+x)^3*d/dx) and A035469 (D = (1+x)^4*d/dx).

(End)

New name from Peter Luschny, Jan 30 2016

A111596 The matrix inverse of the unsigned Lah numbers A271703.

Original entry on oeis.org

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

Offset: 0

Wolfdieter Lang, Aug 23 2005

Also the associated Sheffer triangle to Sheffer triangle A111595.
Coefficients of Laguerre polynomials (-1)^n * n! * L(n,-1,x), which equals (-1)^n * Lag(n,x,-1) below. Lag(n,Lag(.,x,-1),-1) = x^n evaluated umbrally, i.e., with (Lag(.,x,-1))^k = Lag(k,x,-1). - Tom Copeland, Apr 26 2014
Without row n=0 and column m=0 this is, up to signs, the Lah triangle A008297.
The unsigned column sequences are (with leading zeros): A000142, A001286, A001754, A001755, A001777, A001778, A111597-A111600 for m=1..10.
The row polynomials p(n,x) := Sum_{m=0..n} a(n,m)*x^m, together with the row polynomials s(n,x) of A111595 satisfy the exponential (or binomial) convolution identity s(n,x+y) = Sum_{k=0..n} binomial(n,k)*s(k,x)*p(n-k,y), n>=0.
Exponential Riordan array [1,x/(1+x)]. Inverse of the exponential Riordan array [1,x/(1-x)], which is the unsigned version of A111596. - Paul Barry, Apr 12 2007
For the unsigned subtriangle without column number m=0 and row number n=0, see A105278.
Unsigned triangle also matrix product |S1|*S2 of Stirling number matrices.
The unsigned row polynomials are Lag(n,-x,-1), the associated Laguerre polynomials of order -1 with negated argument. See Gradshteyn and Ryzhik, Abramowitz and Stegun and Rota (Finite Operator Calculus) for extensive formulas. - Tom Copeland, Nov 17 2007, Sep 09 2008
An infinitesimal matrix generator for unsigned A111596 is given by A132792. - Tom Copeland, Nov 22 2007
From the formalism of A132792 and A133314 for n > k, unsigned A111596(n,k) = a(k) * a(k+1)...a(n-1) / (n-k)! = a generalized factorial, where a(n) = A002378(n) = n-th term of first subdiagonal of unsigned A111596. Hence Deutsch's remark in A002378 provides an interpretation of A111596(n,k) in terms of combinations of certain circular binary words. - Tom Copeland, Nov 22 2007
Given T(n,k)= A111596(n,k) and 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 = C where C(n,k) = T(n,k)*[a(.)+b(.)]^(n-k), umbrally. - Tom Copeland, Aug 27 2008
Operationally, the unsigned row polynomials may be expressed as p_n(:xD:) = x*:Dx:^n*x^{-1}=x*D^nx^n*x^{-1}= n!*binomial(xD+n-1,n) = (-1)^n n! binomial(-xD,n) = n!L(n,-1,-:xD:), where, by definition, :AB:^n = A^nB^n for any two operators A and B, D = d/dx, and L(n,-1,x) is the Laguerre polynomial of order -1. A similarity transformation of the operators :Dx:^n generates the higher order Laguerre polynomials, which can also be expressed in terms of rising or falling factorials or Kummer's confluent hypergeometric functions (cf. the Mathoverflow post). -  Tom Copeland, Sep 21 2019

			Binomial convolution of row polynomials: p(3,x) = 6*x-6*x^2+x^3; p(2,x) = -2*x+x^2, p(1,x) = x, p(0,x) = 1,
together with those from A111595: s(3,x) = 9*x-6*x^2+x^3; s(2,x) = 1-2*x+x^2, s(1,x) = x, s(0,x) = 1; therefore
9*(x+y)-6*(x+y)^2+(x+y)^3 = s(3,x+y) = 1*s(0,x)*p(3,y) + 3*s(1,x)*p(2,y) + 3*s(2,x)*p(1,y) +1*s(3,x)*p(0,y) = (6*y-6*y^2+y^3) + 3*x*(-2*y+y^2) + 3*(1-2*x+x^2)*y + 9*x-6*x^2+x^3.
From _Wolfdieter Lang_, Apr 28 2014: (Start)
The triangle a(n,m) begins:
n\m  0     1       2       3      4     5   6  7
0:   1
1:   0     1
2:   0    -2       1
3:   0     6      -6       1
4:   0   -24      36     -12      1
5:   0   120    -240     120    -20     1
6:   0  -720    1800   -1200    300   -30   1
7:   0  5040  -15120   12600  -4200   630 -42  1
...
For more rows see the link.
(End)

G. C. Greubel, Rows n=0..100 of triangle, flattened
Wolfdieter Lang, The first 11 rows of the triangle.
M. Abramowitz and I. A. Stegun, eds., Handbook of Mathematical Functions, National Bureau of Standards, Applied Math. Series 55, Tenth Printing, 1972 [alternative scanned copy].
Paul Barry, The Restricted Toda Chain, Exponential Riordan Arrays, and Hankel Transforms, J. Int. Seq. 13 (2010) # 10.8.4, example 4.
Paul Barry, Exponential Riordan Arrays and Permutation Enumeration, J. Int. Seq. 13 (2010) # 10.9.1, example 6.
Paul Barry, Riordan Arrays, Orthogonal Polynomials as Moments, and Hankel Transforms, J. Int. Seq. 14 (2011) # 11.2.2, example 20.
Paul Barry, Combinatorial polynomials as moments, Hankel transforms and exponential Riordan arrays, arXiv preprint arXiv:1105.3044 [math.CO], 2011, also J. Int. Seq. 14 (2011) 11.6.7.
Tom Copeland, A Class of Differential Operators and the Stirling Numbers; Generators, Inversion, and Matrix, Binomial, and Integral Transforms; Lagrange a la Lah
A. Hennessy and P. Barry, Generalized Stirling Numbers, Exponential Riordan Arrays, and Orthogonal Polynomials, J. Int. Seq. 14 (2011) # 11.8.2.
M. Janjic, Some classes of numbers and derivatives, JIS 12 (2009) 09.8.3.
Mathoverflow, Pochhammer symbol of a differential, and hypergeometric polynomials, a question posed by Emilio Pisanty and answered by Tom Copeland, 2012.
J. Taylor, Counting words with Laguerre polynomials, DMTCS Proc., Vol. AS, 2013, p. 1131-1142. [_Tom Copeland_, Jan 08 2016] [Broken link]
J. Taylor, Formal group laws and hypergraph colorings, doctoral thesis, Univ. of Wash., 2016, p. 96. [_Tom Copeland_, Dec 20 2018]
Jian Zhou, On Some Mathematics Related to the Interpolating Statistics, arXiv:2108.10514 [math-ph], 2021.

Row sums: A111884. Unsigned row sums: A000262.
A002868 gives maximal element (in magnitude) in each row.
Cf. A130561 for a natural refinement.
Cf. A264428, A264429, A271703 (unsigned).
Cf. A008297, A089231, A105278 (variants).
Cf. A002378, A008277, A132792, A133314, A132792, A001263.

# The function BellMatrix is defined in A264428.
BellMatrix(n -> `if`(n::odd, -(n+1)!, (n+1)!), 9); # Peter Luschny, Jan 27 2016

a[0, 0] = 1; a[n_, m_] := ((-1)^(n-m))*(n!/m!)*Binomial[n-1, m-1]; Table[a[n, m], {n, 0, 10}, {m, 0, n}] // Flatten (* Jean-François Alcover, Jul 05 2013 *)
T[ n_, k_] := (-1)^n n! Coefficient[ LaguerreL[ n, -1, x], x, k]; (* Michael Somos, Dec 15 2014 *)
rows = 9;
t = Table[(-1)^(n+1) n!, {n, 1, rows}];
T[n_, k_] := BellY[n, k, t];
Table[T[n, k], {n, 0, rows}, {k, 0, n}]  // Flatten (* Jean-François Alcover, Jun 22 2018, after Peter Luschny *)

{T(n, k) = if( n<1 || k<1, n==0 && k==0, (-1)^n * n! * polcoeff( sum(k=1, n, binomial( n-1, k-1) * (-x)^k / k!), k))}; /* Michael Somos, Dec 15 2014 */

lah_number = lambda n, k: factorial(n-k)*binomial(n,n-k)*binomial(n-1,n-k)
A111596_row = lambda n: [(-1)^(n-k)*lah_number(n, k) for k in (0..n)]
for n in range(10): print(A111596_row(n)) # Peter Luschny, Oct 05 2014

# uses[inverse_bell_transform from A264429]
def A111596_matrix(dim):
    fact = [factorial(n) for n in (1..dim)]
    return inverse_bell_transform(dim, fact)
A111596_matrix(10) # Peter Luschny, Dec 20 2015

E.g.f. m-th column: ((x/(1+x))^m)/m!, m>=0.
E.g.f. for row polynomials p(n, x) is exp(x*y/(1+y)).
a(n, m) = ((-1)^(n-m))*|A008297(n, m)| = ((-1)^(n-m))*(n!/m!)*binomial(n-1, m-1), n>=m>=1; a(0, 0)=1; else 0.
a(n, m) = -(n-1+m)*a(n-1, m) + a(n-1, m-1), n>=m>=0, a(n, -1):=0, a(0, 0)=1; a(n, m)=0 if n|a(n,m)| = Sum_{k=m..n} |S1(n,k)|*S2(k,m), n>=0. S2(n,m):=A048993. S1(n,m):=A048994. - Wolfdieter Lang, May 04 2007
From Tom Copeland, Nov 21 2011: (Start)
For this Lah triangle, the n-th row polynomial is given umbrally by
  (-1)^n n! binomial(-Bell.(-x),n), where Bell_n(-x)= exp(x)(xd/dx)^n exp(-x), the n-th Bell / Touchard / exponential polynomial with neg. arg., (cf. A008277). E.g., 2! binomial(-Bell.(-x),2) = -Bell.(-x)*(-Bell.(-x)-1) = Bell_2(-x)+Bell_1(-x) = -2x+x^2.
A Dobinski relation is (-1)^n n! binomial(-Bell.(-x),n)= (-1)^n n! e^x Sum_{j>=0} (-1)^j binomial(-j,n)x^j/j!= n! e^x Sum_{j>=0} (-1)^j binomial(j-1+n,n)x^j/j!. See the Copeland link for the relation to inverse Mellin transform. (End)
The n-th row polynomial is (-1/x)^n e^x (x^2*D_x)^n e^(-x). - Tom Copeland, Oct 29 2012
Let f(.,x)^n = f(n,x) = x!/(x-n)!, the falling factorial,and r(.,x)^n = r(n,x) = (x-1+n)!/(x-1)!, the rising factorial, then the Lah polynomials, Lah(n,t)= n!*Sum{k=1..n} binomial(n-1,k-1)(-t)^k/k! (extra sign factor on odd rows), give the transform Lah(n,-f(.,x))= r(n,x), and Lah(n,r(.,x))= (-1)^n * f(n,x). - Tom Copeland, Oct 04 2014
|T(n,k)| = Sum_{j=0..2*(n-k)} A254881(n-k,j)*k^j/(n-k)!. Note that A254883 is constructed analogously from A254882. - Peter Luschny, Feb 10 2015
The T(n,k) are the inverse Bell transform of [1!,2!,3!,...] and |T(n,k)| are the Bell transform of [1!,2!,3!,...]. See A264428 for the definition of the Bell transform and A264429 for the definition of the inverse Bell transform. - Peter Luschny, Dec 20 2015
Dividing each n-th diagonal by n!, where the main diagonal is n=1, generates a shifted, signed Narayana matrix A001263. - Tom Copeland, Sep 23 2020

New name using a comment from Wolfdieter Lang by Peter Luschny, May 10 2021

A035469 Triangle read by rows, the Bell transform of the triple factorial numbers A007559(n+1) without column 0.

Original entry on oeis.org

1, 4, 1, 28, 12, 1, 280, 160, 24, 1, 3640, 2520, 520, 40, 1, 58240, 46480, 11880, 1280, 60, 1, 1106560, 987840, 295960, 40040, 2660, 84, 1, 24344320, 23826880, 8090880, 1296960, 109200, 4928, 112, 1, 608608000, 643843200

Offset: 1

Wolfdieter Lang

Previous name was: Triangle of numbers related to triangle A035529; generalization of Stirling numbers of second kind A008277, Lah-numbers A008297 and A035342.
a(n,m) enumerates unordered n-vertex m-forests composed of m plane increasing quartic (4-ary) trees. Proof based on the a(n,m) recurrence. See a D. Callan comment on the m=1 case A007559. See also the F. Bergeron et al. reference, especially Table 1, first row and Example 1 for the e.g.f. for m=1. - Wolfdieter Lang, Sep 14 2007
For the definition of the Bell transform see A264428. - Peter Luschny, Jan 19 2016

			Triangle starts:
     {1}
     {4,    1}
    {28,   12,    1}
   {280,  160,   24,    1}
  {3640, 2520,  520,   40,    1}

F. Bergeron, Ph. Flajolet and B. Salvy, Varieties of Increasing Trees, in Lecture Notes in Computer Science vol. 581, ed. J.-C. Raoult, Springer 1922, pp. 24-48.

Peter Bala, Generalized Dobinski formulas
P. Blasiak, K. A. Penson and A. I. Solomon, The general boson normal ordering problem, arXiv:quant-ph/0402027, 2004.
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.
Tom Copeland, A Class of Differential Operators and the Stirling Numbers
Tom Copeland, Mathemagical Forests
Tom Copeland, Addendum to Mathemagical Forests
Milan Janjic, Some classes of numbers and derivatives, JIS 12 (2009) 09.8.3.
Wolfdieter Lang, On generalizations of Stirling number triangles, J. Integer Seqs., Vol. 3 (2000), #00.2.4.
Wolfdieter Lang, First 10 rows.
Shi-Mei Ma, Some combinatorial sequences associated with context-free grammars, arXiv:1208.3104 [math.CO], 2012. - From _N. J. A. Sloane_, Aug 21 2012
E. Neuwirth, Recursively defined combinatorial functions: Extending Galton's board, Discrete Math. 239 (2001) 33-51.
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.

a(n, m)=: S2(4, n, m) is the fourth triangle of numbers in the sequence S2(1, n, m) := A008277(n, m) (Stirling 2nd kind), S2(2, n, m) := A008297(n, m) (Lah), S2(3, n, m) := A035342(n, m). a(n, 1)= A007559(n).
Row sums: A049119(n), n >= 1.
Cf. A094638.

a[n_, m_] /; n >= m >= 1 := a[n, m] = (3(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 *)
rows = 9;
a[n_, m_] := BellY[n, m, Table[Product[3k+1, {k, 0, j}], {j, 0, rows}]];
Table[a[n, m], {n, 1, rows}, {m, 1, n}] // Flatten (* Jean-François Alcover, Jun 22 2018 *)

# uses[bell_matrix from A264428]
# Adds a column 1,0,0,0, ... at the left side of the triangle.
bell_matrix(lambda n: A007559(n+1) , 9) # Peter Luschny, Jan 19 2016

a(n, m) = Sum_{j=m..n} |A051141(n, j)|*S2(j, m) (matrix product), with S2(j, m):=A008277(j, m) (Stirling2 triangle). Priv. comm. to Wolfdieter Lang by E. Neuwirth, Feb 15 2001; see also the 2001 Neuwirth reference. See the general comment on products of Jabotinsky matrices given under A035342.
a(n, m) = n!*A035529(n, m)/(m!*3^(n-m)); a(n+1, m) = (3*n+m)*a(n, m) + a(n, m-1), n >= m >= 1; a(n, m) := 0, n < m; a(n, 0) := 0, a(1, 1)=1;
E.g.f. of m-th column: ((-1+(1-3*x)^(-1/3))^m)/m!.
From Peter Bala, Nov 25 2011: (Start)
E.g.f.: G(x,t) = exp(t*A(x)) = 1 + t*x + (4*t+t^2)*x^2/2! + (28*t + 12*t^2 + t^3)*x^3/3! + ..., where A(x) = -1 + (1-3*x)^(-1/3) satisfies the autonomous differential equation A'(x) = (1+A(x))^4.
The generating function G(x,t) satisfies the partial differential equation t*(dG/dt+G) = (1-3*x)*dG/dx, from which follows the recurrence given above.
The row polynomials are given by D^n(exp(x*t)) evaluated at x = 0, where D is the operator (1+x)^4*d/dx. Cf. A008277 (D = (1+x)*d/dx), A105278 (D = (1+x)^2*d/dx), A035342 (D = (1+x)^3*d/dx) and A049029 (D = (1+x)^5*d/dx).
(End)
Dobinski-type formula for the row polynomials: R(n,x) = exp(-x)*Sum_{k>=0} k*(k+3)*(k+6)*...*(k+3*(n-1))*x^k/k!. - Peter Bala, Jun 23 2014

New name from Peter Luschny, Jan 19 2016

Previous Showing 81-90 of 644 results. Next

cp's OEIS Frontend

A019536 Number of length n necklaces with integer entries that cover an initial interval of positive integers.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Mathematica

PARI

Formula

Extensions

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

Views

Author

Keywords

Comments

References

Links

Crossrefs

Programs

GAP

Magma

Maple

Mathematica

PARI

Python

Formula

A026898 a(n) = Sum_{k=0..n} (n-k+1)^k.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Haskell

Magma

Maple

Mathematica

PARI

PARI

PARI

Sage

Formula

Extensions

A319193 Irregular triangle where T(n,k) is the number of permutations of the integer partition with Heinz number A215366(n,k).

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Maple

Mathematica

Formula

Extensions

A001497 Triangle of coefficients of Bessel polynomials (exponents in decreasing order).

Views

Author

Keywords

Comments

Examples

References

Links

Crossrefs

Programs

Haskell

Magma

Maple

Mathematica

PARI

PARI

Sage

Formula