A003470 -id:A003470 - cp's OEIS frontend

A094587 Triangle of permutation coefficients arranged with 1's on the diagonal. Also, triangle of permutations on n letters with exactly k+1 cycles and with the first k+1 letters in separate cycles.

1, 1, 1, 2, 2, 1, 6, 6, 3, 1, 24, 24, 12, 4, 1, 120, 120, 60, 20, 5, 1, 720, 720, 360, 120, 30, 6, 1, 5040, 5040, 2520, 840, 210, 42, 7, 1, 40320, 40320, 20160, 6720, 1680, 336, 56, 8, 1, 362880, 362880, 181440, 60480, 15120, 3024, 504, 72, 9, 1, 3628800, 3628800

Offset: 0

Paul Barry, May 13 2004

Also, table of Pochhammer sequences read by antidiagonals (see Rudolph-Lilith, 2015). - N. J. A. Sloane, Mar 31 2016
Reverse of A008279. Row sums are A000522. Diagonal sums are A003470. Rows of inverse matrix begin {1}, {-1,1}, {0,-2,1}, {0,0,-3,1}, {0,0,0,-4,1} ... The signed lower triangular matrix (-1)^(n+k)n!/k! has as row sums the signed rencontres numbers Sum_{k=0..n} (-1)^(n+k)n!/k!. (See A000166). It has matrix inverse 1 1,1 0,2,1 0,0,3,1 0,0,0,4,1,...
Exponential Riordan array [1/(1-x),x]; column k has e.g.f. x^k/(1-x). - Paul Barry, Mar 27 2007
From Tom Copeland, Nov 01 2007: (Start)
T is the umbral extension of n!*Lag[n,(.)!*Lag[.,x,-1],0] = (1-D)^(-1) x^n = (-1)^n * n! * Lag(n,x,-1-n) = Sum_{j=0..n} binomial(n,j) * j! * x^(n-j) = Sum_{j=0..n} (n!/j!) x^j. The inverse operator is A132013 with generalizations discussed in A132014.
b = T*a can be characterized several ways in terms of a(n) and b(n) or their o.g.f.'s A(x) and B(x).
1) b(n) = n! Lag[n,(.)!*Lag[.,a(.),-1],0], umbrally,
2) b(n) = (-1)^n n! Lag(n,a(.),-1-n)
3) b(n) = Sum_{j=0..n} (n!/j!) a(j)
4) B(x) = (1-xDx)^(-1) A(x), formally
5) B(x) = Sum_{j=0,1,...} (xDx)^j A(x)
6) B(x) = Sum_{j=0,1,...} x^j * D^j * x^j A(x)
7) B(x) = Sum_{j=0,1,...} j! * x^j * L(j,-:xD:,0) A(x) where Lag(n,x,m) are the Laguerre polynomials of order m, D the derivative w.r.t. x and (:xD:)^j = x^j * D^j. Truncating the operator series at the j = n term gives an o.g.f. for b(0) through b(n).
c = (0!,1!,2!,3!,4!,...) is the sequence associated to T under the list partition transform and the associated operations described in A133314 so T(n,k) = binomial(n,k)*c(n-k). The reciprocal sequence is d = (1,-1,0,0,0,...). (End)
From Peter Bala, Jul 10 2008: (Start)
This array is the particular case P(1,1) of the generalized Pascal triangle P(a,b), a lower unit triangular matrix, shown below:
n\k|0.....................1...............2.......3......4
----------------------------------------------------------
0..|1.....................................................
1..|a....................1................................
2..|a(a+b)...............2a..............1................
3..|a(a+b)(a+2b).........3a(a+b).........3a........1......
4..|a(a+b)(a+2b)(a+3b)...4a(a+b)(a+2b)...6a(a+b)...4a....1
...
The entries A(n,k) of this array satisfy the recursion A(n,k) = (a+b*(n-k-1))*A(n-1,k) + A(n-1,k-1), which reduces to the Pascal formula when a = 1, b = 0.
Various cases are recorded in the database, including: P(1,0) = Pascal's triangle A007318, P(2,0) = A038207, P(3,0) = A027465, P(2,1) = A132159, P(1,3) = A136215 and P(2,3) = A136216.
When b <> 0 the array P(a,b) has e.g.f. exp(x*y)/(1-b*y)^(a/b) = 1 + (a+x)*y + (a*(a+b)+2a*x+x^2)*y^2/2! + (a*(a+b)*(a+2b) + 3a*(a+b)*x + 3a*x^2+x^3)*y^3/3! + ...; the array P(a,0) has e.g.f. exp((x+a)*y).
We have the matrix identities P(a,b)*P(a',b) = P(a+a',b); P(a,b)^-1 = P(-a,b).
An analog of the binomial expansion for the row entries of P(a,b) has been proved by [Echi]. Introduce a (generally noncommutative and nonassociative) product ** on the ring of polynomials in two variables by defining F(x,y)**G(x,y) = F(x,y)G(x,y) + by^2*d/dy(G(x,y)).
Define the iterated product F^(n)(x,y) of a polynomial F(x,y) by setting F^(1) = F(x,y) and F^(n)(x,y) = F(x,y)**F^(n-1)(x,y) for n >= 2. Then (x+a*y)^(n) = x^n + C(n,1)*a*x^(n-1)*y + C(n,2)*a*(a+b)*x^(n-2)*y^2 + ... + C(n,n)*a*(a+b)*(a+2b)*...*(a+(n-1)b)*y^n. (End)
(n+1) * n-th row = reversal of triangle A068424: (1; 2,2; 6,6,3; ...) - Gary W. Adamson, May 03 2009
Let G(m, k, p) = (-p)^k*Product_{j=0..k-1}(j - m - 1/p) and T(n,k,p) = G(n-1,n-k,p) then T(n, k, 1) is this sequence, T(n, k, 2) = A112292(n, k) and T(n, k, 3) = A136214. - Peter Luschny, Jun 01 2009, revised Jun 18 2019
The higher order exponential integrals E(x,m,n) are defined in A163931. For a discussion of the asymptotic expansions of the E(x,m=1,n) ~ (exp(-x)/x)*(1 - n/x + (n^2+n)/x^2 - (2*n+3*n^2+n^3)/x^3 + (6*n+11*n^2+6*n^3+n^4)/x^3 - ...) see A130534. The asymptotic expansion of E(x,m=1,n) leads for n >= 1 to the left hand columns of the triangle given above. Triangle A165674 is generated by the asymptotic expansions of E(x,m=2,n). - Johannes W. Meijer, Oct 07 2009
T(n,k) = n!/k! = number of permutations of [n+1] with exactly k+1 cycles and with elements 1,2,...,k+1 in separate cycles. See link and example below. - Dennis P. Walsh, Jan 24 2011
T(n,k) is the number of n permutations that leave some size k subset of {1,2,...,n} fixed. Sum_{k=0..n}(-1)^k*T(n,k) = A000166(n) (the derangements). - Geoffrey Critzer, Dec 11 2011
T(n,k) = A162995(n-1,k-1), 2 <= k <= n; T(n,k) = A173333(n,k), 1 <= k <= n. - Reinhard Zumkeller, Jul 05 2012
The row polynomials form an Appell sequence. The matrix is a special case of a group of general matrices sketched in A132382. - Tom Copeland, Dec 03 2013
For interpretations in terms of colored necklaces, see A213936 and A173333. - Tom Copeland, Aug 18 2016
See A008279 for a relation of this entry to the e.g.f.s enumerating the faces of permutahedra and stellahedra. - Tom Copeland, Nov 14 2016
Also, T(n,k) is the number of ways to arrange n-k nonattacking rooks on the n X (n-k) chessboard. - Andrey Zabolotskiy, Dec 16 2016
The infinitesimal generator of this triangle is the generalized exponential Riordan array [-log(1-x), x] and equals the unsigned version of A238363. - Peter Bala, Feb 13 2017
Formulas for exponential and power series infinitesimal generators for this triangle T are given in Copeland's 2012 and 2014 formulas as T = unsigned exp[(I-A238385)] =  1/(I - A132440), where I is the identity matrix. - Tom Copeland, Jul 03 2017
If A(0) = 1/(1-x), and A(n) = d/dx(A(n-1)), then A(n) = n!/(1-x)^(n+1) = Sum_{k>=0} (n+k)!/k!*x^k = Sum_{k>=0} T(n+k, k)*x^k. - Michael Somos, Sep 19 2021

			Rows begin {1}, {1,1}, {2,2,1}, {6,6,3,1}, ...
