A000914 -id:A000914 - cp's OEIS frontend

A008275 Triangle read by rows of Stirling numbers of first kind, s(n,k), n >= 1, 1 <= k <= n.

1, -1, 1, 2, -3, 1, -6, 11, -6, 1, 24, -50, 35, -10, 1, -120, 274, -225, 85, -15, 1, 720, -1764, 1624, -735, 175, -21, 1, -5040, 13068, -13132, 6769, -1960, 322, -28, 1, 40320, -109584, 118124, -67284, 22449, -4536, 546, -36, 1, -362880, 1026576, -1172700, 723680, -269325, 63273, -9450, 870, -45, 1

Offset: 1

N. J. A. Sloane

The unsigned numbers are also called Stirling cycle numbers: |s(n,k)| = number of permutations of n objects with exactly k cycles.
The unsigned numbers (read from right to left) also give the number of permutations of 1..n with complexity k, where the complexity of a permutation is defined to be the sum of the lengths of the cycles minus the number of cycles. In other words, the complexity equals the sum of (length of cycle)-1 over all cycles. For n=5, the numbers of permutations with complexity 0,1,2,3,4 are 1, 10, 35, 50, 24. - N. J. A. Sloane, Feb 08 2019
The unsigned numbers are also the number of permutations of 1..n with k left to right maxima (see Khovanova and Lewis, Smith).
With P(n) = the number of integer partitions of n, T(i,n) = the number of parts of the i-th partition of n, D(i,n) = the number of different parts of the i-th partition of n, p(j,i,n) = the j-th part of the i-th partition of n, m(j,i,n) = multiplicity of the j-th part of the i-th partition of n, Sum_[T(i,n)=k]{i=1}^{P(n)} = sum running from i=1 to i=p(n) but taking only partitions with T(i,n)=k parts into account, Product{j=1..T(i,n)} = product running from j=1 to j=T(i,n), Product_{j=1..D(i,n)} = product running from j=1 to j=D(i,n) one has S1(n,k) = Sum_[T(i,n)=k]{i=1}^{P(n)} (n!/Product{j=1..T(i,n)} p(j,i,n))* (1/Product_{j=1..D(i,n)} m(j,i,n)!). For example, S1(6,3) = 225 because n=6 has the following partitions with k=3 parts: (114), (123), (222). Their complexions are: (114): (6!/1*1*4)*(1/2!*1!) = 90, (123): (6!/1*2*3)*(1/1!*1!*1!) = 120, (222): (6!/2*2*2)*(1/3!) = 15. The sum of the complexions is 90+120+15 = 225 = S1(6,3). - Thomas Wieder, Aug 04 2005
Row sums equal 0. - Jon Perry, Nov 14 2005
|s(n,k)| enumerates unordered n-vertex forests composed of k increasing non-plane (unordered) trees. Proof from the e.g.f. of the first column and the F. Bergeron et al. reference, especially Table 1, last row (non-plane "recursive"), given in A049029. - Wolfdieter Lang, Oct 12 2007
|s(n,k)| enumerates unordered increasing n-vertex k-forests composed of k unary trees (out-degree r from {0,1}) whose vertices of depth (distance from the root) j >= 0 come in j+1 colors (j=0 for the k roots). - Wolfdieter Lang, Oct 12 2007, Feb 22 2008
A refinement of the unsigned array is A036039. For an association to forests of "naturally grown" rooted non-planar trees, dispositions of flags on flagpoles, and colorings of the vertices of the complete graphs K_n, see A130534. - Tom Copeland, Mar 30 and Apr 05 2014
The Stirling numbers of the first kind were related to the falling factorial and the convolved, or generalized, Bernoulli numbers B_n by Norlund in 1924 through Sum_{k=1..n+1} T(n+1,k) * x^(k-1) = (x-1)!/(x-1-n)! = (x + B.(0))^n = B_n(x), umbrally evaluated with (B.(0))^k = B_k(0) and the associated Appell polynomial B_n(x) defined by the e.g.f. (t/(exp(t) - 1))^(n+1) * exp(x*t) = exp(B.(x)t). - Tom Copeland, Sep 29 2015
With x = e^z, D_x = d/dx, D_z = d/dz, and p_n(x) the row polynomials of this entry, x^n (D_x)^n = p_n(D_z) = (D_z)! / (D_z - n)! = (xD_x)! / (xD_x - n)!. - Tom Copeland, Nov 27 2015
From the operator relation z + Psi(1) + sum_{n > 0} (-1)^n (-1/n) binomial(D,n) = z + Psi(1+D) with D = d/dz and Psi the digamma function, Zeta(n+1) = Sum_{k > n-1} (1/k) |S(k,n)| / k! for n > 0 and Zeta the Riemann zeta function. - Tom Copeland, Aug 12 2016
Let X_1,...,X_n be i.i.d. random variables with exponential distribution having mean = 1.  Let Y = max{X_1,...,X_n}.  Then (-1)^n*n!/( Sum_{k=1..n+1} a(n+1,k) t^(k-1) ) is the moment generating function of Y.  The expectation of Y is the n-th harmonic number. - Geoffrey Critzer, Dec 25 2018
In the Ewens sampling theory describing the multivariate probability distribution of the sizes of the allelic classes in a sample of size n under the Infinite Alleles Model, |s(n,k)| gives the coefficient in the formula for the probability that a sample of n alleles has exactly k distinct types. - Noah A Rosenberg, Feb 10 2019
Named by Nielson (1906) after the Scottish mathematician James Stirling (1692-1770). - Amiram Eldar, Jun 11 2021 and Oct 02 2023
The first few row polynomials along with a recursion formula are found in a manuscript by Newton written in 1664 or 1665 (p. 169 of Turnbull) giving a geometric presentation of the binomial theorem for rational powers. - Tom Copeland, Dec 10 2022

			|s(3,2)| = 3 for the three unordered 2-forest with 3 vertices and two increasing (nonplane) trees: ((1),(2,3)), ((2),(1,3)), ((3),(1,2)).
Triangle begins:
                                      1
                                 -1,      1
                               2,    -3,      1
                          -6,    11,     -6,     1
                      24,    -50,    35,    -10,    1
                -120,    274,  -225,     85,   -15,    1
             720,  -1764,   1624,  -735,    175,  -21,   1
       -5040,  13068, -13132,  6769,  -1960,   322,  -28,  1
  40320, -109584, 118124, -67284, 22449,  -4536,  546, -36,  1
Another version of the same triangle, from _Joerg Arndt_, Oct 05 2009: (Start)
s(n,k) := number of permutations of n elements with exactly k cycles ("Stirling cycle numbers")
  n|  total   m=1      2      3     4     5    6   7  8 9
  -+-----------------------------------------------------
  1|      1     1
  2|      2     1      1
  3|      6     2      3      1
  4|     24     6     11      6     1
  5|    120    24     50     35    10     1
  6|    720   120    274    225    85    15    1
  7|   5040   720   1764   1624   735   175   21   1
  8|  40320  5040  13068  13132  6769  1960  322  28  1
  9| 362880 40320 109584 118124 67284 22449 4536 546 36 1
(End)
|s(4,2)| = 11 for the eleven unordered 2-forest with 4 vertices, composed of two increasing (nonplane) trees: ((1),((23)(24))), ((2),((13)(14))), ((3),((12)(14))), ((4),((12)(13))); ((1),(2,3,4)),((2),(1,2,3)), ((3), (1,2,4)), ((4),(1,2,3)); ((1,2),(3,4)), ((1,3),(2,4)), ((1,4),(2,3)). - _Wolfdieter Lang_, Feb 22 2008

M. Abramowitz and I. A. Stegun, eds., Handbook of Mathematical Functions, National Bureau of Standards Applied Math. Series 55, 1964 (and various reprintings), p. 833.
Arthur T. Benjamin and Jennifer Quinn, Proofs that really count: the art of combinatorial proof, M.A.A. 2003, p. 93ff.
Boris A. Bondarenko, Generalized Pascal Triangles and Pyramids (in Russian), FAN, Tashkent, 1990, ISBN 5-648-00738-8.
George Boole, Finite Differences, 5th ed. New York, NY: Chelsea, 1970.
Louis Comtet, Advanced Combinatorics, Reidel, 1974; Chapter V, also p. 310.
John H. Conway and Richard K. Guy, The Book of Numbers, Copernicus Press, NY, 1996, p. 93.
Florence Nightingale David, Maurice George Kendall and David Elliot Barton, Symmetric Function and Allied Tables, Cambridge, 1966, p. 226.
Saber N. Elaydi, An Introduction to Difference Equations, 3rd ed. Springer, 2005.
Herman H. Goldstine, A History of Numerical Analysis, Springer-Verlag, 1977; Section 2.7.
Ronald L. Graham, Donald E. Knuth and Oren Patashnik, Concrete Mathematics. Addison-Wesley, Reading, MA, 1990, p. 245. In the second edition, see Chapter 6, especially p. 259.
M. Miyata and J. W. Son, On the complexity of permutations and the metric space of bijections, Tensor, 60 (1998), No. 1, 109-116 (MR1768839).
Isaac Newton, A Method whereby to find ye areas of Those Lines wch can be squared, pp. 168-171 of Turnbull below.
John Riordan, An Introduction to Combinatorial Analysis, p. 48.
Robert Sedgewick and Phillipe Flajolet, An Introduction to the Analysis of Algorithms, Addison-Wesley, Reading, MA, 1996.
H. Turnbull (editor), The Correspondence of Isaac Newton Vol. II 1676-1687, Cambridge Univ. Press, 1960.

