A094645 -id:A094645 - cp's OEIS frontend

A164555 Numerators of the "original" Bernoulli numbers; also the numerators of the Bernoulli polynomials at x=1.

1, 1, 1, 0, -1, 0, 1, 0, -1, 0, 5, 0, -691, 0, 7, 0, -3617, 0, 43867, 0, -174611, 0, 854513, 0, -236364091, 0, 8553103, 0, -23749461029, 0, 8615841276005, 0, -7709321041217, 0, 2577687858367, 0, -26315271553053477373, 0, 2929993913841559, 0, -261082718496449122051

Offset: 0

Paul Curtz, Aug 15 2009

Apart from a sign flip in a(1), the same as A027641.
a(n) is also the numerator of the n-th term of the binomial transform of the sequence of Bernoulli numbers, i.e., of the sequence of fractions A027641(n)/A027642(n).
a(n)/A027642(n) with e.g.f. x/(1-exp(-x)) is the a-sequence for the Sheffer matrix A094645, see the W. Lang link under A006232 for Sheffer a- and z-sequences. - Wolfdieter Lang, Jun 20 2011
a(n)/A027642(n) also give the row sums of the rational triangle of the coefficients of the Bernoulli polynomials A053382/A053383 (falling powers) or A196838/A196839 (rising powers). - Wolfdieter Lang, Oct 25 2011
Given M = the beheaded Pascal's triangle, A074909; with B_n as a vector V, with numerators shown: (1, 1, 1, ...). Then M*V = [1, 2, 3, 4, 5, ...]. If the sign in a(1) is negative in V, then M*V = [1, 0, 0, 0, ...]. - Gary W. Adamson, Mar 09 2012
One might interpret the term "'original' Bernoulli numbers" as numbers given by the e.g.f. x/(1-exp(-x)). - Peter Luschny, Jun 17 2012
Let B(n) = a(n)/A027642(n) then B(n) = Integral_{x=0..1} F_n(x) where F_n(x) are the signed Fubini polynomials F_n(x) = Sum_{k=0..n} (-1)^n*Stirling2(n,k)*k!*(-x)^k (see illustration). - Peter Luschny, Jan 09 2017

			From _Peter Luschny_, Aug 13 2017: (Start)
Integral_{x=0..1} 1 = 1,
Integral_{x=0..1} x = 1/2,
Integral_{x=0..1} 2*x^2 - x = 1/6,
Integral_{x=0..1} 6*x^3 - 6*x^2 + x = 0,
Integral_{x=0..1} 24*x^4 - 36*x^3 + 14*x^2 - x = -1/30,
Integral_{x=0..1} 120*x^5 - 240*x^4 + 150*x^3 - 30*x^2 + x = 0,
...
Integral_{x=0..1} Sum_{k=0..n} (-1)^n*Stirling2(n,k)*k!*(-x)^k = Bernoulli(n). (End)

Jacob Bernoulli, Ars Conjectandi, Basel: Thurneysen Brothers, 1713. See page 97.
John H. Conway and Richard K. Guy, The Book of Numbers, New York: Springer-Verlag, 1996. See pp. 106-108.

Seiichi Manyama, Table of n, a(n) for n = 0..629
Peter Luschny, Illustration of the first terms.
Peter Luschny, The Bernoulli Manifesto, 2013.
Tom Rike, Sums of powers and Bernoulli numbers.

Cf. A027641, A027642, A006232, A053382, A053383, A074909, A094645, A196838, A196839.

Cf. A242179, A106831, A000142.

a164555 n = a164555_list !! n
a164555_list = 1 : map (numerator . sum) (zipWith (zipWith (%))
   (zipWith (map . (*)) (tail a000142_list) a242179_tabf) a106831_tabf)
-- Reinhard Zumkeller, Jul 04 2014

A164555 := proc(n) if n <= 2 then 1; else numer(bernoulli(n)) ; fi; end: # R. J. Mathar, Aug 26 2009
seq(numer(n!*coeff(series(t/(1-exp(-t)),t,n+2),t,n)),n=0..40); # Peter Luschny, Jun 17 2012

CoefficientList[ Series[ x/(1 - Exp[-x]), {x, 0, 40}], x]*Range[0, 40]! // Numerator (* Jean-François Alcover, Mar 04 2013 *)
Table[Numerator[BernoulliB[n,1]], {n, 0, 40}] (* Vaclav Kotesovec, Jan 03 2021 *)

a = lambda n: bernoulli_polynomial(1,n).numerator()
[a(n) for n in (0..40)] # Peter Luschny, Jan 09 2017

a(n) = numerator(B(n)) with B(n) = Sum_{k=0..n} (-1)^(n-k) * C(n+1, k+1) * S(n+k, k) / C(n+k, k) and S the Stirling set numbers. - Peter Luschny, Jun 25 2016
a(n) = numerator(n*EulerPolynomial(n-1, 1)/(2*(2^n-1))) for n>=1. - Peter Luschny, Sep 01 2017
From Artur Jasinski, Jan 01 2021: (Start)
a(n) = numerator(-2*cos(Pi*n/2)*Gamma(n+1)*zeta(n)/(2*Pi)^n) for n != 1.
a(n) = numerator(-n*zeta(1-n)) for n >= 1. In the case n = 0 take the limit. (End)

Edited and extended by R. J. Mathar, Sep 03 2009

Name extended by Peter Luschny, Jan 09 2017

A132393 Triangle of unsigned Stirling numbers of the first kind (see A048994), read by rows, T(n,k) for 0 <= k <= n.

Original entry on oeis.org

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

Offset: 0

Philippe Deléham, Nov 10 2007, Oct 15 2008, Oct 17 2008

Another name: Triangle of signless Stirling numbers of the first kind.
Triangle T(n,k), 0<=k<=n, read by rows given by [0,1,1,2,2,3,3,4,4,5,5,...] DELTA [1,0,1,0,1,0,1,0,1,...] where DELTA is the operator defined in A084938.
A094645*A007318 as infinite lower triangular matrices.
Row sums are the factorial numbers. - Roger L. Bagula, Apr 18 2008
Exponential Riordan array [1/(1-x), log(1/(1-x))]. - Ralf Stephan, Feb 07 2014
Also the Bell transform of the factorial numbers (A000142). For the definition of the Bell transform see A264428 and for cross-references A265606. - Peter Luschny, Dec 31 2015
This is the lower triagonal Sheffer matrix of the associated or Jabotinsky type |S1| = (1, -log(1-x)) (see the W. Lang link under A006232 for the notation and references). This implies the e.g.f.s given below. |S1| is the transition matrix from the monomial basis {x^n} to the rising factorial basis {risefac(x,n)}, n >= 0. - Wolfdieter Lang, Feb 21 2017
T(n, k), for n >= k >= 1, is also the total volume of the n-k dimensional cell (polytope) built from the n-k orthogonal vectors of pairwise different lengths chosen from the set {1, 2, ..., n-1}. See the elementary symmetric function formula for T(n, k) and an example below. - Wolfdieter Lang, May 28 2017
From Wolfdieter Lang, Jul 20 2017: (Start)
The compositional inverse w.r.t. x of y = y(t;x) = x*(1 - t(-log(1-x)/x)) = x + t*log(1-x) is x = x(t;y) = ED(y,t) := Sum_{d>=0} D(d,t)*y^(d+1)/(d+1)!, the e.g.f. of the o.g.f.s D(d,t) = Sum_{m>=0} T(d+m, m)*t^m of the diagonal sequences of the present triangle. See the P. Bala link for a proof (there d = n-1, n >= 1, is the label for the diagonals).
This inversion gives D(d,t) =  P(d, t)/(1-t)^(2*d+1), with the numerator polynomials P(d, t) =  Sum_{m=0..d} A288874(d, m)*t^m. See an example below. See also the P. Bala formula in A112007. (End)
For n > 0, T(n,k) is the number of permutations of the integers from 1 to n which have k visible digits when viewed from a specific end, in the sense that a higher value hides a lower one in a subsequent position. - Ian Duff, Jul 12 2019

			Triangle T(n,k) begins:
  1;
  0,    1;
  0,    1,     1;
  0,    2,     3,     1;
  0,    6,    11,     6,    1;
  0,   24,    50,    35,   10,    1;
  0,  120,   274,   225,   85,   15,   1;
  0,  720,  1764,  1624,  735,  175,  21,  1;
  0, 5040, 13068, 13132, 6769, 1960, 322, 28, 1;
  ...