For n=3 and k=1, T(3,1)=6 since there are exactly 6 permutations of {1,2,3,4} with exactly 2 cycles and with 1 and 2 in separate cycles. The permutations are (1)(2 3 4), (1)(2 4 3), (1 3)(2 4), (1 4)(2 3), (1 3 4)(2), and (1 4 3)(2). - _Dennis P. Walsh_, Jan 24 2011
Triangle begins:
     1,
     1,    1,
     2,    2,    1,
     6,    6,    3,    1,
    24,   24,   12,    4,    1,
   120,  120,   60,   20,    5,    1,
   720,  720,  360,  120,   30,    6,    1,
  5040, 5040, 2520,  840,  210,   42,    7,    1
The production matrix is:
      1,     1,
      1,     1,     1,
      2,     2,     1,    1,
      6,     6,     3,    1,    1,
     24,    24,    12,    4,    1,   1,
    120,   120,    60,   20,    5,   1,   1,
    720,   720,   360,  120,   30,   6,   1,   1,
   5040,  5040,  2520,  840,  210,  42,   7,   1,   1,
  40320, 40320, 20160, 6720, 1680, 336,  56,   8,   1,   1
which is the exponential Riordan array A094587, or [1/(1-x),x], with an extra superdiagonal of 1's.
Inverse begins:
   1,
  -1,  1,
   0, -2,  1,
   0,  0, -3,  1,
   0,  0,  0, -4,  1,
   0,  0,  0,  0, -5,  1,
   0,  0,  0,  0,  0, -6,  1,
   0,  0,  0,  0,  0,  0, -7,  1

