A130850 -id:A130850 - cp's OEIS frontend

A123125 Triangle of Eulerian numbers T(n,k), 0 <= k <= n, read by rows.

1, 0, 1, 0, 1, 1, 0, 1, 4, 1, 0, 1, 11, 11, 1, 0, 1, 26, 66, 26, 1, 0, 1, 57, 302, 302, 57, 1, 0, 1, 120, 1191, 2416, 1191, 120, 1, 0, 1, 247, 4293, 15619, 15619, 4293, 247, 1, 0, 1, 502, 14608, 88234, 156190, 88234, 14608, 502, 1, 0, 1, 1013, 47840, 455192, 1310354, 1310354, 455192, 47840, 1013, 1

Offset: 0

Philippe Deléham, Sep 30 2006

The beginning of this sequence does not quite agree with the usual version, which is A173018. - N. J. A. Sloane, Nov 21 2010
Each row of A123125 is the reverse of the corresponding row in A173018. - Michael Somos, Mar 17 2011
A008292 (subtriangle for k>=1 and n>=1) is the main entry for these numbers.
Triangle T(n,k), 0 <= k <= n, read by rows given by [0,1,0,2,0,3,0,4,0,5,0,...] DELTA [1,0,2,0,3,0,4,0,5,0,6,...] where DELTA is the operator defined in A084938.
Row sums are the factorials. - Roger L. Bagula and Gary W. Adamson, Aug 14 2008
If the initial zero column is deleted, the result is A008292. - Roger L. Bagula and Gary W. Adamson, Aug 14 2008
This result gives an alternative method of calculating the Eulerian numbers by an Umbral Calculus expansion from Comtet. - Roger L. Bagula, Nov 21 2009
This function seems to be equivalent to the PolyLog expansion. - Roger L. Bagula, Nov 21 2009
A raising operator formed from the e.g.f. of this entry is the generator of a sequence of polynomials p(n,x;t) defined in A046802 that specialize to those for A119879 as p(n,x;-1), A007318 as p(n,x;0), A073107 as p(n,x;1), and A046802 as p(n,0;t). See Copeland link for more associations. - Tom Copeland, Oct 20 2015
The Eulerian numbers in this setup count the permutation trees of power n and width k (see the Luschny link). For the associated combinatorial statistic over permutations see the Sage program below and the example section. - Peter Luschny, Dec 09 2015 [See Elder et al. link. Peter Luschny, Jul 13 2022]
From Wolfdieter Lang, Apr 03 2017: (Start)
The row polynomials R(n, x) = Sum_{k=0..n} T(n, k)*x^k are the numerator polynomials of the o.g.f. G(n, x) of n-powers {m^n}_{m>=0} (with 0^0 = 1): G(n, x) = R(n, x)/(1-x)^(n+1). See the Aug 14 2008 formula, where f(x,n) = R(n, x). The e.g.f. of R(n, t) is given in Copeland's Oct 14 2015 formula below.
The first nine column sequences are A000007, A000012, A000295, A000460, A000498, A000505, A000514, A001243, A001244. (End)
With all offsets 0, let A_n(x;y) = (y + E.(x))^n, an Appell sequence in y where E.(x)^k = E_k(x) are the Eulerian polynomials of this entry, A123125. Then the row polynomials of A046802 (the h-polynomials of the stellahedra) are given by h_n(x) = A_n(x;1); the row polynomials of A248727 (the face polynomials of the stellahedra), by f_n(x) = A_n(1 + x;1); the Swiss-knife polynomials of A119879, by Sw_n(x) = A_n(-1;1 + x); and the row polynomials of the Worpitsky triangle (A130850), by w_n(x) = A(1 + x;0). Other specializations of A_n(x;y) give A090582 (the f-polynomials of the permutohedra, cf. also A019538) and A028246 (another version of the Worpitsky triangle). - Tom Copeland, Jan 24 2020
Let b(n) = (1/(n+1))*Sum_{k=0..n-1} (-1)^(n-k+1)*T(n, k+1) / binomial(n, k+1). Then b(n) = Bernoulli(n, 1) = -n*Zeta(1 - n) = Integral_{x=0..1} F_n(x) for n >= 1. Here F_n(x) are the signed Fubini polynomials (A278075). (See also Rzadkowski and Urlinska, example 1.) - Peter Luschny, Feb 15 2021
Patrick J. Burchell (see link) describes the following method: To get the k-th row of the triangle write the nonnegative integers with a fixed exponent k as a sequence, 0^k, 1^k, 2^k, ..., and then apply the first differences to them k + 1 times. - Peter Luschny, Apr 02 2023

			The triangle T(n, k) begins:
  n\k 0 1    2     3      4       5       6      7     8    9 10...
  0:  1
  1:  0 1
  2:  0 1    1
  3:  0 1    4     1
  4:  0 1   11    11      1
  5:  0 1   26    66     26       1
  6:  0 1   57   302    302      57       1
  7:  0 1  120  1191   2416    1191     120      1
  8:  0 1  247  4293  15619   15619    4293    247     1
  9:  0 1  502 14608  88234  156190   88234  14608   502    1
 10:  0 1 1013 47840 455192 1310354 1310354 455192 47840 1013  1
...  Reformatted. - _Wolfdieter Lang_, Feb 14 2015
------------------------------------------------------------------
The width statistic over permutations, n=4.
  [1, 2, 3, 4] => 3; [1, 2, 4, 3] => 2; [1, 3, 2, 4] => 2; [1, 3, 4, 2] => 2;
  [1, 4, 2, 3] => 2; [1, 4, 3, 2] => 1; [2, 1, 3, 4] => 3; [2, 1, 4, 3] => 2;
  [2, 3, 1, 4] => 2; [2, 3, 4, 1] => 3; [2, 4, 1, 3] => 2; [2, 4, 3, 1] => 2;
  [3, 1, 2, 4] => 3; [3, 1, 4, 2] => 3; [3, 2, 1, 4] => 2; [3, 2, 4, 1] => 3;
  [3, 4, 1, 2] => 3; [3, 4, 2, 1] => 2; [4, 1, 2, 3] => 4; [4, 1, 3, 2] => 3;
  [4, 2, 1, 3] => 3; [4, 2, 3, 1] => 3; [4, 3, 1, 2] => 3; [4, 3, 2, 1] => 2;
Gives row(4) = [0, 1, 11, 11, 1]. - _Peter Luschny_, Dec 09 2015
------------------------------------------------------------------
From _Wolfdieter Lang_, Apr 03 2017: (Start)
Recurrence: T(5, 3) = (6-3)*T(4, 2) + 3*T(4, 3) = 3*11 + 3*11 = 66.
O.g.f. column k=2: (x/(1 - 2*x))*E_x*(x/(1-x)) = (x/(1-x))^2/(1-2*x).
E.g.f. column k=2: A(2, x) = x*A(1, x) + x*E(1, x) = x*1 + x*(exp(x)-1) = x*exp(x), hence E(2, x) = (1 + int(x*exp(-x),x ))*exp(2*x) = exp(x)*(exp(x) - (1+x)). See A000295. (End)

L. Comtet, Advanced Combinatorics, Reidel, Holland, 1978, page 245. [Roger L. Bagula, Nov 21 2009]
Ronald L. Graham, Donald E. Knuth and Oren Patashnik, Concrete Mathematics, 2nd ed.; Addison-Wesley, 1994, p. 268, Row reversed table 268. - Wolfdieter Lang, Apr 03 2017
Douglas C. Montgomery and Lynwood A. Johnson, Forecasting and Time Series Analysis, MaGraw-Hill, New York, 1976, page 91. - Roger L. Bagula and Gary W. Adamson, Aug 14 2008

Reinhard Zumkeller, Rows n = 0..125 of triangle, flattened
Paul Barry, Eulerian polynomials as moments, via exponential Riordan arrays, arXiv preprint arXiv:1105.3043 [math.CO], 2011, J. Int. Seq. 14 (2011) # 11.9.5
Paul Barry, Generalized Stirling Numbers, Exponential Riordan Arrays, and Toda Chain Equations, Journal of Integer Sequences, 17 (2014), #14.2.3.
Paul Barry, On a transformation of Riordan moment sequences, arXiv:1802.03443 [math.CO], 2018.
Paul Barry, Generalized Eulerian Triangles and Some Special Production Matrices, arXiv:1803.10297 [math.CO], 2018.
V. Batyrev and M. Blume, The functor of toric varieties associated with Weyl chambers and Losev-Manin moduli spaces, p. 11, arXiv:/0911.3607 [math.AG], 2009. [_Tom Copeland_, Oct 16 2015]
Anna Borowiec and Wojciech Mlotkowski, New Eulerian numbers of type D, arXiv:1509.03758 [math.CO], 2015.
Patrick J. Burchell, A Generalisation of Ramanujan's (back of the envelope) Method for Divergent Series, arXiv:2303.14045 math.NT, 2023.
A. Cohen, Eulerian polynomials of spherical type, Münster J. of Math. 1 (2008). [_Tom Copeland_, Oct 16 2015]
Tom Copeland, The Elliptic Lie Triad: Ricatti and KdV Equations, Infinigens, and Elliptic Genera
Jennifer Elder, Nadia Lafrenière, Erin McNicholas, Jessica Striker and Amanda Welch, Homomesies on permutations -- an analysis of maps and statistics in the FindStat database, arXiv:2206.13409 [math.CO], 2022-2023. (Def. 4.20 and Prop. 4.22.)
FindStat - Combinatorial Statistic Finder, The number of descents of a permutation.
F. Hirzebruch, Eulerian polynomials, Münster J. of Math. 1 (2008), pp. 9-12.
P. Hitczenko and S. Janson, Weighted random staircase tableaux, arXiv preprint arXiv:1212.5498 [math.CO], 2012.
Hsien-Kuei Hwang, Hua-Huai Chern, and Guan-Huei Duh, An asymptotic distribution theory for Eulerian recurrences with applications, arXiv:1807.01412 [math.CO], 2018.
Svante Janson, Euler-Frobenius numbers and rounding, arXiv preprint arXiv:1305.3512 [math.PR], 2013.
Katarzyna Kril and Wojciech Mlotkowski, Permutations of Type B with Fixed Number of Descents and Minus Signs, Volume 26(1) of The Electronic Journal of Combinatorics, 2019.
Wolfdieter Lang, On Sums of Powers of Arithmetic Progressions, and Generalized Stirling, Eulerian and Bernoulli numbers, arXiv:1707.04451 [math.NT], 2017.
Huyile Liang, Yanni Pei, and Yi Wang, Analytic combinatorics of coordination numbers of cubic lattices, arXiv:2302.11856 [math.CO], 2023. See p. 22.
A. Losev and Y. Manin, New moduli spaces of pointed curves and pencils of flat connections, arXiv preprint arXiv:math/0001003 [math.AG], 2000 (p. 8). - _Tom Copeland_, Oct 16 2015
Peter Luschny, Permutation Trees
G. Rzadkowski and M. Urlinska, A Generalization of the Eulerian Numbers, arXiv:1612.06635 [math.CO], 2016.

See A008292 (subtriangle for k>=1 and n>=1), which is the main entry for these numbers. Another version has the zeros at the ends of the rows, as in Concrete Mathematics: see A173018.
Cf. A007318, A130850, A028246, A046802, A009006, A019538, A090582, A119879, A248727.
T(2n,n) gives A180056.

a123125 n k = a123125_tabl !! n !! k
a123125_row n = a123125_tabl !! n
a123125_tabl = [1] : zipWith (:) [0, 0 ..] a008292_tabl
-- Reinhard Zumkeller, Nov 06 2013

gf := 1/(1 - t*exp(x)): ser := series(gf, x, 12):
cx := n -> (-1)^(n + 1)*factor(n!*coeff(ser, x, n)*(t - 1)^(n + 1)):
seq(print(seq(coeff(cx(n), t, k), k = 0..n)), n = 0..9); # Peter Luschny, Feb 11 2021
A123125 := proc(n, k) option remember; if k = n then 1 elif k <= 0 or k > n then 0 else k*procname(n-1, k) + (n-k+1)*procname(n-1, k-1) fi end:
seq(print(seq(A123125(n, k), k=0..n)), n=0..10); # Peter Luschny, Mar 28 2021
# Alternative (Patrick J. Burchell):
t := a -> Statistics:-Difference([0, a]): Trow := k -> (t@@(k+1))([seq(n^k, n = 0..k)]):
seq(print(Trow(n)), n = 0..6); # Peter Luschny, Apr 02 2023