---------------------------------------------------
Production matrix is
  0, 1
  0, 1, 1
  0, 1, 2,  1
  0, 1, 3,  3,  1
  0, 1, 4,  6,  4,  1
  0, 1, 5, 10, 10,  5,  1
  0, 1, 6, 15, 20, 15,  6, 1
  0, 1, 7, 21, 35, 35, 21, 7, 1
  ...
From _Wolfdieter Lang_, May 09 2017: (Start)
Three term recurrence: 50 = T(5, 2) = 1*6 + (5-1)*11 = 50.
Recurrence from the Sheffer a-sequence [1, 1/2, 1/6, 0, ...]: 50 = T(5, 2) = (5/2)*(binomial(1, 1)*1*6 + binomial(2, 1)*(1/2)*11 + binomial(3, 1)*(1/6)*6 + 0) = 50. The vanishing z-sequence produces the k=0 column from T(0, 0) = 1. (End)
Elementary symmetric function T(4, 2) = sigma^{(3)}_2 = 1*2 + 1*3 + 2*3 = 11. Here the cells (polytopes) are 3 rectangles with total area 11. - _Wolfdieter Lang_, May 28 2017
O.g.f.s of diagonals: d=2 (third diagonal) [0, 6, 50, ...] has D(2,t) = P(2, t)/(1-t)^5, with P(2, t) = 2 + t, the n = 2 row of A288874. - _Wolfdieter Lang_, Jul 20 2017
Boas-Buck recurrence for column k = 2 and n = 5: T(5, 2) = (5!*2/3)*((3/8)*T(2,2)/2! + (5/12)*T(3,2)/3! + (1/2)*T(4,2)/4!) = (5!*2/3)*(3/16 + (5/12)*3/3! + (1/2)*11/4!) = 50. The beta sequence begins: {1/2, 5/12, 3/8, ...}. - _Wolfdieter Lang_, Aug 11 2017

Miklos Bona, editor, Handbook of Enumerative Combinatorics, CRC Press, 2015, pages 31, 187, 441, 996.
R. L. Graham, D. E. Knuth, and O. Patashnik, Concrete Mathematics. Addison-Wesley, Reading, MA, 2nd. ed., Table 259, p. 259.
Steve Roman, The Umbral Calculus, Dover Publications, New York (1984), pp. 149-150

Reinhard Zumkeller, Rows n = 0..125 of triangle, flattened
Roland Bacher and P. De La Harpe, Conjugacy growth series of some infinitely generated groups, hal-01285685v2, 2016.
Eli Bagno and David Garber, Combinatorics of q,r-analogues of Stirling numbers of type B, arXiv:2401.08365 [math.CO], 2024. See page 5.
Peter Bala, Diagonals of triangles with generating function exp(t*F(x)).
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.
Jean-Luc Baril and Sergey Kirgizov, The pure descent statistic on permutations, Preprint, 2016.
Jean-Luc Baril and Sergey Kirgizov, Transformation à la Foata for special kinds of descents and excedances, arXiv:2101.01928 [math.CO], 2021.
Jean-Luc Baril and José L. Ramírez, Some distributions in increasing and flattened permutations, arXiv:2410.15434 [math.CO], 2024. See p. 9.
Ricky X. F. Chen, A Note on the Generating Function for the Stirling Numbers of the First Kind, Journal of Integer Sequences, 18 (2015), #15.3.8.
Bishal Deb and Alan D. Sokal, Higher-order Stirling cycle and subset triangles: Total positivity, continued fractions and real-rootedness, arXiv:2507.18959 [math.CO], 2025. See p. 5.
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
W. Steven Gray and Makhin Thitsa, System Interconnections and Combinatorial Integer Sequences, in: System Theory (SSST), 2013 45th Southeastern Symposium on, Date of Conference: 11-11 March 2013, Digital Object Identifier: 10.1109/SSST.2013.6524939.
John M. Holte, Carries, Combinatorics and an Amazing Matrix, The American Mathematical Monthly, Vol. 104, No. 2 (Feb., 1997), pp. 138-149.
Tanya Khovanova and J. B. Lewis, Skyscraper Numbers, J. Int. Seq. 16 (2013) #13.7.2.
Sergey Kitaev and Philip B. Zhang, Distributions of mesh patterns of short lengths, arXiv:1811.07679 [math.CO], 2018.
Wolfdieter Lang, On Sums of Powers of Arithmetic Progressions, and Generalized Stirling, Eulerian and Bernoulli numbers, arXiv:1707.04451 [math.NT], 2017.
Shuzhen Lv and Philip B. Zhang, Joint equidistributions of mesh patterns 123 and 321 with symmetric and antipodal shadings, arXiv:2501.00357 [math.CO], 2024. See p. 11.
Shuzhen Lv and Philip B. Zhang, Joint equidistributions of mesh patterns 123 and 132 with antipodal shadings, arXiv:2506.23148 [math.CO], 2025. See p. 6.
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
Mathematics Stack Exchange, Symmetric (under the swapping) recursions for Stirling numbers of both kinds, Aug 10 2025.
Emanuele Munarini, Shifting Property for Riordan, Sheffer and Connection Constants Matrices, Journal of Integer Sequences, Vol. 20 (2017), Article 17.8.2.
Emanuele Munarini, Combinatorial identities involving the central coefficients of a Sheffer matrix, Applicable Analysis and Discrete Mathematics (2019) Vol. 13, 495-517.
Fedor Petrov, Recursive algorithm for columns of Stirling numbers of the first kind, answer to question on MathOverflow (2025).
E. G. Santos, Counting non-attacking chess pieces placements: Bishops and Anassas, arXiv:2411.16492 [math.CO], 2024. See p. 2.
Umesh Shankar, Log-concavity of rows of triangular arrays satisfying a certain super-recurrence, arXiv:2508.12467 [math.CO], 2025. See p. 4.
Xun-Tuan Su, Deng-Yun Yang, and Wei-Wei Zhang, A note on the generalized factorial, Australasian Journal of Combinatorics, Volume 56 (2013), Pages 133-137.
Benjamin Testart, Completing the enumeration of inversion sequences avoiding one or two patterns of length 3, arXiv:2407.07701 [math.CO], 2024. See p. 37.

Essentially a duplicate of A048994. Cf. A008275, A008277, A112007, A130534, A288874, A354795.

a132393 n k = a132393_tabl !! n !! k
a132393_row n = a132393_tabl !! n
a132393_tabl = map (map abs) a048994_tabl
-- Reinhard Zumkeller, Nov 06 2013

a132393_row := proc(n) local k; seq(coeff(expand(pochhammer (x,n)),x,k),k=0..n) end: # Peter Luschny, Nov 28 2010

p[t_] = 1/(1 - t)^x; Table[ ExpandAll[(n!)SeriesCoefficient[ Series[p[t], {t, 0, 30}], n]], {n, 0, 10}]; a = Table[(n!)* CoefficientList[SeriesCoefficient[ Series[p[t], {t, 0, 30}], n], x], {n, 0, 10}]; Flatten[a] (* Roger L. Bagula, Apr 18 2008 *)
Flatten[Table[Abs[StirlingS1[n,i]],{n,0,10},{i,0,n}]] (* Harvey P. Dale, Feb 04 2014 *)

create_list(abs(stirling1(n,k)),n,0,12,k,0,n); /* Emanuele Munarini, Mar 11 2011 */