Reinhard Zumkeller, Rows n = 0..149 of triangle, flattened
J. Fernando Barbero G., Jesús Salas, and Eduardo J. S. Villaseñor, Bivariate Generating Functions for a Class of Linear Recurrences. I. General Structure, arXiv:1307.2010 [math.CO], 2013.
Paul Barry, The Restricted Toda Chain, Exponential Riordan Arrays, and Hankel Transforms, J. Int. Seq. 13 (2010) # 10.8.4, example 3.
Paul Barry, Exponential Riordan Arrays and Permutation Enumeration, J. Int. Seq. 13 (2010) # 10.9.1, example 5.
Paul Barry, Riordan Arrays, Orthogonal Polynomials as Moments, and Hankel Transforms, J. Int. Seq. 14 (2011) # 11.2.2, example 17.
Paul Barry, Combinatorial polynomials as moments, Hankel transforms and exponential Riordan arrays, arXiv:1105.3044 [math.CO], 2011, also J. Int. Seq. 14 (2011) 11.6.7.
Paul Barry, A note on number triangles that are almost their own production matrix, arXiv:1804.06801 [math.CO], 2018.
Paul Barry, On the inversion of Riordan arrays, arXiv:2101.06713 [math.CO], 2021.
Tom Copeland, Goin' with the Flow: Logarithm of the Derivative Operator Part V, 2014.
T. Copeland, Compositional inverse operators and Sheffer sequences, 2016.
E. Deutsch, L. Ferrari and S. Rinaldi, Production Matrices, Advances in Mathematics, 34 (2005) pp. 101-122.
Othman Echi, Binomial coefficients and Nasir al-Din al-Tusi, Scientific Research and Essays Vol.1 (2), 28-32 November 2006.
H. W. Gould, ed. J. Quaintance, Combinatorial Identities, May 2010 (eqn. 10.35, p.49).
A. Hennessy and P. Barry, Generalized Stirling Numbers, Exponential Riordan Arrays, and Orthogonal Polynomials, J. Int. Seq. 14 (2011) # 11.8.2.
Milan Janjic, Some classes of numbers and derivatives, JIS 12 (2009) 09.8.3.
Peter Luschny, Variants of Variations.
Michelle Rudolph-Lilith, On the Product Representation of Number Sequences, with Application to the Fibonacci Family, arXiv preprint arXiv:1508.07894 [math.NT], 2015.
M. Z. Spivey, On Solutions to a General Combinatorial Recurrence, J. Int. Seq. 14 (2011) # 11.9.7.
Dennis Walsh, A note on permutations with cyclic constraints
Wikipedia, Sheffer sequence

Cf. A000166 (alt. row sums), A000522 (row sums).
Cf. A068424, A036039, A173333, A213936, A263916.
Cf. A000670, A008279, A021009, A048994, A132013, A132014, A132159, A132440, A133314, A218234, A235706, A238385.

a094587 n k = a094587_tabl !! n !! k
a094587_row n = a094587_tabl !! n
a094587_tabl = map fst $ iterate f ([1], 1)
   where f (row, i) = (map (* i) row ++ [1], i + 1)
-- Reinhard Zumkeller, Jul 04 2012

T := proc(n, m): n!/m! end: seq(seq(T(n, m), m=0..n), n=0..9);  # Johannes W. Meijer, Oct 07 2009, revised Nov 25 2012
# Alternative: Note that if you leave out 'abs' you get A021009.
T := proc(n, k) option remember; if n = 0 and k = 0 then 1 elif k < 0 or k > n then 0 else abs((n + k)*T(n-1, k) - T(n-1, k-1)) fi end: #  Peter Luschny, Dec 30 2021

Flatten[Table[Table[n!/k!, {k,0,n}], {n,0,10}]] (* Geoffrey Critzer, Dec 11 2011 *)

def A094587_row(n): return (factorial(n)*exp(x).taylor(x,0,n)).list()
for n in (0..7): print(A094587_row(n)) # Peter Luschny, Sep 28 2017