T. D. Noe, Rows 1 to 100 of triangle, flattened.
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].
Nikita Alexeev and Peter Zograf, Hultman numbers, polygon gluings and matrix integrals, arXiv:1111.3061 [math.PR], 2011.
Joerg Arndt, Matters Computational (The Fxtbook), p. 277.
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.
Elena Barcucci, Alberto Del Lungo and Renzo Pinzani, "Deco" polyominoes, permutations and random generation, Theoretical Computer Science, Vol. 159 (1996), pp. 29-42.
Nicolas Behr, Vincent Danos, and Ilias Garnier, Combinatorial Conversion and Moment Bisimulation for Stochastic Rewriting Systems, arXiv:1904.07313 [cs.LO], 2019.
Hacene Belbachir and Mourad Rahmani, Alternating Sums of the Reciprocals of Binomial Coefficients, Journal of Integer Sequences, Vol. 15 (2012), Article 12.2.8.
Edward A. Bender, Central and local limit theorems applied to asymptotic enumeration Journal of Combinatorial Theory, Series A, Vol. 15, No. 1 (1973), pp. 91-111. See Example 5.1.
Boris A. Bondarenko, Generalized Pascal Triangles and Pyramids, English translation published by Fibonacci Association, Santa Clara Univ., Santa Clara, CA, 1993; see p. 32.
Federico Castillo, Damian de la Fuente, Nicolas Libedinsky, and David Plaza, On the size of Bruhat intervals, arXiv:2309.08539 [math.CO], 2023.
José Luis Cereceda,, Generalized Akiyama-Tanigawa Algorithm for Hypersums of Powers of Integers, J. Int. Seq., Vol. 16 (2013), Article 13.3.2.
José Luis Cereceda, Iterative Procedure for Hypersums of Powers of Integers, Journal of Integer Sequences, Vol. 17 (2014), Article 14.5.3.
Ricky X. F. Chen, A Note on the Generating Function for the Stirling Numbers of the First Kind, Journal of Integer Sequences, Vol. 18 (2015), Article 15.3.8.
Tom Copeland, A Class of Differential Operators and the Stirling Numbers, Generators, Inversion, and Matrix, Binomial, and Integral Transforms, Lagrange a la Lah.
Harry Crane, The ubiquitous Ewens sampling formula, Statistical Science, Vol. 31 (2016), pp. 1-19.
Thierry Dana-Picard and David G. Zeitoun, Sequences of definite integrals, infinite series and Stirling numbers, International Journal of Mathematical Education in Science and Technology, Vol. 43, No. 2 (2012), pp. 219-230.
Sajal K. Das, Joydeep Ghosh and Narsingh Deo, Stirling networks: a versatile combinatorial topology for multiprocessor systems, Discrete Applied Mathematics, Vol. 37 (1992), pp. 119-146.
Robert M. Dickau, Stirling numbers of the first kind.
Askar Dzhumadil’daev and Damir Yeliussizov, Walks, partitions, and normal ordering, Electronic Journal of Combinatorics, Vol. 22, No. 4 (2015), #P4.10.
B. S. El-Desouky, Nenad P. Cakić and Toufik Mansour, Modified approach to generalized Stirling numbers via differential operators, Appl. Math. Lett., Vol. 23, No. 1 (2010), pp. 115-120.
FindStat - Combinatorial Statistic Finder, The number of saliances of a permutation, The number of cycles in the cycle decomposition of a permutation.
Bill Gosper, Colored illustrations of triangle of Stirling numbers of first kind read mod 2, 3, 4, 5, 6, 7.
Daniel B. Grünberg, On asymptotics, Stirling numbers, gamma function and polylogs, Results in Mathematics, Vol. 49, No. 1 (2006), pp. 89-125; arXiv preprint, arXiv:math/0607514 [math.CO], 2006.
Jerome Hines, A generalization of the S-Stirling numbers, Math. Mag., Vol. 29, No. 4 (1956), pp. 200-203.
Yoshinari Inaba, Hyper-Sums of Powers of Integers and the Akiyama-Tanigawa Matrix, Journal of Integer Sequences, 8 (2005), Article 05.2.7.
Charles Knessl and Joseph B. Keller, Stirling number asymptotics from recursion equations using the ray method, Stud. Appl. Math., Vol. 84, No. 1 (1991), pp. 43-56.
Tanya Khovanova and Joel Brewster Lewis, Skyscrapers.
Donald E. Knuth, Convolution polynomials, arXiv:math/9207221 [math.CA], 1993; The Mathematica J., 2 (1992), 67-78.
Wolfdieter Lang, On generalizations of Stirling number triangles, J. Integer Seqs., Vol. 3 (2000), Article 00.2.4.
Elliott H. Lieb, Concavity properties and a generating function for Stirling numbers, Journal of Combinatorial Theory, Vol. 5, No. 2 (1968), pp. 203-206.
Daniel E. Loeb, A generalization of Stirling numbers, arXiv:math/9502217 [math.CO], 1995.
Mahid M. Mangontarum and Jacob Katriel, On q-Boson Operators and q-Analogues of the r-Whitney and r-Dowling Numbers, J. Int. Seq., Vol. 18 (2015), Article 15.9.8.
Toufik Mansour, Augustine Munagi and Mark Shattuck, Recurrence Relations and Two-Dimensional Set Partitions, J. Int. Seq., Vol. 14 (2011), Article 11.4.1.
Toufik Mansour, Matthias Schork and Mark Shattuck, The Generalized Stirling and Bell Numbers Revisited, Journal of Integer Sequences, Vol. 15 (2012), Article 12.8.3.
B. H. Margolius, Transient and periodic solution to the time-inhomogeneous quasi-birth death process, Queueing Systems, Volume 56, Numbers 3-4 / August, 2007. [From _N. J. A. Sloane_, Jul 09 2009]
Dragoslav S. Mitrinović, Sur une relation de récurrence relative aux nombres de Bernoulli, Comptes rendus de l'Académie des sciences de Paris, Vol. 250 (1960), pp. 4266-4267.
Dragoslav S. Mitrinović and Ružica S. Mitrinović, Sur les nombres de Stirling et les nombres de Bernoulli d'ordre supérieur, Univ. Beograd. Publ. Elektrotehn. Fak. Ser. Mat. Fiz., No. 43 (1960), pp. 1-63.
Dragoslav S. Mitrinović and Ružica S. Mitrinović, Sur une classe de nombres se rattachant aux nombres de Stirling--Appendice: Table des nombres de Stirling, Univ. Beograd. Publ. Elektrotehn. Fak. Ser. Mat. Fiz., No. 60 (1961), pp. 1-15 and 17-62.
Dragoslav S. Mitrinović and Ružica S. Mitrinović, Tableaux d'une classe de nombres reliés aux nombres de Stirling, Univ. Beograd. Publ. Elektrotehn. Fak. Ser. Mat. Fiz., No. 77 (1962), pp. 1-77.
Dragoslav S. Mitrinović and Ružica S. Mitrinović, Tableaux d'une classe de nombres reliés aux nombres de Stirling, Univ. Beograd. Publ. Elektrotehn. Fak. Ser. Mat. Fiz., No. 77 (1962), pp. 1-77 [jstor stable version].
Emanuele Munarini, Combinatorial identities involving the central coefficients of a Sheffer matrix, Applicable Analysis and Discrete Mathematics, Vol. 13 (2019), pp. 495-517.
Niels Nielson, Handbuch der theorie der gammafunktion, Teubner, Leipzig, 1906, see p. 67.
Niels Nörlund, Vorlesungen über Differenzenrechnung, Springer, Berlin, 1924.
Niels Nörlund, Vorlesungen über Differenzenrechnung, Chelsea Pub. Co., New York, 1954. [archive.org]
K. A. Penson, P. Blasiak, A. Horzela, G. H. E. Duchamp and A. I. Solomon, Laguerre-type derivatives: Dobinski relations and combinatorial identities, J. Math. Phys. Vol. 50 (2009), 083512.
Grzegorz Rządkowski, Two formulas for Successive Derivatives and Their Applications, J. Int. Seq., Vol. 12 (2009), Article 09.8.2.
Maxie D. Schmidt, Generalized j-Factorial Functions, Polynomials, and Applications, J. Int. Seq., Vol. 13 (2010), Article 10.6.7.
Benjamin Schreyer, Rigged Horse Numbers and their Modular Periodicity, arXiv:2409.03799 [math.CO], 2024. See p. 12.
Raymond Scurr and Gloria Olive, Stirling numbers revisited, Discrete Math., Vol. 189, No. 1-3 (1998), pp. 209-219. MR1637761 (99d:11019).
Mark Shattuck, Convolution identities for Stirling numbers of the first kind via involution, INTEGERS, Vol. 12 (2012), #A59. - From _N. J. A. Sloane_, Feb 04 2013
Mark Shattuck, Combinatorial Proofs of Some Stirling Number Convolution Formulas, J. Int. Seq., Vol. 25 (2022), Article 22.2.2.
Warren D. Smith, Puzzle: Counting new records.
Michael Z. Spivey, A generalized recurrence for Bell Numbers, JIS, Vol. 11 (2008), Article 08.2.5.
Michael Z. Spivey, On Solutions to a General Combinatorial Recurrence, J. Int. Seq., Vol. 14 (2011), Article 11.9.7.
Richard P. Stanley, Ordering events in Minkowski space, arXiv:math/0501256 [math.CO], 2005.
James Stirling, The Differential Method, London, 1749; see p. 10.
Nico M. Temme, Asymptotic Estimates of Stirling Numbers, Studies in Applied Mathematics, Vol. 89, No. 3 (1993), pp. 233-243; Summary by Helmut Prodinger.
A. N. Timashev, On asymptotic expansions of Stirling numbers of the first and second kinds, (Russian) Diskret. Mat. 10(3) (1998), 148-159; translation in Discrete Math. Appl., Vol. 8, No. 5 (1998), pp. 533-544.
Eric Weisstein's World of Mathematics, Permutation Cycle.
Eric Weisstein's World of Mathematics, Stirling Number of the First Kind.
Thomas Wieder, Comments on A008275.
OEIS Wiki, Factorial polynomials.
Index entries for "core" sequences

Diagonals: A000217, A000914, A001303, A000915, A053567, etc.
Cf. A048994, A008277 (Stirling numbers of second kind), A039814, A039815, A039816, A039817, A048993, A087748.
Cf. A084938, A094216, A008276 (row reversed), A008277, A008278, A094262, A121632, A130534 (unsigned version), A087755 (triangle mod 2), A000142 (row sums of absolute values).

a008275 n k = a008275_tabl !! (n-1) !! (k-1)
a008275_row n = a008275_tabl !! (n-1)
a008275_tabl = map tail $ tail a048994_tabl
-- Reinhard Zumkeller, Mar 18 2013

with (combinat):seq(seq(stirling1(n, k), k=1..n), n=1..10); # Zerinvary Lajos, Jun 03 2007
for i from 0 to 9 do seq(stirling1(i, j), j = 1 .. i) od; # Zerinvary Lajos, Nov 29 2007

Flatten[Table[StirlingS1[n, k], {n, 1, 10}, {k, 1, n}]] (* Jean-François Alcover, May 18 2011 *)
Flatten@Table[Coefficient[Product[x-k, {k, 0, n-1}], x, Range[n]], {n, Range[10]}] (* Oliver Seipel, Jun 11 2024 *)
a[n_, n_] := 1; a[n_, 0] := 0; a[0, k_] := 0;
a[n_, k_] := a[n, k] = a[n-1, k-1] + (n-1) a[n-1, k];
Flatten@Table[(-1)^(n-k) a[n, k], {n, 1, 10}, {k, 1, n}] (* Oliver Seipel, Jun 11 2024 *)

create_list(stirling1(n+1,k+1),n,0,30,k,0,n); /* Emanuele Munarini, Jun 01 2012 */

T(n,k)=if(n<1,0,n!*polcoeff(binomial(x,n),k))

T(n,k)=if(n<1,0,n!*polcoeff(polcoeff((1+x+x*O(x^n))^y,n),k))