column(n,k) = my(v1, v2); v1 = vector(n-1, i, 0); v2 = vector(n, i, 0); v2[1] = 1; for(i=1, n-1, v1[i] = (i+k)*(i+k-1)/2*v2[i]; for(j=1, i-1, v1[j] *= (i-j)*(i+k)/(i-j+2)); v2[i+1] = vecsum(v1)/i); v2 \\ generates n first elements of the k-th column starting from the first nonzero element. - Mikhail Kurkov, Mar 05 2025

T(n,k) = T(n-1,k-1)+(n-1)*T(n-1,k), n,k>=1; T(n,0)=T(0,k); T(0,0)=1.
Sum_{k=0..n} T(n,k)*x^(n-k) = A000012(n), A000142(n), A001147(n), A007559(n), A007696(n), A008548(n), A008542(n), A045754(n), A045755(n) for x = 0, 1, 2, 3, 4, 5, 6, 7, 8 respectively. - Philippe Deléham, Nov 13 2007
Expand 1/(1-t)^x = Sum_{n>=0}p(x,n)*t^n/n!; then the coefficients of the p(x,n) produce the triangle. - Roger L. Bagula, Apr 18 2008
Sum_{k=0..n} T(n,k)*2^k*x^(n-k) = A000142(n+1), A000165(n), A008544(n), A001813(n), A047055(n), A047657(n), A084947(n), A084948(n), A084949(n) for x = 1, 2, 3, 4, 5, 6, 7, 8, 9 respectively. - Philippe Deléham, Sep 18 2008
a(n) = Sum_{k=0..n} T(n,k)*3^k*x^(n-k) = A001710(n+2), A001147(n+1), A032031(n), A008545(n), A047056(n), A011781(n), A144739(n), A144756(n), A144758(n) for x=1,2,3,4,5,6,7,8,9,respectively. - Philippe Deléham, Sep 20 2008
Sum_{k=0..n} T(n,k)*4^k*x^(n-k) = A001715(n+3), A002866(n+1), A007559(n+1), A047053(n), A008546(n), A049308(n), A144827(n), A144828(n), A144829(n) for x=1,2,3,4,5,6,7,8,9 respectively. - Philippe Deléham, Sep 21 2008
Sum_{k=0..n} x^k*T(n,k) = x*(1+x)*(2+x)*...*(n-1+x), n>=1. - Philippe Deléham, Oct 17 2008
From Wolfdieter Lang, Feb 21 2017: (Start)
E.g.f. k-th column: (-log(1 - x))^k, k >= 0.
E.g.f. triangle (see the Apr 18 2008 Baluga comment): exp(-x*log(1-z)).
E.g.f. a-sequence: x/(1 - exp(-x)). See A164555/A027642. The e.g.f. for the z-sequence is 0. (End)
From Wolfdieter Lang, May 28 2017: (Start)
The row polynomials R(n, x) = Sum_{k=0..n} T(n, k)*x^k, for n >= 0, are R(n, x) = risefac(x,n-1) := Product_{j=0..n-1} x+j, with the empty product for n=0 put to 1. See the Feb 21 2017 comment above. This implies:
T(n, k) = sigma^{(n-1)}_(n-k), for n >= k >= 1, with the elementary symmetric functions sigma^{(n-1)}_m of degree m in the n-1 symbols 1, 2, ..., n-1, with binomial(n-1, m) terms. See an example below.(End)
Boas-Buck type recurrence for column sequence k: T(n, k) = (n!*k/(n - k)) * Sum_{p=k..n-1} beta(n-1-p)*T(p, k)/p!, for n > k >= 0, with input T(k, k) = 1, and beta(k) = A002208(k+1)/A002209(k+1). See a comment and references in A286718. - Wolfdieter Lang, Aug 11 2017
T(n,k) = Sum_{j=k..n} j^(j-k)*binomial(j-1, k-1)*A354795(n,j) for n > 0. - Mélika Tebni, Mar 02 2023
n-th row polynomial: n!*Sum_{k = 0..2*n} (-1)^k*binomial(-x, k)*binomial(-x, 2*n-k) = n!*Sum_{k = 0..2*n} (-1)^k*binomial(1-x, k)*binomial(-x, 2*n-k). - Peter Bala, Mar 31 2024
From Mikhail Kurkov, Mar 05 2025: (Start)
For a general proof of the formulas below via generating functions, see Mathematics Stack Exchange link.
Recursion for the n-th row (independently of other rows): T(n,k) = 1/(n-k)*Sum_{j=2..n-k+1} binomial(-k,j)*T(n,k+j-1)*(-1)^j for 1 <= k < n with T(n,n) = 1.
Recursion for the k-th column (independently of other columns): T(n,k) = 1/(n-k)*Sum_{j=2..n-k+1} (j-2)!*binomial(n,j)*T(n-j+1,k) for 1 <= k < n with T(n,n) = 1 (see Fedor Petrov link). (End)

A049444 Generalized Stirling number triangle of first kind.

Original entry on oeis.org

1, -2, 1, 6, -5, 1, -24, 26, -9, 1, 120, -154, 71, -14, 1, -720, 1044, -580, 155, -20, 1, 5040, -8028, 5104, -1665, 295, -27, 1, -40320, 69264, -48860, 18424, -4025, 511, -35, 1, 362880, -663696, 509004, -214676, 54649, -8624, 826, -44, 1, -3628800, 6999840, -5753736

Offset: 0

Wolfdieter Lang

T(n, k) = ^2P_n^k in the notation of the given reference with T(0, 0) := 1. The monic row polynomials s(n,x) := Sum_{m=0..n} T(n, k)*x^k which are s(n, x) = Product_{j=0..n-1} (x-(2+j)), n >= 1 and s(0, x)=1 satisfy s(n, x+y) = Sum_{k=0..n} binomial(n, k)*s(k,x)*S1(n-k, y), with the Stirling1 polynomials S1(n, x) = Sum_{m=1..n} (A008275(n, m)*x^m) and S1(0, x)=1.
In the umbral calculus (see the S. Roman reference given in A048854) the s(n, x) polynomials are called Sheffer polynomials for (exp(2*t), exp(t)-1). This translates to the usual exponential Riordan (Sheffer) notation (1/(1+x)^2, log(1+x)).
See A143491 for the unsigned version of this array and A143494 for the inverse. - Peter Bala, Aug 25 2008
Corresponding to the generalized Stirling number triangle of second kind A137650. - Peter Luschny, Sep 18 2011
Unsigned, reversed rows (cf. A145324, A136124) are the dimensions of the cohomology of a complex manifold with a symmetric group (S_n) action. See p. 17 of the Hyde and Lagarias link. See also the Murri link for an interpretation as the Betti numbers of the moduli space M(0,n) of smooth Riemann surfaces. - Tom Copeland, Dec 09 2016
The row polynomials s(n, x) = (-1)^n*risingfactorial(2 - x, n) are related to the column sequences of the unsigned Abel triangle A137452(n, k), for k >= 2. See the formula there. - Wolfdieter Lang, Nov 21 2022

			The Triangle  begins:
n\k       0       1        2       3       4      5      6    7   8 9 ...
0:        1
1:       -2       1
2:        6      -5        1
3:      -24      26       -9       1
4:      120    -154       71     -14       1
5      -720    1044     -580     155     -20      1
6:     5040   -8028     5104   -1665     295    -27      1
7:   -40320   69264   -48860   18424   -4025    511    -35    1
8:   362880 -663696   509004 -214676   54649  -8624    826  -44
9: -3628800 6999840 -5753736 2655764 -761166 140889 -16884 1266 -54 1
...  [reformatted by _Wolfdieter Lang_, Nov 21 2022]

Y. Manin, Frobenius Manifolds, Quantum Cohomology and Moduli Spaces, American Math. Soc. Colloquium Publications Vol. 47, 1999. [From Tom Copeland, Jun 29 2008]
S. Roman, The Umbral Calculus, Academic Press, 1984 (also Dover Publications, 2005).