T(n, k) = n!/k! if n >= k >= 0, otherwise 0.
T(n, k) = Sum_{i=k..n} |S1(n+1, i+1)*S2(i, k)| * (-1)^i, with S1, S2 the Stirling numbers.
T(n,k) = (n-k)*T(n-1,k) + T(n-1,k-1). E.g.f.: exp(x*y)/(1-y) = 1 + (1+x)*y + (2+2*x+x^2)*y^2/2! + (6+6*x+3*x^2+x^3)*y^3/3!+ ... . - Peter Bala, Jul 10 2008
A094587 = 1 / ((-1)*A129184 * A127648 + I), I = Identity matrix. - Gary W. Adamson, May 03 2009
From Johannes W. Meijer, Oct 07 2009: (Start)
The o.g.f. of right hand column k is Gf(z;k) = (k-1)!/(1-z)^k, k => 1.
The recurrence relations of the right hand columns lead to Pascal's triangle A007318. (End)
Let f(x) = (1/x)*exp(-x). The n-th row polynomial is R(n,x) = (-x)^n/f(x)*(d/dx)^n(f(x)), and satisfies the recurrence equation R(n+1,x) = (x+n+1)*R(n,x)-x*R'(n,x). Cf. A132159. - Peter Bala, Oct 28 2011
A padded shifted version of this lower triangular matrix with zeros in the first column and row except for a one in the diagonal position is given by integral(t=0 to t=infinity) exp[-t(I-P)] = 1/(I-P) = I + P^2 + P^3 + ... where P is the infinitesimal generator matrix A218234 and I the identity matrix. The non-padded version is given by P replaced by A132440. - Tom Copeland, Oct 25 2012
From Peter Bala, Aug 28 2013: (Start)
The row polynomials R(n,x) form a Sheffer sequence of polynomials with associated delta operator equal to d/dx. Thus d/dx(R(n,x)) = n*R(n-1,x). The Sheffer identity is R(n,x + y) = Sum_{k=0..n} binomial(n,k)*y^(n-k)*R(k,x).
Let P(n,x) = Product_{k=0..n-1} (x + k) denote the rising factorial polynomial sequence with the convention that P(0,x) = 1. Then this is triangle of connection constants when expressing the basis polynomials P(n,x + 1) in terms of the basis P(n,x). For example, row 3 is (6, 6, 3, 1) so P(3,x + 1) = (x + 1)*(x + 2)*(x + 3) = 6 + 6*x + 3*x*(x + 1) + x*(x + 1)*(x + 2). (End)
From Tom Copeland, Apr 21 & 26, and Aug 13 2014: (Start)
T-I = M = -A021009*A132440*A021009 with e.g.f. y*exp(x*y)/(1-y). Cf. A132440. Dividing the n-th row of M by n generates the (n-1)th row of T.
T = 1/(I - A132440) = {2*I - exp[(A238385-I)]}^(-1) = unsigned exp[(I-A238385)] = exp[A000670(.)*(A238385-I)] = , umbrally, where I = identity matrix.
The e.g.f. is exp(x*y)/(1-y), so the row polynomials form an Appell sequence with lowering operator d/dx and raising operator x + 1/(1-D).
With L(n,m,x)= Laguerre polynomials of order m, the row polynomials are (-1)^n*n!*L(n,-1-n,x) = (-1)^n*(-1!/(-1-n)!)*K(-n,-1-n+1,x) = n!* K(-n,-n,x) where K is Kummer's confluent hypergeometric function (as a limit of n+s as s tends to zero).
Operationally, (-1)^n*n!*L(n,-1-n,-:xD:) = (-1)^n*x^(n+1)*:Dx:^n*x^(-1-n) = (-1)^n*x*:xD:^n*x^(-1) = (-1)^n*n!*binomial(xD-1,n) = n!*K(-n,-n,-:xD:) where :AB:^n = A^n*B^n for any two operators. Cf. A235706 and A132159.
The n-th row of signed M has the coefficients of d[(-:xD:)^n]/d(:Dx:)= f[d/d(-:xD:)](-:xD:)^n with f(y)=y/(y-1), :Dx:^n= n!L(n,0,-:xD:), and (-:xD:)^n = n!L(n,0,:Dx:). M has the coefficients of [D/(1-D)]x^n. (End)
From Tom Copeland, Nov 18 2015: (Start)
Coefficients of the row polynomials of the e.g.f. Sum_{n>=0} P_n(b1,b2,..,bn;t) x^n/n! = e^(P.(..;t) x) = e^(xt) / (1-b.x) = (1 + b1 x + b2 x^2 + b3 x^3 + ...) e^(xt) = 1 + (b1 + t) x + (2 b2 + 2 b1 t + t^2) x^2/2! + (6 b3 + 6 b2 t + 3 b1 t^2 + t^3) x^3/3! + ... , with lowering operator L = d/dt, i.e., L P_n(..;t) = n * P_(n-1)(..;t), and raising operator R = t + d[log(1 + b1 D + b2 D^2 + ...)]/dD = t - Sum_{n>=1} F(n,b1,..,bn) D^(n-1), i.e., R P_n(..,;t) = P_(n+1)(..;t), where D = d/dt and F(n,b1,..,bn) are the Faber polynomials of A263916.
Also P_n(b1,..,bn;t) = CIP_n(t-F(1,b1),-F(2,b1,b2),..,-F(n,b1,..,bn)), the cycle index polynomials A036039.
(End)
The raising operator R = x + 1/(1-D) = x + 1 + D + D^2 + ... in matrix form acting on an o.g.f. (formal power series) is the transpose of the production matrix M below. The linear term x is the diagonal of ones after transposition. The other transposed diagonals come from D^m x^n = n! / (n-m)! x^(n-m). Then P(n,x) = (1,x,x^2,..) M^n (1,0,0,..)^T is a matrix representation of R P(n-1,x) = P(n,x). - Tom Copeland, Aug 17 2016
The row polynomials have e.g.f. e^(xt)/(1-t) = exp(t*q.(x)), umbrally. With p_n(x) the row polynomials of A132013, q_n(x) = v_n(p.(u.(x))), umbrally, where u_n(x) = (-1)^n v_n(-x) = (-1)^n Lah_n(x), the Lah polynomials with e.g.f. exp[x*t/(t-1)]. This has the matrix form [T] = [q] = [v]*[p]*[u]. Conversely, p_n(x) = u_n (q.(v.(x))). - Tom Copeland, Nov 10 2016
From the Appell sequence formalism, 1/(1-b.D) t^n = P_n(b1,b2,..,bn;t), the generalized row polynomials noted in the Nov 18 2015 formulas, consistent with the 2007 comments. - Tom Copeland, Nov 22 2016
From Peter Bala, Feb 18 2017: (Start)
G.f.: Sum_{n >= 1} (n*x)^(n-1)/(1 + (n - t)*x)^n = 1 + (1 + t)*x + (2 + 2*t + t^2)*x^2 + ....
n-th row polynomial R(n,t) = Sum_{k = 0..n} (-1)^(n-k)*binomial(n,k)*(x + k)^k*(x + k - t)^(n-k) = Sum_{k = 0..n} (-1)^(n-k)*binomial(n,k)*(x + k)^(n-k)*(x + k + t)^k, for arbitrary x. The particular case of the latter sum when x = 0 and t = 1 is identity 10.35 in Gould, Vol.4. (End)
Rodrigues-type formula for the row polynomials: R(n, x) = -exp(x)*Int(exp(-x)* x^n, x), for n >= 0. Recurrence: R(n, x) =  x^n + n*R(n-1, x), for  n >= 1, and R(0, x) = 1. d/dx(R(n, x)) = R(n, x) - x^n, for n >= 0 (compare with the formula from Peter Bala, Aug 28 2013). - Wolfdieter Lang, Dec 23 2019
T(n, k) = Sum_{i=0..n-k} A048994(n-k, i) * n^i for 0 <= k <= n. - Werner Schulte, Jul 26 2022