f[x_, n_] := f[x, n] = (1 - x)^(n + 1)*Sum[k^n*x^k, {k, 0, Infinity}];
Table[CoefficientList[f[x, n], x], {n,0,9}] // Flatten (* Roger L. Bagula, Aug 14 2008 *)
t[n_ /; n >= 0, 0] = 1; t[n_, k_] /; k<0 || k>n = 0; t[n_, k_] := t[n, k] = (n-k) t[n-1, k-1] + (k+1) t[n-1, k]; T[n_, k_] := t[n, n-k];
Table[T[n, k], {n,0,10}, {k, 0, n}] // Flatten (* Jean-François Alcover, May 26 2019 *)
A123125[n_, k_] := Sum[(-1)^j*(n - j - k + 1)^n * Binomial[n + 1, j], {j, 0, n - k}];
Table[A123125[n, k], {n, 0, 9}, {k, 0, n}] // TableForm  (* Peter Luschny, Aug 12 2022 *)

from math import isqrt, comb
def A123125(n):
    a = (m:=isqrt(k:=n+1<<1))+(k>m*(m+1))
    b = comb(a+1,2)-n
    return sum(-(b-j)**(a-1)*comb(a,j) if j&1 else (b-j)**(a-1)*comb(a,j) for j in range(b)) # Chai Wah Wu, Nov 13 2024

def statistic_eulerian(pi):
    if not pi: return 0
    h, i, branch, next = 0, len(pi), [0], pi[0]
    while True:
        while next < branch[len(branch)-1]:
            del(branch[len(branch)-1])
        current = 0
        h += 1
        while next > current:
            i -= 1
            if i == 0: return h
            branch.append(next)
            current, next = next, pi[i]
def A123125_row(n):
    L = [0]*(n+1)
    for p in Permutations(n):
        L[statistic_eulerian(p)] += 1
    return L
[A123125_row(n) for n in range(7)] # Peter Luschny, Dec 09 2015

Sum_{k=0..n} T(n,k) = n! = A000142(n).
Sum_{k=0..n} 2^k*T(n,k) = A000629(n).
Sum_{k=0..n} 3^k*T(n,k) = abs(A009362(n+1)).
Sum_{k=0..n} 2^(n-k)*T(n,k) = A000670(n).
Sum_{k=0..n} T(n,k)*3^(n-k) = A122704(n). - Philippe Deléham, Nov 07 2007
G.f.: f(x,n) = (1 - x)^(n + 1)*Sum_{k>=0} k^n*x^k. - Roger L. Bagula and Gary W. Adamson, Aug 14 2008. f is not the g.f. of the triangle, it is the polynomial of row n. See an Apr 03 2017 comment above - Wolfdieter Lang, Apr 03 2017
Sum_{k=0..n} T(n,k)*x^k = A000007(n), A000142(n), A000629(n), A123227(n), A201355(n), A201368(n) for x = 0, 1, 2, 3, 4, 5 respectively. - Philippe Deléham, Dec 01 2011
E.g.f. (1-t)/(1-t*exp((1-t)x)). A123125 * A007318 = A130850 = unsigned A075263, related to reversed A028246. A007318 * A123125 = A046802. Evaluating the row polynomials at -1, giving the alternating-sign row sum, generates A009006. - Tom Copeland, Oct 14 2015
From Wolfdieter Lang, Apr 03 2017: (Start)
T(n, k) = A173018(n, n-k), 0 <= k <= n. Row reversed Euler's triangle. See Graham et al., p. 268.
Recurrence (from A173018): T(n, 0) = 1 if n=0 else 0; T(n, k) = 0 if n < k and T(n, k) = (n+1-k)*T(n-1, k-1) + k*T(n-1, k) else.
T(n, k) = Sum_{j=0..k} (-1)^(k-j)*binomial(n-j, k-j)*S2(n, j)*j!, 0 <= k <= n, else 0. For S2(n, k)*k! see A131689.
The recurrence for the o.g.f. of the sequence of column k is
  G(k, x) = (x/(1 - k*x))*(E_x - (k-2))*G(k-1, x), with the Euler operator E_x  = x*d_x, for k >= 1, with G(0, x) = 1. (Proof from the recurrence of T(n, k)).
The e.g.f of the sequence of column k is found from E(k, x) = (1 + int(A(k, x),x)*exp(-k*x))*exp(k*x), k >= 1, with the recurrence
  A(k, x) = x*A(k-1, x) +(1 + (1-k)*(1-x))*E(k-1, x) for k >= 1, with  A(0,x)= 0. (Proof from the recurrence of T(n, k)). (End)
T(n, k) = Sum_{j=0..n-k} (-1)^j*(n-j-k+1)^n*binomial(n + 1, j). - Peter Luschny, Aug 12 2022
G.f.: Sum_{m >= 0} x^m/(1/(1-x)-m*t). - Mamuka Jibladze, Mar 12 2025

A028246 Triangular array a(n,k) = (1/k)Sum_{i=0..k} (-1)^(k-i)binomial(k,i)*i^n; n >= 1, 1 <= k <= n, read by rows.

Original entry on oeis.org

1, 1, 1, 1, 3, 2, 1, 7, 12, 6, 1, 15, 50, 60, 24, 1, 31, 180, 390, 360, 120, 1, 63, 602, 2100, 3360, 2520, 720, 1, 127, 1932, 10206, 25200, 31920, 20160, 5040, 1, 255, 6050, 46620, 166824, 317520, 332640, 181440, 40320, 1, 511, 18660, 204630, 1020600, 2739240, 4233600, 3780000, 1814400, 362880

Offset: 1

N. J. A. Sloane, Doug McKenzie (mckfam4(AT)aol.com)

Let M = n X n matrix with (i,j)-th entry a(n+1-j, n+1-i), e.g., if n = 3, M = [1 1 1; 3 1 0; 2 0 0]. Given a sequence s = [s(0)..s(n-1)], let b = [b(0)..b(n-1)] be its inverse binomial transform and let c = [c(0)..c(n-1)] = M^(-1)*transpose(b). Then s(k) = Sum_{i=0..n-1} b(i)*binomial(k,i) = Sum_{i=0..n-1} c(i)*k^i, k=0..n-1. - Gary W. Adamson, Nov 11 2001
From Gary W. Adamson, Aug 09 2008: (Start)
Julius Worpitzky's 1883 algorithm generates Bernoulli numbers.
By way of example [Wikipedia]:
  B0 =  1;
  B1 = 1/1 -  1/2;
  B2 = 1/1 -  3/2 +   2/3;
  B3 = 1/1 -  7/2 +  12/3 -    6/4;
  B4 = 1/1 - 15/2 +  50/3 -   60/4 +   24/5;
  B5 = 1/1 - 31/2 + 180/3 -  390/4 +  360/5 -  120/6;
  B6 = 1/1 - 63/2 + 602/3 - 2100/4 + 3360/5 - 2520/6 + 720/7;
  ...
Note that in this algorithm, odd n's for the Bernoulli numbers sum to 0, not 1, and the sum for B1 = 1/2 = (1/1 - 1/2). B3 = 0 = (1 - 7/2 + 13/3 - 6/4) = 0. The summation for B4 = -1/30. (End)
Pursuant to Worpitzky's algorithm and given M = A028246 as an infinite lower triangular matrix, M * [1/1, -1/2, 1/3, ...] (i.e., the Harmonic series with alternate signs) = the Bernoulli numbers starting [1/1, 1/2, 1/6, ...]. - Gary W. Adamson, Mar 22 2012
From Tom Copeland, Oct 23 2008: (Start)
G(x,t) = 1/(1 + (1-exp(x*t))/t) = 1 + 1 x + (2 + t)*x^2/2! + (6 + 6t + t^2)*x^3/3! + ... gives row polynomials for A090582, the f-polynomials for the permutohedra (see A019538).
G(x,t-1) = 1 + 1*x + (1 + t)*x^2 / 2! + (1 + 4t + t^2)*x^3 / 3! + ... gives row polynomials for A008292, the h-polynomials for permutohedra.
G[(t+1)x,-1/(t+1)] = 1 + (1+ t) x + (1 + 3t + 2 t^2) x^2 / 2! + ... gives row polynomials for the present triangle. (End)
The Worpitzky triangle seems to be an apt name for this triangle. - Johannes W. Meijer, Jun 18 2009
If Pascal's triangle is written as a lower triangular matrix and multiplied by A028246 written as an upper triangular matrix, the product is a matrix where the (i,j)-th term is (i+1)^j. For example,
  1,0,0,0     1,1,1, 1      1,1, 1, 1
  1,1,0,0  *  0,1,3, 7   =  1,2, 4, 8
  1,2,1,0     0,0,2,12      1,3, 9,27
  1,3,3,1     0,0,0, 6      1,4,16,64
  So, numbering all three matrices' rows and columns starting at 0, the (i,j) term of the product is (i+1)^j. - Jack A. Cohen (ProfCohen(AT)comcast.net), Aug 03 2010
The Fi1 and Fi2 triangle sums are both given by sequence A000670. For the definition of these triangle sums see A180662. The mirror image of the Worpitzky triangle is A130850. - Johannes W. Meijer, Apr 20 2011
Let S_n(m) = 1^m + 2^m + ... + n^m. Then, for n >= 0, we have the following representation of S_n(m) as a linear combination of the binomial coefficients:
  S_n(m) = Sum_{i=1..n+1} a(i+n*(n+1)/2)*C(m,i). E.g., S_2(m) = a(4)*C(m,1) + a(5)*C(m,2) + a(6)*C(m,3) = C(m,1) + 3*C(m,2) + 2*C(m,3). - Vladimir Shevelev, Dec 21 2011