Reinhard Zumkeller, Rows n = 0..125 of triangle, flattened
E. Getzler, Operads and moduli spaces of genus 0 Riemann surfaces, arXiv:alg-geom/9411004, 1994, (see p. 23, g(x,t)). [From _Tom Copeland_, Dec 11 2011]
T. Hyde and J. Lagarias Polynomial splitting measures and cohomology of the pure braid group, arXiv preprint arXiv:1604.05359 [math.RT], 2016.
Y. Manin, Generating functions in algebraic geometry and sums over trees, arXiv:alg-geom/9407005, 1994, (Eqn. 0.7 and 1.7). [From _Tom Copeland_, Dec 10 2011]
D. S. Mitrinovic and M. S. Mitrinovic, Tableaux d'une classe de nombres reliés aux nombres de Stirling, Univ. Beograd. Publ. Elektrotehn. Fak. Ser. Mat. Fiz. No. 77 (1962), 1-77.
R. Murri, Fatgraph Algorithms and the Homology of the Kontsevich Complex, arXiv:1202.1820 [math.AG], 2012, (see Table 1, p. 3). [From _Tom Copeland_, Sep 18 2012]

Unsigned column sequences are A000142(n+1), A001705-A001709. Row sums (signed triangle): n!*(-1)^n, row sums (unsigned triangle): A001710(n-2). Cf. A008275 (Stirling1 triangle).
Cf. A000035, A084938, A094645, A094646.
Cf. A143491, A143494.
Cf. A136124, A137452, A137650.

a049444 n k = a049444_tabl !! n !! k
a049444_row n = a049444_tabl !! n
a049444_tabl = map fst $ iterate (\(row, i) ->
   (zipWith (-) ([0] ++ row) $ map (* i) (row ++ [0]), i + 1)) ([1], 2)
-- Reinhard Zumkeller, Mar 11 2014

A049444_row := proc(n) local k,i;
add(add(Stirling1(n, n-i), i=0..k)*x^(n-k-1),k=0..n-1);
seq(coeff(%,x,k),k=1..n-1) end:
seq(print(A049444_row(n)),n=1..7); # Peter Luschny, Sep 18 2011
A049444:= (n, k)-> add((-1)^(n-j)*(n-j+1)!*binomial(n, j)*Stirling1(j, k), j=0..n):
seq(print(seq(A049444(n, k), k=0..n)), n=0..11);  # Mélika Tebni, May 02 2022

t[n_, i_] = Sum[(-1)^k*Binomial[n, k]*(k+1)!*StirlingS1[n-k, i], {k, 0, n-i}]; Flatten[Table[t[n, i], {n, 0, 9}, {i, 0, n}]] [[1 ;; 48]]
(* Jean-François Alcover, Apr 29 2011, after Milan Janjic *)

T(n, k) = T(n-1, k-1) - (n+1)*T(n-1, k), n >= k >= 0; T(n, k) = 0, n < k; T(n, -1) = 0, T(0, 0) = 1.
E.g.f. for k-th column of signed triangle: ((log(1+x))^k)/(k!*(1+x)^2).
Triangle (signed) = [-2, -1, -3, -2, -4, -3, -5, -4, -6, -5, ...] DELTA [1, 0, 1, 0, 1, 0, 1, 0, ...]; triangle (unsigned) = [2, 1, 3, 2, 4, 3, 5, 4, 6, 5, ...] DELTA [1, 0, 1, 0, 1, 0, 1, 0, 1, 0, ...], where DELTA is Deléham's operator defined in A084938 (unsigned version in A143491).
E.g.f.: (1 + x)^(y-2). - Vladeta Jovovic, May 17 2004 [For row polynomials s(n, y)]
With P(n, t) = Sum_{j=0..n-2} T(n-2,j) * t^j and P(1, t) = -1 and P(0, t) = 1, then G(x, t) = -1 + exp[P(.,t)*x] = [(1+x)^t - 1 - t^2 * x] / [t(t-1)], whose compositional inverse in x about 0 is given in A074060. G(x, 0) = -log(1+x) and G(x, 1) = (1+x) log(1+x) - 2x. G(x, q^2) occurs in formulas on pages 194-196 of the Manin reference. - Tom Copeland, Feb 17 2008
If we define f(n,i,a) = Sum_{k=0..n-i} binomial(n,k)*Stirling1(n-k,i)*Product_{j=0..k-1} (-a-j), then T(n,i) = f(n,i,2), for n=1,2,...; i=0..n. - Milan Janjic, Dec 21 2008
T(n, k) = Sum_{j=0..n} (-1)^(n-j)*(n-j+1)!*binomial(n, j)*Stirling1(j, k). - Mélika Tebni, May 02 2022
From Wolfdieter Lang, Nov 24 2022: (Start)
Recurrence for row polynomials {s(n, x)}_{n>=0}: s(0, x) = 1, s(n, x) = (x - 2)*exp(-(d/dx)) s(n-1, x), for n >= 1. This is adapted from the general Sheffer result given by S. Roman, Corollary 3.7.2., p. 50.
Recurrence for column sequence {T(n, k)}{n>=k}: T(n, n) = 1, T(n, k) = (n!/(n-k))*Sum{j=k..n-1} (1/j!)*(a(n-1-j) + k*beta(n-1-j))*T(n-1, k), for k >= 0, where alpha = repeat(-2, 2) and beta(n) = [x^n] (d/dx)log(log(x)/x) = (-1)^(n+1)*A002208(n+1)/A002209(n+1), for n >= 0. This is the adapted Boas-Buck recurrence, also given in Rainville, Theorem 50., p. 141, For the references and a comment see A046521. (End)

Second formula corrected by Philippe Deléham, Nov 09 2008

A105794 Inverse of a generalized Stirling number triangle of first kind.

Original entry on oeis.org

1, -1, 1, 1, -1, 1, -1, 1, 0, 1, 1, -1, 1, 2, 1, -1, 1, 0, 5, 5, 1, 1, -1, 1, 10, 20, 9, 1, -1, 1, 0, 21, 70, 56, 14, 1, 1, -1, 1, 42, 231, 294, 126, 20, 1, -1, 1, 0, 85, 735, 1407, 924, 246, 27, 1, 1, -1, 1, 170, 2290, 6363, 6027, 2400, 435, 35, 1

Offset: 0

Paul Barry, Apr 20 2005

Inverse of number triangle A105793. Row sums are the generalized Bell numbers A000296.
T(n,k)*(-1)^(n-k) gives the inverse Sheffer matrix of A094645. In the umbral notation (cf. Roman, p. 17, quoted in A094645) this is Sheffer for (1-t,-log(1-t)). - Wolfdieter Lang, Jun 20 2011
From Peter Bala, Jul 10 2013: (Start)
For a graph G and a positive integer k, the graphical Stirling number S(G;k) is the number of partitions of the vertex set of G into k nonempty independent sets (an independent set in G is a subset of the vertices, no two elements of which are adjacent). Omitting the first two rows and columns of this triangle produces the triangle of graphical Stirling numbers of cycles on n vertices (interpreting a 2-cycle as a single edge). In comparison, the classical Stirling numbers of the second kind, A008277, are the graphical Stirling numbers of path graphs on n vertices. See Galvin and Than. See also A227341.
This is the triangle of connection constants for expressing the polynomial sequence (x-1)^n in terms of the falling factorial polynomials, that is, (x-1)^n = Sum_{k = 0..n} T(n,k)*x_(k), where x_(k) = x*(x-1)*...*(x-k+1) denotes the falling factorial.
The row polynomials are a particular case of the Actuarial polynomials - see Roman 4.3.4.  (End)

			The triangle starts with
  n=0:  1;
  n=1: -1,  1;
  n=2:  1, -1, 1;
  n=3: -1,  1, 0, 1;
  n=4:  1, -1, 1, 2, 1;
  n=5: -1,  1, 0, 5, 5, 1;
  ... - _Wolfdieter Lang_, Jun 20 2011