Edited by Johannes W. Meijer, Oct 07 2009

New description from Dennis P. Walsh, Jan 24 2011

A086764 Triangle T(n, k), read by row, related to Euler's difference table A068106 (divide column k of A068106 by k!).

Original entry on oeis.org

1, 0, 1, 1, 1, 1, 2, 3, 2, 1, 9, 11, 7, 3, 1, 44, 53, 32, 13, 4, 1, 265, 309, 181, 71, 21, 5, 1, 1854, 2119, 1214, 465, 134, 31, 6, 1, 14833, 16687, 9403, 3539, 1001, 227, 43, 7, 1, 133496, 148329, 82508, 30637, 8544, 1909, 356, 57, 8, 1

Offset: 0

Philippe Deléham, Aug 02 2003

The k-th column sequence, k >= 0, without leading zeros, enumerates the ways to distribute n beads, n >= 1, labeled differently from 1 to n, over a set of (unordered) necklaces, excluding necklaces with exactly one bead, and k+1 indistinguishable, ordered, fixed cords, each allowed to have any number of beads. Beadless necklaces as well as beadless cords each contribute a factor 1, hence for n=0 one has 1. See A000255 for the description of a fixed cord with beads. This comment derives from a family of recurrences found by Malin Sjodahl for a combinatorial problem for certain quark and gluon diagrams (Feb 27 2010). - Wolfdieter Lang, Jun 02 2010

			Formatted as a square array:
      1      3     7    13   21   31  43 57 ... A002061;
      2     11    32    71  134  227 356    ... A094792;
      9     53   181   465 1001 1909        ... A094793;
     44    309  1214  3539 8544             ... A094794;
    265   2119  9403 30637                  ... A023043;
   1854  16687 82508                        ... A023044;
  14833 148329                              ... A023045;
Formatted as a triangular array (mirror of A076731):
       1;
       0      1;
       1      1     1;
       2      3     2     1;
       9     11     7     3    1;
      44     53    32    13    4    1;
     265    309   181    71   21    5    1;
    1854   2119  1214   465  134   31    6   1;
   14833  16687  9403  3539 1001  227   43   7   1;
  133496 148329 82508 30637 8544 1909  356  57   8   1;

Indranil Ghosh, Rows 0..50, flattened
W. Y. C. Chen et al., Higher-order log-concavity in Euler's difference table, Discrete Math., 311 (2011), 2128-2134. (These are the numbers d^k_n.)
Fanja Rakotondrajao, k-Fixed-Points-Permutations, Integers: Electronic journal of combinatorial number theory 7 (2007) A36.

Columns: A000166, A000155, A000153, A000261, A001909, A001910, A176732, A176733, A176734, A176735, A176736.
Mirror image of A076731.
Cf. A002061, A003470, A008288, A023043, A023044, A023045, A068106.
Cf. A076731, A094792, A094793, A094794, A095000, A108625, A143007.
Cf. A143409, A143410, A143411, A143413.

A086764:= func< n,k | (&+[(-1)^j*Binomial(n-k,j)*Factorial(n-j): j in [0..n]])/Factorial(k) >;
[A086764(n,k): k in [0..n], n in [0..12]]; // G. C. Greubel, Oct 05 2023

T[n_,k_]:=(1/k!)*Sum[(-1)^j*Binomial[n-k,j]*(n-j)!,{j,0,n}];Flatten[Table[T[n,k],{n,0,11},{k,0,n}]] (* Indranil Ghosh, Feb 20 2017 *)
T[n_, k_] := (n!/k!) HypergeometricPFQ[{k-n},{-n},-1];
Table[T[n,k], {n,0,9}, {k,0,n}] // Flatten (* Peter Luschny, Oct 05 2017 *)

def A086764(n,k): return sum((-1)^j*binomial(n-k,j)*factorial(n-j) for j in range(n+1))//factorial(k)
flatten([[A086764(n,k) for k in range(n+1)] for n in range(13)]) # G. C. Greubel, Oct 05 2023