Given the set X = [1..n] and 1 <= k <= n, then a(n,k) is the number of sets T of size k of subset S of X such that S is either empty or else contains 1 and another element of X and such that any two elemements of T are either comparable or disjoint. - Michael Somos, Apr 20 2013
Working with the row and column indexing starting at -1, a(n,k) gives the number of k-dimensional faces in the first barycentric subdivision of the standard n-dimensional simplex (apply Brenti and Welker, Lemma 2.1). For example, the barycentric subdivision of the 2-simplex (a triangle) has 1 empty face, 7 vertices, 12 edges and 6 triangular faces giving row 4 of this triangle as (1,7,12,6). Cf. A053440. - Peter Bala, Jul 14 2014
See A074909 and above g.f.s for associations among this array and the Bernoulli polynomials and their umbral compositional inverses. - Tom Copeland, Nov 14 2014
An e.g.f. G(x,t) = exp[P(.,t)x] = 1/t - 1/[t+(1-t)(1-e^(-xt^2))] = (1-t) * x + (-2t + 3t^2 - t^3) * x^2/2! + (6t^2 - 12t^3 + 7t^4 - t^5) * x^3/3! + ... for the shifted, reverse, signed polynomials with the first element nulled, is generated by the infinitesimal generator g(u,t)d/du = [(1-u*t)(1-(1+u)t)]d/du, i.e., exp[x * g(u,t)d/du] u eval. at u=0 generates the polynomials. See A019538 and the G. Rzadkowski link below for connections to the Bernoulli and Eulerian numbers, a Ricatti differential equation, and a soliton solution to the KdV equation. The inverse in x of this e.g.f. is Ginv(x,t) = (-1/t^2)*log{[1-t(1+x)]/[(1-t)(1-tx)]} = [1/(1-t)]x + [(2t-t^2)/(1-t)^2]x^2/2 + [(3t^2-3t^3+t^4)/(1-t)^3]x^3/3  + [(4t^3-6t^4+4t^5-t^6)/(1-t)^4]x^4/4 + ... . The numerators are signed, shifted A135278 (reversed A074909), and the rational functions are the columns of A074909. Also, dG(x,t)/dx = g(G(x,t),t) (cf. A145271). (Analytic G(x,t) added, and Ginv corrected and expanded on Dec 28 2015.) - Tom Copeland, Nov 21 2014
The operator R = x + (1 + t) + t e^{-D} / [1 + t(1-e^(-D))] = x + (1+t) + t - (t+t^2) D + (t+3t^2+2t^3) D^2/2! - ... contains an e.g.f. of the reverse row polynomials of the present triangle, i.e., A123125 * A007318 (with row and column offset 1 and 1). Umbrally, R^n 1 = q_n(x;t) = (q.(0;t)+x)^n, with q_m(0;t) = (t+1)^(m+1) - t^(m+1), the row polynomials of A074909, and D = d/dx. In other words, R generates the Appell polynomials associated with the base sequence A074909. For example, R 1 = q_1(x;t) = (q.(0;t)+x) = q_1(0;t) + q__0(0;t)x = (1+2t) + x, and R^2 1 = q_2(x;t) = (q.(0;t)+x)^2 = q_2(0:t) + 2q_1(0;t)x + q_0(0;t)x^2 = 1+3t+3t^2 + 2(1+2t)x + x^2. Evaluating the polynomials at x=0 regenerates the base sequence. With a simple sign change in R, R generates the Appell polynomials associated with A248727. - Tom Copeland, Jan 23 2015
For a natural refinement of this array, see A263634. - Tom Copeland, Nov 06 2015
From Wolfdieter Lang, Mar 13 2017: (Start)
The e.g.f. E(n, x) for {S(n, m)}{m>=0} with S(n, m) = Sum{k=1..m} k^n, n >= 0, (with undefined sum put to 0) is exp(x)*R(n+1, x) with the exponential row polynomials R(n, x) = Sum_{k=1..n} a(n, k)*x^k/k!. E.g., e.g.f. for n = 2, A000330: exp(x)*(1*x/1!+3*x^2/2!+2*x^3/3!).
The o.g.f. G(n, x) for {S(n, m)}{m >=0} is then found by Laplace transform to be G(n, 1/p) = p*Sum{k=1..n} a(n+1, k)/(p-1)^(2+k).
  Hence G(n, x) = x/(1 - x)^(n+2)*Sum_{k=1..n} A008292(n,k)*x^(k-1).
  E.g., n=2: G(2, 1/p) = p*(1/(p-1)^2 + 3/(p-1)^3 + 2/(p-1)^4) = p^2*(1+p)/(p-1)^4; hence G(2, x) = x*(1+x)/(1-x)^4.
This works also backwards: from the o.g.f. to the e.g.f. of {S(n, m)}_{m>=0}. (End)
a(n,k) is the number of k-tuples of pairwise disjoint and nonempty subsets of a set of size n. - Dorian Guyot, May 21 2019
From Rajesh Kumar Mohapatra, Mar 16 2020: (Start)
a(n-1,k) is the number of chains of length k in a partially ordered set formed from subsets of an n-element set ordered by inclusion such that the first term of the chains is either the empty set or an n-element set.
Also, a(n-1,k) is the number of distinct k-level rooted fuzzy subsets of an n-set ordered by set inclusion. (End)
The relations on p. 34 of Hasan (also p. 17 of Franco and Hasan) agree with the relation between A019538 and this entry given in the formula section. - Tom Copeland, May 14 2020
T(n,k) is the size of the Green's L-classes in the D-classes of rank (k-1) in the semigroup of partial transformations on an (n-1)-set. - Geoffrey Critzer, Jan 09 2023
T(n,k) is the number of strongly connected binary relations on [n] that have period k (A367948) and index 1.  See Theorem 5.4.25(6) in Ki Hang Kim reference. - Geoffrey Critzer, Dec 07 2023

			The triangle a(n, k) starts:
n\k 1   2    3     4      5      6      7      8     9
1:  1
2:  1   1
3:  1   3    2
4:  1   7   12     6
5:  1  15   50    60     24
6:  1  31  180   390    360    120
7:  1  63  602  2100   3360   2520    720
8:  1 127 1932 10206  25200  31920  20160   5040
9:  1 255 6050 46620 166824 317520 332640 181440 40320
... [Reformatted by _Wolfdieter Lang_, Mar 26 2015]
-----------------------------------------------------
Row 5 of triangle is {1,15,50,60,24}, which is {1,15,25,10,1} times {0!,1!,2!,3!,4!}.
From _Vladimir Shevelev_, Dec 22 2011: (Start)
Also, for power sums, we have
S_0(n) = C(n,1);
S_1(n) = C(n,1) +    C(n,2);
S_2(n) = C(n,1) +  3*C(n,2) +  2*C(n,3);
S_3(n) = C(n,1) +  7*C(n,2) + 12*C(n,3) +  6*C(n,4);
S_4(n) = C(n,1) + 15*C(n,2) + 50*C(n,3) + 60*C(n,4) + 24*C(n,5); etc.
(End)
For X = [1,2,3], the sets T are {{}}, {{},{1,2}}, {{},{1,3}}, {{},{1,2,3}}, {{},{1,2},{1,2,3}}, {{},{1,3},{1,2,3}} and so a(3,1)=1, a(3,2)=3, a(3,3)=2. - _Michael Somos_, Apr 20 2013

Ki Hang Kim, Boolean Matrix Theory and Applications, Marcel Dekker, New York and Basel (1982).

Seiichi Manyama, Table of n, a(n) for n = 1..10000
V. S. Abramovich, Power sums of natural numbers, Kvant, no. 5 (1973), 22-25. (in Russian)
Peter Bala, Deformations of the Hadamard product of power series
Paul Barry, Three Études on a sequence transformation pipeline, arXiv:1803.06408 [math.CO], 2018.
Melike Baykal-Gürsoy, Marcelo Figueroa-Candia, and Zhe Duan, Completion times of jobs on two-state service processes and their asymptotic behavior, Probability in the Engineering and Information Sciences (2024), 1-29. See p. 23.
H. Belbachir, M. Rahmani, and B. Sury, Sums Involving Moments of Reciprocals of Binomial Coefficients, J. Int. Seq. 14 (2011) #11.6.6.
Hacene Belbachir and Mourad Rahmani, Alternating Sums of the Reciprocals of Binomial Coefficients, Journal of Integer Sequences, Vol. 15 (2012), #12.2.8.
F. Brenti and V. Welker, f-vectors of barycentric subdivisions, arXiv:math/0606356v1 [math.CO], Math. Z., 259(4), 849-865, 2008.
Patibandla Chanakya and Putla Harsha, Generalized Nested Summation of Powers of Natural Numbers, arXiv:1808.08699 [math.NT], 2018. See Table 1.
Tom Copeland, Generators, Inversion, and Matrix, Binomial, and Integral Transforms
Colin Defant, Troupes, Cumulants, and Stack-Sorting, arXiv:2004.11367 [math.CO], 2020.
E. Delucchi, A. Pixton and L. Sabalka. Face vectors of subdivided simplicial complexes arXiv:1002.3201v3 [math.CO], Discrete Mathematics, Volume 312, Issue 2, January 2012, Pages 248-257.
G. H. E Duchamp, N. Hoang-Nghia, and A. Tanasa, A word Hopf algebra based on the selection/quotient principle, arXiv:1207.6522 [math.CO], 2012-2013; Séminaire Lotharingien de Combinatoire 68 (2013), Article B68c.
M. Dukes and C. D. White, Web Matrices: Structural Properties and Generating Combinatorial Identities, arXiv:1603.01589 [math.CO], 2016.
Nick Early, Honeycomb tessellations and canonical bases for permutohedral blades, arXiv:1810.03246 [math.CO], 2018.
S. Franco and A. Hasan, Graded Quivers, Generalized Dimer Models and Toric Geometry , arXiv preprint arXiv:1904.07954 [hep-th], 2019
A. Hasan, Physics and Mathematics of Graded Quivers, dissertation, Graduate Center, City University of New York, 2019.
H. Hasse, Ein Summierungsverfahren für die Riemannsche Zeta-Reihe, Math. Z. 32, 458-464 (1930).
Guy Louchard, Werner Schachinger, and Mark Daniel Ward, The number of distinct adjacent pairs in geometrically distributed words: a probabilistic and combinatorial analysis, arXiv:2203.14773 [math.PR], 2022. See p. 5.
Shi-Mei Ma, A family of two-variable derivative polynomials for tangent and secant, El J. Combinat. 20(1) (2013), P11.
Richard J. Mathar, Integrals Associated with the Digamma Integral Representation, arXiv:2308.14154 [math.GM], 2023. See p. 5.
Rajesh Kumar Mohapatra and Tzung-Pei Hong, On the Number of Finite Fuzzy Subsets with Analysis of Integer Sequences, Mathematics (2022) Vol. 10, No. 7, 1161.
A. Riskin and D. Beckwith, Problem 10231, Amer. Math. Monthly, 102 (1995), 175-176.
G. Rzadkowski, Bernoulli numbers and solitons revisited, Journal of Nonlinear Mathematical Physics, Volume 17, Issue 1, 2010.
John K. Sikora, On Calculating the Coefficients of a Polynomial Generated Sequence Using the Worpitzky Number Triangles, arXiv:1806.00887 [math.NT], 2018.
G. J. Simmons, A combinatorial problem associated with a family of combination locks, Math. Mag., 37 (1964), 127-132 (but there are errors). The triangle is on page 129.
N. J. A. Sloane, Transforms
Sam Vandervelde, The Worpitzky Numbers Revisited, Amer. Math. Monthly, 125:3 (2018), 198-206.
Wikipedia, Bernoulli number.
Wikipedia, Barycentric subdivision
David C. Wood, The computation of polylogarithms (2014).

Dropping the column of 1's gives A053440.
Without the k in the denominator (in the definition), we get A019538. See also the Stirling number triangle A008277.
Cf. A087127, A087107, A087108, A087109, A087110, A087111, A084938 A075263.
Row sums give A000629(n-1) for n >= 1.
Cf. A027642, A002445. - Gary W. Adamson, Aug 09 2008
Appears in A161739 (RSEG2 triangle), A161742 and A161743. - Johannes W. Meijer, Jun 18 2009
Binomial transform is A038719. Cf. A131689.
Cf. A007318, A008292, A046802, A074909, A090582, A123125, A130850, A135278, A141618, A145271, A163626, A248727, A263634.
Cf. A119879.
From Rajesh Kumar Mohapatra, Mar 29 2020: (Start)
A000007(n-1) (column k=1), A000225(n-1) (column k=2), A028243(n-1) (column k=3), A028244(n-1) (column k=4), A028245(n-1) (column k=5), for n > 0.
Diagonal gives A000142(n-1), for n >=1.
Next-to-last diagonal gives A001710,
Third, fourth, fifth, sixth, seventh external diagonal respectively give A005460, A005461, A005462, A005463, A005464. (End)

Flat(List([1..10], n-> List([1..n], k-> Stirling2(n,k)* Factorial(k-1) ))); # G. C. Greubel, May 30 2019

[[StirlingSecond(n,k)*Factorial(k-1): k in [1..n]]: n in [1..10]]; // G. C. Greubel, May 30 2019

a := (n,k) -> add((-1)^(k-i)*binomial(k,i)*i^n, i=0..k)/k;
seq(print(seq(a(n,k),k=1..n)),n=1..10);
T := (n,k) -> add(eulerian1(n,j)*binomial(n-j,n-k), j=0..n):
seq(print(seq(T(n,k),k=0..n)),n=0..9); # Peter Luschny, Jul 12 2013

a[n_, k_] = Sum[(-1)^(k-i) Binomial[k,i]*i^n, {i,0,k}]/k; Flatten[Table[a[n, k], {n, 10}, {k, n}]] (* Jean-François Alcover, May 02 2011 *)

{T(n, k) = if( k<0 || k>n, 0, n! * polcoeff( (x / log(1 + x + x^2 * O(x^n) ))^(n+1), n-k))}; /* Michael Somos, Oct 02 2002 */

{T(n,k) = stirling(n,k,2)*(k-1)!}; \\ G. C. Greubel, May 31 2019