S. Roman, The umbral calculus, Pure and Applied Mathematics 111, Academic Press Inc., New York, 1984. Reprinted by Dover in 2005.

Robert Israel, Table of n, a(n) for n = 0..9869
Peter Bala, Notes on A105794
Peter Bala, A 3 parameter family of generalized Stirling numbers
B. Duncan and R. Peele, Bell and Stirling Numbers for Graphs, Journal of Integer Sequences 12 (2009), article 09.7.1.
D. Galvin and D. T. Thanh, Stirling numbers of forests and cycles, Electr. J. Comb. Vol. 20(1): P73 (2013).
Sophie Morier-Genoud, Counting Coxeter's friezes over a finite field via moduli spaces, arXiv:1907.12790 [math.CO], 2019.

Cf. A000296 (row sums), A105793 (matrix inverse), A227341.

Cf. A007318, A008277, A048993.

B:= Matrix(12,12,shape=triangular[lower],(n,k) -> combinat:-stirling1(n-1,k-1)+(n-1)*combinat:-stirling1(n-2,k-1)):
A:= B^(-1):
seq(seq(A[i,j],j=1..i),i=1..12); # Robert Israel, Jan 19 2015
T := (n, k) -> add((-1)^(n - i)*binomial(n, i)*Stirling2(i, k), i=0..n):
seq(seq(T(n, k), k=0..n), n=0..9);  # Peter Luschny, Feb 15 2025

Table[Sum[(-1)^(n - i)*Binomial[n, i] StirlingS2[i, k], {i, 0, n}], {n, 0, 10}, {k, 0, n}] // Flatten (* Michael De Vlieger, Oct 14 2019 *)

Term k in row n is given by {(-1)^(k+n) * (Sum_{j=0..k} (-1)^j * binomial(k,j) * (1-j)^n) / k! }; i.e., a finite difference. - Tom Copeland, Jun 05 2006
O.g.f. for row n = n-th finite difference of the Touchard (Bell) polynomials, T(x,j) and so the e.g.f. for these finite differences and therefore the sequence = exp{x*[exp(t)-1]-t} = exp{t*[T(x,.)-1]} umbrally. - Tom Copeland, Jun 05 2006
The e.g.f. A(x,t) = exp(x*(exp(t)-1)-t) satisfies the partial differential equation x*dA/dx - dA/dt = (1-x)*A.
Recurrence relation: T(n+1,k) = T(n,k-1) + (k-1)*T(n,k).
Let f(x) = ((x-1)/x)*exp(x). For n >= 1, the n-th row polynomial R(n,x) = x*exp(-x)*(x*d/dx)^(n-1)(f(x)) and satisfies the recurrence equation R(n+1,x) = (x-1)*R(n,x) + x*R'(n,x). - Peter Bala, Oct 28 2011
Let f(x) = exp(exp(x)-x). Then R(n,exp(x)) = 1/f(x)*(d/dx)^n(f(x)). Similar formulas hold for A008277, A039755, A111577, A143494 and A154537. - Peter Bala, Mar 01 2012
From Peter Bala, Jul 10 2013: (Start)
T(n,k) = Sum_{i = 0..n-1} (-1)^i*Stirling2(n-1-i,k-1), for n >= 1, k >= 1.
The k-th column o.g.f. is (1/(1+x))*x^k/((1-x)*(1-2*x)*...*(1-(k-1)*x)) (the empty product occurring in the denominator when k = 0 and k = 1 is taken equal to 1).
Dobinski-type formula for the row polynomials: R(n,x) = exp(-x)*Sum_{k >= 0} (k-1)^n*x^k/k!.
Sum_{k = 0..n} binomial(n,k)*R(k,x) = n-th Bell polynomial (n-th row polynomial of A048993). (End)
From Peter Bala, Jan 13 2014: (Start)
T(n,k) = Sum_{i = 0..n} (-1)^(n-i)*binomial(n,i)*Stirling2(i,k).
T(n,k) = Sum_{i = 0..n} (-2)^(n-i)*binomial(n,i)*Stirling2(i+1,k+1).
Matrix product P^(-1) * S = P^(-2) * S1, where P = A007318, S = A048993 and S1 = A008277. (End)
From Werner Schulte, Feb 15 2025: (Start)
Conjecture 1: Sum_{k=0..n} (-1)^(n-k) * T(n, k) * (k+m)! / m! = (m+2)^n for m >= 0.
Conjecture 2: (-1)^n - B(n) = Sum_{k=1..n} (-1)^k * T(n, k) * (k-1)! / (k+1) where B(n) = B(n, 0) is n-th Bernoulli number. (End)

A049458 Generalized Stirling number triangle of first kind.

Original entry on oeis.org

1, -3, 1, 12, -7, 1, -60, 47, -12, 1, 360, -342, 119, -18, 1, -2520, 2754, -1175, 245, -25, 1, 20160, -24552, 12154, -3135, 445, -33, 1, -181440, 241128, -133938, 40369, -7140, 742, -42, 1, 1814400, -2592720, 1580508, -537628

Offset: 0

Wolfdieter Lang

a(n,m)= ^3P_n^m in the notation of the given reference with a(0,0) := 1. The monic row polynomials s(n,x) := sum(a(n,m)*x^m,m=0..n) which are s(n,x)= product(x-(3+k),k=0..n-1), n >= 1 and s(0,x)=1 satisfy s(n,x+y) = sum(binomial(n,k)*s(k,x)*S1(n-k,y),k=0..n), with the Stirling1 polynomials S1(n,x)=sum(A008275(n,m)*x^m, m=1..n) and S1(0,x)=1.
In the umbral calculus (see the S. Roman reference given in A048854) the s(n,x) polynomials are called Sheffer polynomials for (exp(3*t),exp(t)-1).
See A143492 for the unsigned version of this array and A143495 for the inverse. - Peter Bala, Aug 25 2008

			1;
-3, 1;
12, -7, 1;
-60, 47, -12, 1;
360, -342, 119, -18, 1;
s(2,x) = 12-7*x+x^2. S1(2,x) = -x+x^2 (Stirling1 polynomial).

Mitrinovic, D. S.; Mitrinovic, R. S.; Tableaux d'une classe de nombres relies aux nombres de Stirling. Univ. Beograd. Pubi. Elektrotehn. Fak. Ser. Mat. Fiz. No. 77 1962, 77 pp.

Reinhard Zumkeller, Rows n = 0..125 of triangle, flattened

Unsigned column sequences are: A001710-A001714. Row sums (signed triangle): (n+1)!*(-1)^n. Row sums (unsigned triangle): A001715(n+3).
Cf. A000035 A084938.
Cf. A094645, A094646.
A143492, A143495. - Peter Bala, Aug 25 2008

a049458 n k = a049458_tabl !! n !! k
a049458_row n = a049458_tabl !! n
a049458_tabl = map fst $ iterate (\(row, i) ->
   (zipWith (-) ([0] ++ row) $ map (* i) (row ++ [0]), i + 1)) ([1], 3)
-- Reinhard Zumkeller, Mar 11 2014

A049458_row := n -> seq((-1)^(n-k)*coeff(expand(pochhammer(x+3, n)), x, k), k=0..n): seq(print(A049458_row(n)),n=0..8); # Peter Luschny, May 16 2013

t[n_, k_] := (-1)^(n - k)*Coefficient[ Pochhammer[x + 3, n], x, k]; Table[t[n, k], {n, 0, 8}, {k, 0, n}] // Flatten (* Jean-François Alcover, Jul 17 2013, after Peter Luschny *)

a(n, m)= a(n-1, m-1) - (n+2)*a(n-1, m), n >= m >= 0; a(n, m) := 0, n