T(n, n) = 1; T(n+1, n) = n.
T(n+2, n) = A002061(n+1) = n^2 + n + 1; T(n+3, n) = n^3 + 3*n^2 + 5*n + 2.
T(n, k) = (k + 1)*T(n, k + 1) - T(n-1, k); T(n, n) = 1; T(n, k) = 0, if k > n.
T(n, k) = (n-1)*T(n-1, k) + (n-k-1)*T(n-2, k).
k!*T(n, k) = A068106(n, k). [corrected by Georg Fischer, Aug 13 2022]
Sum_{k>=0} T(n, k) = A003470(n+1).
T(n, k) = (1/k!) * Sum_{j>=0} (-1)^j*binomial(n-k, j)*(n-j)!. - Philippe Deléham, Jun 13 2005
From Peter Bala, Aug 14 2008: (Start)
The following remarks all relate to the array read as a square array: e.g.f for column k: exp(-y)/(1-y)^(k+1); e.g.f. for array: exp(-y)/(1-x-y) = (1 + x + x^2 + x^3 + ...) + (x + 2*x^2 + 3*x^3 + 4*x^4 + ...)*y + (1 + 3*x + 7*x^2 + 13*x^3 + ...)*y^2/2! + ... .
This table is closely connected to the constant e. The row, column and diagonal entries of this table occur in series formulas for e.
Row n for n >= 2: e = n!*(1/T(n,0) + (-1)^n*[1/(1!*T(n,0)*T(n,1)) + 1/(2!*T(n,1)*T(n,2)) + 1/(3!*T(n,2)*T(n,3)) + ...]). For example, row 3 gives e = 6*(1/2 - 1/(1!*2*11) - 1/(2!*11*32) - 1/(3!*32*71) - ...). See A095000.
Column 0: e = 2 + Sum_{n>=2} (-1)^n*n!/(T(n,0)*T(n+1,0)) = 2 + 2!/(1*2) - 3 !/(2*9) + 4!/(9*44) - ... .
Column k, k >= 1: e = (1 + 1/1! + 1/2! + ... + 1/k!) + 1/k!*Sum_{n >= 0} (-1)^n*n!/(T(n,k)*T(n+1,k)). For example, column 3 gives e = 8/3 + 1/6*(1/(1*3) - 1/(3*13) + 2/(13*71) - 6/(71*465) + ...).
Main diagonal: e = 1 + 2*(1/(1*1) - 1/(1*7) + 1/(7*71) - 1/(71*1001) + ...).
First subdiagonal: e = 8/3 + 5/(3*32) - 7/(32*465) + 9/(465*8544) - ... .
Second subdiagonal: e = 2*(1 + 2^2/(1*11) - 3^2/(11*181) + 4^2/(181*3539) - ...). See A143413.
Third subdiagonal: e = 3 - (2*3*5)/(2*53) + (3*4*7)/(53*1214) - (4*5*9)/(1214*30637) + ... .
For the corresponding results for the constants 1/e, sqrt(e) and 1/sqrt(e) see A143409, A143410 and A143411 respectively. For other arrays similarly related to constants see A008288 (for log(2)), A108625 (for zeta(2)) and A143007 (for zeta(3)). (End)
G.f. for column k is hypergeom([1,k+1],[],x/(x+1))/(x+1). - Mark van Hoeij, Nov 07 2011
T(n, k) = (n!/k!)*hypergeom([k-n], [-n], -1). - Peter Luschny, Oct 05 2017

More terms from David Wasserman, Mar 28 2005
Additional comments from Zerinvary Lajos, Mar 30 2006
Edited by N. J. A. Sloane, Sep 24 2011

A143409 Square array read by antidiagonals: form the Euler-Seidel matrix for the sequence {k!} and then divide column k by k!.

Original entry on oeis.org

1, 2, 1, 5, 3, 1, 16, 11, 4, 1, 65, 49, 19, 5, 1, 326, 261, 106, 29, 6, 1, 1957, 1631, 685, 193, 41, 7, 1, 13700, 11743, 5056, 1457, 316, 55, 8, 1, 109601, 95901, 42079, 12341, 2721, 481, 71, 9, 1, 986410, 876809, 390454, 116125, 25946, 4645, 694, 89, 10, 1

Offset: 0

Peter Bala, Aug 14 2008

The Euler-Seidel matrix for the sequence {k!} is array A076571 read as a square, whose k-th column entries have a common factor of k!. Removing these common factors gives the current table.
This table is closely connected to the constant 1/e. The row, column and diagonal entries of this table occur in series acceleration formulas for 1/e.
For a similar table based on the differences of the sequence {k!} and related to the constant e, see A086764. For other arrays similarly related to constants see A143410 (for sqrt(e)), A143411 (for 1/sqrt(e)), A008288 (for log(2)), A108625 (for zeta(2)) and A143007 (for zeta(3)).

			The Euler-Seidel matrix for the sequence {k!} begins