# Assuming offset (n, k) = (0, 0).
def T(n, k):
    if k >  n: return 0
    if k == 0: return 1
    return k*T(n - 1, k - 1) + (k + 1)*T(n - 1, k)
for n in range(9):
    print([T(n, k) for k in range(n + 1)])  # Peter Luschny, Apr 26 2022

def A163626_row(n) :
    x = polygen(ZZ,'x')
    A = []
    for m in range(0, n, 1) :
        A.append((-x)^m)
        for j in range(m, 0, -1):
            A[j - 1] = j * (A[j - 1] - A[j])
    return list(A[0])
for i in (1..7) : print(A163626_row(i))  # Peter Luschny, Jan 25 2012

[[stirling_number2(n,k)*factorial(k-1) for k in (1..n)] for n in (1..10)] # G. C. Greubel, May 30 2019

E.g.f.: -log(1-y*(exp(x)-1)). - Vladeta Jovovic, Sep 28 2003
a(n, k) = S2(n, k)*(k-1)! where S2(n, k) is a Stirling number of the second kind (cf. A008277). Also a(n,k) = T(n,k)/k, where T(n, k) = A019538.
Essentially same triangle as triangle [1, 0, 2, 0, 3, 0, 4, 0, 5, 0, 6, 0, 7, ...] DELTA [1, 1, 2, 2, 3, 3, 4, 4, 5, 5, ...] where DELTA is Deléham's operator defined in A084938, but the notation is different.
Sum of terms in n-th row = A000629(n) - Gary W. Adamson, May 30 2005
The row generating polynomials P(n, t) are given by P(1, t)=t, P(n+1, t) = t(t+1)(d/dt)P(n, t) for n >= 1 (see the Riskin and Beckwith reference). - Emeric Deutsch, Aug 09 2005
From Gottfried Helms, Jul 12 2006: (Start)
Delta-matrix as can be read from H. Hasse's proof of a connection between the zeta-function and Bernoulli numbers (see link below).
Let P = lower triangular matrix with entries P[row,col] = binomial(row,col).
Let J = unit matrix with alternating signs J[r,r]=(-1)^r.
Let N(m) = column matrix with N(m)(r) = (r+1)^m, N(1)--> natural numbers.
Let V = Vandermonde matrix with V[r,c] = (r+1)^c.
V is then also N(0)||N(1)||N(2)||N(3)... (indices r,c always beginning at 0).
Then Delta = P*J * V and B' = N(-1)' * Delta, where B is the column matrix of Bernoulli numbers and ' means transpose, or for the single k-th Bernoulli number B_k with the appropriate column of Delta,
B_k = N(-1)' * Delta[ *,k ] = N(-1)' * P*J * N(k).
Using a single column instead of V and assuming infinite dimension, H. Hasse showed that in x = N(-1) * P*J * N(s), where s can be any complex number and s*zeta(1-s) = x.
His theorem reads: s*zeta(1-s) = Sum_{n>=0..inf} (n+1)^-1*delta(n,s), where delta(n,s) = Sum_{j=0..n} (-1)^j * binomial(n,j) * (j+1)^s.
(End)
a(n,k) = k*a(n-1,k) + (k-1)*a(n-1,k-1) with a(n,1) = 1 and a(n,n) = (n-1)!. - Johannes W. Meijer, Jun 18 2009
Rephrasing the Meijer recurrence above: Let M be the (n+1)X(n+1) bidiagonal matrix with M(r,r) = M(r,r+1) = r, r >= 1, in the two diagonals and the rest zeros. The row a(n+1,.) of the triangle is row 1 of M^n. - Gary W. Adamson, Jun 24 2011
From Tom Copeland, Oct 11 2011: (Start)
With e.g.f.. A(x,t) = G[(t+1)x,-1/(t+1)]-1 (from 2008 comment) = -1 + 1/[1-(1+t)(1-e^(-x))] = (1+t)x + (1+3t+2t^2)x^2/2! + ..., the comp. inverse in x is
B(x,t)= -log(t/(1+t)+1/((1+t)(1+x))) = (1/(1+t))x - ((1+2t)/(1+t)^2)x^2/2 + ((1+3t+3t^2)/(1+t)^3)x^3/3 + .... The numerators are the row polynomials of A074909, and the rational functions are (omitting the initial constants) signed columns of the re-indexed Pascal triangle A007318.
Let h(x,t)= 1/(dB/dx) = (1+x)(1+t(1+x)), then the row polynomial P(n,t) = (1/n!)(h(x,t)*d/dx)^n x, evaluated at x=0, A=exp(x*h(y,t)*d/dy) y, eval. at y=0, and dA/dx = h(A(x,t),t), with P(1,t)=1+t. (Series added Dec 29 2015.)(End)
Let  denote the Eulerian numbers A173018(n,k), then T(n,k) = Sum_{j=0..n} *binomial(n-j,n-k). - Peter Luschny, Jul 12 2013
Matrix product A007318 * A131689. The n-th row polynomial R(n,x) = Sum_{k >= 1} k^(n-1)*(x/(1 + x))^k, valid for x in the open interval (-1/2, inf). Cf A038719. R(n,-1/2) = (-1)^(n-1)*(2^n - 1)*Bernoulli(n)/n. - Peter Bala, Jul 14 2014
a(n,k) = A141618(n,k) / C(n,k-1). - Tom Copeland, Oct 25 2014
For the row polynomials, A028246(n,x) = A019538(n-1,x) * (1+x). - Tom Copeland, Dec 28 2015
A248727 = A007318*(reversed A028246) = A007318*A130850 = A007318*A123125*A007318 = A046802*A007318. - Tom Copeland, Nov 14 2016
n-th row polynomial R(n,x) = (1+x) o (1+x) o ... o (1+x) (n factors), where o denotes the black diamond multiplication operator of Dukes and White. See example E11 in the Bala link. - Peter Bala, Jan 12 2018
From Dorian Guyot, May 21 2019: (Start)
Sum_{i=0..k} binomial(k,i) * a(n,i) = (k+1)^n.
Sum_{k=0..n} a(n,k) = 2*A000670(n).
(End)
With all offsets 0, let A_n(x;y) = (y + E.(x))^n, an Appell sequence in y where E.(x)^k = E_k(x) are the Eulerian polynomials of A123125. Then the row polynomials of this entry, A028246, are given by x^n * A_n(1 + 1/x;0). Other specializations of A_n(x;y) give A046802, A090582, A119879, A130850, and A248727. - Tom Copeland, Jan 24 2020
The row generating polynomials R(n,x) = Sum_{i=1..n} a(n,i) * x^i satisfy the recurrence equation R(n+1,x) = R(n,x) + Sum_{k=0..n-1} binomial(n-1,k) * R(k+1,x) * R(n-k,x) for n >= 1 with initial value R(1,x) = x. - Werner Schulte, Jun 17 2021

Definition corrected by Li Guo, Dec 16 2006
Typo in link corrected by Johannes W. Meijer, Oct 17 2009
Error in title corrected by Johannes W. Meijer, Sep 24 2010
Edited by M. F. Hasler, Oct 29 2014

A046802 T(n, k) = Sum_{j=k..n} binomial(n, j)*E1(j, j-k), where E1 are the Eulerian numbers A173018. Triangle read by rows, T(n, k) for 0 <= k <= n.

Original entry on oeis.org

1, 1, 1, 1, 3, 1, 1, 7, 7, 1, 1, 15, 33, 15, 1, 1, 31, 131, 131, 31, 1, 1, 63, 473, 883, 473, 63, 1, 1, 127, 1611, 5111, 5111, 1611, 127, 1, 1, 255, 5281, 26799, 44929, 26799, 5281, 255, 1, 1, 511, 16867, 131275, 344551, 344551, 131275, 16867, 511, 1, 1, 1023, 52905

Offset: 0

N. J. A. Sloane

T(n,k) is the number of positroid cells of the totally nonnegative Grassmannian G+(k,n) (cf. Postnikov/Williams). It is the triangle of the h-vectors of the stellahedra. - Tom Copeland, Oct 10 2014
See A248727 for a simple transformation of the row polynomials of this entry that produces the umbral compositional inverses of the polynomials of A074909, related to the face polynomials of the simplices. - Tom Copeland, Jan 21 2015
From Tom Copeland, Jan 24 2015: (Start)
The reciprocal of this entry's e.g.f. is [t e^(-xt) - e^(-x)] / (t-1) = 1 - (1+t) x + (1+t+t^2) x^2/2! - (1+t+t^2+t^3) x^3/3! + ... = e^(q.(0;t)x), giving the base sequence (q.(0;t))^n = q_n(0;t) = (-1)^n [1-t^(n+1)] / (1-t) for the umbral compositional inverses (q.(0;t)+z)^n = q_n(z;t) of the Appell polynomials associated with this entry, p_n(z;t) below, i.e., q_n(p.(z;t)) = z^n = p_n(q.(z;t)), in umbral notation. The relations in A133314 then apply between the two sets of base polynomials. (Inserted missing index in a formula - Mar 03 2016.)
The associated o.g.f. for the umbral inverses is Og(x) = x / (1-x q.(0:t)) = x / [(1+x)(1+tx)] = x / [1+(1+t)x+tx^2]. Applying A134264 to h(x) = x / Og(x) = 1 + (1+t) x + t x^2 leads to an o.g.f. for the Narayana polynomials A001263 as the comp. inverse Oginv(x) = [1-(1+t)x-sqrt[1-2(1+t)x+((t-1)x)^2]] / (2xt). Note that Og(x) gives the signed h-polynomials of the simplices and that Oginv(x) gives the h-polynomials of the simplicial duals of the Stasheff polynomials, or type A associahedra. Contrast this with A248727 = A046802 * A007318, which has o.g.f.s related to the corresponding f-polynomials. (End)
The Appell polynomials p_n(x;t) in the formulas below specialize to the Swiss-knife polynomials of A119879 for t = -1, so the Springer numbers A001586 are given by 2^n p_n(1/2;-1). - Tom Copeland, Oct 14 2015
The row polynomials are the h-polynomials associated to the stellahedra, whose f-polynomials are the row polynomials of A248727. Cf. page 60 of the Buchstaber and Panov link. - Tom Copeland, Nov 08 2016
The row polynomials are the h-polynomials of the stellohedra, which enumerate partial permutations according to descents. Cf. Section 10.4 of the Postnikov-Reiner-Williams reference. - Lauren Williams, Jul 05 2022
From p. 60 of the Buchstaber and Panov link, S = P * C / T where S, P, C, and T are the bivariate e.g.f.s of the h vectors of the stellahedra, permutahedra, hypercubes, and (n-1)-simplices, respectively. - Tom Copeland, Jan 09 2017
The number of Le-diagrams of type (k, n) this means the diagram uses the bounding box size k x (n-k), equivalently the number of Grassmann necklaces of type (k, n) and also the number of decorated permutations with k anti-exceedances. - Thomas Scheuerle, Dec 29 2024

			The triangle T(n, k) begins:
n\k 0   1     2      3      4      5      6     7
0:  1
1:  1   1
2:  1   3     1
3:  1   7     7      1
4:  1  15    33     15      1
5:  1  31   131    131     31      1
6:  1  63   473    883    473     63      1
7:  1 127  1611   5111   5111   1611    127     1
... Reformatted. - _Wolfdieter Lang_, Feb 14 2015

L. Comtet, Advanced Combinatorics, Reidel, Holland, 1974, page 245 [From Roger L. Bagula, Nov 21 2009]
D. Singh, The numbers L(m,n) and their relations with prepared Bernoulli and Eulerian numbers, Math. Student, 20 (1952), 66-70.