Triangle (signed) = [ -3, -1, -4, -2, -5, -3, -6, -4, -7, ...] DELTA [1, 0, 1, 0, 1, 0, 1, 0, ...]; triangle (unsigned) = [3, 1, 4, 2, 5, 3, 6, 4, 7, 5, ...] DELTA [1, 0, 1, 0, 1, 0, 1, 0, ...]; where DELTA is Deléham's operator defined in A084938 (unsigned version in A143492).

E.g.f.: (1+y)^(x-3). - Vladeta Jovovic, May 17 2004

If we define f(n,i,a)=sum(binomial(n,k)*stirling1(n-k,i)*product(-a-j,j=0..k-1),k=0..n-i), then T(n,i) = f(n,i,3), for n=1,2,...;i=0...n. - Milan Janjic, Dec 21 2008

Second formula corrected by Philippe Deléham, Nov 09 2008

A105954 Array read by descending antidiagonals: A(n, k) = (n + 1)! * H(k, n + 1), where H(n, k) is a higher-order harmonic number, H(0, k) = 1/k and H(n, k) = Sum_{j=1..k} H(n-1, j), for 0 <= k <= n.

Original entry on oeis.org

1, 1, 1, 1, 3, 2, 1, 5, 11, 6, 1, 7, 26, 50, 24, 1, 9, 47, 154, 274, 120, 1, 11, 74, 342, 1044, 1764, 720, 1, 13, 107, 638, 2754, 8028, 13068, 5040, 1, 15, 146, 1066, 5944, 24552, 69264, 109584, 40320, 1, 17, 191, 1650, 11274, 60216, 241128, 663696, 1026576, 362880

Offset: 0

Leroy Quet, Jun 26 2005

Antidiagonal sums are A093345 (n! * (1 + Sum_{i=1..n}((1/i)*Sum_{j=0..i-1} 1/j!))). - Gerald McGarvey, Aug 27 2005
A recasting of A093905 and A067176. - R. J. Mathar, Mar 01 2009
The triangular array of this sequence is the reversal of A165675 which is related to the asymptotic expansion of the higher order exponential integral E(x,m=2,n); see also A165674. - Johannes W. Meijer, Oct 16 2009

			A(2, 2) = (1 + (1 + 1/2) + (1 + 1/2 + 1/3))*6 = 26.
Array A(n, k) begins:
  [n\k]  0       1       2        3        4        5          6
  -------------------------------------------------------------------
  [0]    1,      1,      1,       1,       1,       1,         1, ...
  [1]    1,      3,      5,       7,       9,       11,       13, ...
  [2]    2,     11,     26,      47,      74,      107,      146, ...
  [3]    6,     50,    154,     342,     638,     1066,     1650, ...
  [4]   24,    274,   1044,    2754,    5944,    11274,    19524, ...
  [5]  120,   1764,   8028,   24552,   60216,   127860,   245004, ...
  [6]  720,  13068,  69264,  241128,  662640,  1557660,  3272688, ...
  [7] 5040, 109584, 663696, 2592720, 7893840, 20355120, 46536624, ...

G. C. Greubel, Table of n, a(n) for the first 27 rows, flattened
Arthur T. Benjamin, David Gaebler and Robert Gaebler, A Combinatorial Approach to Hyperharmonic Numbers, INTEGERS, Electronic Journal of Combinatorial Number Theory, Volum 3, #A15, 2003.

Column 0 = A000142 (factorial numbers).
Column 1 = A000254 (Stirling numbers of first kind s(n, 2)) starting at n=1.
Column 2 = A001705 (Generalized Stirling numbers: a(n) = n!*Sum_{k=0..n-1}(k+1)/(n-k)), starting at n=1.
Column 3 = A001711 (Generalized Stirling numbers: a(n) = Sum_{k=0..n}(-1)^(n+k)*(k+1)*3^k*stirling1(n+1, k+1)).
Column 4 = A001716 (Generalized Stirling numbers: a(n) = Sum_{k=0..n}(-1)^(n+k)*(k+1)*4^k*stirling1(n+1, k+1)).
Column 5 = A001721 (Generalized Stirling numbers: a(n) = Sum_{k=0..n}(-1)^(n+k)*binomial(k+1, 1)*5^k*stirling1(n+1, k+1)).
Column 6 = A051524 (2nd unsigned column of A051338) starting at n=1.
Column 7 = A051545 (2nd unsigned column of A051339) starting at n=1.
Column 8 = A051560 (2nd unsigned column of A051379) starting at n=1.
Column 9 = A051562 (2nd unsigned column of A051380) starting at n=1.
Column 10= A051564 (2nd unsigned column of A051523) starting at n=1.
2nd row is A005408 (2n - 1, starting at n=1).
3rd row is A080663 (3n^2 - 1, starting at n=1).
Main diagonal gives A384024.
Cf. A000254, A165674 and A165675, A028421 and A126671.

H := proc(n, k) option remember; if n = 0 then 1/k else add(H(n - 1, j), j = 1..k) fi end: A := (n, k) -> (n + 1)!*H(k, n + 1):
# Alternative with standard harmonic number:
A := (n, k) -> if k = 0 then n! else (harmonic(n + k) - harmonic(k - 1))*(n + k)! / (k - 1)! fi:
for n from 0 to 7 do seq(A(n, k), k = 0..6) od;
# Alternative with hypergeometric formula:
A := (n, k) -> (n+1)*((n + k)! / k!)*hypergeom([-n, 1, 1], [2, k+1], 1):
seq(print(seq(simplify(A(n, k)), k = 0..6)), n=0..7); # Peter Luschny, Jul 01 2022

H[0, m_] := 1/m; H[n_, m_] := Sum[H[n - 1, k], {k, m}]; a[n_, m_] := m!H[n, m]; Flatten[ Table[ a[i, n - i], {n, 10}, {i, n - 1, 0, -1}]]
Table[ a[n, m], {m, 8}, {n, 0, m + 1}] // TableForm (* to view the table *)
(* Robert G. Wilson v, Jun 27 2005 *)

a(n, k) = polcoef(prod(j=0, n, 1+(j+k)*x), n); \\ Seiichi Manyama, May 19 2025

A(n, k) = (Harmonic(n + k) - Harmonic(k - 1))*(n + k)!/(k - 1)! if k > 0, otherwise n!.
From Gerald McGarvey, Aug 27 2005, edited by Peter Luschny, Jul 02 2022: (Start)
E.g.f. for column k: -log(1 - x)/(x*(1 - x)^k).
Row 3 is r(n) = 4*n^3 + 18*n^2 + 22*n + 6.
Row 4 is r(n) = 5*n^4 + 40*n^3 + 105*n^2 + 100*n + 24.
Row 5 is r(n) = 6*n^5 + 75*n^4 + 340*n^3 + 675*n^2 + 548*n + 120.
Row 6 is r(n) = 7*n^6 + 126*n^5 + 875*n^4 + 2940*n^3 + 4872*n^2 + 3528*n + 720.
Row 7 is r(n) = 8*n^7 + 196*n^6 + 1932*n^5 + 9800*n^4 + 27076*n^3 + 39396*n^2 + 26136*n + 5040.
The sum of the polynomial coefficients for the n-th row is |S1(n, 2)|, which are the unsigned Stirling1 numbers which appear in column 1.
A(m, n) = Sum_{k=1..m} n*A094645(m, n)*(n+1)^(k-1). (A094645 is Generalized Stirling number triangle of first kind, e.g.f.: (1-y)^(1-x).) (End)
In Gerard McGarvey's formulas for the row coefficients we find Wiggen's triangle A028421 and their o.g.f.s lead to Wood's polynomials A126671; see A165674. - Johannes W. Meijer, Oct 16 2009
A(n, k) = (n + 1)*((n + k)! / k!)*hypergeom([-n, 1, 1], [2, k + 1], 1). - Peter Luschny, Jul 01 2022
A(n,k) = [x^n] Product_{j=0..n} (1 + (j+k)*x). - Seiichi Manyama, May 19 2025

More terms from Robert G. Wilson v, Jun 27 2005