vecstirling(n)=Vec(factorback(vector(n-1,i,1-i*'x))) /* (A function that returns all the s(n,k) as a vector) */ \\ Bill Allombert (Bill.Allombert(AT)math.u-bordeaux1.fr), Mar 16 2009

s(n, k) = s(n-1, k-1) - (n-1)*s(n-1, k), n, k >= 1; s(n, 0) = s(0, k) = 0; s(0, 0) = 1.
The unsigned numbers a(n, k)=|s(n, k)| satisfy a(n, k) = a(n-1, k-1) + (n-1)*a(n-1, k), n, k >= 1; a(n, 0) = a(0, k) = 0; a(0, 0) = 1.
E.g.f.: for m-th column (unsigned): ((-log(1-x))^m)/m!.
s(n, k) = T(n-1, k-1), n>1 and k>1, where T(n, k) is the triangle [ -1, -1, -2, -2, -3, -3, -4, -4, -5, -5, -6, -6, ...] DELTA [1, 0, 1, 0, 1, 0, 1, 0, 1, ...] and DELTA is Deléham's operator defined in A084938. The unsigned numbers are also |s(n, k)| = T(n-1, k-1), for n>0 and k>0, where T(n, k) = [1, 1, 2, 2, 3, 3, 4, 4, 5, 5, ...] DELTA [1, 0, 1, 0, 1, 0, 1, 0, ...].
Sum_{i=0..n} (-1)^(n-i) * StirlingS1(n, i) * binomial(i, k) = (-1)^(n-k) *  StirlingS1(n+1, k+1). - Carlo Wood (carlo(AT)alinoe.com), Feb 13 2007
G.f. for row n: Product_{j=1..n} (x-j) (e.g., (x-1)*(x-2)*(x-3) = x^3 - 6*x^2 + 11*x - 6). - Jon Perry, Nov 14 2005
s(n,k) = A048994(n,k), for k=1..n. - Reinhard Zumkeller, Mar 18 2013 (Corrected by N. J. A. Sloane, May 07 2025 at the suggestion of Manfred Boergens, May 07 2025)
As lower triangular matrices A008277*A008275 = I, the identity matrix. - Tom Copeland, Apr 25 2014
a(n,k) = s(n,k) = lim_{y -> 0} Sum_{j=0..k} (-1)^j*binomial(k,j)*((-j*y)!/(-j*y-n)!)*y^(-k)/k! = Sum_{j=0..k} (-1)^(n-j)*binomial(k,j)*((j*y - 1 + n)!/(j*y-1)!)*y^(-k)/k!. - Tom Copeland, Aug 28 2015
From Daniel Forgues Jan 16 2016: (Start)
Let x_(0) := 1 (empty product), and for n >= 1:
  x_(n) := Product_{k=0..n-1} (x-k), called a factorial term (Boole, 1970) or a factorial polynomial (Elaydi, 2005: p. 60), and also x_(-n) := 1 / [Product_{k=0..n-1} (x+k)].
Then, for n >= 1: x_(n) = Sum_{k=1..n} T(n,k) * x^k, 1 / [x_(-n)] = Sum_{k=1..n} |T(n,k)| * x^k, x^n = Sum_{k=1..n} A008277(n,k) * x_(k), where A008277(n,k) are Stirling numbers of the second kind.
The row sums (of either signed or absolute values) are Sum_{k=1..n} T(n,k) = 0^(n-1), Sum_{k=1..n} |T(n,k)| = T(n+1,1) = n!. (End)
s(n,m) = ((-1)^(n-m)/n)*Sum_{i=0..m-1} C(2*n-m-i, m-i-1)*A008517(n-m+1,n-m-i+1). - Vladimir Kruchinin, Feb 14 2018
Orthogonal relation: Sum_{i=0..n} i^p*Sum_{j=k..n} (-1)^(i+j) * binomial(j,i) * Stirling1(j,k)/j! = delta(p,k), i,k,p <= n, n >= 1. - Leonid Bedratyuk, Jul 27 2020
From Zizheng Fang, Dec 28 2020: (Start)
Sum_{k=1..n} (-1)^k * k * T(n, k) = -T(n+1, 2).
Sum_{k=1..n} k * T(n, k) = (-1)^n * (n-2)! = T(n-1, 1) for n>=2. (End)
n-th row polynomial = n!*Sum_{k = 0..2*n} (-1)^(n+k)*binomial(x, k)*binomial(x-1, 2*n-k) = n!*Sum_{k = 0..2*n+1} (-1)^(n+k+1)*binomial(x, k)*binomial(x-1, 2*n+1-k). - Peter Bala, Mar 29 2024

A211795 Number of (w,x,y,z) with all terms in {1,...,n} and wx < 2y*z.

Original entry on oeis.org

0, 1, 11, 58, 177, 437, 894, 1659, 2813, 4502, 6836, 10008, 14121, 19449, 26117, 34372, 44422, 56597, 71044, 88160, 108115, 131328, 158074, 188773, 223604, 263172, 307719, 357715, 413493, 475690, 544480, 620632, 704381, 796413

Offset: 0

Clark Kimberling, Apr 27 2012

nonn

Each sequence in the following guide counts 4-tuples
(w,x,y,z) such that the indicated relation holds and the four numbers w,x,y,z are in {1,...,n}.  The notation "m div" means that m divides every term of the sequence.
A211058 ... wx <= yz
A211787 ... wx <= 2yz
A211795 ... wx < 2yz
A211797 ... wx > 2yz
A211809 ... wx >= 2yz
A211812 ... wx <= 3yz
A211917 ... wx < 3yz
A211918 ... wx > 3yz
A211919 ... wx >= 3yz
A211920 ... 2wx < 3yz
A211921 ... 2wx <= 3yz
A211922 ... 2wx > 3yz
A211923 ... 2wx >= 3yz
A212019 ... wx = 2yz ..... 2 div
A212020 ... wx = 3yz ..... 2 div
A212021 ... 2wx = 3yz .... 2 div
A212047 ... wx = 4yz
A212048 ... 3wx = 4yz .... 2 div
A212049 ... wx = 5yz ..... 2 div
A212050 ... 2wx = 5yz .... 2 div
A212051 ... 3wx = 5yz .... 2 div
A212052 ... 4wx = 5yz .... 2 div
A209978 ... wx = yz + 1 .. 2 div
A212053 ... wx <= yz + 1
A212054 ... wx > yz + 1
A212055 ... wx <= yz + 2
A212056 ... wx > yz + 2
A197168 ... wx = yz + 2 .. 2 div
A061201 ... w = xyz
A212057 ... w < xyz
A212058 ... w >= xyz
A212059 ... w = xyz - 1
A212060 ... w = xyz - 2
A212061 ... wx = (yz)^2
A212062 ... w^2 = xyz
A212063 ... w^2 < xyz
A212064 ... w^2 >= xyz
A212065 ... w^2 <= xyz
A212066 ... w^2 > xyz
A212067 ... w^3 = xyz
A002623 ... w = 2x + y + z
A006918 ... w = 2x + 2y + z
A000601 ... w = x + 2y + 3z (except for initial 0's)
A212068 ... 2w = x + y + z
A212069 ... 3w = x + y + z (w = average{x,y,z})
A212088 ... 3w < x + y + z
A212089 ... 3w >= x + y + z
A212090 ... w < x + y + z
A000332 ... w >= x + y + z
A212145 ... w < 2x + y + z
A001752 ... w >= 2x + y + z
A001400 ... w = 2x +3y + 4z
A005900 ... w = -x + y + z
A192023 ... w = -x + y + z + 2
A212091 ... w^2 = x^2 + y^2 + z^2 ... 3 div
A212087 ... w^2 + x^2 = y^2 + z^2
A212092 ... w^2 < x^2 + y^2 + z^2
A212093 ... w^2 <= x^2 + y^2 + z^2
A212094 ... w^2 > x^2 + y^2 + z^2
A212095 ... w^2 >= x^2 + y^2 + z^2
A212096 ... w^3 = x^3 + y^3 + z^3 ... 6 div
A212097 ... w^3 < x^3 + y^3 + z^3
A212098 ... w^3 <= x^3 + y^3 + z^3
A212099 ... w^3 > x^3 + y^3 + z^3
A212100 ... w^3 >= x^3 + y^3 + z^3
A212101 ... wx^2 = yz^2
A212102 ... 1/w = 1/x + 1/y + 1/z
A212103 ... 3/w = 1/x + 1/y + 1/z; w = h.m. of {x,y,z}
A212104 ... 3/w >= 1/x + 1/y + 1/z; w >= h.m.
A212105 ... 3/w < 1/x + 1/y + 1/z; w < h.m.
A212106 ... 3/w > 1/x + 1/y + 1/z; w > h.m.
A212107 ... 3/w <= 1/x + 1/y + 1/z; w <= h.m.
A212133 ... median(w,x,y,z) = mean(w,x,y,z)
A212134 ... median(w,x,y,z) <= mean(w,x,y,z)
A212135 ... median(w,x,y,z) > mean(w,x,y,z)
A212241 ... wx + yz > n
A212243 ... 2wx + yz = n
A212244 ... w = xyz - n
A212245 ... w = xyz - 2n
A212246 ... 2w = x + y + z - n
A212247 ... 3w = x + y + z + n
A212249 ... 3w < x + y + z + n
A212250 ... 3w >= x + y + z + n
A212251 ... 3w = x + y + z + n + 1
A212252 ... 3w = x + y + z + n + 2
A212254 ... w = x + 2y + 3z - n
A212255 ... w^2 = mean(x^2, y^2, z^2)
A212256 ... 4/w = 1/x + 1/y +1/z + 1/n
In the list above, if the relation in the second column is of the form "w rel ax + by + cz" then the sequence is linearly recurrent.  In the list below, the same is true for expressions involving more than one relation.
A000332 ... w < x <= y < z .... C(n,4)
A000914 ... w < x <= y < z .... Stirling 1st kind
A000914 ... w < x <= y >= z ... Stirling 1st kind
A050534 ... w < x < y >= z .... tritriangular
A001296 ... w <= x <= y >= z .. 4-dim pyramidal
A006322 ... x < x > y >= z
A002418 ... w < x >= y < z
A050534 ... w < x >=y >= z
A212415 ... w < x >= y <= z
A001296 ... w < x >= y <= z
A212246 ... w <= x > y <= z
A006322 ... w <= x >= y <= z
A212501 ... w > x < y >= z
A212503 ... w < 2x and y < 2z ..... A (note below)
A212504 ... w < 2x and y > 2z ..... A
A212505 ... w < 2x and y >= 2z .... A
A212506 ... w <= 2x and y <= 2z ... A
A212507 ... w < 2x and y <= 2z .... B
A212508 ... w < 2x and y < 3z ..... C
A212509 ... w < 2x and y <= 3z .... C
A212510 ... w < 2x and y > 3z ..... C
A212511 ... w < 2x and y >= 3z .... C
A212512 ... w <= 2x and y < 3z .... C
A212513 ... w <= 2x and y <= 3z ... C
A212514 ... w <= 2x and y > 3z .... C
A212515 ... w <= 2x and y >= 3z ... C
A212516 ... w > 2x and y < 3z ..... C
A212517 ... w > 2x and y <= 3z .... C
A212518 ... w > 2x and y > 3z ..... C
A212519 ... w > 2x and y >= 3z .... C
A212520 ... w >= 2x and y < 3z .... C
A212521 ... w >= 2x and y <= 3z ... C
A212522 ... w >= 2x and y > 3z .... C
A212523 ... w + x < y + z
A212560 ... w + x <= y + z
A212561 ... w + x = 2y + 2z
A212562 ... w + x < 2y + 2z ....... B
A212563 ... w + x <= 2y + 2z ...... B
A212564 ... w + x > 2y + 2z ....... B
A212565 ... w + x >= 2y + 2z ...... B
A212566 ... w + x = 3y + 3z
A212567 ... 2w + 2x = 3y + 3z
A212570 ... |w - x| = |x - y| + |y - z|
A212571 ... |w - x| < |x - y| + |y - z| ... B ... 4 div
A212572 ... |w - x| <= |x - y| + |y - z| .. B
A212573 ... |w - x| > |x - y| + |y - z| ... B ... 2 div
A212574 ... |w - x| >= |x - y| + |y - z| .. B
A212575 ... 2|w - x| = |x - y| + |y - z|
A212576 ... |w - x| = 2|x - y| + 2|y - z|
A212577 ... |w - x| = 2|x - y| - |y - z|
A212578 ... 2|w - x| = |x - y| - |y - z|
A212579 ... min{|w-x|,|w-y|} = min{|x-y|,|x-z|}
A212692 ... w = |x - y| + |y - z| ............... 2 div
A212568 ... w < |x - y| + |y - z| ............... 2 div
A212573 ... w <= |x - y| + |y - z| .............. 2 div
A212574 ... w > |x - y| + |y - z|
A212575 ... w >= |x - y| + |y - z|
A212676 ... w + x = |x - y| + |y - z| ......... H
A212677 ... w + y = |x - y| + |y - z|
A212678 ... w + x + y = |x - y| + |y - z|
A006918 ... w + x + y + z = |x - y| + |y - z| . H
A212679 ... |x - y| = |y - z| ................. H
A212680 ... |x - y| = |y - z| + 1 ..............H 2 div
A212681 ... |x - y| < |y - z| ................... 2 div
A212682 ... |x - y| >= |y - z|
A212683 ... |x - y| = w + |y - z| ............... 2 div
A212684 ... |x - y| = n - w + |y - z|
A212685 ... |w - x| = w + |y - z|
A186707 ... |w - x| < w + |y - z| ... (Note D)
A212714 ... |w - x| >= w + |y - z| .......... H . 2 div
A212686 ... 2*|w - x| = n + |y - z| ............. 4 div
A212687 ... 2*|w - x| < n + |y - z| ......... B
A212688 ... 2*|w - x| < n + |y - z| ......... B . 2 div
A212689 ... 2*|w - x| > n + |y - z| ......... B . 2 div
A212690 ... 2*|w - x| <= n + |y - z| ........ B
A212691 ... w + |x - y| = |x - z| + |y - z| . E . 2 div
...
In the above lists, all the terms of (w,x,y,z) are in {1,...,n}, but in the next lists they are all in {0,...,n}, and sequences are all linearly recurrent.
  R=range{w,x,y,z}=max{w,x,y,z}-min{w,x,y,z}.
A212740 ... max{w,x,y,z} < 2*min{w,x,y,z} .... A
A212741 ... max{w,x,y,z} >= 2*min{w,x,y,z} ... A
A212742 ... max{w,x,y,z} <= 2*min{w,x,y,z} ... A
A212743 ... max{w,x,y,z} > 2*min{w,x,y,z} .... A . 2 div
A212744 ... w=range (=max-min) ............... E
A212745 ... w=max{w,x,y,z} - 2*min{w,x,y,z}
A212746 ... R is in {w,x,y,z} ................ E
A212569 ... R is not in {w,x,y,z} ............ E
A212749 ... w=R or xA212750 ... w=R or x=R or yA212751 ... w=R or x=R or yA212752 ... wR ......... A
A212753 ... wR or z>R ......... D
A212754 ... wR or y>R or z>R ......... D
A002415 ... w = x + R ........................ D
A212755 ... |w - x| = R ...................... D
A212756 ... 2w = x + R
A212757 ... 2w = R
A212758 ... w = floor(R/2)
A002413 ... w = floor((x+y+z/2))
A212759 ... w, x, y are even
A212760 ... w is even and x = y + z .......... E
A212761 ... w is odd and x and y are even .... F . 2 div
A212762 ... w and x are odd y is even ........ F . 2 div
A212763 ... w, x, y are odd .................. F
A212764 ... w, x, y are even and z is odd .... F
A030179 ... w and x are even and y and z odd
A212765 ... w is even and x,y,z are odd ...... F
A212766 ... w is even and x is odd ........... A . 2 div
A212767 ... w and x are even and w+x=y+z ..... E
A212889 ... R is even ........................ A
A212890 ... R is odd ......................... A . 2 div
A212742 ... w-x, x-y, y-z are all even ....... A
A212892 ... w-x, x-y, y-z are all odd ........ A
A212893 ... w-x, x-y, y-z have same parity ... A
A005915 ... min{|w-x|, |x-y|, |y-z|} = 0
A212894 ... min{|w-x|, |x-y|, |y-z|} = 1
A212895 ... min{|w-x|, |x-y|, |y-z|} = 2
A179824 ... min{|w-x|, |x-y|, |y-z|} > 0
A212896 ... min{|w-x|, |x-y|, |y-z|} <= 1
A212897 ... min{|w-x|, |x-y|, |y-z|} > 1
A212898 ... min{|w-x|, |x-y|, |y-z|} <= 2
A212899 ... min{|w-x|, |x-y|, |y-z|} > 2
A212901 ... |w-x| = |x-y| = |y-z|
A212900 ... |w-x|, |x-y|, |y-z| are distinct . G
A212902 ... |w-x| < |x-y| < |y-z| ............ G
A212903 ... |w-x| <= |x-y| <= |y-z| .......... G
A212904 ... |w-x| + |x-y| + |y-z| = n ........ H
A212905 ... |w-x| + |x-y| + |y-z| = 2n ....... H
...
Note A: A212503-A212506 (and others) have these recurrence coefficients: 2,2,-6,0,6,-2,-2,1.
B: 3,-1,-5,5,1,-3,1
C: 0,2,2,-1,-4,0,2,0,-2,0,4,1,-2,-2,0,1
D: 4,-5,0,5,-4,1
E: 1,3,-3,-3,3,1,-1
F: 1,4,-4,-6,6,4,-4,-1,1
G: 2,1,-3,-1,1,3,-1,-2,1
H: 2,1,-4,1,2,-1

		Examples
					a(2)=11 counts these (w,x,y,z): (1,1,1,1), (1,1,1,2), (1,1,2,1), (2,1,2,1), (2,1,1,2), (1,2,2,1), (1,2,1,2), (1,1,2,2), (1,2,2,2), (2,1,2,2), (2,2,2,2).
		
	
References
			A. Barvinok, Lattice Points and Lattice Polytopes, Chapter 7 in Handbook of Discrete and Computational Geometry, CRC Press, 1997, 133-152.
P. Gritzmann and J. M. Wills, Lattice Points, Chapter 3.2 in Handbook of Convex Geometry, vol. B, North-Holland, 1993, 765-797.

		
Links
			Bo Gyu Jeong, Table of n, a(n) for n = 0..200

		
CrossrefsCf. A000583 (n^4), A210000, A211809, A212959.
Programs
			
				Mathematica
				t = Compile[{{n, _Integer}}, Module[{s = 0},
    (Do[If[w*x < 2 y*z, s = s + 1], {w, 1, #},
      {x, 1, #}, {y, 1, #}, {z, 1, #}] &[n]; s)]];
Map[t[#] &, Range[0, 40]] (* A211795 *)
(* Peter J. C. Moses, Apr 13 2012 *)
				

		
Formulaa(n) = n^4 - A211809(n).

A264428 Triangle read by rows, Bell transform of Bell numbers.

Original entry on oeis.org

1, 0, 1, 0, 1, 1, 0, 2, 3, 1, 0, 5, 11, 6, 1, 0, 15, 45, 35, 10, 1, 0, 52, 205, 210, 85, 15, 1, 0, 203, 1029, 1330, 700, 175, 21, 1, 0, 877, 5635, 8946, 5845, 1890, 322, 28, 1, 0, 4140, 33387, 63917, 50358, 20055, 4410, 546, 36, 1, 0, 21147, 212535, 484140, 450905, 214515, 57855, 9240, 870, 45, 1

Offset: 0

Peter Luschny, Nov 13 2015

Consider the sequence S0 -> T0 -> S1 -> T1 -> S2 -> T2 -> ... Here Sn -> Tn indicates the Bell transform mapping a sequence Sn to a triangle Tn as defined in the link and Tn -> S{n+1} the operator associating a triangle with the sequence of its row sums. If
S0 = A000012 = <1,1,1,...> then
T0 = A048993 # Stirling subset numbers,
S1 = A000110 # Bell numbers,
T1 = A264428 # Bell transform of Bell numbers,
S2 = A187761 # second-order Bell numbers,
T2 = A264430 # Bell transform of second-order Bell numbers,
S3 = A264432 # third-order Bell numbers.
This construction is closely related to permutations trees and A179455. Sn is A179455_col(n+1) prepended by A179455_diag(k) = k! for k <= n. In other words, Sn 'converges' to n! for n -> oo.
Given a sequence (s(n))n>=0 with s(0) = 0 and with e.g.f. B(x) = Sum_{n >= 1} s(n)*x^n/n!, then the Bell matrix associated with s(n) equals the exponential Riordan array [1, B(x)] belonging to the Lagrange subgroup of the exponential Riordan group. Omitting the first row and column from the Bell matrix produces the exponential Riordan array [d/dx(B(x)), B(x)] belonging to the Derivative subgroup of the exponential Riordan group. - Peter Bala, Jun 07 2016

			Triangle starts:
[1]
[0,   1]
[0,   1,    1]
[0,   2,    3,    1]
[0,   5,   11,    6,    1]
[0,  15,   45,   35,   10,    1]
[0,  52,  205,  210,   85,   15,   1]
[0, 203, 1029, 1330,  700,  175,  21,  1]
[0, 877, 5635, 8946, 5845, 1890, 322, 28, 1]

G. C. Greubel, Table of n, a(n) for n = 0..1325
Peter Luschny, The Bell transform
Peter Luschny, Permutation Trees

Cf. A000012, A000110, A000217, A000914, A027801, A048993, A051836, A179455, A187761 (row sums), A264430, A264432, A265312.

# Computes sequence in matrix form.
BellMatrix := proc(f, len) local T, A; A := [seq(f(n), n=0..len-2)];
T := proc(n, k) option remember; if k=0 then k^n else
add(binomial(n-1,j-1)*T(n-j,k-1)*A[j], j=1..n-k+1) fi end;
Matrix(len, (n,k)->T(n-1,k-1), shape=triangular[lower]) end:
BellMatrix(n -> combinat:-bell(n), 9); # Peter Luschny, Jan 21 2016
# Alternative, using the recurrence of Peter Bala:
R := proc(n) option remember; if n = 0 then 1 else
t*add(binomial(n-1,k)*combinat:-bell(k)*R(n-k-1,t),k=0..n-1) fi end:
T_row := n-> seq(coeff(R(n), t, k), k=0..n):
seq(print(T_row(n)),n=0..8); # Peter Luschny, Jun 09 2016

BellMatrix[f_Function|f_Symbol, len_] := With[{t = Array[f, len, 0]}, Table[BellY[n, k, t], {n, 0, len-1}, {k, 0, len-1}]];
rows = 11;
M = BellMatrix[BellB, rows];
Table[M[[n, k]], {n, 1, rows}, {k, 1, n}] // Flatten (* Jean-François Alcover, Jan 21 2016, updated Jul 14 2018 *)
With[{r = 8}, Flatten[Table[BellY[n, k, BellB[Range[0, r]]], {n, 0, r}, {k, 0, n}]]] (* Jan Mangaldan, May 22 2016 *)

bell_matrix(f, len) = { my( m = matrix(len, len) );  m[1, 1] = 1;
  for( n = 1, len-1, m[n+1, 2] = f(n-1) );
  for( n = 0, len-1, for( k = 1, n,
     m[n+1, k+1] = sum(j = 1, n-k+1, binomial(n-1,j-1)*m[n-j+1,k]*m[j+1,2]) ));
  return( m )
}
f(n) = polcoeff( sum( k=0, n, prod( i=1, k, x / (1 - i*x)), x^n * O(x)), n);
bell_matrix(f, 9) \\ Peter Luschny, Jan 24 2016

from functools import cache
from math import comb as binomial
def BellMatrix(f, size):
    A = [f(n) for n in range(size - 1)]
    @cache
    def T(n, k):
        if k == 0: return k ** n
        return sum(
            binomial(n - 1, j) * T(n - j - 1, k - 1) * A[j]
            for j in range(n - k + 1) )
    return [[T(n, k) for k in range(n + 1)] for n in range(size)]
@cache
def b(n, k=0): return n < 1 or k*b(n-1, k) + b(n-1, k+1)
print(BellMatrix(b, 9))  # Peter Luschny, Jun 14 2022

# The functions below are referenced in various other sequences.
def bell_transform(n, a): # partition_based
    row = []
    fn = factorial(n)
    for k in (0..n):
        result = 0
        for p in Partitions(n, length=k):
            factorial_product = 1
            power_factorial_product = 1
            for part, count in p.to_exp_dict().items():
                factorial_product *= factorial(count)
                power_factorial_product *= factorial(part)**count
            coefficient = fn//(factorial_product*power_factorial_product)
            result += coefficient*prod([a[i-1] for i in p])
        row.append(result)
    return row
def bell_matrix(generator, dim):
    G = [generator(k) for k in srange(dim)]
    row = lambda n: bell_transform(n, G)
    return matrix(ZZ, [row(n)+[0]*(dim-n-1) for n in srange(dim)])
def inverse_bell_matrix(generator, dim):
    G = [generator(k) for k in srange(dim)]
    row = lambda n: bell_transform(n, G)
    M = matrix(ZZ, [row(n)+[0]*(dim-n-1) for n in srange(dim)]).inverse()
    return matrix(ZZ, dim, lambda n,k: (-1)^(n-k)*M[n,k])
bell_numbers = [sum(bell_transform(n, [1]*10)) for n in range(11)]
for n in range(11): print(bell_transform(n, bell_numbers))

From Peter Bala, Jun 07 2016: (Start)
E.g.f.: exp(t*B(x)), where B(x) = Integral_{u = 0..x} exp(exp(u) - 1) du =  x + x^2/2! + 2*x^3/3! + 5*x^4/4! + 15*x^5/5! + 52*x^6/6! + ....
Row polynomial recurrence: R(n+1,t) = t*Sum_{k = 0 ..n} binomial(n,k)*Bell(k)* R(n-k,t) with R(0,t) = 1. (End)

A130534 Triangle T(n,k), 0 <= k <= n, read by rows, giving coefficients of the polynomial (x+1)(x+2)...(x+n), expanded in increasing powers of x. T(n,k) is also the unsigned Stirling number |s(n+1, k+1)|, denoting the number of permutations on n+1 elements that contain exactly k+1 cycles.

Original entry on oeis.org

1, 1, 1, 2, 3, 1, 6, 11, 6, 1, 24, 50, 35, 10, 1, 120, 274, 225, 85, 15, 1, 720, 1764, 1624, 735, 175, 21, 1, 5040, 13068, 13132, 6769, 1960, 322, 28, 1, 40320, 109584, 118124, 67284, 22449, 4536, 546, 36, 1, 362880, 1026576, 1172700, 723680, 269325, 63273, 9450, 870, 45, 1

Offset: 0

Philippe Deléham, Aug 09 2007

This triangle is an unsigned version of the triangle of Stirling numbers of the first kind, A008275, which is the main entry for these numbers. - N. J. A. Sloane, Jan 25 2011
Or, triangle T(n,k), 0 <= k <= n, read by rows given by [1,1,2,2,3,3,4,4,5,5,6,6,...] DELTA [1,0,1,0,1,0,1,0,1,0,1,0,...] where DELTA is the operator defined in A084938.
Reversal of A094638.
Equals A132393*A007318, as infinite lower triangular matrices. - Philippe Deléham, Nov 13 2007
From Johannes W. Meijer, Oct 07 2009: (Start)
The higher order exponential integrals E(x,m,n) are defined in A163931. The asymptotic expansion of the exponential integrals E(x,m=1,n) ~ (exp(-x)/x)*(1 - n/x + n*(n+1)/x^2 - n*(n+1)*(n+2)/x^3 + ...), see Abramowitz and Stegun. This formula follows from the general formula for the asymptotic expansion, see A163932. We rewrite E(x,m=1,n) ~ (exp(-x)/x)*(1 - n/x + (n^2+n)/x^2 - (2*n+3*n^2+n^3)/x^3 + (6*n+11*n^2+6*n^3+n^4)/x^3 - ...) and observe that the T(n,m) are the polynomials coefficients in the denominators. Looking at the a(n,m) formula of A028421, A163932 and A163934, and shifting the offset given above to 1, we can write T(n-1,m-1) = a(n,m) = (-1)^(n+m)*Stirling1(n,m), see the Maple program.
The asymptotic expansion leads for values of n from one to eleven to known sequences, see the cross-references. With these sequences one can form the triangles A008279 (right-hand columns) and A094587 (left-hand columns).
See A163936 for information about the o.g.f.s. of the right-hand columns of this triangle.
(End)
The number of elements greater than i to the left of i in a permutation gives the i-th element of the inversion vector. (Skiena-Pemmaraju 2003, p. 69.) T(n,k) is the number of n-permutations that have exactly k 0's in their inversion vector. See evidence in Mathematica code below. - Geoffrey Critzer, May 07 2010
T(n,k) counts the rooted trees with k+1 trunks in forests of "naturally grown" rooted trees with n+2 nodes. This corresponds to sums of coefficients of iterated derivatives representing vectors, Lie derivatives, or infinitesimal generators for flow fields and formal group laws. Cf. links in A139605. - Tom Copeland, Mar 23 2014
A refinement is A036039. - Tom Copeland, Mar 30 2014
From Tom Copeland, Apr 05 2014: (Start)
With initial n=1 and row polynomials of T as p(n,x)=x(x+1)...(x+n-1), the powers of x correspond to the number of trunks of the rooted trees of the "naturally-grown" forest referred to above. With each trunk allowed m colors, p(n,m) gives the number of such non-plane colored trees for the forest with each tree having n+1 vertices.
p(2,m) = m + m^2 = A002378(m) = 2*A000217(m) = 2*(first subdiag of |A238363|).
p(3,m) = 2m + 3m^2 + m^3 = A007531(m+2) = 3*A007290(m+2) = 3*(second subdiag A238363).
p(4,m) = 6m + 11m^2 + 6m^3 + m^4 = A052762(m+3) = 4*A033487(m) = 4*(third subdiag).
From the Joni et al. link, p(n,m) also represents the disposition of n distinguishable flags on m distinguishable flagpoles.
The chromatic polynomial for the complete graph K_n is the falling factorial, which encodes the colorings of the n vertices of K_n and gives a shifted version of p(n,m).
E.g.f. for the row polynomials: (1-y)^(-x).
(End)
A relation to derivatives of the determinant |V(n)| of the n X n Vandermonde matrix V(n) in the indeterminates c(1) thru c(n):
  |V(n)| = Product_{1<=jTom Copeland, Apr 10 2014
From Peter Bala, Jul 21 2014: (Start)
Let M denote the lower unit triangular array A094587 and for k = 0,1,2,... define M(k) to be the lower unit triangular block array
  /I_k 0\
  \ 0  M/
having the k X k identity matrix I_k as the upper left block; in particular, M(0) = M. Then the present triangle equals the infinite matrix product M(0)*M(1)*M(2)*... (which is clearly well defined). See the Example section. (End)
For the relation of this rising factorial to the moments of Viennot's Laguerre stories, see the Hetyei link, p. 4. - Tom Copeland, Oct 01 2015
Can also be seen as the Bell transform of n! without column 0 (and shifted enumeration). For the definition of the Bell transform see A264428. - Peter Luschny, Jan 27 2016

			Triangle  T(n,k) begins:
n\k         0        1        2       3       4      5      6     7    8  9 10
n=0:        1
n=1:        1        1
n=2:        2        3        1
n=3:        6       11        6       1
n=4:       24       50       35      10       1
n=5:      120      274      225      85      15      1
n=6:      720     1764     1624     735     175     21      1
n=7:     5040    13068    13132    6769    1960    322     28     1
n=8:    40320   109584   118124   67284   22449   4536    546    36    1
n=9:   362880  1026576  1172700  723680  269325  63273   9450   870   45  1
n=10: 3628800 10628640 12753576 8409500 3416930 902055 157773 18150 1320 55  1
[Reformatted and extended by _Wolfdieter Lang_, Feb 05 2013]
T(3,2) = 6 because there are 6 permutations of {1,2,3,4} that have exactly 2 0's in their inversion vector: {1, 2, 4, 3}, {1, 3, 2, 4}, {1, 3, 4, 2}, {2, 1, 3, 4},{2, 3, 1, 4}, {2, 3, 4, 1}. The respective inversion vectors are {0, 0, 1}, {0, 1, 0}, {0, 2, 0}, {1, 0, 0}, {2, 0, 0}, {3, 0, 0}. - _Geoffrey Critzer_, May 07 2010
T(3,1)=11 since there are exactly 11 permutations of {1,2,3,4} with exactly 2 cycles, namely, (1)(234), (1)(243), (2)(134), (2)(143), (3)(124), (3)(142), (4)(123), (4)(143), (12)(34), (13)(24), and (14)(23). - _Dennis P. Walsh_, Jan 25 2011
From _Peter Bala_, Jul 21 2014: (Start)
With the arrays M(k) as defined in the Comments section, the infinite product M(0*)M(1)*M(2)*... begins
  / 1          \/1        \/1        \      / 1           \
  | 1  1       ||0 1      ||0 1      |      | 1  1        |
  | 2  2  1    ||0 1 1    ||0 0 1    |... = | 2  3  1     |
  | 6  6  3 1  ||0 2 2 1  ||0 0 1 1  |      | 6 11  6  1  |
  |24 24 12 4 1||0 6 6 3 1||0 0 2 2 1|      |24 50 35 10 1|
  |...         ||...      ||...      |      |...          |
(End)

John H. Conway and Richard K. Guy, The Book of Numbers, New York: Springer-Verlag, 1996. See pp. 93-94.
Sriram Pemmaraju and Steven Skiena, Computational Discrete Mathematics, Cambridge University Press, 2003, pp. 69-71. [Geoffrey Critzer, May 07 2010]

T. D. Noe, Rows n = 0..50 of triangle, flattened
M. Abramowitz and I. A. Stegun, eds., Handbook of Mathematical Functions, National Bureau of Standards, Applied Math. Series 55, Tenth Printing, 1972, Chapter 5, pp. 227-251. [From _Johannes W. Meijer_, Oct 07 2009]
A. Chervov, Decomplexification of the Capelli identities and holomorphic factorization, arxiv 1203.5759 [math.QA], Mar 2012. [_Tom Copeland_, Apr 10 2014]
FindStat - Combinatorial Statistic Finder, The number of saliances of the permutation, The number of cycles in the cycle decomposition of a permutation.
Martin Griffiths, Generating Functions for Extended Stirling Numbers of the First Kind, Journal of Integer Sequences, 17 (2014), #14.6.4.
G. Hetyei, Meixner polynomials of the second kind and quantum algebras representing su(1,1), arXiv preprint arXiv:0909.4352 [math.QA], 2009.
S. Joni, G. Rota, and B. Sagan, From Sets to Functions: Three Elementary Examples, Discrete Mathematics, vol. 37, no. 2-3, pp. 193-202, 1981. [_Tom Copeland_, Apr 05 2014]
Matthieu Josuat-Verges, A q-analog of Schläfli and Gould identities on Stirling numbers, Preprint, arXiv:1610.02965 [math.CO], 2016.
Marin Knežević, Vedran Krčadinac, and Lucija Relić, Matrix products of binomial coefficients and unsigned Stirling numbers, arXiv:2012.15307 [math.CO], 2020.
Lucas Sá and Antonio M. García-García, The Wishart-Sachdev-Ye-Kitaev model: Q-Laguerre spectral density and quantum chaos, arXiv:2104.07647 [hep-th], 2021.
Igor Victorovich Statsenko, On the ordinal numbers of triangles of generalized special numbers, Innovation science No 2-2, State Ufa, Aeterna Publishing House, 2024, pp. 15-19. In Russian.
Dennis Walsh, A short note on unsigned Stirling numbers

See A008275, which is the main entry for these numbers; A094638 (reversed rows).
Diagonals: A000012 A000217 A000914 A001303 A000915 A053567 A112002 A191685. Columns A000142 A000254 A000399 A000454 A000482 A001233 A001234.
From Johannes W. Meijer, Oct 07 2009: (Start)
Row sums equal A000142.
The asymptotic expansions lead to A000142 (n=1), A000142(n=2; minus a(0)), A001710 (n=3), A001715 (n=4), A001720 (n=5), A001725 (n=6), A001730 (n=7), A049388 (n=8), A049389 (n=9), A049398 (n=10), A051431 (n=11), A008279 and A094587.
Cf. A163931 (E(x,m,n)), A028421 (m=2), A163932 (m=3), A163934 (m=4), A163936.
(End)
Cf. A136662.

a130534 n k = a130534_tabl !! n !! k
a130534_row n = a130534_tabl !! n
a130534_tabl = map (map abs) a008275_tabl
-- Reinhard Zumkeller, Mar 18 2013

with(combinat): A130534 := proc(n,m): (-1)^(n+m)*stirling1(n+1,m+1) end proc: seq(seq(A130534(n,m), m=0..n), n=0..10); # Johannes W. Meijer, Oct 07 2009, revised Sep 11 2012
# The function BellMatrix is defined in A264428.
# Adds (1,0,0,0, ..) as column 0 (and shifts the enumeration).
BellMatrix(n -> n!, 9); # Peter Luschny, Jan 27 2016

Table[Table[ Length[Select[Map[ToInversionVector, Permutations[m]], Count[ #, 0] == n &]], {n, 0, m - 1}], {m, 0, 8}] // Grid (* Geoffrey Critzer, May 07 2010 *)
rows = 10;
t = Range[0, rows]!;
T[n_, k_] := BellY[n, k, t];
Table[T[n, k], {n, 1, rows}, {k, 1, n}] // Flatten (* Jean-François Alcover, Jun 22 2018, after Peter Luschny *)

T(0,0) = 1, T(n,k) = 0 if k > n or if n < 0, T(n,k) = T(n-1,k-1) + n*T(n-1,k). T(n,0) = n! = A000142(n). T(2*n,n) = A129505(n+1). Sum_{k=0..n} T(n,k) = (n+1)! = A000142(n+1). Sum_{k=0..n} T(n,k)^2 = A047796(n+1). T(n,k) = |Stirling1(n+1,k+1)|, see A008275. (x+1)(x+2)...(x+n) = Sum_{k=0..n} T(n,k)*x^k. [Corrected by Arie Bos, Jul 11 2008]
Sum_{k=0..n} T(n,k)*x^k = A000007(n), A000142(n), A000142(n+1), A001710(n+2), A001715(n+3), A001720(n+4), A001725(n+5), A001730(n+6), A049388(n), A049389(n), A049398(n), A051431(n) for x = -1, 0, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, respectively. - Philippe Deléham, Nov 13 2007
For k=1..n, let A={a_1,a_2,...,a_k} denote a size-k subset of {1,2,...,n}. Then T(n,n-k) = Sum(Product_{i=1..k} a_i) where the sum is over all subsets A. For example, T(4,1)=50 since 1*2*3 + 1*2*4 + 1*3*4 + 2*3*4 = 50. - Dennis P. Walsh, Jan 25 2011
The preceding formula means T(n,k) = sigma_{n-k}(1,2,3,..,n) with the (n-k)-th elementary symmetric function sigma with the indeterminates chosen as 1,2,...,n. See the Oct 24 2011 comment in A094638 with sigma called there a. - Wolfdieter Lang, Feb 06 2013
From Gary W. Adamson, Jul 08 2011: (Start)
n-th row of the triangle = top row of M^n, where M is the production matrix:
  1, 1;
  1, 2, 1;
  1, 3, 3, 1;
  1, 4, 6, 4, 1;
  ... (End)
Exponential Riordan array [1/(1 - x), log(1/(1 - x))]. Recurrence: T(n+1,k+1) = Sum_{i=0..n-k} (n + 1)!/(n + 1 - i)!*T(n-i,k). - Peter Bala, Jul 21 2014

A094638 Triangle read by rows: T(n,k) = |s(n,n+1-k)|, where s(n,k) are the signed Stirling numbers of the first kind A008276 (1 <= k <= n; in other words, the unsigned Stirling numbers of the first kind in reverse order).

Original entry on oeis.org

1, 1, 1, 1, 3, 2, 1, 6, 11, 6, 1, 10, 35, 50, 24, 1, 15, 85, 225, 274, 120, 1, 21, 175, 735, 1624, 1764, 720, 1, 28, 322, 1960, 6769, 13132, 13068, 5040, 1, 36, 546, 4536, 22449, 67284, 118124, 109584, 40320, 1, 45, 870, 9450, 63273, 269325, 723680, 1172700, 1026576, 362880

Offset: 1

André F. Labossière, May 17 2004

Triangle of coefficients of the polynomial (x+1)(x+2)...(x+n), expanded in decreasing powers of x. - T. D. Noe, Feb 22 2008
Row n also gives the number of permutation of 1..n with complexity 0,1,...,n-1. See the comments in A008275. - N. J. A. Sloane, Feb 08 2019
T(n,k) is the number of deco polyominoes of height n and having k columns. A deco polyomino is a directed column-convex polyomino in which the height, measured along the diagonal, is attained only in the last column. Example: T(2,1)=1 and T(2,2)=1 because the deco polyominoes of height 2 are the vertical and horizontal dominoes, having, respectively, 1 and 2 columns. - Emeric Deutsch, Aug 14 2006
Sum_{k=1..n} k*T(n,k) = A121586. - Emeric Deutsch, Aug 14 2006
Let the triangle U(n,k), 0 <= k <= n, read by rows, be given by [1,0,1,0,1,0,1,0,1,0,1,...] DELTA [1,1,2,2,3,3,4,4,5,5,6,...] where DELTA is the operator defined in A084938; then T(n,k) = U(n-1,k-1). - Philippe Deléham, Jan 06 2007
From Tom Copeland, Dec 15 2007: (Start)
Consider c(t) = column vector(1, t, t^2, t^3, t^4, t^5, ...).
Starting at 1 and sampling every integer to the right, we obtain (1,2,3,4,5,...). And T * c(1) = (1, 1*2, 1*2*3, 1*2*3*4,...), giving n! for n > 0. Call this sequence the right factorial (n+)!.
Starting at 1 and sampling every integer to the left, we obtain (1,0,-1,-2,-3,-4,-5,...). And T * c(-1) = (1, 1*0, 1*0*-1, 1*0*-1*-2,...) = (1, 0, 0, 0, ...), the left factorial (n-)!.
Sampling every other integer to the right, we obtain (1,3,5,7,9,...). T * c(2) = (1, 1*3, 1*3*5, ...) = (1,3,15,105,945,...), giving A001147 for n > 0, the right double factorial, (n+)!!.
Sampling every other integer to the left, we obtain (1,-1,-3,-5,-7,...). T * c(-2) = (1, 1*-1, 1*-1*-3, 1*-1*-3*-5,...) = (1,-1,3,-15,105,-945,...) = signed A001147, the left double factorial, (n-)!!.
Sampling every 3 steps to the right, we obtain (1,4,7,10,...). T * c(3) = (1, 1*4, 1*4*7,...) = (1,4,28,280,...), giving A007559 for n > 0, the right triple factorial, (n+)!!!.
Sampling every 3 steps to the left, we obtain (1,-2,-5,-8,-11,...), giving T * c(-3) = (1, 1*-2, 1*-2*-5, 1*-2*-5*-8,...) = (1,-2,10,-80,880,...) = signed A008544, the left triple factorial, (n-)!!!.
The list partition transform A133314 of [1,T * c(t)] gives [1,T * c(-t)] with all odd terms negated; e.g., LPT[1,T*c(2)] = (1,-1,-1,-3,-15,-105,-945,...) = (1,-A001147). And e.g.f. for [1,T * c(t)] = (1-xt)^(-1/t).
The above results hold for t any real or complex number. (End)
Let R_n(x) be the real and I_n(x) the imaginary part of Product_{k=0..n} (x + I*k). Then, for n=1,2,..., we have R_n(x) = Sum_{k=0..floor((n+1)/2)}(-1)^k*Stirling1(n+1,n+1-2*k)*x^(n+1-2*k), I_n(x) = Sum_{k=0..floor(n/2)}(-1)^(k+1)*Stirling1(n+1,n-2*k)*x^(n-2*k). - Milan Janjic, May 11 2008
T(n,k) is also the number of permutations of n with "reflection length" k (i.e., obtained from 12..n by k not necessarily adjacent transpositions). For example, when n=3, 132, 213, 321 are obtained by one transposition, while 231 and 312 require two transpositions. - Kyle Petersen, Oct 15 2008
From Tom Copeland, Nov 02 2010: (Start)
[x^(y+1) D]^n = x^(n*y) [T(n,1)(xD)^n + T(n,2)y (xD)^(n-1) + ... + T(n,n)y^(n-1)(xD)], with D the derivative w.r.t. x.
E.g., [x^(y+1) D]^4 = x^(4*y) [(xD)^4 + 6 y(xD)^3 + 11 y^2(xD)^2 + 6 y^3(xD)].
(xD)^m can be further expanded in terms of the Stirling numbers of the second kind and operators of the form x^j D^j. (End)
With offset 0, 0 <= k <= n: T(n,k) is the sum of products of each size k subset of {1,2,...,n}. For example, T(3,2) = 11 because there are three subsets of size two: {1,2},{1,3},{2,3}. 1*2 + 1*3 + 2*3 = 11. - Geoffrey Critzer, Feb 04 2011
The Kn11, Fi1 and Fi2 triangle sums link this triangle with two sequences, see the crossrefs. For the definitions of these triangle sums see A180662. The mirror image of this triangle is A130534. - Johannes W. Meijer, Apr 20 2011
T(n+1,k+1) is the elementary symmetric function a_k(1,2,...,n), n >= 0, k >= 0, (a_0(0):=1). See the T. D. Noe and Geoffrey Critzer comments given above. For a proof see the Stanley reference, p. 19, Second Proof. - Wolfdieter Lang, Oct 24 2011
Let g(t) = 1/d(log(P(j+1,-t)))/dt (see Tom Copeland's 2007 formulas). The Mellin transform (t to s) of t*Dirac[g(t)] gives Sum_{n=1..j} n^(-s), which as j tends to infinity gives the Riemann zeta function for Re(s) > 1. Dirac(x) is the Dirac delta function. The complex contour integral along a circle of radius 1 centered at z=1 of z^s/g(z) gives the same result. - Tom Copeland, Dec 02 2011
Rows are coefficients of the polynomial expansions of the Pochhammer symbol, or rising factorial, Pch(n,x) = (x+n-1)!/(x-1)!. Expansion of Pch(n,xD) = Pch(n,Bell(.,:xD:)) in a polynomial with terms :xD:^k=x^k*D^k gives the Lah numbers A008297. Bell(n,x) are the unsigned Bell polynomials or Stirling polynomials of the second kind A008277. - Tom Copeland, Mar 01 2014
From Tom Copeland, Dec 09 2016: (Start)
The Betti numbers, or dimension, of the pure braid group cohomology. See pp. 12 and 13 of the Hyde and Lagarias link.
Row polynomials and their products appear in presentation of the Jack symmetric functions of R. Stanley. See Copeland link on the Witt differential generator.
(End)
From Tom Copeland, Dec 16 2019: (Start)
The e.g.f. given by Copeland in the formula section appears in a combinatorial Dyson-Schwinger equation of quantum field theory in Yeats in Thm. 2 on p. 62 related to a Hopf algebra of rooted trees. See also the Green function on p. 70.
Per comments above, this array contains the coefficients in the expansion in polynomials of the Euler, or state number, operator xD of the rising factorials Pch(n,xD) = (xD+n-1)!/(xD-1)! = x [:Dx:^n/n!]x^{-1} = L_n^{-1}(-:xD:), where :Dx:^n = D^n x^n and :xD:^n = x^n D^n. The polynomials L_n^{-1} are the Laguerre polynomials of order -1, i.e., normalized Lah polynomials.
The Witt differential operators L_n = x^(n+1) D and the row e.g.f.s appear in Hopf and dual Hopf algebra relations presented by Foissy. The Witt operators satisfy L_n L_k - L_k L_n = (k-n) L_(n+k), as for the dual Hopf algebra. (End)

			Triangle starts:
  1;
  1,  1;
  1,  3,  2;
  1,  6, 11,  6;
  1, 10, 35, 50, 24;
...

M. Miyata and J. W. Son, On the complexity of permutations and the metric space of bijections, Tensor, 60 (1998), No. 1, 109-116 (MR1768839).
Jerome Spanier and Keith B. Oldham, "Atlas of Functions", Hemisphere Publishing Corp., 1987, chapter 18, equations 18:4:2 - 18:4:8 at page 151.
R. P. Stanley, Enumerative Combinatorics, Vol. 1, Cambridge University Press, 1997.

T. D. Noe, Rows n = 1..51 of triangle, flattened
E. Barcucci, A. Del Lungo and R. Pinzani, "Deco" polyominoes, permutations and random generation, Theoretical Computer Science, 159, 1996, 29-42.
F. Bergeron, Ph. Flajolet and B. Salvy, Varieties of Increasing Trees, Lecture Notes in Computer Science vol. 581, ed. J.-C. Raoult, Springer 1992, pp. 24-48. Added Mar 01 2014
F. Bergeron, Philippe Flajolet and Bruno Salvy, Varieties of increasing trees, HAL, Rapport De Recherche Inria. Added Mar 01 2014
Tom Copeland, Addendum to Mathemagical Forests
Tom Copeland, A Class of Differential Operators and the Stirling Numbers
Tom Copeland, Mellin Interpolation of Differential Ops and Associated Infinigens and Appell Polynomials: The Ordered, Laguerre, and Scherk-Witt-Lie Diff Ops
Tom Copeland, Witt differential generator for special Jack symmetric functions / polynomials
FindStat - Combinatorial Statistic Finder, The absolute length of a permutation
L. Foissy, Faa di Bruno subalgebras of the Hopf algebra of planar trees from combinatorial Dyson-Schwinger equations, arXiv:0707.1204 [math.RA], (2007).
O. Furdui, T. Trif, On the Summation of Certain Iterated Series, J. Int. Seq. 14 (2011) #11.6.1
F. Hivert, J.-C. Novelli and J.-Y. Thibon, The Algebra of Binary Search Trees, Theoretical Computer Science, 339 (2005), 129-165.
T. Hyde and J. Lagarias Polynomial splitting measures and cohomology of the pure braid group, arXiv preprint arXiv:1604.05359 [math.RT], 2016.
MathOverflow, Motivation of Virasoro algebra, an answer by Tom Copeland to an MO question posed in 2012.
Romeo Mestrovic, Lucas' theorem: its generalizations, extensions and applications (1878--2014), arXiv preprint arXiv:1409.3820 [math.NT], 2014.
Robert E. Moritz, On the sum of products of n consecutive integers, Univ. Washington Publications in Math., 1 (No. 3, 1926), 44-49 [Annotated scanned copy]
M. D. Schmidt, Generalized j-Factorial Functions, Polynomials, and Applications , J. Int. Seq. 13 (2010), 10.6.7.
M. Z. Spivey, On Solutions to a General Combinatorial Recurrence, J. Int. Seq. 14 (2011) # 11.9.7.
K. Yeats, A Combinatorial Perspective on Quantum Field Theory, SpringerBriefs in Mathematical Physics, Vol. 15, 2017.

A008276 gives the (signed) Stirling numbers of the first kind.
Cf. A000108, A014137, A001246, A033536, A000984, A094639, A006134, A082894, A002897, A079727, A000217 (2nd column), A000914 (3rd column), A001303 (4th column), A000915 (5th column), A053567 (6th column), A000142 (row sums).
Triangle sums (see the comments): A124380 (Kn11), A001710 (Fi1, Fi2). - Johannes W. Meijer, Apr 20 2011
Cf. A121586, A130534, A143491.

Flat(List([1..10], n-> List([1..n], k-> Stirling1(n,n-k+1) ))); # G. C. Greubel, Dec 29 2019

a094638 n k = a094638_tabl !! (n-1) !! (k-1)
a094638_row n = a094638_tabl !! (n-1)
a094638_tabl = map reverse a130534_tabl
-- Reinhard Zumkeller, Aug 01 2014

[(-1)^(k+1)*StirlingFirst(n,n-k+1): k in [1..n], n in [1..10]]; // G. C. Greubel, Dec 29 2019

T:=(n,k)->abs(Stirling1(n,n+1-k)): for n from 1 to 10 do seq(T(n,k),k=1..n) od; # yields sequence in triangular form. # Emeric Deutsch, Aug 14 2006

Table[CoefficientList[Series[Product[1 + i x, {i,n}], {x,0,20}], x], {n,0,6}] (* Geoffrey Critzer, Feb 04 2011 *)
Table[Abs@StirlingS1[n, n-k+1], {n, 10}, {k, n}]//Flatten (* Michael De Vlieger, Aug 29 2015 *)

create_list(abs(stirling1(n+1,n-k+1)),n,0,10,k,0,n); /* Emanuele Munarini, Jun 01 2012 */

{T(n,k)=if(n<1 || k>n,0,(n-1)!*polcoeff(polcoeff(x*y/(1 - x*y+x*O(x^n))^(1 + 1/y),n,x),k,y))} /* Paul D. Hanna, Jul 21 2011 */

[[stirling_number1(n, n-k+1) for k in (1..n)] for n in (1..10)] # G. C. Greubel, Dec 29 2019

With P(n,t) = Sum_{k=0..n-1} T(n,k+1) * t^k = 1*(1+t)*(1+2t)...(1+(n-1)*t) and P(0,t)=1, exp[P(.,t)*x] = (1-tx)^(-1/t). T(n,k+1) = (1/k!) (D_t)^k (D_x)^n [ (1-tx)^(-1/t) - 1 ] evaluated at t=x=0. (1-tx)^(-1/t) - 1 is the e.g.f. for a plane m-ary tree when t=m-1. See Bergeron et al. in "Varieties of Increasing Trees". - Tom Copeland, Dec 09 2007
First comment and formula above rephrased as o.g.f. for row n: Product_{i=0...n} (1+i*x). - Geoffrey Critzer, Feb 04 2011
n-th row polynomials with alternate signs are the characteristic polynomials of the (n-1)x(n-1) matrices with 1's in the superdiagonal, (1,2,3,...) in the main diagonal, and the rest zeros. For example, the characteristic polynomial of [1,1,0; 0,2,1; 0,0,3] is x^3 - 6*x^2 + 11*x - 6. - Gary W. Adamson, Jun 28 2011
E.g.f.: A(x,y) = x*y/(1 - x*y)^(1 + 1/y) = Sum_{n>=1, k=1..n} T(n,k)*x^n*y^k/(n-1)!. - Paul D. Hanna, Jul 21 2011
With F(x,t) = (1-t*x)^(-1/t) - 1 an e.g.f. for the row polynomials P(n,t) of A094638 with P(0,t)=0, G(x,t)= [1-(1+x)^(-t)]/t is the comp. inverse in x. Consequently, with H(x,t) = 1/(dG(x,t)/dx) = (1+x)^(t+1),
  P(n,t) = [(H(x,t)*d/dx)^n] x evaluated at x=0; i.e.,
  F(x,t) = exp[x*P(.,t)] = exp[x*H(u,t)*d/du] u, evaluated at u = 0.
  Also, dF(x,t)/dx = H(F(x,t),t). - Tom Copeland, Sep 20 2011
T(n,k) = |A008276(n,k)|. - R. J. Mathar, May 19 2016
The row polynomials of this entry are the reversed row polynomials of A143491 multiplied by (1+x). E.g., (1+x)(1 + 5x + 6x^2) = (1 + 6x + 11x^2 + 6x^3). - Tom Copeland, Dec 11 2016
Regarding the row e.g.f.s in Copeland's 2007 formulas, e.g.f.s for A001710, A001715, and A001720 give the compositional inverses of the e.g.f. here for t = 2, 3, and 4 respectively. - Tom Copeland, Dec 28 2019

Edited by Emeric Deutsch, Aug 14 2006

A001296 4-dimensional pyramidal numbers: a(n) = (3n+1)binomial(n+2, 3)/4. Also Stirling2(n+2, n).

Original entry on oeis.org

0, 1, 7, 25, 65, 140, 266, 462, 750, 1155, 1705, 2431, 3367, 4550, 6020, 7820, 9996, 12597, 15675, 19285, 23485, 28336, 33902, 40250, 47450, 55575, 64701, 74907, 86275, 98890, 112840, 128216, 145112, 163625, 183855, 205905, 229881, 255892, 284050, 314470

Offset: 0

N. J. A. Sloane

Permutations avoiding 12-3 that contain the pattern 31-2 exactly once.
Kekulé numbers for certain benzenoids. - Emeric Deutsch, Nov 18 2005
Partial sums of A002411. - Jonathan Vos Post, Mar 16 2006
If Y is a 3-subset of an n-set X then, for n>=6, a(n-5) is the number of 6-subsets of X having at least two elements in common with Y. - Milan Janjic, Nov 23 2007
Starting with 1 = binomial transform of [1, 6, 12, 10, 3, 0, 0, 0, ...]. Equals row sums of triangle A143037. - Gary W. Adamson, Jul 18 2008
Rephrasing the Perry formula of 2003: a(n) is the sum of all products of all two  numbers less than or equal to n, including the squares. Example: for n=3 the sum of these products is 1*1 + 1*2 + 1*3 + 2*2 + 2*3 + 3*3 = 25. - J. M. Bergot, Jul 16 2011
Half of the partial sums of A011379. [Jolley, Summation of Series, Dover (1961), page 12 eq (66).] - R. J. Mathar, Oct 03 2011
Also the number of (w,x,y,z) with all terms in {1,...,n+1} and w < x >= y > z (see A211795). - Clark Kimberling, May 19 2012
Convolution of A000027 with A000326. - Bruno Berselli, Dec 06 2012
This sequence is related to A000292 by a(n) = n*A000292(n) - Sum_{i=0..n-1} A000292(i) for n>0. - Bruno Berselli, Nov 23 2017
a(n-2) is the maximum number of intersections made from the perpendicular bisectors of all pair combinations of n points. - Ian Tam, Dec 22 2020

			G.f. = x + 7*x^2 + 25*x^3 + 65*x^4 + 140*x^5 + 266*x^6 + 462*x^7 + 750*x^8 + 1155*x^9 + ...

M. Abramowitz and I. A. Stegun, eds., Handbook of Mathematical Functions, National Bureau of Standards Applied Math. Series 55, 1964 (and various reprintings), p. 835.
A. H. Beiler, Recreations in the Theory of Numbers, Dover, NY, 1964, p. 195.
L. Comtet, Advanced Combinatorics, Reidel, 1974, p. 227, #16.
S. J. Cyvin and I. Gutman, Kekulé structures in benzenoid hydrocarbons, Lecture Notes in Chemistry, No. 46, Springer, New York, 1988 (see p. 166, Table 10.4/I/3).
F. N. David, M. G. Kendall and D. E. Barton, Symmetric Function and Allied Tables, Cambridge, 1966, p. 223.
N. J. A. Sloane, A Handbook of Integer Sequences, Academic Press, 1973 (includes this sequence).
N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).

T. D. Noe, Table of n, a(n) for n = 0..1000
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].
P. Aluffi, Degrees of projections of rank loci, arXiv:1408.1702 [math.AG], 2014. ["After compiling the results of many explicit computations, we noticed that many of the numbers d_{n,r,S} appear in the existing literature in contexts far removed from the enumerative geometry of rank conditions; we owe this surprising (to us) observation to perusal of [Slo14]."]
S. Butler and P. Karasik, A note on nested sums, J. Int. Seq. 13 (2010), 10.4.4, page 5.
M. Griffiths, Remodified Bessel Functions via Coincidences and Near Coincidences, Journal of Integer Sequences, Vol. 14 (2011), Article 11.7.1.
L. Hogben, Choice and Chance by Cardpack and Chessboard, Vol. 1, Chanticleer Press, NY, 1950, p. 36.
C. Krishnamachaki, The operator (xD)^n, J. Indian Math. Soc., 15 (1923),3-4. [Annotated scanned copy]
Feihu Liu, Guoce Xin, and Chen Zhang, Ehrhart Polynomials of Order Polytopes: Interpreting Combinatorial Sequences on the OEIS, arXiv:2412.18744 [math.CO], 2024. See p. 29.
T. Mansour, Restricted permutations by patterns of type 2-1, arXiv:math/0202219 [math.CO], 2002.
Simon Plouffe, Approximations de séries génératrices et quelques conjectures, Dissertation, Université du Québec à Montréal, 1992; arXiv:0911.4975 [math.NT], 2009.
Simon Plouffe, 1031 Generating Functions, Appendix to Thesis, Montreal, 1992
Eric Weisstein's World of Mathematics, Stirling numbers of the 2nd kind.
Index to sequences related to pyramidal numbers
Index entries for linear recurrences with constant coefficients, signature (5,-10,10,-5,1).

a(n)=f(n, 2) where f is given in A034261.
Cf. A000217, A000326, A001297, A001298, A002411, A008277, A008517, A094262.
a(n)= A093560(n+3, 4), (3, 1)-Pascal column.
Cf. A220212 for a list of sequences produced by the convolution of the natural numbers with the k-gonal numbers.
Cf. similar sequences listed in A241765 and A254142.
Cf. A000914.

/* A000027 convolved with A000326: */ A000326:=func; [&+[(n-i+1)*A000326(i): i in [0..n]]: n in [0..40]]; // Bruno Berselli, Dec 06 2012

[(3*n+1)*Binomial(n+2,3)/4: n in [0..40]]; // Vincenzo Librandi, Jul 30 2014

A001296:=-(1+2*z)/(z-1)**5; # Simon Plouffe in his 1992 dissertation for sequence without the leading zero

Table[n*(1+n)*(2+n)*(1+3*n)/24, {n, 0, 100}]
CoefficientList[Series[x (1 + 2 x)/(1 - x)^5, {x, 0, 40}], x] (* Vincenzo Librandi, Jul 30 2014 *)
Table[StirlingS2[n+2, n], {n, 0, 40}] (* Jean-François Alcover, Jun 24 2015 *)
Table[ListCorrelate[Accumulate[Range[n]],Range[n]],{n,0,40}]//Flatten (* or *) LinearRecurrence[{5,-10,10,-5,1},{0,1,7,25,65},40] (* Harvey P. Dale, Aug 14 2017 *)

t(n)=n*(n+1)/2
for(i=1,30,print1(","sum(j=1,i,j*t(j))))

{a(n) = n * (n+1) * (n+2) * (3*n+1) / 24}; /* Michael Somos, Sep 04 2017 */

[stirling_number2(n+2,n) for n in range(0,38)] # Zerinvary Lajos, Mar 14 2009

a(n) = n*(1+n)*(2+n)*(1+3*n)/24. - T. D. Noe, Jan 21 2008
G.f.: x*(1+2*x)/(1-x)^5. - Paul Barry, Jul 23 2003
a(n) = Sum_{j=0..n} j*A000217(j). - Jon Perry, Jul 28 2003
E.g.f. with offset -1: exp(x)*(1*(x^2)/2! + 4*(x^3)/3! + 3*(x^4)/4!). For the coefficients [1, 4, 3] see triangle A112493.
E.g.f. x*exp(x)*(24 + 60*x + 28*x^2 + 3*x^3)/24 (above e.g.f. differentiated).
a(n) = 4*a(n-1) - 6*a(n-2) + 4*a(n-3) - a(n-4) + 3. - Kieren MacMillan, Sep 29 2008
a(n) = 5*a(n-1) - 10*a(n-2) + 10*a(n-3) - 5*a(n-4) + a(n-5). - Jaume Oliver Lafont, Nov 23 2008
O.g.f. is D^2(x/(1-x)) = D^3(x), where D is the operator x/(1-x)*d/dx. - Peter Bala, Jul 02 2012
a(n) = A153978(n)/2. - J. M. Bergot, Aug 09 2013
a(n) = A002817(n) + A000292(n-1). - J. M. Bergot, Aug 29 2013; [corrected by Cyril Damamme, Feb 26 2018]
a(n) = A000914(n+1) - 2 * A000330(n+1). - Antal Pinter, Dec 31 2015
a(n) = A080852(3,n-1). - R. J. Mathar, Jul 28 2016
a(n) = 1*(1+2+...+n) + 2*(2+3+...+n) + ... + n*n. For example, a(6) = 266 = 1(1+2+3+4+5+6) + 2*(2+3+4+5+6) + 3*(3+4+5+6) + 4*(4+5+6) + 5*(5+6) + 6*(6).- J. M. Bergot, Apr 20 2017
a(n) = A000914(-2-n) for all n in Z. - Michael Somos, Sep 04 2017
a(n) = A000292(n) + A050534(n+1). - Cyril Damamme, Feb 26 2018
From Amiram Eldar, Jul 02 2020: (Start)
Sum_{n>=1} 1/a(n) = (6/5) * (47 - 3*sqrt(3)*Pi - 27*log(3)).
Sum_{n>=1} (-1)^(n+1)/a(n) = (6/5) * (16*log(2) + 6*sqrt(3)*Pi - 43). (End)

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

Original entry on oeis.org

0, 2, 9, 24, 50, 90, 147, 224, 324, 450, 605, 792, 1014, 1274, 1575, 1920, 2312, 2754, 3249, 3800, 4410, 5082, 5819, 6624, 7500, 8450, 9477, 10584, 11774, 13050, 14415, 15872, 17424, 19074, 20825, 22680, 24642, 26714, 28899, 31200, 33620, 36162, 38829, 41624

Offset: 0

N. J. A. Sloane

a(n) is the largest number that is not the sum of distinct numbers of form kn+1, k >= 0. - David W. Wilson, Dec 11 1999
Sum of the nontriangular numbers between successive triangular numbers. 1, (2), 3, (4, 5), 6, (7, 8, 9), 10, (11, 12, 13, 14), 15, ... Sum of the terms in brackets. Or sum of n consecutive integers beginning with T(n) + 1, where T(n) = n(n+1)/2. - Amarnath Murthy, Aug 27 2005
Apparently this is also the splittance (as defined by Hammer & Simeone, 1977) of the Kneser graphs of the form K(n+3,2). - Felix Goldberg (felixg(AT)tx.technion.ac.il), Jul 13 2009
Row sums of triangle A159797. - Omar E. Pol, Jul 24 2009
The same results occur when one plots the points (1,3), (3,6), (6,10), (10,15), and so on, for all the triangular numbers and finds the area beneath. Take three consecutive triangular numbers and label them a, b, c; the area created is simply (b-a)*(b+c)/2. Thus for 6,10,15 the area beneath the line defined by the points (6,10) and (10,15) is (10-6)*(10+15)/2 = 50. - J. M. Bergot, Jun 28 2011
Let P = ab where a and b are nonequal prime numbers > 1. Let Q be the product of all divisors of P^n. Q can be expressed as P^k, where k = n*(n+1)^2/2. This follows from the fact that all divisors are of the form a^i*b^j, for i,j from 0 to n. An example is given below. In the more general case, where P is the product of m nonequal prime numbers, k = n*(n+1)^m/2. When m=3, the sequence is the same as A092364. - James A. Raymond & Douglas Raymond, Dec 04 2011
For n > 0: sum of n-th row in A014132, seen as a triangle read by rows. - Reinhard Zumkeller, Dec 12 2012
Partial sums of A005449. - Omar E. Pol, Jan 12 2013
a(n) is the sum of x (or y) coordinates for an n X n square lattice in the upper right quadrant of Z^2 whose corner points are (0, 0), (0, n), (n, 0), and (n, n). - Joseph Wheat, Feb 03 2018
a(n) is the number of permutations of [n+2] that contain exactly 2 elements which are not left-to-right minimal. E.g., the 9 permutations comprising a(2) are 2134, 2143, 3124, 3142, 4123, 4132, 2314, 2413, 3412. - Andy Niedermaier, May 07 2022

			Let P^n=6^2. The product of the divisors of 36 = 10077796 = 6^9, i.e., for n=2, k=9. - _James A. Raymond_ & Douglas Raymond, Dec 04 2011

N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).

T. D. Noe, Table of n, a(n) for n = 0..1000
Quang T. Bach, Roshil Paudyal and Jeffrey B. Remmel, A Fibonacci analogue of Stirling numbers, arXiv preprint arXiv:1510.04310 [math.CO], 2015.
Paul Barry, On the Gap-sum and Gap-product Sequences of Integer Sequences, arXiv:2104.05593 [math.CO], 2021.
Dexter Jane L. Indong and Gilbert R. Peralta, Inversions of permutations in Symmetric, Alternating, and Dihedral Gropus, JIS, Vol. 11 (2008), Article 08.4.3.
S. M. Losanitsch, Die Isomerie-Arten bei den Homologen der Paraffin-Reihe, Chem. Ber., Vol. 30 (1897), pp. 1917-1926.
S. M. Losanitsch, Die Isomerie-Arten bei den Homologen der Paraffin-Reihe, Chem. Ber., Vol. 30 (1897), pp. 1917-1926. (Annotated scanned copy)
Luis Manuel Rivera, Integer sequences and k-commuting permutations, arXiv preprint arXiv:1406.3081 [math.CO], 2014-2015.
Index entries for linear recurrences with constant coefficients, signature (4,-6,4,-1).
Index entries for two-way infinite sequences.

Cf. A002411: -a(-1-n).
Cf. A000217, A000292, A000330, A014132, A035006, A045991, A092364, A159797, A163274.
Cf. A000914 (partial sums), A005449 (first differences).
Cf. similar sequences of the type n*(n+1)*(n+k)/2 listed in A267370.
Cf. A002378, A027480.
A bisection of A330298.

List([0..10^3],n->n*(n+1)^2/2); # Muniru A Asiru, Feb 04 2018

a006002 n = n * (n + 1) ^ 2 `div` 2  -- Reinhard Zumkeller, Dec 12 2012

[n*(n+1)^2/2 : n in [0..50]]; // Wesley Ivan Hurt, Nov 17 2014

seq(binomial(n+1,2)*(n+1), n=0..36); # Zerinvary Lajos, Apr 25 2007

Table[(n^3-n^2)/2, {n, 41}] (* Zerinvary Lajos, Mar 21 2007 *)
LinearRecurrence[{4,-6,4,-1}, {0,2,9,24}, 40] (* Harvey P. Dale, Aug 14 2012 *)
Accumulate @ # (# + 1) & [Range[0,50]] (* Waldemar Puszkarz, Jan 24 2015 *)

PARI
```
a(n)=n*(n+1)^2/2
```

G.f.: x*(x + 2)/(1 - x)^4. - Michael Somos, Jan 30 2004
a(n) = (n + 1) * binomial(n+1, 2). - Zerinvary Lajos, Jan 10 2006
a(n) = A035006(n+1)/4. - Johannes W. Meijer, Feb 04 2010
a(n) = 2*binomial(n+1, 2) + 3*binomial(n+1, 3). - Gary Detlefs, Jun 06 2010
a(n) = 4*a(n-1) - 6*a(n-2) + 4*a(n-3) - a(n-4). - Harvey P. Dale, Aug 14 2012
a(n) = A000292(n) + A000330(n). - Omar E. Pol, Jan 11 2013
a(n) = A045991(n+1)/2. - J. M. Bergot, Aug 10 2013
a(n) = Sum_{j=1..n} Sum_{i=1..j} (2*j - i + 1). - Wesley Ivan Hurt, Nov 17 2014
a(n) = Sum_{i=0..n} n*(n - i) + i. - Bruno Berselli, Jan 13 2016
a(n) = t(n, A000217(n)), where t(h,k) = A000217(h) + h*k. - Bruno Berselli, Feb 28 2017
Sum_{n>0} 1/a(n) = 4 - Pi^2/3. - Jaume Oliver Lafont, Jul 11 2017 [corrected by Amiram Eldar, Jan 28 2022]
E.g.f.: exp(x)*x*(4 + 5*x + x^2)/2. - Stefano Spezia, May 21 2021
Sum_{n>=1} (-1)^(n+1)/a(n) = Pi^2/6 + 4*log(2) - 4. - Amiram Eldar, Jan 28 2022
From J.S. Seneschal, Jun 27 2024: (Start)
a(n) = (A002378(n)^2/2)/n = (n+1)/2 * A002378(n).
a(n) = A027480(n) - A000217(n). (End)

A112007 Coefficient triangle for polynomials used for o.g.f.s for unsigned Stirling1 diagonals.

Original entry on oeis.org

1, 2, 1, 6, 8, 1, 24, 58, 22, 1, 120, 444, 328, 52, 1, 720, 3708, 4400, 1452, 114, 1, 5040, 33984, 58140, 32120, 5610, 240, 1, 40320, 341136, 785304, 644020, 195800, 19950, 494, 1, 362880, 3733920, 11026296, 12440064, 5765500, 1062500, 67260, 1004, 1

Offset: 0

Wolfdieter Lang, Sep 12 2005

This is the row reversed second-order Eulerian triangle A008517(k+1,k+1-m). For references see A008517.
The o.g.f. for the k-th diagonal, k >= 1, of the unsigned Stirling1 triangle |A008275| is G1(1,x)=1/(1-x) if k=1 and G1(k,x) = g1(k-2,x)/(1-x)^(2*k-1), if k >= 2, with the row polynomials g1(k;x):=Sum_{m=0..k} a(k,m)*x^m.
The recurrence eq. for the row polynomials is g1(k,x)=((k+1)+k*x)*g1(k-1,x) + x*(1-x)*(d/dx)g1(k-1,x), k >= 1, with input g1(0,x):=1.
The column sequences start with A000142 (factorials), A002538, A002539, A112008, A112485.
This o.g.f. computation was inspired by Bender et al. article where the Stirling polynomials have been rediscussed.
The A163936 triangle is identical to the triangle given above except for an extra right hand column [1, 0, 0, 0, ... ]. The A163936 triangle is related to the higher order exponential integrals E(x,m,n), see A163931 and A163932. - Johannes W. Meijer, Oct 16 2009

			Triangle begins:
    1;
    2,   1;
    6,   8,   1;
   24,  58,  22,   1;
  120, 444, 328,  52,   1;
  ...
G.f. for k=3 sequence A000914(n-1), [2,11,35,85,175,322,546,...], is G1(3,x)= g1(1,x)/(1-x)^5= (2+x)/(1-x)^5.

Robert Israel, Table of n, a(n) for n = 0..10010 (rows 0 to 140, flattened)
Peter Bala, Diagonals of triangles with generating function exp(t*F(x))
C. M. Bender, D. C. Brody and B. K. Meister, Bernoulli-like polynomials associated with Stirling Numbers, arXiv:math-ph/0509008 [math-ph], 2005.
E. Burlachenko, Composition polynomials of the RNA matrix and B-composition polynomials of the Riordan pseudo-involution, arXiv:1907.12272 [math.NT], 2019.
Tom Copeland, Generators, Inversion, and Matrix, Binomial, and Integral Transforms
Wolfdieter Lang, First 10 rows.
Wolfdieter Lang, On Generating functions of Diagonals Sequences of Sheffer and Riordan Number Triangles, arXiv:1708.01421 [math.NT], August 2017.
R. Paris, A uniform asymptotic expansion for the incomplete gamma function, Journal of Computational and Applied Mathematics, 148 (2002), p. 223-239. See 234. [_Tom Copeland_, Jan 03 2016]
Andrew Elvey Price, Alan D. Sokal, Phylogenetic trees, augmented perfect matchings, and a Thron-type continued fraction (T-fraction) for the Ward polynomials, arXiv:2001.01468 [math.CO], 2020.

Row sums give A001147(k+1) = (2*k+1)!!, k>=0.

a:= proc(k,m) option remember; if m >= 0 and k >= 0 then (k+m+1)*procname(k-1,m)+(k-m+1)*procname(k-1,m-1) else 0 fi end proc:
a(0,0):= 1:
seq(seq(a(k,m),m=0..k),k=0..10); # Robert Israel, Jul 20 2017

a[k_, m_] = Sum[(-1)^(k + n + 1)*Binomial[2k + 3, n]*StirlingS1[m + k - n + 2, m + 1 - n], {n, 0, m}]; Flatten[Table[a[k, m], {k, 0, 8}, {m, 0, k}]][[1 ;; 45]] (* Jean-François Alcover, Jun 01 2011, after Johannes W. Meijer *)

a(k, m)=sum(n=0, m, (-1)^(k + n + 1)*binomial(2*k + 3, n)*stirling(m + k - n + 2, m + 1 - n, 1));
for(k=0, 10, for(m=0, k, print1(a(k, m),", "))) \\ Indranil Ghosh, Jul 21 2017

a(k, m) = (k+m+1)*a(k-1, m) + (k-m+1)*a(k-1, m-1), if k >= m >= 0, a(0, 0)=1; a(k, -1):=0, otherwise 0.
a(k,m) = Sum_{n=0..m} (-1)^(k+n+1)*C(2*k+3,n)*Stirling1(m+k-n+2,m+1-n). - Johannes W. Meijer, Oct 16 2009
The compositional inverse (with respect to x) of y = y(t,x) = (x+t*log(1-x)) is  x = x(t,y) = 1/(1-t)*y + t/(1-t)^3*y^2/2! + (2*t+t^2)/(1-t)^5*y^3/3! + (6*t+8*t^2+t^3)/(1-t)^7*y^4/4! + .... The numerator polynomials of the rational functions in t are the row polynomials of this triangle. As observed above, the rational functions in t are the generating functions for the diagonals of |A008275|. See the Bala link for a proof. Cf. A008517. - Peter Bala, Dec 02 2011

A047659 Number of ways to place 3 nonattacking queens on an n X n board.

Original entry on oeis.org

0, 0, 0, 0, 24, 204, 1024, 3628, 10320, 25096, 54400, 107880, 199400, 348020, 579264, 926324, 1431584, 2148048, 3141120, 4490256, 6291000, 8656860, 11721600, 15641340, 20597104, 26797144, 34479744, 43915768, 55411720, 69312516, 86004800, 105919940

Offset: 0

Paul Zimmermann

Lucas mentions that the number of ways of placing p <= n non-attacking queens on an n X n chessboard is given by a polynomial in n of degree 2p and attribute the result to Mantel, professor in Delft. Cf. Stanley, exercise 15.

E. Landau, Naturwissenschaftliche Wochenschrift (Aug. 2 1896).
R. P. Stanley, Enumerative Combinatorics, vol. I, exercise 15 in chapter 4 (and its solution) asks one to show the existence of a rational generating function for the number of ways of placing k non-attacking queens on an n X n chessboard.

Vincenzo Librandi, Table of n, a(n) for n = 0..1000
S. Chaiken, C. R. H. Hanusa and T. Zaslavsky, A q-queens problem I. General theory, Jan 26 2013. - _N. J. A. Sloane_, Feb 16 2013
S. Chaiken, C. R. H. Hanusa and T. Zaslavsky, A q-Queens Problem. IV. Queens, Bishops, Nightriders (and Rooks), arXiv:1609.00853 [math.CO], Sep 03 2016.
Christopher R. H. Hanusa, T Zaslavsky, S Chaiken, A q-Queens Problem. IV. Queens, Bishops, Nightriders (and Rooks), arXiv preprint arXiv:1609.00853 [math.CO], 2016-2020.
V. Kotesovec, Non-attacking chess pieces, 6ed, 2013, p. 11.
Edmund Landau, Ueber das Achtdamenproblem und seine Verallgemeinerung, Naturwissenschaftliche Wochenschrift, Aug 02 1896.
Edouard Lucas, Récréations mathématiques, Gauthier-Villars, Paris, 1882-1894, Vol. I, p. 228.
Antal Pinter, Numerical solution of the k=3 Queens problem, 2011, P(3) at p.8-9.
I. Rivin, I. Vardi and P. Zimmermann, The n-queens problem, Amer. Math. Monthly, 101 (1994), 629-639.
Wenxi Wang, Muhammad Usman, Alyas Almaawi, Kaiyuan Wang, Kuldeep S. Meel, Sarfraz Khurshid, A Study of Symmetry Breaking Predicates and Model Counting, National University of Singapore (2020).
Index entries for linear recurrences with constant coefficients, signature (5,-8,0,14,-14,0,8,-5,1).

Cf. A036464, A061994, A108792, A176186, A178721.

Column k=3 of A348129.

[(3*(2*n-1)*(-1)^n +4*n^6 -40*n^5 +158*n^4 -300*n^3 +264*n^2 -86*n +3)/24: n in [0..35]]; // Vincenzo Librandi, Sep 21 2015

f:=n-> n^6/6 - 5*n^5/3 + 79*n^4/12 - 25*n^3/2 + 11*n^2 - 43*n/12 + 1/8 + (-1)^n*(n/4 - 1/8); [seq(f(n),n=1..40)]; # N. J. A. Sloane, Feb 16 2013

Table[If[EvenQ[n],n (n-2)^2 (2n^3-12n^2+23n-10)/12,(n-1)(n-3) (2n^4- 12n^3+25n^2-14n+1)/12],{n,0,30}] (* or *) LinearRecurrence[ {5,-8,0,14,-14,0,8,-5,1},{0,0,0,0,24,204,1024,3628,10320},30] (* Harvey P. Dale, Nov 06 2011 *)

a(n)=if(n%2, (n - 1)*(n - 3)*(2*n^4 - 12*n^3 + 25*n^2 - 14*n + 1), n*(n - 2)^2*(2*n^3 - 12*n^2 + 23*n - 10))/12 \\ Charles R Greathouse IV, Feb 09 2017

a(n) = n(n - 2)^2(2n^3 - 12n^2 + 23n - 10)/12 if n is even and (n - 1)(n - 3)(2n^4 - 12n^3 + 25n^2 - 14n + 1)/12 if n is odd (Landau, 1896).
a(n) = 5a(n - 1) - 8a(n - 2) + 14a(n - 4) - 14a(n - 5) + 8a(n - 7) - 5a(n - 8) + a(n - 9) for n >= 9.
G.f.: 4(9*x^4 + 35*x^3 + 49*x^2 + 21*x + 6)*x^4/((1 - x)^7*(1 + x)^2).
a(0)=0, a(1)=0, a(2)=0, a(3)=0, a(4)=24, a(5)=204, a(6)=1024, a(7)=3628, a(8)=10320, a(n) = 5*a(n-1)-8*a(n-2)+14*a(n-4)-14*a(n-5)+8*a(n-7)- 5*a(n-8)+ a(n-9). - Harvey P. Dale, Nov 06 2011
a(n) = n^6/6 - 5*n^5/3 + 79*n^4/12 - 25*n^3/2 + 11*n^2 - 43*n/12 + 1/8 + (-1)^n*(n/4 - 1/8) [Chaiken et al.]. - N. J. A. Sloane, Feb 16 2013
a(n) = (3*(2*n-1)*(-1)^n +4*n^6 -40*n^5 +158*n^4 -300*n^3 +264*n^2 -86*n +3)/24. - Antal Pinter, Oct 03 2014
E.g.f.: (exp(2*x)*(3 - 6*x^2 + 8*x^3 + 18*x^4 + 20*x^5 + 4*x^6) -3 - 6*x) / (24*exp(x)). - Vaclav Kotesovec, Feb 15 2015
For n>3, a(n) = A179058(n) -4*(n-2)*A000914(n-2) -2*(n-2)*A002415(n-1) + 2*A008911(n-1) +8*(A001752(n-4) +A007009(n-3)). - Antal Pinter, Sep 20 2015
In general, for m <= n, n >= 3, the number of ways to place 3 nonattacking queens on an m X n board is n^3/6*(m^3 - 3*m^2 + 2*m) - n^2/2*(3*m^3 - 9*m^2 + 6*m) + n/6*(2*m^4 + 20*m^3 - 77*m^2 + 58*m) - 1/24*(39*m^4 - 82*m^3 - 36*m^2 + 88*m) + 1/16*(2*m - 4*n + 1)*(1 + (-1)^(m+1)) + 1/2*(1 + abs(n - 2*m + 3) - abs(n - 2*m + 4))*(1/24*((n - 2*m + 11)^4 - 42*(n - 2*m + 11)^3 + 656*(n - 2*m + 11)^2 - 4518*(n - 2*m + 11) + 11583) - 1/16*(4*m - 2*n - 1)*(1 + (-1)^(n+1))) [Panos Louridas, idee & form 93/2007, pp. 2936-2938]. - Vaclav Kotesovec, Feb 20 2016

The formula given in the Rivin et al. paper is wrong.

Entry improved by comments from Antonio G. Astudillo (afg_astudillo(AT)hotmail.com), May 30 2001

A094216 Triangle read by rows giving the coefficients of formulas generating each variety of S1(n,k) (unsigned Stirling numbers of first kind). The p-th row (p>=1) contains T(i,p) for i=1 to 2p, where T(i,p) satisfies Sum_{i=1..2p} T(i,p) * C(n,i).

Original entry on oeis.org

1, 1, 2, 7, 8, 3, 6, 38, 93, 111, 65, 15, 24, 226, 874, 1821, 2224, 1600, 630, 105, 120, 1524, 8200, 24860, 47185, 58465, 47474, 24430, 7245, 945, 720, 11628, 81080, 326712, 852690, 1522375, 1905168, 1676325, 1018682, 407925, 97020, 10395, 5040

Offset: 1

André F. Labossière, May 27 2004, Feb 21 2007

The formulas S1(n+p,n) obtained are those of S1(n+2,n) { A000914 }, S1(n+3,n) { A001303 }, S1(n+4,n) { A000915 }, S1(n+5,n) { A053567 } and so on.

			Row 5 contains 120,1524,8200,24860,47185,58465,47474,24430,7245,945, so the formula generating S1(n+5,n) numbers { A053567 } will be the following : 120*n +1524*C(n,2) +8200*C(n,3) +24860*C(n,4) +47185*C(n,5) +58465*C(n,6) +47474*C(n,7) +24430*C(n,8) +7245*C(n,9) +945*C(n,10). And then substituting for the 10th number of such a S1(n+p,n) gives S1(15,10) = 37312275.

Milton Abramowitz and Irene A. Stegun, eds., Handbook of Mathematical Functions, National Bureau of Standards Applied Math. Series 55, 1964, 9th Printing (1970), pp. 833-834.

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].
Francis L. Miksa (1901-1975), Stirling numbers of the first kind, "27 leaves reproduced from typewritten manuscript on deposit in the UMT File", Mathematical Tables and Other Aids to Computation, vol. 10, no. 53, January 1956, pp. 37-38 (Reviews and Descriptions of Tables and Books, 7[I]).
Dragoslav S. Mitrinovic (1908-1995), Sur les nombres de Stirling de première espèce et les polynômes de Stirling, AMS 11B73_05A19, Publications de la Faculté d'Electrotechnique de l'Université de Belgrade, Série Mathématiques et Physique (ISSN 0522-8441), no. 23, 1959 (5.V.1959), pp. 1-20.
John J. O'Connor and Edmund F. Robertson, James Stirling (1692-1770), (September 1998).
Eric Weisstein's World of Mathematics, Stirling numbers of the first kind.
Stephen Wolfram, Wolfram Research, Mathematica 5.2, webMathematica 2.

Cf. A000914, A001303, A000915, A053567, A008275, A008276.

Cf. A000012, A000217, A001147, A000142, A094262.

row[m_] := Module[{eq, t}, eq[n_] := Array[t, 2 m].Table[Binomial[n, k], {k, 1, 2 m}] == Abs[StirlingS1[n + m, n]]; Array[t, 2 m] /. Solve[ Array[ eq, 2 m]] // First];
Array[row, 7] // Flatten (* Jean-François Alcover, Nov 14 2019 *)

a(1,k) = k!
...
a(2*k-5,k) = a(2*k,k) * (175000*k^8 -2117500*k^7 +10856650*k^6 -30743377*k^5 +52511770*k^4 -55386931*k^3 +35321832*k^2 -12560580*k+1944000) / (1632960*k^3 -7348320*k^2 +9389520*k -3061800).
a(2*k-4,k) = a(2*k,k) * (2500*k^6 -17400*k^5 +48511*k^4 -69378*k^3 +53929*k^2 -21906*k +3744) / (7776*k^2-15552*k+5832).
a(2*k-3,k) = a(2*k,k) * (1250*k^4-4225*k^3+5023*k^2-2600*k+528) / (1620*k-810).
a(2*k-2,k) = a(2*k,k) * (50*k^3-93*k^2+55*k-12) / (36*k-18).
a(2*k-1,k) = a(2*k,k) * (5*k-2) / 3.
a(2*k,k) = (2*k)! / (k!*2^k).

Showing 1-10 of 54 results. Next

cp's OEIS Frontend

A008275 Triangle read by rows of Stirling numbers of first kind, s(n,k), n >= 1, 1 <= k <= n.

Views

Author

Keywords

Comments

Examples

References

Links

Crossrefs

Programs

Haskell

Maple

Mathematica

Maxima

PARI

PARI

PARI

Formula

A211795 Number of (w,x,y,z) with all terms in {1,...,n} and w*x < 2*y*z.

Views

Author

Keywords

Comments

Examples

References

Links

Crossrefs

Programs

Mathematica

Formula

A264428 Triangle read by rows, Bell transform of Bell numbers.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Maple

Mathematica

PARI

Python

Sage

Formula

A130534 Triangle T(n,k), 0 <= k <= n, read by rows, giving coefficients of the polynomial (x+1)(x+2)...(x+n), expanded in increasing powers of x. T(n,k) is also the unsigned Stirling number |s(n+1, k+1)|, denoting the number of permutations on n+1 elements that contain exactly k+1 cycles.

Views

Author

Keywords

Comments

Examples

References

Links

Crossrefs

Programs

Haskell

Maple

Mathematica

Formula

A094638 Triangle read by rows: T(n,k) = |s(n,n+1-k)|, where s(n,k) are the signed Stirling numbers of the first kind A008276 (1 <= k <= n; in other words, the unsigned Stirling numbers of the first kind in reverse order).

Views

Author

Keywords

Comments

Examples

References

Links

Crossrefs

Programs

GAP

Haskell

Magma

Maple

Mathematica

Maxima

PARI

Sage

Formula

Extensions

A211795 Number of (w,x,y,z) with all terms in {1,...,n} and wx < 2y*z.

A001296 4-dimensional pyramidal numbers: a(n) = (3n+1)binomial(n+2, 3)/4. Also Stirling2(n+2, n).

A094216 Triangle read by rows giving the coefficients of formulas generating each variety of S1(n,k) (unsigned Stirling numbers of first kind). The p-th row (p>=1) contains T(i,p) for i=1 to 2p, where T(i,p) satisfies Sum_{i=1..2p} T(i,p) * C(n,i).