==============================================
n\k|.....0.....1.....2.....3.....4.....5.....6
==============================================
0..|.....1.....1.....2.....6....24...120...720
1..|.....2.....3.....8....30...144...840
2..|.....5....11....38...174...984
3..|....16....49...212..1158
4..|....65...261..1370
5..|...326..1631
6..|..1957
...
Dividing the k-th column by k! gives
==============================================
n\k|.....0.....1.....2.....3.....4.....5.....6
==============================================
0..|.....1.....1.....1.....1.....1.....1.....1
1..|.....2.....3.....4.....5.....6.....7
2..|.....5....11....19....29....41
3..|....16....49...106...193
4..|....65...261...685
5..|...326..1631
6..|..1957
...
Examples of series formula for 1/e:
Row 2: 1/e = 2*(1/5 - 1/(1!*5*11) + 1/(2!*11*19) - 1/(3!*19*29) + ...).
Column 4: 24/e = 9 - (0!/(1*6) + 1!/(6*41) + 2!/(41*316) + ...).
...
Displayed as a triangle:
0 |     1
1 |     2,     1
2 |     5,     3,    1
3 |    16,    11,    4,    1
4 |    65,    49,   19,    5,   1
5 |   326,   261,  106,   29,   6,  1
6 |  1957,  1631,  685,  193,  41,  7, 1
7 | 13700, 11743, 5056, 1457, 316, 55, 8, 1

Robert Israel, Table of n, a(n) for n = 0..10010 (antidiagonals 0 to 140, flattened)
D. Dumont, Matrices d'Euler-Seidel, Sem. Loth. Comb. B05c (1981) 59-78.
Eric Weisstein's World of Mathematics Poisson-Charlier polynomial

Cf. A008288, A076571, A086764, A108625, A143007, A143410, A143411, A143413, A001517 (main diagonal), A028387 (row 2), A000522 (column 0), A001339 (column 1), A082030 (column 2), A095000 (column 3), A095177 (column 4).

Cf. A273596, A003470.

T := (n, k) -> 1/k!*add(binomial(n,j)*(k+j)!, j = 0..n):
for n from 0 to 9 do seq(T(n, k), k = 0..9) end do;
# Alternate:
T:= proc(n,k) option remember;
  if n = 0 then return 1 fi;
  (n+k)*procname(n-1,k) + procname(n-1,k-1);
end proc:
seq(seq(T(s-n,n),n=0..s),s=0..10); # Robert Israel, Jul 07 2017
# Or:
A143409 := (n,k) -> hypergeom([k+1, k-n], [], -1):
seq(seq(simplify(A143409(n,k)),k=0..n),n=0..9); # Peter Luschny, Oct 05 2017

T[n_, k_] := HypergeometricPFQ[{k+1,k-n}, {}, -1];
Table[T[n,k], {n,0,9}, {k,0,n}] // Flatten (* Peter Luschny, Oct 05 2017 *)

T(n,k) = (1/k!)*Sum_{j = 0..n} binomial(n,j)*(k+j)!.
T(n,k) = ((n+k)!/k!)*Num_Pade(n,k), where Num_Pade(n,k) denotes the numerator of the Padé approximation for the function exp(x) of degree (n,k) evaluated at x = 1.
Recurrence relations:
T(n,k) = T(n-1,k) + (k+1)*T(n-1,k+1);
T(n,k) = (n+k)*T(n-1,k) + T(n-1,k-1).
E.g.f. for column k: exp(y)/(1-y)^(k+1).
E.g.f. for array: exp(y)/(1-x-y) = (1 + x + x^2 + ...) + (2 + 3*x + 4*x^2 + ...)*y + (5 + 11*x + 19*x^2 + ...)*y^2/2! + ... .
Row n lists the values of the Poisson-Charlier polynomial x^(n) + C(n,1)*x^(n-1) + C(n,2)*x^(n-2) + ... + C(n,n) for x = 1,2,3,..., where x^(m) denotes the rising factorial x*(x+1)*...*(x+m-1).
Main diagonal is A001517.
Series formulas for 1/e:
Row n: 1/e = n!*[1/T(n,0) - 1/(1!*T(n,0)*T(n,1)) + 1/(2!*T(n,1)*T(n,2)) - 1/(3!*T(n,2)*T(n,3)) + ...].
Column k: k!/e = A000166(k) + (-1)^(k+1)*[0!/(T(0,k)*T(1,k)) + 1!/(T(1,k)*T(2,k)) + 2!/(T(2,k)*T(3,k)) + ...].
Main diagonal: 1/e = 1 - 2*Sum_{n>=0} (-1)^n/(T(n,n)*T(n+1,n+1)) = 1 - 2*[1/(1*3) - 1/(3*19) + 1/(19*193) - ...].
Second subdiagonal: 1/e = 2*(1^2/(1*5) - 2^2/(5*49) + 3^2/(49*685) - ...).
Compare with A143413.
From Peter Luschny, Oct 05 2017: (Start)
T(n, k) = hypergeom([k+1, k-n], [], -1).
When seen as a triangular array then the row sums are A273596 and the alternating row sums are A003470. (End)

A060475 Triangular array formed from successive differences of factorial numbers, then with factorials removed.

Original entry on oeis.org

1, 1, 0, 1, 1, 1, 1, 2, 3, 2, 1, 3, 7, 11, 9, 1, 4, 13, 32, 53, 44, 1, 5, 21, 71, 181, 309, 265, 1, 6, 31, 134, 465, 1214, 2119, 1854, 1, 7, 43, 227, 1001, 3539, 9403, 16687, 14833, 1, 8, 57, 356, 1909, 8544, 30637, 82508, 148329, 133496, 1, 9, 73, 527, 3333, 18089, 81901, 296967, 808393, 1468457, 1334961