Edited by Peter Luschny, Jul 02 2022

A081047 Difference of Stirling numbers of the first kind.

Original entry on oeis.org

1, 0, -1, -5, -26, -154, -1044, -8028, -69264, -663696, -6999840, -80627040, -1007441280, -13575738240, -196287356160, -3031488633600, -49811492505600, -867718162483200, -15974614352793600, -309920046408806400, -6320046028584960000

Offset: 0

Paul Barry, Mar 05 2003

G. C. Greubel, Table of n, a(n) for n = 0..400
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, Volume 43, 2012 - Issue 2.
Motohico Mulase, In Search of a Hidden Curve, arXiv:2501.00716 [math.QA], 2025. See p. 27.

Cf. A001705, A008275, A081046.

With[{nn = 100}, CoefficientList[Series[(1 + Log[1 - x])/(1 - x), {x, 0, nn}], x] Range[0, nn]!] (* G. C. Greubel, Jan 21 2017 *)

E.g.f.: (1+log(1-x))/(1-x). - Paul Barry, Nov 26 2008
a(n) = abs(s(n+1, 1))-abs(s(n+1, 2)), where s(n, m) is a (signed) Stirling number of the first kind (A008275). (corrected by Wolfdieter Lang, Jun 20 2011)
a(n) = A094645(n+2,2), n>=0. - _Wolfdieter Lang, Jun 20 2011

A094646 Generalized Stirling number triangle of first kind.

Original entry on oeis.org

1, -2, 1, 2, -3, 1, 0, 2, -3, 1, 0, 2, -1, -2, 1, 0, 4, 0, -5, 0, 1, 0, 12, 4, -15, -5, 3, 1, 0, 48, 28, -56, -35, 7, 7, 1, 0, 240, 188, -252, -231, 0, 42, 12, 1, 0, 1440, 1368, -1324, -1638, -231, 252, 114, 18, 1, 0, 10080, 11016, -7900, -12790, -3255, 1533, 1050, 240, 25, 1

Offset: 0

Vladeta Jovovic, May 17 2004

Triangle T(n,k), 0 <= k <= n, read by rows, given by [ -2, 1, -1, 2, 0, 3, 1, 4, 2, 5, ...] DELTA [1, 0, 1, 0, 1, 0, 1, 0, 1, 0, ...] where DELTA is the operator defined in A084938. - Philippe Deléham, Aug 23 2006
From Wolfdieter Lang, Jun 23 2011: (Start)
The row polynomials s(n,x):=Sum_{k=0..n} T(n,k)*x^k satisfy risefac(x-2,n)=s(n,x), with the rising factorials risefac(x-2,n):=Product_{j=0..n-1} (x-2+j), n >= 1, risefac(x-2,0)=1. Compare this with the formula risefac(x,n)=|S1|(n,x), with the row polynomials |S1|(n,x) of A132393 (unsigned Stirling1).
This is the third triangle of an a-family of Sheffer arrays, call them |S1|(a), with e.g.f. of the row polynomials |S1|(a;x;z) = ((1-z)^a)*exp(-x*log(1-z)). In the notation showing the column e.g.f.s this is Sheffer ((1-z)^a,-log(1-z)). In the umbral notation (see the Roman reference, given under A094645) this is called Sheffer for (exp(a*t),1-exp(-t)). For a=0 this becomes the unsigned Stirling1 triangle |S1|(0) = A132393 with row polynomials |S1|(0;n,x) =: s1(n,x).
E.g.f. column number k (with leading zeros): ((1-x)^a)*((-log(1-x))^k)/k!, k >= 0.
E.g.f. for row sums is (1-x)^(a-1), and the e.g.f. for the alternating row sums is (1-x)^(a+1).
Row polynomial recurrence:
  |S1|(a;n,x)=(x+(n-1-a))*|S1|(a;n-1,m), |S1|(a;0,x)=1.
Meixner identity (see the reference under A060338):
  |S1|(a;n,x) - |S1|(a;n,x-1) = n*|S1|(a;n-1,x), n >= 1,
Also (from the corollary 3.7.2 on p. 50 of the Roman reference): |S1|(a;n,x) = (x-a)*|S1|(a;n-1,x+1), n >= 1.
Recurrence: |S1|(a;n,k) = |S1|(a;n-1,k-1) + (n-(a+1))*|S1|(a;n-1,k); |S1|(a;n,k)=0 if n < m, |S1|(a;n,-1)=0, |S1|(a;0,0)=1.
Connection to |Stirling1|=|S1|(0):
  |S1|(a;n,k) = Sum_{p=0..a} |S1|(a;a,p)*abs(Stirling1(n-a,k-p)), n >= a.
The exponential convolution identity is
  |S1|(a;n,x+y) = Sum_{k=0..n} binomial(n,k)*|S1|(a;k,y)*s1(n-k,x), n >= 0, with symmetry x <-> y.
The Sheffer a- and z-sequences are (see the W. Lang link under A006232): Sha(a;n)=A164555(n)/A027642(n) (independent of a) with e.g.f. x/(1-exp(-x)), and the z-sequence has e.g.f. (exp(a*x)-1)/(exp(-x)-1).
The inverse Sheffer matrix has e.g.f. exp(a*z)*exp(x*(1-exp(-z))), in short notation (exp(a*z),1-exp(-z)),
  (or in umbral notation ((1-t)^a,-log(1-t))).
(End)

			Triangle begins
   1;
  -2,  1;
   2, -3,  1;
   0,  2, -3,  1;
   0,  2, -1, -2,  1;
   0,  4,  0, -5,  0,  1;
   ...
risefac(x-2,3) = (x-2)*(x-1)*x = 2*x-3*x^2+x^3.
-1 = T(4,2) = T(3,1) + 1*T(3,2) =  2 + (-3).
T(4,3) = 2*abs(S1(2,3)) - 3*abs(S1(2,2)) + 1*abs(S1(2,1)) = 2*0 - 3*1 + 1*1 = -2.

Cf. A049444, A049458, A094645.

A094646_row := n -> seq((-1)^(n-k)*coeff(expand(pochhammer(x-n+3, n)), x, k), k=0..n): seq(print(A094646_row(n)), n = 0..6); # Peter Luschny, May 16 2013

Flatten[ Table[ CoefficientList[ Pochhammer[x-2, n], x], {n, 0, 10}]] (* Jean-François Alcover, Sep 26 2011 *)

E.g.f.: (1-y)^(2-x).
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 = 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13 respectively. - Philippe Deléham, Nov 13 2007
If we define f(n,i,a) = Sum_{k=0..n-i} binomial(n,k)*Stirling1(n-k,i)*Product_{j=0..k-1} (-a-j), then |T(n,i)| = |f(n,i,-2)|, for n=1,2,...; i=0..n. - Milan Janjic, Dec 21 2008
From Wolfdieter Lang, Jun 23 2011: (Start)
risefac(x-2,n) = Sum_{k=0..n} T(n,k)*x^k, n >= 0, with the rising factorials (see a comment above).
Recurrence: T(n,k) = T(n-1,k-1) + (n-3)*T(n-1,k); T(n,k)=0 if n < m, T(n,-1)=0, T(0,0)=1.
T(n,k) = 2*abs(S1(n-2,k)) - 3*abs(S1(n-2,k-1)) + abs(S1(n-2,k-2)), n >= 2, with S1(n,k) = Stirling1(n,k) = A048994(n,k).
E.g.f. column number k (with leading zeros):
((1-x)^2)*((-log(1-x))^k)/k!, k >= 0.
E.g.f. for row sums is 1-x, i.e., [1,-1,0,0,...],
  and the e.g.f. for the alternating row sums is (1-x)^3. i.e., [1,-3,3,1,0,0,...]. (End)

A105793 Expansion of e.g.f. (1 + y)^(1 + x).

Original entry on oeis.org