Michael De Vlieger, Table of n, a(n) for n = 0..11475 (rows 0 <= n <= 150, flattened)
Paul Barry, General Eulerian Polynomials as Moments Using Exponential Riordan Arrays, Journal of Integer Sequences, 16 (2013), #13.9.6.
Paul Barry, Three Études on a sequence transformation pipeline, arXiv:1803.06408 [math.CO], 2018.
Paul Barry, Series reversion with Jacobi and Thron continued fractions, arXiv:2107.14278 [math.NT], 2021.
Carolina Benedetti, Anastasia Chavez, and Daniel Tamayo, Quotients of uniform positroids, arXiv:1912.06873 [math.CO], 2019.
V. Buchstaber and T. Panov Toric Topology, arXiv:1210.2368v3 [math.AT], 2014.
Jan Geuenich, Tilting modules for the Auslander algebra of K[x]/(xn), arXiv:1803.10707 [math.RT], 2018.
Jun Ma, Shi-mei Ma, and Yeong-Nan Yeh, Recurrence relations for binomial-eulerian polynomials, arXiv:1711.09016 [math.CO], 2017, Thm. 2.1.
A. Postnikov, Total positivity, Grassmannians, and networks, arXiv:math/0609764 [math.CO], 2006.
A. Postnikov, V. Reiner, and L. Williams, Faces of generalized permutohedra, arXiv:math/0609184 [math.CO], 2006-2007.
D. Singh, The numbers L(m,n) and their relations with prepared Bernoulli and Eulerian numbers, Math. Student, 20 (1952), 66-70. [Annotated scanned copy]
L. K. Williams, Enumeration of totally positive Grassmann cells, arXiv:math/0307271 [math.CO], 2003-2004.
L. Williams, The Positive Grassmannian (from a mathematician's perspective), 2014

Row sums give A000522.
Cf. A008292, A123125, A248727, A074909, A007318, A000225, A066810, A028246, A001263, A119879, A001586, A019538, A090582, A123125, A130850.
Cf. A122045, A007047, A000522, A026898, A036919.

T := (n, k) -> add(binomial(n, r)*combinat:-eulerian1(r, r-k), r = k .. n):
for n from 0 to 8 do seq(T(n, k), k=0..n) od; # Peter Luschny, Jun 27 2018

t[, 1] = 1; t[n, n_] = 1; t[n_, 2] = 2^(n-1)-1;
t[n_, k_] = Sum[((i-k+1)^i*(k-i)^(n-i-1) - (i-k+2)^i*(k-i-1)^(n-i-1))*Binomial[n-1, i], {i, 0, k-1}];
T[n_, k_] := t[n+1, k+1]; Table[T[n, k], {n, 0, 12}, {k, 0, n}] // Flatten
(* Jean-François Alcover, Jan 22 2015, after Tom Copeland *)
T[ n_, k_] := Coefficient[n! SeriesCoefficient[(1-x) Exp[t] / (1 - x Exp[(1-x) t]), {t, 0, n}] // Simplify, x, k];
Table[T[n, k], {n, 0, 10}, {k, 0, n}] (* Michael Somos, Jan 22 2015 *)

E.g.f.: (y-1)*exp(x*y)/(y-exp((y-1)*x)). - Vladeta Jovovic, Sep 20 2003
p(t,x) = (1 - x)*exp(t)/(1 - x*exp(t*(1 - x))). - Roger L. Bagula, Nov 21 2009
With offset=0, T(n,0)=1 otherwise T(n,k) = sum_{i=0..k-1} C(n,i)((i-k)^i*(k-i+1)^(n-i) - (i-k+1)^i*(k-i)^(n-i)) (cf. Williams). - Tom Copeland, Oct 10 2014
With offset 0, T = A007318 * A123125. Second column is A000225; 3rd, appears to be A066810. - Tom Copeland, Jan 23 2015
A raising operator (with D = d/dx) associated with this entry's row polynomials is R = x + t + (1-t) / [1-t e^{(1-t)D}]  = x + t + 1 + t D + (t+t^2) D^2/2! + (t+4t^2+t^3) D^3/3! + ... , containing the e.g.f. for the Eulerian polynomials of A123125. Then R^n 1 = (p.(0;t)+x)^n = p_n(x;t) are the Appell polynomials with this entry's row polynomials p_n(0;t) as the base sequence. Examples of this formalism are given in A028246 and A248727. - Tom Copeland, Jan 24 2015
With offset 0, T = A007318*(padded A090582)*(inverse of A097805) = A007318*(padded A090582)*(padded A130595) = A007318*A123125 = A007318*(padded A090582)*A007318*A097808*A130595, where padded matrices are of the form of padded A007318, which is A097805. Inverses of padded matrices are just the padded versions of inverses of the unpadded matrices. This relates the face vectors, or f-vectors, and h-vectors of the permutahedra / permutohedra to those of the stellahedra / stellohedra. - Tom Copeland, Nov 13 2016
Umbrally, the row polynomials (offset 0) are r_n(x) = (1 + q.(x))^n, where (q.(x))^k = q_k(x) are the row polynomials of A123125. - Tom Copeland, Nov 16 2016
From the previous umbral statement, OP(x,d/dy) y^n =  (y + q.(x))^n, where OP(x,y) = exp[y * q.(x)] = (1-x)/(1-x*exp((1-x)y)), the e.g.f. of A123125, so OP(x,d/dy) y^n evaluated at y = 1 is  r_n(x), the n-th row polynomial of this entry, with offset 0. - Tom Copeland, Jun 25 2018
Consolidating some formulas in this entry and A248727, in umbral notation for concision, with all offsets 0: Let A_n(x;y) = (y + E.(x))^n, an Appell sequence in y where E.(x)^k = E_k(x) are the Eulerian polynomials of A123125. Then the row polynomials of this entry (A046802, the h-polynomials of the stellahedra) are given by h_n(x) = A_n(x;1); the row polynomials of A248727 (the face polynomials of the stellahedra), by f_n(x) = A_n(1 + x;1); the Swiss-knife polynomials of A119879, by Sw_n(x) = A_n(-1;1 + x); and the row polynomials of the Worpitsky triangle (A130850), by w_n(x) = A(1 + x;0). Other specializations of A_n(x;y) give A090582 (the f-polynomials of the permutohedra, cf. also A019538) and A028246 (another version of the Worpitsky triangle). - Tom Copeland, Jan 24 2020
From Peter Luschny, Apr 30 2021: (Start)
Sum_{k=0..n} (-1)^k*T(n, k) = A122045(n).
Sum_{k=0..n} 2^(n-k)*T(n,k) = A007047(n).
Sum_{k=0..n} T(n, n-k) = A000522(n).
Sum_{k=0..n} T(n-k, k) = Sum_{k=0..n} (n - k)^k = A026898(n-1) for n >= 1.
Sum_{k=0..n} k*T(n, k) = A036919(n) = floor(n*n!*e/2).
(End)

More terms from Vladeta Jovovic, Sep 20 2003
First formula corrected by Wolfdieter Lang, Feb 14 2015
Offset set to 0 and edited by Peter Luschny, Apr 30 2021

A090582 T(n, k) = Sum_{j=0..n-k} (-1)^jbinomial(n - k + 1, j)(n - k + 1 - j)^n. Triangle read by rows, T(n, k) for 1 <= k <= n.

Original entry on oeis.org

1, 2, 1, 6, 6, 1, 24, 36, 14, 1, 120, 240, 150, 30, 1, 720, 1800, 1560, 540, 62, 1, 5040, 15120, 16800, 8400, 1806, 126, 1, 40320, 141120, 191520, 126000, 40824, 5796, 254, 1, 362880, 1451520, 2328480, 1905120, 834120, 186480, 18150, 510, 1, 3628800, 16329600, 30240000, 29635200, 16435440, 5103000, 818520, 55980, 1022, 1

Offset: 1

Hugo Pfoertner, Jan 11 2004

Let Q(m, n) = Sum_(k=0..n-1) (-1)^k * binomial(n, k) * (n-k)^m. Then Q(m,n) is the numerator of the probability P(m,n) = Q(m,n)/n^m of seeing each card at least once if m >= n cards are drawn with replacement from a deck of n cards, written in a two-dimensional array read by antidiagonals with Q(m,m) as first row and Q(m,1) as first column.
The sequence is given as a matrix with the first row containing the cases #draws = size_of_deck. The second row contains #draws = 1 + size_of_deck. If "mn" indicates m cards drawn from a deck with n cards then the locations in the matrix are:
  11 22 33 44 55 66 77 ...
  21 32 43 54 65 76 87 ...
  31 42 53 64 75 86 97 ...
  41 52 63 74 85 .. .. ...
read by antidiagonals ->:
  11, 22, 21, 33, 32, 31, 44, 43, 42, 41, 55, 54, 53, 52, ....
The probabilities are given by Q(m,n)/n^m:
.(m,n):.....11..22..21..33..32..31..44..43..42..41...55...54..53..52..51
.....Q:......1...2...1...6...6...1..24..36..14...1..120..240.150..30...1
...n^m:......1...4...1..27...8...1.256..81..16...1.3125.1024.243..32...1
%.Success:.100..50.100..22..75.100...9..44..88.100....4...23..62..94.100
P(n,n) = n!/n^n which can be approximated by sqrt(Pi*(2n+1/3))/e^n (Gosper's approximation to n!).
Let P[n] be the set of all n-permutations. Build a superset Q[n] of P[n] composed of n-permutations in which some (possibly all or none) ascents have been designated. An ascent in a permutation s[1]s[2]...s[n] is a pair of consecutive elements s[i],s[i+1] such that s[i] < s[i+1]. As a triangular array read by rows T(n,k) is the number of elements in Q[n] that have exactly k distinguished ascents, n >= 1, 0 <= k <= n-1. Row sums are A000670. E.g.f. is y/(1+y-exp(y*x)). For example, T(3,1)=6 because there are four 3-permutations with one ascent, with these we would also count 1->2 3, and 1 2->3 where exactly one ascent is designated by "->". (After Flajolet and Sedgewick.) - Geoffrey Critzer, Nov 13 2012
Sum_(k=1..n) Q(n, k)*binomial(x, k) = x^n such that Sum_{k=1..n} Q(n, i)*binomial(x+1,i+1) = Sum_{k=1..x} k^n. - David A. Corneth, Feb 17 2014
A141618(n,n-k+1) = a(n,k) * C(n,k-1) / k. - Tom Copeland, Oct 25 2014
See A074909 and above g.f.s below for associations among this array and the Bernoulli polynomials and their umbral compositional inverses. - Tom Copeland, Nov 14 2014
For connections to toric varieties and Eulerian polynomials (in addition to those noted below), see the Dolgachev and Lunts and the Stembridge links in A019538. - Tom Copeland, Dec 31 2015
See A008279 for a relation between the e.g.f.s enumerating the faces of permutahedra (this entry) and stellahedra. - Tom Copeland, Nov 14 2016
From the Hasan and Franco and Hasan papers: The m-permutohedra for m=1,2,3,4 are the line segment, hexagon, truncated octahedron and omnitruncated 5-cell. The first three are well-known from the study of elliptic models, brane tilings and brane brick models. The m+1 torus can be tiled by a single (m+2)-permutohedron. Relations to toric Calabi-Yau Kahler manifolds are also discussed. - Tom Copeland, May 14 2020

			For m = 4, n = 2, we draw 4 times from a deck of two cards. Call the cards "a" and "b" - of the 16 possible combinations of draws (each of which is equally likely to occur), only two do not contain both a and b: a, a, a, a and b, b, b, b. So the probability of seeing both a and b is 14/16. Therefore Q(m, n) = 14.
Table starts:
  [1] 1;
  [2] 2,      1;
  [3] 6,      6,       1;
  [4] 24,     36,      14,      1;
  [5] 120,    240,     150,     30,      1;
  [6] 720,    1800,    1560,    540,     62,     1;
  [7] 5040,   15120,   16800,   8400,    1806,   126,    1;
  [8] 40320,  141120,  191520,  126000,  40824,  5796,   254,   1;
  [9] 362880, 1451520, 2328480, 1905120, 834120, 186480, 18150, 510, 1.

Reinhard Zumkeller, Rows n = 1..125 of triangle, flattened
Paul Barry, Three Études on a sequence transformation pipeline, arXiv:1803.06408 [math.CO], 2018.
Kevin Buhr, Michael Press and DMB, Collecting a deck of cards on the street. Thread in NG sci.math.
José L. Cereceda, Figurate numbers and sums of powers of integers, arXiv:2001.03208 [math.NT], 2020. See triangle Table 1 p. 5.
Tom Copeland, Generators, Inversion, and Matrix, Binomial, and Integral Transforms
P. Flajolet and R. Sedgewick, Analytic Combinatorics, 2009, page 209.
S. Franco and A. Hasan, Graded Quivers, Generalized Dimer Models and Toric Geometry, arXiv preprint arXiv:1904.07954 [hep-th], 2019.
A. Hasan, Physics and Mathematics of Graded Quivers, dissertation, Graduate Center, City University of New York, 2019.
Hugo Pfoertner, "Hit all" probabilities: Program and results.

Cf. A073593 first m >= n giving at least 50% probability, A085813 ditto for 95%, A055775 n^n/n!, A090583 Gosper's approximation to n!.
Reflected version of A019538.
Cf. A233734 (central terms).
Cf. A074909, A008292, A028246.
Cf. A046802, A119879, A123125, A130850, and A248727.

a090582 n k = a090582_tabl !! (n-1) !! (k-1)
a090582_row n = a090582_tabl !! (n-1)
a090582_tabl = map reverse a019538_tabl
-- Reinhard Zumkeller, Dec 15 2013

T := (n, k) -> add((-1)^j*binomial(n - k + 1, j)*(n - k + 1 - j)^n, j = 0..n-k):
# Or:
T := (n, k) -> (n - k + 1)!*Stirling2(n, n - k + 1):
for n from 1 to 9 do seq( T(n, k), k = 1..n) od; # Peter Luschny, May 21 2021

In[1]:= Table[Table[k! StirlingS2[n, k], {k, n, 1, -1}], {n, 1, 6}] (* Victor Adamchik, Oct 05 2005 *)
nn=6; a=y/(1+y-Exp[y x]); Range[0,nn]! CoefficientList[Series[a, {x,0,nn}], {x,y}]//Grid (* Geoffrey Critzer, Nov 10 2012 *)

a(n)={m=ceil((-1+sqrt(1+8*n))/2);k=m+1+binomial(m,2)-n;k*sum(i=1,k,(-1)^(i+k)*i^(m-1)*binomial(k-1,i-1))} \\ David A. Corneth, Feb 17 2014

T(n, k) = (n - k + 1)!*Stirling2(n, n - k + 1), generated by Stirling numbers of the second kind multiplied by a factorial. - Victor Adamchik, Oct 05 2005
Triangle T(n,k), 1 <= k <= n, read by rows given by [1, 1, 2, 2, 3, 3, 4, 4, 5, 5, ...] DELTA [0, 1, 0, 2, 0, 3, 0, 4, 0, 5, 0, 6, 0, ...] where DELTA is the operator defined in A084938. - Philippe Deléham, Nov 10 2006
From Tom Copeland, Oct 07 2008: (Start)
G(x,t) = 1/ (1 + (1-exp(x*t))/t) = 1 + 1*x + (2 + t)*x^2/2! + (6 + 6*t + t^2)*x^3/3! + ... gives row polynomials of A090582, the f-polynomials for the permutohedra (see A019538).
G(x,t-1) = 1 + 1*x + (1 + t)*x^2/2! + (1 + 4*t + t^2)*x^3/3! + ... gives row polynomials for A008292, the h-polynomials of permutohedra.
G[(t+1)x,-1/(t+1)] = 1 + (1 + t)*x + (1 + 3*t + 2*t^2)*x^2/2! + ... gives row polynomials of A028246. (End)
From Tom Copeland, Oct 11 2011: (Start)
With e.g.f. A(x,t) = G(x,t) - 1, the compositional inverse in x is
B(x,t) = log((t+1)-t/(1+x))/t. Let h(x,t) = 1/(dB/dx) = (1+x)*(1+(1+t)x), then the row polynomial P(n,t) is given by (1/n!)*(h(x,t)*d/dx)^n x, evaluated at x=0, A=exp(x*h(y,t)*d/dy) y, eval. at y=0, and dA/dx = h(A(x,t),t). (End)
k <= 0 or k > n yields Q(n, k) = 0; Q(1,1) = 1; Q(n, k) = k * (Q(n-1, k) + Q(n-1, k-1)). - David A. Corneth, Feb 17 2014
T = A008292*A007318. - Tom Copeland, Nov 13 2016
With all offsets 0 for this entry, let A_n(x;y) = (y + E.(x))^n, an Appell sequence in y where E.(x)^k = E_k(x) are the Eulerian polynomials of A123125 with offsets -1 so that the array becomes A008292; i.e., we ignore the first row and first column of A123125. Then the row polynomials of this entry, A090582, are given by A_n(1 + x;0). Other specializations of A_n(x;y) give A028246, A046802, A119879, A130850, and A248727. - Tom Copeland, Jan 24 2020

A119879 Exponential Riordan array (sech(x),x).

Original entry on oeis.org

1, 0, 1, -1, 0, 1, 0, -3, 0, 1, 5, 0, -6, 0, 1, 0, 25, 0, -10, 0, 1, -61, 0, 75, 0, -15, 0, 1, 0, -427, 0, 175, 0, -21, 0, 1, 1385, 0, -1708, 0, 350, 0, -28, 0, 1, 0, 12465, 0, -5124, 0, 630, 0, -36, 0, 1, -50521, 0, 62325, 0, -12810, 0, 1050, 0, -45, 0, 1

Offset: 0

Paul Barry, May 26 2006

Row sums have e.g.f. exp(x)*sech(x) (signed version of A009006). Inverse of masked Pascal triangle A119467. Transforms the sequence with e.g.f. g(x) to the sequence with e.g.f. g(x)*sech(x).
Coefficients of the Swiss-Knife polynomials for the computation of Euler, tangent and Bernoulli number (triangle read by rows). Another version in A153641. - Philippe Deléham, Oct 26 2013
Relations to Green functions and raising/creation and lowering/annihilation/destruction operators are presented in Hodges and Sukumar and in Copeland's discussion of this sequence and 2020 pdf. - Tom Copeland, Jul 24 2020

			Triangle begins:
     1;
     0,    1;
    -1,    0,     1;
     0,   -3,     0,   1;
     5,    0,    -6,   0,   1;
     0,   25,     0, -10,   0,   1;
   -61,    0,    75,   0, -15,   0,   1;
     0, -427,     0, 175,   0, -21,   0,  1;
  1385,    0, -1708,   0, 350,   0, -28,  0,  1;

G. C. Greubel, Rows n=0..100 of triangle, flattened
Paul Barry, Riordan Arrays, Orthogonal Polynomials as Moments, and Hankel Transforms, J. Int. Seq. 14 (2011) # 11.2.2, example 28.
Tom Copeland, The Elliptic Lie Triad: Ricatti and KdV Equations, Infinigens, and Elliptic Genera, 2015.
Tom Copeland, Discussion of this sequence
Tom Copeland, SKipping over Dimensions, Juggling Zeros in the Matrix, 2020.
A. Hodges and C. V. Sukumar, Bernoulli, Euler, permutations and quantum algebras, Proc. R. Soc. A Oct. 2007 vol 463 no. 463 2086 2401-2414.
Peter Luschny, Additive decompositions of classical numbers via the Swiss-Knife polynomials [_Tom Copeland_, Oct 20 2015]
Miguel Méndez and Rafael Sánchez, On the combinatorics of Riordan arrays and Sheffer polynomials: monoids, operads and monops, arXiv:1707.00336 [math.CO], 2017, Section 4.3, Example 4.
Miguel A. Méndez and Rafael Sánchez Lamoneda, Monops, Monoids and Operads: The Combinatorics of Sheffer Polynomials, The Electronic Journal of Combinatorics 25(3) (2018), #P3.25.

Row sums are A155585. - Johannes W. Meijer, Apr 20 2011
Rows reversed: A081658.
Cf. A109449, A153641, A162660. - Philippe Deléham, Oct 26 2013
Cf. A000182, A046802, A119467, A133314.
Cf. A019538, A028246, A090582, A123125, A130850, A248727.

T := (n,k) -> binomial(n,k)*2^(n-k)*euler(n-k,1/2): # Peter Luschny, Jan 25 2009

T[n_, k_] := Binomial[n, k] 2^(n-k) EulerE[n-k, 1/2];
Table[T[n, k], {n, 0, 10}, {k, 0, n}] // Flatten (* Jean-François Alcover, Jun 20 2018, after Peter Luschny *)

{T(n,k) = binomial(n,k)*2^(n-k)*(2/(n-k+1))*(subst(bernpol(n-k+1, x), x, 1/2) - 2^(n-k+1)*subst(bernpol(n-k+1, x), x, 1/4))};
for(n=0,5, for(k=0,n, print1(T(n,k), ", "))) \\ G. C. Greubel, Feb 25 2019

@CachedFunction
def A119879_poly(n, x) :
    return 1 if n == 0  else add(A119879_poly(k, 0)*binomial(n, k)*(x^(n-k)-1+n%2) for k in range(n)[::2])
def A119879_row(n) :
    R = PolynomialRing(ZZ, 'x')
    return R(A119879_poly(n,x)).coeffs()  # Peter Luschny, Jul 16 2012
# Alternatively:

# uses[riordan_array from A256893]
riordan_array(sech(x), x, 9, exp=true) # Peter Luschny, Apr 19 2015

Number triangle whose k-th column has e.g.f. sech(x)*x^k/k!.
T(n,k) = C(n,k)*2^(n-k)*E_{n-k}(1/2) where C(n,k) is the binomial coefficient and E_{m}(x) are the Euler polynomials. - Peter Luschny, Jan 25 2009
The coefficients in ascending order of x^i of the polynomials p{0}(x) = 1 and p{n}(x) = Sum_{k=0..n-1; k even} binomial(n,k)*p{k}(0)*((n mod 2) - 1 + x^(n-k)). - Peter Luschny, Jul 16 2012
E.g.f.: exp(x*z)/cosh(x). - Peter Luschny, Aug 01 2012
Sum_{k=0..n} T(n,k)*x^k = A122045(n), A155585(n), A119880(n), A119881(n) for x = 0, 1, 2, 3 respectively. - Philippe Deléham, Oct 27 2013
With all offsets 0, let A_n(x;y) = (y + E.(x))^n, an Appell sequence in y where E.(x)^k = E_k(x) are the Eulerian polynomials of A123125. Then the row polynomials of A046802 (the h-polynomials of the stellahedra) are given by h_n(x) = A_n(x;1); the row polynomials of A248727 (the face polynomials of the stellahedra), by f_n(x) = A_n(1 + x;1); the Swiss-knife polynomials of this entry, A119879, by Sw_n(x) = A_n(-1;1 + x); and the row polynomials of the Worpitsky triangle (A130850), by w_n(x) = A(1 + x;0). Other specializations of A_n(x;y) give A090582 (the f-polynomials of the permutohedra, cf. also A019538) and A028246 (another version of the Worpitsky triangle). - Tom Copeland, Jan 24 2020
Triangle equals P*((I + P^2)/2)^(-1), where P denotes Pascal's triangle A007318. - Peter Bala, Mar 07 2024

A263634 Irregular triangle read by rows: row n gives coefficients of n-th logarithmic polynomial L_n(x_1, x_2, ...) with monomials sorted into standard order.

Original entry on oeis.org

1, -1, 1, 2, -3, 1, -6, 12, -4, -3, 1, 24, -60, 20, 30, -5, -10, 1, -120, 360, -120, -270, 30, 120, 30, -6, -15, -10, 1, 720, -2520, 840, 2520, -210, -1260, -630, 42, 210, 140, 210, -7, -21, -35, 1

Offset: 1

N. J. A. Sloane, Oct 29 2015

"Standard order" here means as produced by Maple's "sort" command.
From Petros Hadjicostas, May 27 2020: (Start)
According to the Maple help files for the "sort" command, polynomials in multiple variables are "sorted in total degree with ties broken by lexicographic order (this is called graded lexicographic order)."
Thus for example, x_1^2*x_3 = x_1*x_1*x_3 > x_1*x_2*x_2 = x_1*x_2^2, while x_1^2*x_4 = x_1*x_1*x_4 > x_1*x_2*x_3. (End)
Row sums are 0 (for n > 1). Numbers of terms in rows are partition numbers A000041.
From Tom Copeland, Nov 06 2015: (Start)
With the formal Taylor series f(x) = 1 + x[1] x + x[2] x^2/2! + ... , the partition polynomials of this entry give d[log(f(x))]/dx = L_1(x[1]) + L_2(x[1], x[2]) x + L_3(...) x^2/2! + ..., and the coefficients of the reduced polynomials with x[n] = t are signed A028246.
The raising operator R = x + d[log(f(D)]/dD = x + L_1(x[1]) + L_2[x[1], x[2]) D + L_3(x[1], x[2], x[3]) D^2/2! + ... with D = d/dx generates an Appell sequence of polynomials, given umbrally by P_n(x[1], ..., x[n]; x) = (x[.] + x)^n = Sum_{k=0..n} binomial(n,k) x[k] * x^(n-k) = R^n 1 with the e.g.f. f(t)*e^(x*t) = exp[t P.(x[1], ..., x[.]; x)]. P_0 = x[0] = 1.
The umbral compositional inverse Appell sequence is generated by R = x - d[log(f(D))]/dD with e.g.f. e^(x*t)/f(t) = exp[t IP.(x[1], ..., x[.]; x)], so umbrally IP_n(x[1], ..., x[n]; P.(x[1], ..., x[n]; x)) = x^n = P_n(x[1], ..., x[n]; IP.(x[1], ..., x[n]; x)). An unsigned array for the reduced IP_n(x[1], ..., x[n]; x) polynomials with IP_0 = x[0] = 1 and x[n] = -1 for n > 0 is A154921, for which f(t) = 2 - e^t. (End)
From Tom Copeland, Sep 08 2016: (Start)
The Appell formalism allows a matrix representation in the power basis x^n of the raising operator R that incorporates this array's partition polynomials L_n(x[1], ..., x[n]):
VP_(n+1) = VP_n * R = VP_n * XPS^(-1) * MX * XPS, where XPS is the matrix formed from multiplying the n-th diagonal of the Pascal matrix PS of A007318 by the indeterminate x[n], with x[0] = 1 for the main diagonal of ones, i.e., XPS[n,k] = PS[n,k] * x[n-k]; the matrix MX is A129185; the matrix XPS^(-1) is the inverse of XPS, which can be formed by multiplying the diagonals of the Pascal matrix by the partition polynomials IPT(n, x[1], ..., x[n]) of A133314, i.e., XPS^(-1)[n,k] = PS[n,k] * IPT(n-k, x[1], ...); and VP_n is the row vector in the power basis representing the Appell polynomial P_n(x) formed from the basic sequence of moments 1, x[1], x[2], ..., i.e., umbrally P_n(x) = (x[.] + x)^n = Sum_{k=0..n} binomial(n,k) * x[k] * x^(n-k).
Then R = XPS^(-1) * MX * XPS is the Pascal matrix PS with an additional first superdiagonal of ones and the other lower diagonals multiplied by the partition polynomials of this array, i.e., R[n,k] = PS[n,k] * L_{n+1-k}(x[1], ..., x[n+1-k]) except for the first superdiagonal of ones.
Consistently, VP_n = (1, 0, 0, ...) * R^n = (1, 0, 0, ...) * XPS^(-1) * MX^n * XPS = (1, 0, 0, ...) * MX^n * XPS = the n-th row vector of XPS, which is the vector representation of P_n(x) = (x[.] + x)^n with x[0] = 1.
See the Copeland link for the umbral representation R = exp[g.*D] * x * exp[h.*D] that reflects the matrix representations.
The Stirling partition polynomials of the first kind St1_n(a[1], a[2], ..., a[n]) of A036039, the Stirling partition polynomials of the second kind St2_n(b[1], b[2], ..., b[n]) of A036040, and the refined Lah polynomials Lah_n[c[1], c[2], ..., c[n]) of A130561 are Appell sequences in the respective distinguished indeterminates a[1], b[1], and c[1]. Comparing the formulas for their raising operators with that in this entry, L_n(x[1], x[2], ..., x[n]) evaluates to
A) (n-1)! * a[n] for x[n] = St1_n(a[1], a[2], ..., a[n]);
B) b[n] for x[n] = St2_n(b[1], b[2], ..., b[n]);
C) n! * c[n] for x[n] = Lah_n(c[1], c[2], ..., c[n]).
Conversely, from the respective e.g.f.s (added Sep 12 2016)
D) x[n] = St1_n(L_1(x[1])/0!, ..., L_n(x[1], ..., x[n])/(n-1)!);
E) x[n] = St2_n(L_1(x[1]), ..., L_n(x[1], ..., x[n]));
F) x[n] = Lah_n(L_1(x[1])/1!, ..., L_n(x[1], ..., x[n])/n!).
Given only the Appell sequence with no closed form for the e.g.f., the raising operator can be generated using this formalism, as has been partially done for A134264. (End)
For the Appell sequences above, the raising operator is related to the recursion P_(n+1)(x) = x * P_n(x) + Sum_{k=0..n} binomial(n,k) * L_(n-k+1)(x[1], ..., x[n+k-1]) * P_k(x). For a derivation and connections to formal cumulants (c_n = L_n(x[1], ...)) and moments (m_n = x[n]), see the Copeland link on noncrossing partitions. With x = 0, the recursion reduces to x[n+1] =  Sum_{k = 0..n} binomial(n,k) * L_(n-k+1)(x[1], ..., x[n+k-1]) * x[k] with x[0] = 1. This array is a differently ordered version of A127671. - Tom Copeland, Sep 13 2016
With x[n] = x^(n-1), a signed version of A130850 is obtained. - Tom Copeland, Nov 14 2016
See p. 2 of Getzler for a relation to stable graphs called necklaces used in computations for Deligne-Mumford-Knudsen moduli spaces of stable curves of genus 1. - Tom Copeland, Nov 15 2019
For a relation to a combinatorial Faa di Bruno Hopf algebra related to functional composition, as presented by Connes and Moscovici, see Figueroa et al. - Tom Copeland, Jan 17 2020
From Tom Copeland, May 17 2020: (Start)
The e.g.f. of an Appell sequence is f(t) e^(x*t) with f(0) = 1. Given the Laguerre-Polya class function f(t) = e^(-a*t^2 + b*t) Product_m (1 - t/z_m) e^(t/z_m) with a = 0 for simplicity (more generally a >= 0) and b real and where the product runs over all the zeros z_m of f(t) with all zeros real and Sum_m 1/(z_m)^2 convergent, the raising operator of the Appell polynomials is R = x + b - Sum_{k > 0} c_(k+1) D^k with c_k = Sum_m (1/(z_m)^k), i.e., traces of powers of the reciprocals of the zeros. From R in earlier comments, b = L_1(x_1) and otherwise c_k = -L_k(x_1, ..., x_k).
The Laguerre / Turan / de Gua inequalities (Csordas and Williamson and Skovgaard) imply that all the zeros of each Appell polynomial are real and simple and its extrema are local maxima above the x-axis and local minima below and are located above or below the zeros of the next lower degree Appell polynomial. (End)
From Tom Copeland, Oct 15 2020: (Start)
With a_n = n! * b_n = (n-1)! * c_n for n  > 0, represent a function with f(0) = a_0 = b_0 = 1 as an
A) exponential generating function (e.g.f), or formal Taylor series: f(x) = e^{a.x} = 1 + Sum_{n > 0} a_n * x^n/n!
B) ordinary generating function (o.g.f.), or formal power series: f(x)  = 1/(1-b.x) = 1 + Sum_{n > 0}  b_n * x^n
C) logarithmic generating function (l.g.f): f(x) = 1 - log(1 - c.x) = 1 + Sum_{n > 0}  c_n * x^n /n.
Expansions of log(f(x)) are given in
I) A127671 and A263634 for the e.g.f: log[ e^{a.*x} ] =  e^{L.(a_1,a_2,...)x} = Sum_{n > 0} L_n(a_1,...,a_n) * x^n/n!, the logarithmic polynomials, cumulant expansion polynomials
II) A263916 for the o.g.f.: log[ 1/(1-b.x) ] =  log[ 1 - F.(b_1,b_2,...)x ] = -Sum_{n > 0} F_n(b_1,...,b_n) * x^n/n, the Faber polynomials.
Expansions of exp(f(x)-1) are given in
III) A036040 for an e.g.f: exp[ e^{a.x} - 1 ] = e^{BELL.(a_1,...)x}, the Bell/Touchard/exponential partition polynomials, a.k.a. the Stirling partition polynomials of the second kind
IV) A130561 for an o.g.f.: exp[ b.x/(1-b.x) ] = e^{LAH.(b.,...)x}, the Lah partition polynomials
V) A036039 for an l.g.f.: exp[ -log(1-c.x) ] =  e^{CIP.(c_1,...)x}, the cycle index polynomials of the symmetric groups S_n, a.k.a. the Stirling partition polynomials of the first kind.
Since exp and log are a compositional inverse pair, one can extract the indeterminates of the log set of partition polynomials from the exp set and vice versa. For a discussion of the relations among these polynomials and the combinatorics of connected and disconnected graphs/maps, see Novak and LaCroix on classical moments and cumulants and the two books on statistical mechanics referenced in A036040. (End)
Ignoring signs, these polynomials appear in Schröder in the set of equations (II) on p. 343 and in Stewart's translation on p. 31. - Tom Copeland, Aug 25 2021

			The first few polynomials are:
(1) x[1].
(2) -x[1]^2 + x[2].
(3) 2*x[1]^3 - 3*x[1]*x[2] + x[3].
(4) -6*x[1]^4 + 12*x[1]^2*x[2] - 4*x[1]*x[3] - 3*x[2]^2 + x[4].
(5) 24*x[1]^5 - 60*x[1]^3*x[2] + 20*x[1]^2*x[3] + 30*x[1]*x[2]^2 - 5*x[1]*x[4] - 10*x[2]*x[3] + x[5].
(6) -120*x[1]^6 + 360*x[1]^4*x[2] - 120*x[1]^3*x[3] - 270*x[1]^2*x[2]^2 + 30*x[1]^2*x[4] + 120*x[1]*x[2]*x[3] + 30*x[2]^3 - 6*x[1]*x[5] - 15*x[2]*x[4] - 10*x[3]^2 + x[6].
...
[1]    1
[2]   -1,    1
[3]    2,   -3,     1
[4]   -6,   12,    -4,    -3,   1
[5]   24,  -60,    20,    30,  -5,  -10,   1
[6] -120,  360,  -120,  -270,  30,  120,  30, -6, -15, -10, 1

L. Comtet, Advanced Combinatorics, Reidel, 1974, pp. 140, 156, 308.

Peter Luschny, Row n for n = 1..20.
Tom Copeland, Appell polynomials, cumulants, noncrossing partitions, Dyck paths, and inversion, 2014.
Tom Copeland, The creation / raising operators for Appell sequences, 2015.
G. Csordas and J. Williamson, The zeros of the Jensen polynomials are simple, Proceed. of the AMS, 49(1) (1975), 263-264.
H. Figueroa, J. Gracia-Bondia, and J. Varilly, Faa di Bruno Hopf algebras, arXiv:0508337 [math.CO], 2005; see p. 3.
E. Getzler, The semi-classical approximation for modular operads, arXiv:alg-geom/9612005, 1996; see p. 2.
J. Novak and M. LaCroix, Three lectures on free probability, arXiv:1205.2097 [math.CO], 2012.
E. Schröder, Ueber unendlich viele Algorithmen zur Auflösung der Gleichungen, Mathematische Annalen vol. 2, 317-365, 1870.
H. Skovgaard, On inequalities of the Turan type, Math. Scand. 2 (1954), 65-73.
G. Stewart, On infinitely many algorithms for solving equations, 1993, (translation into English of Schröder's paper above)

Cf. A000041, A178867, A263633, A028246, A007318, A129185, A133314, A154921.
Cf. A036039, A036040, A127671, A130561, A130850, A134264.
Cf. A263916.

triangle := proc(numrows) local E, s, Q;
E := add(x[i]*t^i/i!, i=1..numrows);
s := series(log(1 + E), t, numrows+1);
Q := k -> sort(expand(k!*coeff(s, t, k)));
seq(print(coeffs(Q(k))), k=1..numrows) end:
triangle(6); # updated by Peter Luschny, May 27 2020

G.f.: Log(1 + Sum_{i >= 1} x_i*t^i/i!) = Sum_{n >= 1} L_n(x_1, x_2, ...)*t^n/n!. [Comtet, p. 140, Eq. [5a]. - corrected by Tom Copeland, Sep 08 2016]

Conjecture: row polynomials are R(n,1) for n > 0 where R(n,k) = R(n-1,k+1) - Sum_{j=1..n-1} binomial(n-2,j-1)*R(j,k)*R(n-j,1) for n > 1, k > 0 with R(1,k) = x_k for k > 0. - Mikhail Kurkov, Mar 30 2025

A248727 A046802(x,y) --> A046802(x,y+1), transform of e.g.f. for the graded number of positroids of the totally nonnegative Grassmannians G+(k,n); enumerates faces of the stellahedra.

Original entry on oeis.org

1, 2, 1, 5, 5, 1, 16, 24, 10, 1, 65, 130, 84, 19, 1, 326, 815, 720, 265, 36, 1, 1957, 5871, 6605, 3425, 803, 69, 1, 13700, 47950, 65646, 44240, 15106, 2394, 134, 1, 109601, 438404, 707840, 589106, 267134, 63896, 7094, 263, 1

Offset: 0

Tom Copeland, Oct 12 2014

This is a transform of A046802 treating it as an array of h-vectors, so y is replaced by (y+1) in the e.g.f. for A046802.
An e.g.f. for the reversed row polynomials with signs is given by exp(a.(0;t)x) = [e^{(1+t)x} [1+t(1-e^(-x))]]^(-1) = 1 - (1+2t)x + (1+5t+5t^2)x^2/2! + ... . The reciprocal is an e.g.f. for the reversed face polynomials of the simplices A074909, i.e., exp(b.(0;t)x) = e^{(1+t)x} [1+t(1-e^(-x))] = 1 + (1+2t)x +(1+3t+3t^2) x^2/2! + ... , so the relations of A133314 apply between the two sets of polynomials. In particular, umbrally [a.(0;t)+b.(0;t)]^n vanishes except for n=0 for which it's unity, implying the two sets of Appell polynomials formed from the two bases, a_n(z;t) = (a.(0;t)+z)^n and b_n(z;t) = (b.(0;t) + z)^n, are an umbral compositional inverse pair, i.e., b_n(a.(x;t);t)= x^n = a_n(b.(x;t);t). Raising operators for these Appell polynomials are related to the polynomials of A028246, whose reverse polynomials are given by A123125 * A007318. Compare: A248727 = A007318 * A123125 * A007318 and A046802 = A007318 * A123125. See A074909 for definitions and related links. - Tom Copeland, Jan 21 2015
The o.g.f. for the umbral inverses is Og(x) = x / (1 - x b.(0;t)) = x / [(1-tx)(1-(1+t)x)] = x + (1+2t) x^2 + (1+3t+3t^2) x^3 + ... . Its compositional inverse is an o.g.f for signed A033282, the reverse f-polynomials for the simplicial duals of the Stasheff polytopes, or associahedra of type A, Oginv(x) =[1+(1+2t)x-sqrt[1+2(1+2t)x+x^2]] / (2t(1+t)x) = x - (1+2t) x^2 + (1+5t+5t^2) x^3 + ... . Contrast this with the o.g.f.s related to the corresponding h-polynomials in A046802. - Tom Copeland, Jan 24 2015
Face vectors, or coefficients of the face polynomials, of the stellahedra, or stellohedra. See p. 59 of Buchstaber and Panov. - Tom Copeland, Nov 08 2016
See A008279 for a relation between the e.g.f.s enumerating the faces of permutahedra and stellahedra. - Tom Copeland, Nov 14 2016

			The triangle T(n, k) starts:
n\k    0     1     2     3     4    5   6  7 ...
1:     1
2:     2     1
3:     5     5     1
4:    16    24    10     1
5:    65   130    84    19     1
6:   326   815   720   265    36    1
7:  1957  5871  6605  3425   803   69   1
8: 13700 47950 65646 44240 15106 2394 134  1
... reformatted, _Wolfdieter Lang_, Mar 27 2015

P. Barry, Three Études on a sequence transformation pipeline, arXiv:1803.06408 [math.CO], 2018.
L. Berry, S. Forcey, M. Ronco, and P. Showers, Polytopes and Hopf algebras of painted trees: Fan graphs and Stellohedra, arXiv:1608.08546 [math.CO], 2018.
L. Berry, S. Forcey, M. Ronco, and P. Showers, Species substitution, graph suspension, and graded Hopf algebras of painted tree polytopes, arXiv:1608.08546 [math.CO], 2019.
V. Buchstaber and T. Panov, Toric Topology, arXiv:1210.2368v3 [math.AT], 2014.
R. Da Rosa, D. Jensen, and D. Ranganathan, Toric graph associahedra and compactifications of M_(0,n), arXiv:1411.0537 [math.AG], 2015.
S. Forcey, M. Ronco, and P. Showers, Polytopes and algebras of grafted trees: Stellohedra, arXiv:1608.08546v2 [math.CO], 2016.
Stefan Forcey, The Hedra Zoo
I. Limonchenko, Moment-angle manifolds, 2-truncated cubes and Massey operations, arXiv:1510.07778 [math.AT], 2017.
M. Lin, Graph Cohomology, 2016, (Fig. 2.5 is a stellahedron).
T. Manneville and V. Pilaud, Compatibility fans for graphical nested complexes, arXiv:1501.07152v3 [math.CO], 2015-2016.
MathOverflow, Analogue of conic sections for the permutohedra, associahedra, and noncrossing partitions, an MO question posed by T. Copeland, 2017. (See Buchstaber references therein.)
V. Pilaud, The Associahedron and its Friends, presentation for Séminaire Lotharingien de Combinatoire, April 4-6, 2016. [From _Tom Copeland_, Jun 26 2018]

Cf. A046802, A007047, A122045, A000522, A036918, A001339, A052944.
Cf. A074909, A028246, A133314, A007318, A123125, A033282.
Cf. A008279, A008292, A090582, A097808, A130850.
Cf. A019538, A119879.

(* t = A046802 *) t[, 1] = 1; t[n, n_] = 1; t[n_, 2] = 2^(n - 1) - 1; t[n_, k_] = Sum[((i - k + 1)^i*(k - i)^(n - i - 1) - (i - k + 2)^i*(k - i - 1)^(n - i - 1))*Binomial[n - 1, i], {i, 0, k - 1}]; T[n_, j_] := Sum[Binomial[k, j]*t[n + 1, k + 1], {k, j, n}]; Table[T[n, k], {n, 0, 10}, {k, 0, n}] // Flatten (* Jean-François Alcover, Jan 23 2015, after Tom Copeland *)

Let M(n,k)= sum{i=0,..,k-1, C(n,i)[(i-k)^i*(k-i+1)^(n-i)- (i-k+1)^i*(k-i)^(n-i)]} with M(n,0)=1. Then M(n,k)= A046802(n,k), and T(n,j)= sum(k=j,..,n, C(k,j)*M(n,k)) for j>0 with T(n,0)= 1 + sum(k=1,..,n, M(n,k)) for n>0 and T(0,0)=1.
E.g.f: y * exp[x*(y+1)]/[y+1-exp(x*y)].
Row sums are A007047. Row polynomials evaluated at -1 are unity. Row polynomials evaluated at -2 are A122045.
First column is A000522. Second column appears to be A036918/2 = (A001339-1)/2 = n*A000522(n)/2.
Second diagonal is A052944. (Changed from conjecture to fact on Nov 08 2016.)
The raising operator for the reverse row polynomials with row signs is R = x - (1+t) - t e^(-D) / [1 + t(1-e^(-D))] evaluated at x = 0, with D = d/dx. Also R = x - d/dD log[exp(a.(0;t)D], or R = - d/dz log[e^(-xz) exp(a.(0;t)z)] = - d/dz log[exp(a.(-x;t)z)] with the e.g.f. defined in the comments and z replaced by D. Note that t e ^(-D) / [1+t(1-e^(-D))] = t - (t+t^2) D + (t+3t^2+2t^3) D^2/2! - ... is an e.g.f. for the signed reverse row polynomials of A028246. - Tom Copeland, Jan 23 2015
Equals A007318*(padded A090582)*A007318*A097808 = A007318*(padded (A008292*A007318))*A007318*A097808 = A007318*A130850 = A007318*(mirror of A028246). Padded means in the same way that A097805 is padded A007318. - Tom Copeland, Nov 14 2016
Umbrally, the row polynomials are p_n(x) = (1 + q.(x))^n, where (q.(x))^k = q_k(x) are the row polynomials of A130850. - Tom Copeland, Nov 16 2016
From the previous umbral statement, OP(x,d/dy) y^n =  (y + q.(x))^n, where OP(x,y) = exp[y * q.(x)] = x/((1+x)*exp(-x*y) - 1), the e.g.f. of A130850, so OP(x,d/dy) y^n evaluated at y = 1 is  p_n(x), the n-th row polynomial of this entry, with offset 0. - Tom Copeland, Jun 25 2018
Consolidating some formulas in this entry and A046082, in umbral notation for concision, with all offsets 0: Let A_n(x;y) = (y + E.(x))^n, an Appell sequence in y where E.(x)^k = E_k(x) are the Eulerian polynomials of A123125. Then the row polynomials of A046802 (the h-polynomials of the stellahedra) are given by h_n(x) = A_n(x;1); the row polynomials of this entry (A248727, the face polynomials of the stellahedra), by f_n(x) = A_n(1 + x;1); the Swiss-knife polynomials of A119879, by Sw_n(x) = A_n(-1;1 + x); and the row polynomials of the Worpitsky triangle (A130850), by w_n(x) = A(1 + x;0). Other specializations of A_n(x;y) give A090582 (the f-polynomials of the permutohedra, cf. also A019538) and A028246 (another version of the Worpitsky triangle). - Tom Copeland, Jan 24 2020

Title expanded by Tom Copeland, Nov 08 2016

cp's OEIS Frontend

A123125 Triangle of Eulerian numbers T(n,k), 0 <= k <= n, read by rows.

Views

Author

Keywords

Comments

Examples

References

Links

Crossrefs

Programs

Haskell

Maple

Mathematica

Python

Sage

Formula

A028246 Triangular array a(n,k) = (1/k)*Sum_{i=0..k} (-1)^(k-i)*binomial(k,i)*i^n; n >= 1, 1 <= k <= n, read by rows.

Views

Author

Keywords

Comments

Examples

References

Links

Crossrefs

Programs

GAP

Magma

Maple

Mathematica

PARI

PARI

Python

Sage

Sage

Formula

Extensions

A046802 T(n, k) = Sum_{j=k..n} binomial(n, j)*E1(j, j-k), where E1 are the Eulerian numbers A173018. Triangle read by rows, T(n, k) for 0 <= k <= n.

Views

Author

Keywords

Comments

Examples

References

Links

Crossrefs

Programs

Maple

Mathematica

Formula

Extensions

A090582 T(n, k) = Sum_{j=0..n-k} (-1)^j*binomial(n - k + 1, j)*(n - k + 1 - j)^n. Triangle read by rows, T(n, k) for 1 <= k <= n.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Haskell

Maple

Mathematica

PARI

Formula

A119879 Exponential Riordan array (sech(x),x).

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Maple

Mathematica

PARI

Sage

Sage

A028246 Triangular array a(n,k) = (1/k)Sum_{i=0..k} (-1)^(k-i)binomial(k,i)*i^n; n >= 1, 1 <= k <= n, read by rows.

A090582 T(n, k) = Sum_{j=0..n-k} (-1)^jbinomial(n - k + 1, j)(n - k + 1 - j)^n. Triangle read by rows, T(n, k) for 1 <= k <= n.