Offset: 0

Henry Bottomley, Mar 16 2001

T(n,k) is also the number of partial bijections (of an n-element set) with a fixed domain of size k and without fixed points. Equivalently, T(n,k) is the number of partial derangements with a fixed domain of size k in the symmetric inverse semigroup (monoid), I sub n. - Abdullahi Umar, Sep 14 2008

			Triangle begins
  1,
  1,  0,
  1,  1,  1,
  1,  2,  3,  2,
  1,  3,  7, 11,  9,
  1,  4, 13, 32, 53, 44,
  ...

G. C. Greubel, Rows n=0..100 of triangle, flattened
A. Laradji, and A. Umar, Combinatorial results for the symmetric inverse semigroup, Semigroup Forum 75 (2007), 221-236. - _Abdullahi Umar_, Sep 14 2008
L. Takacs, The Problem of Coincidences, Archive for History of Exact Sciences, Volume 21, No. 3, Sept. 1980. pp 229-244, paragraph 10 (Catalan).

Columns include A000012, A001477, A002061.
Diagonals include A000166, A000255, A000153, A000261, A001909, A001910.
Main diagonal is abs of A002119.
Similar to A076731.
Row sums equal A003470. - Johannes W. Meijer, Jul 27 2011
Cf. A047920, A008290.

[[Factorial(k)*(&+[(-1)^j*Binomial(n-j, k-j)/Factorial(j): j in [0..k]]): k in [0..n]]: n in [0..10]]; // G. C. Greubel, Mar 04 2019

A060475 := proc(n,k): k! * add(binomial(n-j,k-j)*(-1)^j/j!, j=0..k) end:
seq(seq(A060475(n,k), k=0..n), n=0..7); # Johannes W. Meijer, Jul 27 2011
T := (n,k) -> KummerU(-k, -n, -1):
seq(seq(simplify(T(n, k)), k = 0..n), n = 0..10); # Peter Luschny, Jul 07 2022

t[n_, k_] := k!*Sum[Binomial[n - j, k - j]*(-1)^j/j!, {j, 0, k}]; Table[t[n, k], {n, 0, 10}, {k, 0, n}] // Flatten (* Robert G. Wilson v, Aug 08 2011 *)

{T(n,k) = k!*sum(j=0,k, (-1)^j*binomial(n-j, k-j)/j!)};
for(n=0,10, for(k=0,n, print1(T(n,k), ", "))) \\ G. C. Greubel, Mar 04 2019

[[factorial(k)*sum((-1)^j*binomial(n-j, k-j)/factorial(j) for j in (0..k)) for k in (0..n)] for n in (0..10)] # G. C. Greubel, Mar 04 2019

T(n,k) = A047920(n,k)/(n-k)! = (n-1)*T(n-1,k-1) + (k-1)*T(n-2,k-2) = (n-k+1)*T(n, k-1) - T(n-1,k-1).
From Abdullahi Umar, Sep 14 2008: (Start)
T(n,k) = k! * Sum_{j=0..k} C(n-j,k-j)*(-1)^j/j!.
C(n,k)*T(n,k) = A144089(n, k). (End)
T(n,k) = A076732(n+1,k+1)/(k+1). - Johannes W. Meijer, Jul 27 2011
E.g.f. as a square array: A(x,y) = exp(-x)/(1 - x - y) = (1 + y + y^2 + y^3 + ...) + (y + 2*y^2 + 3*y^3 + 4*y^4 + ...)*x + (1 + 3*y + 7*y^2 + 13*y^3 + ...)*x^2/2! + (2 + 11*y + 32*y^2 + 71*y^3 + ...)*x^3/3! + .... Observe that (1 - y)*A(x*(1 - y),y) = exp(x*(y - 1))/(1 - x) is the e.g.f. for A008290. - Peter Bala, Sep 25 2013
T(n, k) = KummerU(-k, -n, -1). - Peter Luschny, Jul 07 2022

cp's OEIS Frontend

A259859 a(0)=0; thereafter A003470(n-1) + A003470(n) - 1.

Views

Author

Keywords

Links

Crossrefs

Programs

Maple

Mathematica

Formula

Extensions

A094587 Triangle of permutation coefficients arranged with 1's on the diagonal. Also, triangle of permutations on n letters with exactly k+1 cycles and with the first k+1 letters in separate cycles.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Haskell

Maple

Mathematica

Sage

Formula

Extensions

A086764 Triangle T(n, k), read by row, related to Euler's difference table A068106 (divide column k of A068106 by k!).

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Magma

Mathematica

SageMath

Formula

Extensions

A143409 Square array read by antidiagonals: form the Euler-Seidel matrix for the sequence {k!} and then divide column k by k!.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Maple

Mathematica

Formula

A060475 Triangular array formed from successive differences of factorial numbers, then with factorials removed.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Magma

Maple

Mathematica

PARI

Sage

Formula

A076732 Table T(n,k) giving number of ways of obtaining exactly one correct answer on an (n,k)-matching problem (1 <= k <= n).

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Maple

Mathematica

Formula