1, 1, 1, 0, 1, 1, 0, -1, 0, 1, 0, 2, -1, -2, 1, 0, -6, 5, 5, -5, 1, 0, 24, -26, -15, 25, -9, 1, 0, -120, 154, 49, -140, 70, -14, 1, 0, 720, -1044, -140, 889, -560, 154, -20, 1, 0, -5040, 8028, -64, -6363, 4809, -1638, 294, -27, 1, 0, 40320, -69264, 8540, 50840, -44835, 17913, -3990, 510, -35, 1

Offset: 0

Paul Barry, Apr 20 2005

Generalized Stirling number triangle of first kind. Row sums are (1,2,2,0,0,0,...) = 2C(2,n) - 2C(1,n) + C(0,n). Inverse is A105794.

Triangle T(n,k), 0 <= k <= n, read by rows, given by [1, -1, 0, -2, -1, -3, -2, -4, -3, -5, -4, -6, ...] DELTA [1, 0, 1, 0, 1, 0, 1, 0, 1, 0, 1, ...] where DELTA is the operator defined in A084938. - Philippe Deléham, Aug 23 2006

			From _Wolfdieter Lang_, Jun 19 2017: (Start)
The triangle T(n, k) starts
  n\k  0     1      2    3     4      5     6     7   8   9 10 ...
  0:   1
  1:   1     1
  2:   0     1      1
  3:   0    -1      0    1
  4:   0     2     -1   -2     1
  5:   0    -6      5    5    -5      1
  6:   0    24    -26  -15    25     -9     1
  7:   0  -120    154   49  -140     70   -14     1
  8:   0   720  -1044 -140   889   -560   154   -20   1
  9:   0 -5040   8028  -64 -6363   4809 -1638   294 -27   1
  10:  0 40320 -69264 8540 50840 -44835 17913 -3990 510 -35  1
  ... reformatted
Recurrence from a- and z-sequence (see above): T(1, 0) = T(0, 0) = 1; T(1, 1) = (1/1)*(1*T(0, 0)) = 1, T(2, 0) = 2*(T(1, 0) - T(1, 1)) = 0, T(2, 1) = (2/1)*(T(1,0) + (-1/2)*T(1, 1)) = 1. T(3, 1) = (3/1)*(0 + (-1/2)*T(2, 1) + (1/6)*T(2, 2)) = -1. (End)

Cf. A007318, A027641/A027642, A033999, A048994, A094645, A269953.

t[0, 0] = 1; t[n_, k_] := StirlingS1[n, k] + n*StirlingS1[n-1, k]; Table[t[n, k], {n, 0, 10}, {k, 0, n}] // Flatten (* Jean-François Alcover, Dec 04 2013, after Giuliano Cabrele *)

E.g.f.: (1+y)^(1+x); rows have g.f. k!*binomial(x+1, k); Columns have g.f. (1+x)*log(1+x)^k.
If we define f(n,i,a) = Sum_{k=0..n-i} binomial(n,k)*Stirling1(n-k,i)*Product_{j=0..k-1} (-a-j), then T(n,i) = f(n,i,-1), for n=1,2,...; i=0...n. - Milan Janjic, Dec 21 2008
So T(n,k) = Stirling1(n,k) + n*Stirling1(n-1,k), Stirling1 being the (signed) Stirling numbers of first kind A048994. In terms of lower triangular matrices, 0<= k <= n, T is also the product [Stirling1] * [Pascal] = [A048994] * [A007318], i.e., T(n,k) = Sum_{j=0..n} Stirling1(n,j) * binomial(j,k). - Giuliano Cabrele, Jan 19 2009
This is the triangle of connection constants for expressing the basis of falling factorial polynomials x_(k) := x*(x-1)*...*(x-k+1) in terms of the polynomial sequence (x-1)^n, that is, x_(n) = Sum_{k = 0..n} T(n,k)*(x-1)^k. - Peter Bala, Jul 10 2013
From Wolfdieter Lang, Jun 19 2017: (Start)
Triangle T is the (infinite) matrix product of A048994 (Stirling1) and A007318 (Pascal): T(n,k) = Sum_{m=k..n} Stirling1(n, m)*Pascal(m, k), n >= k >= 0, and 0 for n < k. Note that the Pascal matrix is Sheffer (e^t, t) of the Appell type.
T is the Sheffer (aka exponential Riordan) matrix (1+t, log(1+t)).
E.g.f. column k: (1+x)*(log(1+x))^k/k!, k >= 0.
The a-sequence for T is A027641/A027642 (Bernoulli), and the z-sequence is A033999 (repeat(1,-1)) (see a W. Lang link under A006232 for a- and z-sequences for Sheffer matrices, also for references).
Therefore the combined recurrence is: T(n, 0) = n*Sum_{j=0..n-1} (-1)^j*T(n-1, j), n >= 1, T(0, 0) = 1, and T(n, m) = (n/m)*Sum_{j=0..n-m} binomial(m-1+j, m-1)*Bernoulli(j)*T(n-1, m-1+j), n >= m >= 1. (End)

A189924 a(n) = abs(Stirling1(n+2,2)) - abs(Stirling1(n+2,3)).

Original entry on oeis.org

1, 2, 5, 15, 49, 140, -64, -8540, -146124, -2124936, -30374136, -445116672, -6793958016, -108691150464, -1826654613120, -32257962443520, -598196854045440, -11635261535301120, -237044583523514880, -5050811716879104000

Offset: 0

Wolfdieter Lang, Jun 21 2011

This is the fourth (k=3) column sequence in triangle A094645 without leading zeros.

G. C. Greubel, Table of n, a(n) for n = 0..445

Cf. A081047 (column k=2).

[Abs(StirlingFirst(n+2,2)) - Abs(StirlingFirst(n+2,3)): n in [0..30]]; // G. C. Greubel, Jan 13 2018

Table[Abs[StirlingS1[n+2,2]]-Abs[StirlingS1[n+2,3]],{n,0,20}] (* Harvey P. Dale, May 21 2015 *)

a(n)=abs(stirling(n+2,2))-abs(stirling(n+2,3)) \\ Charles R Greathouse IV, Jun 27 2011

a(n) = abs(Stirling1(n+2,2)) - abs(Stirling1(n+2,3)), with the unsigned Stirling1 numbers abs(Stirling1(n,k)) = A132393(n,k).

E.g.f.: (1/2)*(2-log(1-x)^2)/(1-x)^2 (from differentiating three times (1-x)*((-log(1-x))^3)/3!).

Showing 1-10 of 11 results. Next

cp's OEIS Frontend

A164555 Numerators of the "original" Bernoulli numbers; also the numerators of the Bernoulli polynomials at x=1.

Views

Author

Keywords

Comments

Examples

References

Links

Crossrefs

Programs

Haskell

Maple

Mathematica

Sage

Formula

Extensions

A132393 Triangle of unsigned Stirling numbers of the first kind (see A048994), read by rows, T(n,k) for 0 <= k <= n.

Views

Author

Keywords

Comments

Examples

References

Links

Crossrefs

Programs

Haskell

Maple

Mathematica

Maxima

PARI

Formula

A049444 Generalized Stirling number triangle of first kind.

Views

Author

Keywords

Comments

Examples

References

Links

Crossrefs

Programs

Haskell

Maple

Mathematica

Formula

Extensions

A105794 Inverse of a generalized Stirling number triangle of first kind.

Views

Author

Keywords

Comments

Examples

References

Links

Crossrefs

Programs

Maple

Mathematica

Formula

A049458 Generalized Stirling number triangle of first kind.

Views

Author

Keywords

Comments

Examples

References

Links

Crossrefs

Programs

Haskell

Maple

Mathematica

Formula

Extensions

A105954 Array read by descending antidiagonals: A(n, k) = (n + 1)! * H(k, n + 1), where H(n, k) is a higher-order harmonic number, H(0, k) = 1/k and H(n, k) = Sum_{j=1..k} H(n-1, j), for 0 <= k <= n.

Views

Author

Keywords