cp's OEIS Frontend

This is a front-end for the Online Encyclopedia of Integer Sequences, made by Christian Perfect. The idea is to provide OEIS entries in non-ancient HTML, and then to think about how they're presented visually. The source code is on GitHub.

Showing 1-6 of 6 results.

A256890 Triangle T(n,k) = t(n-k, k); t(n,m) = f(m)*t(n-1,m) + f(n)*t(n,m-1), where f(x) = x + 2.

Original entry on oeis.org

1, 2, 2, 4, 12, 4, 8, 52, 52, 8, 16, 196, 416, 196, 16, 32, 684, 2644, 2644, 684, 32, 64, 2276, 14680, 26440, 14680, 2276, 64, 128, 7340, 74652, 220280, 220280, 74652, 7340, 128, 256, 23172, 357328, 1623964, 2643360, 1623964, 357328, 23172, 256, 512, 72076, 1637860, 10978444, 27227908, 27227908, 10978444, 1637860, 72076, 512
Offset: 0

Views

Author

Dale Gerdemann, Apr 12 2015

Keywords

Comments

Related triangles may be found by varying the function f(x). If f(x) is a linear function, it can be parameterized as f(x) = a*x + b. With different values for a and b, the following triangles are obtained:
a\b 1.......2.......3.......4.......5.......6
The row sums of these, and similarly constructed number triangles, are shown in the following table:
a\b 1.......2.......3.......4.......5.......6.......7.......8.......9
The formula can be further generalized to: t(n,m) = f(m+s)*t(n-1,m) + f(n-s)*t(n,m-1), where f(x) = a*x + b. The following table specifies triangles with nonzero values for s (given after the slash).
a\b 0 1 2 3
-2 A130595/1
-1
0
With the absolute value, f(x) = |x|, one obtains A038221/3, A038234/4, A038247/5, A038260/6, A038273/7, A038286/8, A038299/9 (with value for s after the slash).
If f(x) = A000045(x) (Fibonacci) and s = 1, the result is A010048 (Fibonomial).
In the notation of Carlitz and Scoville, this is the triangle of generalized Eulerian numbers A(r, s | alpha, beta) with alpha = beta = 2. Also the array A(2,1,4) in the notation of Hwang et al. (see page 31). - Peter Bala, Dec 27 2019

Examples

			Array, t(n, k), begins as:
   1,    2,      4,        8,        16,         32,          64, ...;
   2,   12,     52,      196,       684,       2276,        7340, ...;
   4,   52,    416,     2644,     14680,      74652,      357328, ...;
   8,  196,   2644,    26440,    220280,    1623964,    10978444, ...;
  16,  684,  14680,   220280,   2643360,   27227908,   251195000, ...;
  32, 2276,  74652,  1623964,  27227908,  381190712,  4677894984, ...;
  64, 7340, 357328, 10978444, 251195000, 4677894984, 74846319744, ...;
Triangle, T(n, k), begins as:
    1;
    2,     2;
    4,    12,      4;
    8,    52,     52,       8;
   16,   196,    416,     196,      16;
   32,   684,   2644,    2644,     684,      32;
   64,  2276,  14680,   26440,   14680,    2276,     64;
  128,  7340,  74652,  220280,  220280,   74652,   7340,   128;
  256, 23172, 357328, 1623964, 2643360, 1623964, 357328, 23172,   256;
		

Crossrefs

Programs

  • Magma
    A256890:= func< n,k | (&+[(-1)^(k-j)*Binomial(j+3,j)*Binomial(n+4,k-j)*(j+2)^n: j in [0..k]]) >;
    [A256890(n,k): k in [0..n], n in [0..10]]; // G. C. Greubel, Oct 18 2022
    
  • Mathematica
    Table[Sum[(-1)^(k-j)*Binomial[j+3, j] Binomial[n+4, k-j] (j+2)^n, {j,0,k}], {n,0, 9}, {k,0,n}]//Flatten (* Michael De Vlieger, Dec 27 2019 *)
  • PARI
    t(n,m) = if ((n<0) || (m<0), 0, if ((n==0) && (m==0), 1, (m+2)*t(n-1, m) + (n+2)*t(n, m-1)));
    tabl(nn) = {for (n=0, nn, for (k=0, n, print1(t(n-k, k), ", ");); print(););} \\ Michel Marcus, Apr 14 2015
    
  • SageMath
    def A256890(n,k): return sum((-1)^(k-j)*Binomial(j+3,j)*Binomial(n+4,k-j)*(j+2)^n for j in range(k+1))
    flatten([[A256890(n,k) for k in range(n+1)] for n in range(11)]) # G. C. Greubel, Oct 18 2022

Formula

T(n,k) = t(n-k, k); t(0,0) = 1, t(n,m) = 0 if n < 0 or m < 0 else t(n,m) = f(m)*t(n-1,m) + f(n)*t(n,m-1), where f(x) = x + 2.
Sum_{k=0..n} T(n, k) = A001715(n).
T(n,k) = Sum_{j = 0..k} (-1)^(k-j)*binomial(j+3,j)*binomial(n+4,k-j)*(j+2)^n. - Peter Bala, Dec 27 2019
Modified rule of Pascal: T(0,0) = 1, T(n,k) = 0 if k < 0 or k > n else T(n,k) = f(n-k) * T(n-1,k-1) + f(k) * T(n-1,k), where f(x) = x + 2. - Georg Fischer, Nov 11 2021
From G. C. Greubel, Oct 18 2022: (Start)
T(n, n-k) = T(n, k).
T(n, 0) = A000079(n). (End)

A143495 Triangle read by rows: 3-Stirling numbers of the second kind.

Original entry on oeis.org

1, 3, 1, 9, 7, 1, 27, 37, 12, 1, 81, 175, 97, 18, 1, 243, 781, 660, 205, 25, 1, 729, 3367, 4081, 1890, 380, 33, 1, 2187, 14197, 23772, 15421, 4550, 644, 42, 1, 6561, 58975, 133057, 116298, 47271, 9702, 1022, 52, 1, 19683, 242461, 724260, 830845, 447195
Offset: 3

Views

Author

Peter Bala, Aug 20 2008

Keywords

Comments

This is the case r = 3 of the r-Stirling numbers of the second kind. The 3-Stirling numbers of the second kind give the number of ways of partitioning the set {1,2,...,n} into k nonempty disjoint subsets with the restriction that the elements 1, 2 and 3 belong to distinct subsets. For remarks on the general case see A143494 (r = 2). The corresponding array of 3-Stirling numbers of the first kind is A143492. The theory of r-Stirling numbers of both kinds is developed in [Broder]. For 3-Lah numbers refer to A143498.
From Peter Bala, Sep 19 2008: (Start)
Let D be the derivative operator d/dx and E the Euler operator x*d/dx. Then x^(-3)*E^n*x^3 = Sum_{k = 0..n} T(n+3,k+3)*x^k*D^k.
The row generating polynomials R_n(x) := Sum_{k= 3..n} T(n,k)*x^k satisfy the recurrence R_(n+1)(x) = x*R_n(x) + x*d/dx(R_n(x)) with R_3(x) = x^3. It follows that the polynomials R_n(x) have only real zeros (apply Corollary 1.2. of [Liu and Wang]).
Relation with the 3-Eulerian numbers E_3(n,j) := A144697(n,j): T(n,k) = (3!/k!)*Sum_{j = n-k..n-3} E_3(n,j)*binomial(j,n-k) for n >= k >= 3.
(End)
T(n,k) = S(n,k,3), n>=k>=3, in Mikhailov's first paper, eq.(28) or (A3). E.g.f. column k from (A20) with k->3, r->k. Therefore, with offset [0,0], this triangle is the Sheffer triangle (exp(3*x),exp(x)-1) with e.g.f. of column no. m>=0: exp(3*x)*((exp(x)-1)^m)/m!. See one of the formulas given below. For Sheffer matrices see the W. Lang link under A006232 with the S. Roman reference, also found in A132393. - Wolfdieter Lang, Sep 29 2011

Examples

			Triangle begins
  n\k|....3....4....5....6....7....8
  ==================================
  3..|....1
  4..|....3....1
  5..|....9....7....1
  6..|...27...37...12....1
  7..|...81..175...97...18....1
  8..|..243..781..660..205...25....1
  ...
T(5,4) = 7. The set {1,2,3,4,5} can be partitioned into four subsets such that 1, 2 and 3 belong to different subsets in 7 ways: {{1}{2}{3}{4,5}}, {{1}{2}{5}{3,4}}, {{1}{2}{4}{3,5}}, {{1}{3}{4}{2,5}}, {{1}{3}{5}{2,4}}, {{2}{3}{4}{1,5}} and {{2}{3}{5}{1,4}}.
From _Peter Bala_, Feb 23 2025: (Start)
The array factorizes as
/ 1               \       /1             \ /1             \ /1            \
| 3    1           |     | 3   1          ||0  1           ||0  1          |
| 9    7   1       |  =  | 9   4   1      ||0  3   1       ||0  0  1       | ...
|27   37  12   1   |     |27  13   5  1   ||0  9   4  1    ||0  0  3  1    |
|81  175  97  18  1|     |81  40  18  6  1||0 27  13  5  1 ||0  0  9  4  1 |
|...               |     |...             ||...            ||...           |
where, in the infinite product on the right-hand side, the first array is the Riordan array (1/(1 - 3*x), x/(1 - x)). See A106516. (End)
		

Crossrefs

Cf. A005061 (column 4), A005494 (row sums), A008277, A016753 (column 5), A028025 (column 6), A049458 (matrix inverse), A106516, A143492, A143494, A143496, A143498.

Programs

  • Maple
    A143495 := (n, k) -> (1/(k-3)!)*add((-1)^(k-i-1)*binomial(k-3,i)*(i+3)^(n-3), i = 0..k-3): for n from 3 to 12 do seq(A143495(n, k), k = 3..n) end do;
  • Mathematica
    nmax = 12; t[n_, k_] := 1/(k-3)!* Sum[ (-1)^(k-j-1)*Binomial[k-3, j]*(j+3)^(n-3), {j, 0, k-3}]; Flatten[ Table[ t[n, k], {n, 3, nmax}, {k, 3, n}]] (* Jean-François Alcover, Dec 07 2011, after Maple *)
  • Sage
    @CachedFunction
    def stirling2r(n, k, r) :
        if n < r: return 0
        if n == r: return 1 if k == r else 0
        return stirling2r(n-1, k-1, r) + k*stirling2r(n-1, k, r)
    A143495 = lambda n, k: stirling2r(n, k, 3)
    for n in (3..8): [A143495(n, k) for k in (3..n)] # Peter Luschny, Nov 19 2012

Formula

T(n+3,k+3) = (1/k!)*Sum_{i = 0..k} (-1)^(k-i)*binomial(k,i)*(i+3)^n, n,k >= 0.
T(n,k) = Stirling2(n,k) - 3*Stirling2(n-1,k) + 2*Stirling2(n-2,k), n,k >= 3.
Recurrence relation: T(n,k) = T(n-1,k-1) + k*T(n-1,k) for n > 3, with boundary conditions: T(n,2) = T(2,n) = 0 for all n; T(3,3) = 1; T(3,k) = 0 for k > 3.
Special cases: T(n,3) = 3^(n-3); T(n,4) = 4^(n-3) - 3^(n-3).
E.g.f. (k+3) column (with offset 3): (1/k!)*exp(3x)*(exp(x)-1)^k.
O.g.f. k-th column: Sum_{n >= k} T(n,k)*x^n = x^k/((1-3*x)*(1-4*x)*...*(1-k*x)).
E.g.f.: exp(3*t + x*(exp(t)-1)) = Sum_{n >= 0} Sum_{k = 0..n} T(n+3,k+3)*x^k*t^n/n! = Sum_{n >= 0} B_n(3;x)*t^n/n! = 1 + (3+x)*t/1! + (9+7*x+x^2)*t^2/2! + ..., where the row polynomials, B_n(3;x) := Sum_{k = 0..n} T(n+3,k+3)*x^k, may be called the 3-Bell polynomials.
Dobinski-type identities: Row polynomial B_n(3;x) = exp(-x)*Sum_{i >= 0} (i+3)^n*x^i/i!; Sum_{k = 0..n} k!*T(n+3,k+3)*x^k = Sum_{i >= 0} (i+3)^n*x^i/(1+x)^(i+1).
The T(n,k) are the connection coefficients between the falling factorials and the shifted monomials (x+3)^(n-3). For example, 9 + 7*x + x*(x-1) = (x+3)^2 and 27 + 37*x + 12x*(x-1) + x*(x-1)*(x-2) = (x+3)^3.
This array is the matrix product P^2 * S, where P denotes Pascal's triangle, A007318 and S denotes the lower triangular array of Stirling numbers of the second kind, A008277 (apply Theorem 10 of [Neuwirth]). The inverse array is A049458, the signed 3-Stirling numbers of the first kind.

A144696 Triangle of 2-Eulerian numbers.

Original entry on oeis.org

1, 1, 2, 1, 7, 4, 1, 18, 33, 8, 1, 41, 171, 131, 16, 1, 88, 718, 1208, 473, 32, 1, 183, 2682, 8422, 7197, 1611, 64, 1, 374, 9327, 49780, 78095, 38454, 5281, 128, 1, 757, 30973, 264409, 689155, 621199, 190783, 16867, 256
Offset: 2

Views

Author

Peter Bala, Sep 19 2008

Keywords

Comments

Let [n] denote the ordered set {1,2,...,n}. The symmetric group S_n consists of the injective mappings p:[n] -> [n]. Such a permutation p has an excedance at position i, 1 <= i < n, if p(i) > i. One well-known interpretation of the Eulerian numbers A(n,k) is that they count the permutations in the symmetric group S_n with k excedances. The triangle of Eulerian numbers is A008292 (but with an offset of 1 in the column numbering). We generalize this definition to restricted permutations as follows.
Let r be a nonnegative integer and let Permute(n,n-r) denote the set of injective maps p:[n-r] -> [n], which we think of as permutations of n numbers taken n-r at a time. Clearly, |Permute(n,n-r)| = n!/r!. We say that p has an excedance at position i, 1 <= i <= n-r, if p(i) > i. Then the r-Eulerian number, denoted by A(r;n,k), is defined as the number of permutations in Permute(n,n-r) having k excedances. Thus the current array of 2-Eulerian numbers gives the number of permutations in Permute(n,n-2) with k excedances. See the example section below for some numerical examples.
Clearly A(0;n,k) = A(n,k). The case r = 1 also produces the ordinary Eulerian numbers A(n,k). There is an obvious bijection from Permute(n,n) to Permute(n,n-1) that preserves the number of excedances of a permutation. Consequently, the 1-Eulerian numbers are equal to the 0-Eulerian numbers: A(1;n,k) = A(0;n,k) = A(n,k).
For other cases of r-Eulerian numbers see A144697 (r = 3), A144698 (r = 4) and A144699 (r = 5). There is also a concept of r-Stirling numbers of the first and second kinds - see A143491 and A143494. If we multiply the entries of the current array by a factor of 2 and then reverse the rows we obtain A120434.
An alternative interpretation of the current array due to [Strosser] involves the 2-excedance statistic of a permutation (see also [Foata & Schutzenberger, Chapitre 4, Section 3]). We define a permutation p in Permute(n,n-2) to have a 2-excedance at position i (1 <=i <= n-2) if p(i) >= i + 2.
Given a permutation p in Permute(n,n-2), define ~p to be the permutation in Permute(n,n-2) that takes i to n+1 - p(n-i-1). The map ~ is a bijection of Permute(n,n-2) with the property that if p has (resp. does not have) an excedance in position i then ~p does not have (resp. has) a 2-excedance at position n-i-1. Hence ~ gives a bijection between the set of permutations with k excedances and the set of permutations with (n-k) 2-excedances. Thus reading the rows of this array in reverse order gives a triangle whose entries count the permutations in Permute(n,n-2) with k 2-excedances.
Example: Represent a permutation p:[n-2] -> [n] in Permute(n,n-2) by its image vector (p(1),...,p(n-2)). In Permute(10,8) the permutation p = (1,2,4,7,10,6,5,8) does not have an excedance in the first two positions (i = 1 and 2) or in the final three positions (i = 6, 7 and 8). The permutation ~p = (3,6,5,1,4,7,9,10) has 2-excedances only in the first three positions and the final two positions.
From Peter Bala, Dec 27 2019: (Start)
This is the array A(1,1,3) in the notation of Hwang et al. (p. 25), where the authors remark that the r-Eulerian numbers were first studied by Shanlan Li (Duoji Bilei, Ch. 4), who gave the summation formulas
Sum_{i = 2..n+1} (i-1)*C(i,2) = C(n+3,4) + 2*C(n+2,4)
Sum_{i = 2..n+1} (i-1)^2*C(i,2) = C(n+4,5) + 7*C(n+3,5) + 4*C(n+2,5)
Sum_{i = 2..n+1} (i-1)^3*C(i,2) = C(n+5,6) + 18*C(n+4,6) + 33*C(n+3,6) + 8*C(n+2,6). (End)

Examples

			The triangle begins
===========================================
n\k|..0.....1.....2.....3.....4.....5.....6
===========================================
2..|..1
3..|..1.....2
4..|..1.....7.....4
5..|..1....18....33.....8
6..|..1....41...171...131....16
7..|..1....88...718..1208...473....32
8..|..1...183..2682..8422..7197..1611....64
...
Row 4 = [1,7,4]: We represent a permutation p:[n-2] -> [n] in Permute(n,n-2) by its image vector (p(1),...,p(n-2)). Here n = 4. The permutation (1,2) has no excedances; 7 permutations have a single excedance, namely, (1,3), (1,4), (2,1), (3,1), (3,2), (4,1) and (4,2); the remaining 4 permutations, (2,3), (2,4), (3,4) and (4,3) each have two excedances.
		

References

  • J. Riordan. An introduction to combinatorial analysis. New York, J. Wiley, 1958.
  • R. Strosser. Séminaire de théorie combinatoire, I.R.M.A., Université de Strasbourg, 1969-1970.
  • Li, Shanlan (1867). Duoji bilei (Series summation by analogies), 4 scrolls. In Zeguxizhai suanxue (Mathematics from the Studio Devoted to the Imitation of the Ancient Chinese Tradition) (Jinling ed.), Volume 4.
  • Li, Shanlan (2019). Catégories analogues d’accumulations discrètes (Duoji bilei), traduit et commenté par Andrea Bréard. La Bibliothèque Chinoise. Paris: Les Belles Lettres.

Crossrefs

Cf. A000079 (right diagonal), A001710 (row sums).
Cf. A000182 (related to alt. row sums).

Programs

  • Magma
    m:=2; [(&+[(-1)^(k-j)*Binomial(n+1,k-j)*Binomial(j+m,m-1)*(j+1)^(n-m+1): j in [0..k]])/m: k in [0..n-m], n in [m..m+10]]; // G. C. Greubel, Jun 04 2022
    
  • Maple
    with(combinat):
    T:= (n,k) -> 1/2!*add((-1)^(k-j)*binomial(n+1,k-j)*(j+1)^(n-1)*(j+2), j = 0..k):
    for n from 2 to 10 do
    seq(T(n,k),k = 0..n-2)
    end do;
  • Mathematica
    T[n_, k_]:= 1/2!*Sum[(-1)^(k-j)*Binomial[n+1, k-j]*(j+1)^(n-1)*(j+2), {j, 0, k}];
    Table[T[n, k], {n,2,10}, {k,0,n-2}]//Flatten (* Jean-François Alcover, Oct 15 2019 *)
  • SageMath
    m=2 # A144696
    def T(n,k): return (1/m)*sum( (-1)^(k-j)*binomial(n+1,k-j)*binomial(j+m,m-1)*(j+1)^(n-m+1) for j in (0..k) )
    flatten([[T(n,k) for k in (0..n-m)] for n in (m..m+10)]) # G. C. Greubel, Jun 04 2022

Formula

T(n,k) = (1/2!)*Sum_{j = 0..k} (-1)^(k-j)*binomial(n+1,k-j)*(j+1)^(n-1)*(j+2);
T(n,n-k) = (1/2!)*Sum_{j = 2..k} (-1)^(k-j)*binomial(n+1,k-j)*j^(n-1)*(j-1).
Recurrence relations:
T(n,k) = (k+1)*T(n-1,k) + (n-k)*T(n-1,k-1) with boundary conditions T(n,0) = 1 for n >= 2, T(2,k) = 0 for k >= 1. Special cases: T(n,n-2) = 2^(n-2); T(n,n-3) = A066810(n-1).
E.g.f. (with suitable offsets): (1/2)*[(1 - x)/(1 - x*exp(t - t*x))]^2 = 1/2 + x*t + (x + 2*x^2)*t^2/2! + (x + 7*x^2 + 4*x^3)*t^3/3! + ... .
The row generating polynomials R_n(x) satisfy the recurrence R_(n+1)(x) = (n*x+1)*R_n(x) + x*(1-x)*d/dx(R_n(x)) with R_2(x) = 1. It follows that the polynomials R_n(x) for n >= 3 have only real zeros (apply Corollary 1.2. of [Liu and Wang]).
The (n+1)-th row generating polynomial = (1/2!)*Sum_{k = 1..n} (k+1)!*Stirling2(n,k) *x^(k-1)*(1-x)^(n-k).
For n >= 2,
(1/2)*(x*d/dx)^(n-1) (1/(1-x)^2) = x/(1-x)^(n+1) * Sum_{k = 0..n-2} T(n,k)*x^k,
(1/2)*(x*d/dx)^(n-1) (x^2/(1-x)^2) = 1/(1-x)^(n+1) * Sum_{k = 2..n} T(n,n-k)*x^k,
1/(1-x)^(n+1)*Sum_{k = 0..n-2} T(n,k)*x^k = (1/2!) * Sum_{m = 0..inf} (m+1)^(n-1)*(m+2)*x^m,
1/(1-x)^(n+1)*Sum_{k = 2..n} T(n,n-k)*x^k = (1/2!) * Sum_{m = 2..inf} m^(n-1)*(m-1)*x^m.
Worpitzky-type identities:
Sum_{k = 0..n-2} T(n,k)*binomial(x+k,n) = (1/2!)*x^(n-1)*(x - 1);
Sum_{k = 2..n} T(n,n-k)*binomial(x+k,n) = (1/2!)*(x + 1)^(n-1)*(x + 2).
Relation with Stirling numbers (Frobenius-type identities):
T(n+1,k-1) = (1/2!) * Sum_{j = 0..k} (-1)^(k-j)*(j+1)!* binomial(n-j,k-j)*Stirling2(n,j) for n,k >= 1;
T(n+1,k-1) = (1/2!) * Sum_{j = 0..n-k} (-1)^(n-k-j)*(j+1)!* binomial(n-j,k)*S(2;n+2,j+2) for n,k >= 1 and
T(n+2,k) = (1/2!) * Sum_{j = 0..n-k} (-1)^(n-k-j)*(j+2)!* binomial(n-j,k)*S(2;n+2,j+2) for n,k >= 0, where S(2;n,k) denotes the 2-Stirling numbers A143494(n,k).
The row polynomials of this array are related to the Eulerian polynomials. For example, 1/2*x*d/dx [x*(x + 4*x^2 + x^3)/(1-x)^4] = x^2*(1 + 7*x + 4*x^2)/(1-x)^5 and 1/2*x*d/dx [x*(x + 11*x^2 + 11*x^3 + x^4)/(1-x)^5] = x^2*(1 + 18*x + 33*x^2 + 8*x^3)/(1-x)^6.
Row sums A001710. Alternating row sums [1, -1, -2, 8, 16, -136, -272, 3968, 7936, ... ] are alternately (signed) tangent numbers and half tangent numbers - see A000182.
Sum_{k = 0..n-2} 2^k*T(n,k) = A069321(n-1). Sum_{k = 0..n-2} 2^(n-k)*T(n,k) = 4*A083410(n-1).
For n >=2, the shifted row polynomial t*R(n,t) = (1/2)*D^(n-1)(f(x,t)) evaluated at x = 0, where D is the operator (1-t)*(1+x)*d/dx and f(x,t) = (1+x*t/(t-1))^(-2). - Peter Bala, Apr 22 2012

A144699 Triangle of 5-Eulerian numbers.

Original entry on oeis.org

1, 1, 5, 1, 16, 25, 1, 39, 171, 125, 1, 86, 786, 1526, 625, 1, 181, 3046, 11606, 12281, 3125, 1, 372, 10767, 70792, 142647, 92436, 15625, 1, 755, 36021, 380071, 1279571, 1553145, 663991, 78125, 1, 1522, 116368, 1880494, 9818494, 19555438, 15519952, 4608946, 390625
Offset: 5

Views

Author

Peter Bala, Sep 19 2008

Keywords

Comments

This is the case r = 5 of the r-Eulerian numbers, denoted by A(r;n,k), defined as follows. Let [n] denote the ordered set {1,2,...,n} and let r be a nonnegative integer. Let Permute(n,n-r) denote the set of injective maps p:[n-r] -> [n], which we think of as permutations of n numbers taken n-r at a time. Clearly, |Permute(n,n-r)| = n!/r!. We say that the permutation p has an excedance at position i, 1 <= i <= n-r, if p(i) > i. Then the r-Eulerian number A(r;n,k) is defined as the number of permutations in Permute(n,n-r) with k excedances. Thus the 5-Eulerian numbers count the permutations in Permute(n,n-5) with k excedances. For other cases see A008292 (r = 0 and r = 1), A144696 (r = 2), A144697 (r = 3) and A144698 (r = 4).
An alternative interpretation of the current array due to [Strosser] involves the 5-excedance statistic of a permutation (see also [Foata & Schutzenberger, Chapitre 4, Section 3]). We define a permutation p in Permute(n,n-5) to have a 5-excedance at position i (1 <=i <= n-5) if p(i) >= i + 5.
Given a permutation p in Permute(n,n-5), define ~p to be the permutation in Permute(n,n-5) that takes i to n+1 - p(n-i-4). The map ~ is a bijection of Permute(n,n-5) with the property that if p has (resp. does not have) an excedance in position i then ~p does not have (resp. has) a 5-excedance at position n-i-4. Hence ~ gives a bijection between the set of permutations with k excedances and the set of permutations with (n-k) 5-excedances. Thus reading the rows of this array in reverse order gives a triangle whose entries count the permutations in Permute(n,n-5) with k 5-excedances.
Example: Represent a permutation p:[n-5] -> [n] in Permute(n,n-5) by its image vector (p(1),...,p(n-5)). In Permute(12,7) the permutation p = (1,2,4,12,3,6,5) does not have an excedance in the first two positions (i = 1 and 2) or in the final three positions (i = 5, 6 and 7). The permutation ~p = (8,7,10,1,9,11,12) has 5-excedances only in the first three positions and the final two positions.

Examples

			Triangle begins
=======================================================
n\k|..0.......1......2......3.........4.......5.......6
=======================================================
5..|..1
6..|..1.......5
7..|..1......16......25
8..|..1......39.....171.....125
9..|..1......86.....786....1526.....625
10.|..1.....181....3046...11606...12281....3125
11.|..1.....372...10767...70792..142647...92436...15625
...
T(7,1) = 16: We represent a permutation p:[n-5] -> [n] in Permute(n,n-5) by its image vector (p(1),...,p(n-5)). The 16 permutations in Permute(7,2) having 1 excedance are (1,3), (1,4), (1,5), (1,6), (1,7), (3,2), (4,2), (5,2), (6,2), (7,2), (2,1), (3,1), (4,1), (5,1), (6,1) and (7,1).
		

References

  • R. Strosser, Séminaire de théorie combinatoire, I.R.M.A., Université de Strasbourg, 1969-1970.

Crossrefs

Programs

  • Magma
    m:=5; [(&+[(-1)^(k-j)*Binomial(n+1,k-j)*Binomial(j+m,m-1)*(j+1)^(n-m+1): j in [0..k]])/m: k in [0..n-m], n in [m..m+10]]; // G. C. Greubel, Jun 08 2022
    
  • Maple
    with(combinat):
    T:= (n,k) -> 1/5!*add((-1)^(k-j)*binomial(n+1,k-j)*(j+1)^(n-4)*(j+2)*(j+3)*(j+4)*(j+5),j = 0..k):
    for n from 5 to 13 do
    seq(T(n,k),k = 0..n-5)
    end do;
  • Mathematica
    T[n_, k_] /; 0 < k <= n-5 := T[n, k] = (k+1) T[n-1, k] + (n-k) T[n-1, k-1];
    T[, 0] = 1; T[, _] = 0;
    Table[T[n, k], {n, 5, 13}, {k, 0, n-5}] // Flatten (* Jean-François Alcover, Nov 11 2019 *)
  • SageMath
    m=5 # A144699
    def T(n,k): return (1/m)*sum( (-1)^(k-j)*binomial(n+1,k-j)*binomial(j+m,m-1)*(j+1)^(n-m+1) for j in (0..k) )
    flatten([[T(n,k) for k in (0..n-m)] for n in (m..m+10)]) # G. C. Greubel, Jun 08 2022

Formula

T(n,k) = (1/5!)*Sum_{j = 0..k} (-1)^(k-j)*binomial(n+1,k-j)*(j+1)^(n-4)*(j+2)*(j+3)*(j+4)*(j+5).
T(n,n-k) = (1/5!)*Sum_{j = 5..k} (-1)^(k-j)*binomial(n+1,k-j)*j^(n-4)*(j-1)*(j-2)*(j-3)*(j-4).
Recurrence relation:
T(n,k) = (k+1)*T(n-1,k) + (n-k)*T(n-1,k-1) with boundary conditions T(n,0) = 1 for n >= 5, T(5,k) = 0 for k >= 1.
E.g.f. (with suitable offsets): (1/5)*((1 - x)/(1 - x*exp(t - t*x)))^5 = 1/5 + x*t + (x + 5*x^2)*t^2/2! + (x + 16*x^2 + 25*x^3)*t^3/3! + ... .
The row generating polynomials R_n(x) satisfy the recurrence R_(n+1)(x) = (n*x+1)*R_n(x) + x*(1-x)*d/dx(R_n(x)) with R_5(x) = 1. It follows that the polynomials R_n(x) for n >= 6 have only real zeros (apply Corollary 1.2. of [Liu and Wang]).
The (n+4)-th row generating polynomial = (1/5!)*Sum_{k = 1..n} (k+4)!*Stirling2(n,k)*x^(k-1)*(1-x)^(n-k).
For n >= 5,
(1/5)*(x*d/dx)^(n-4) (1/(1-x)^5) = (x/(1-x)^(n+1)) * Sum_{k = 0..n-5} T(n,k)*x^k,
(1/5)*(x*d/dx)^(n-4) (x^5/(1-x)^5) = (1/(1-x)^(n+1)) * Sum_{k = 5..n} T(n,n-k)*x^k,
(1/(1-x)^(n+1)) * Sum_{k = 0..n-5} T(n,k)*x^k = (1/5!) * Sum_{m = 0..infinity} (m+1)^(n-4)*(m+2)*(m+3)*(m+4)*(m+5)*x^m,
(1/(1-x)^(n+1)) * Sum_{k = 5..n} T(n,n-k)*x^k = (1/5!) * Sum_{m = 5..infinity} m^(n-4)*(m-1)*(m-2)*(m-3)*(m-4)*x^m.
Worpitzky-type identities:
Sum_{k = 0..n-5} T(n,k)*binomial(x+k,n) = (1/5!)*x^(n-4)*(x-1)*(x-2)*(x-3)*(x-4).
Sum_{k = 5..n} T(n,n-k)* binomial(x+k,n) = (1/5!)*(x+1)^(n-4)*(x+2)*(x+3)*(x+4)*(x+5).
Relation with Stirling numbers (Frobenius-type identities):
T(n+4,k-1) = (1/5!) * Sum_{j = 0..k} (-1)^(k-j)* (j+4)!* binomial(n-j,k-j)*Stirling2(n,j) for n,k >= 1;
T(n+4,k-1) = (1/5!) * Sum_{j = 0..n-k} (-1)^(n-k-j)*(j+4)!* binomial(n-j,k)*S(5;n+5,j+5) for n,k >= 0 and
T(n+5,k) = (1/5!) * Sum_{j = 0..n-k} (-1)^(n-k-j)*(j+5)!* binomial(n-j,k)*S(5;n+5,j+5) for n,k >= 0, where S(5;n,k) denotes the 5-Stirling numbers of the second kind, which are given by the formula S(5;n+5,j+5) = 1/j!*Sum_{i = 0..j} (-1)^(j-i)*binomial(j,i)*(i+5)^n for n,j >= 0.

A144698 Triangle of 4-Eulerian numbers.

Original entry on oeis.org

1, 1, 4, 1, 13, 16, 1, 32, 113, 64, 1, 71, 531, 821, 256, 1, 150, 2090, 6470, 5385, 1024, 1, 309, 7470, 40510, 65745, 33069, 4096, 1, 628, 25191, 221800, 612295, 592884, 194017, 16384, 1, 1267, 81853, 1113919, 4835875, 7843369, 4915423, 1101157, 65536
Offset: 4

Views

Author

Peter Bala, Sep 19 2008

Keywords

Comments

This is the case r = 4 of the r-Eulerian numbers, denoted by A(r;n,k), defined as follows. Let [n] denote the ordered set {1,2,...,n} and let r be a nonnegative integer. Let Permute(n,n-r) denote the set of injective maps p:[n-r] -> [n], which we think of as permutations of n numbers taken n-r at a time. Clearly, |Permute(n,n-r)| = n!/r!. We say that the permutation p has an excedance at position i, 1 <= i <= n-r, if p(i) > i. Then the r-Eulerian number A(r;n,k) is defined as the number of permutations in Permute(n,n-r) with k excedances. Thus the 4-Eulerian numbers are the number of permutations in Permute(n,n-4) with k excedances. For other cases see A008292 (r = 0 and r = 1), A144696 (r = 2), A144697 (r = 3) and A144699 (r = 5).
An alternative interpretation of the current array due to [Strosser] involves the 4-excedance statistic of a permutation (see also [Foata & Schutzenberger, Chapter 4, Section 3]). We define a permutation p in Permute(n,n-4) to have a 4-excedance at position i (1 <= i <= n-4) if p(i) >= i + 4.
Given a permutation p in Permute(n,n-4), define ~p to be the permutation in Permute(n,n-4) that takes i to n+1 - p(n-i-3). The map ~ is a bijection of Permute(n,n-4) with the property that if p has (resp. does not have) an excedance in position i then ~p does not have (resp. has) a 4-excedance at position n-i-3. Hence ~ gives a bijection between the set of permutations with k excedances and the set of permutations with (n-k) 4-excedances. Thus reading the rows of this array in reverse order gives a triangle whose entries are the number of permutations in Permute(n,n-4) with k 4-excedances.
Example: Represent a permutation p:[n-4] -> [n] in Permute(n,n-4) by its image vector (p(1),...,p(n-4)). In Permute(10,6) the permutation p = (1,2,4,10,3,6) does not have an excedance in the first two positions (i = 1 and 2) or in the final two positions (i = 5 and 6). The permutation ~p = (5,8,1,7,9,10) has 4-excedances only in the first two positions and the final two positions.

Examples

			Triangle begins
  ===+=============================================
  n\k|  0      1      2      3      4      5      6
  ===+=============================================
   4 |  1
   5 |  1      4
   6 |  1     13     16
   7 |  1     32    113     64
   8 |  1     71    531    821    256
   9 |  1    150   2090   6470   5385   1024
  10 |  1    309   7470  40510  65745  33069   4096
  ...
T(6,1) = 13: We represent a permutation p:[n-4] -> [n] in Permute(n,n-4) by its image vector (p(1),...,p(n-4)). The 13 permutations in Permute(6,2) having 1 excedance are (1,3), (1,4), (1,5), (1,6), (3,2), (4,2), (5,2), (6,2), (2,1), (3,1), (4,1), (5,1) and (6,1).
		

References

  • R. Strosser, Séminaire de théorie combinatoire, I.R.M.A., Université de Strasbourg, 1969-1970.

Crossrefs

Cf. A001720 (row sums), A000302 (right diagonal).

Programs

  • Magma
    m:=4; [(&+[(-1)^(k-j)*Binomial(n+1,k-j)*Binomial(j+m,m-1)*(j+1)^(n-m+1): j in [0..k]])/m: k in [0..n-m], n in [m..m+10]]; // G. C. Greubel, Jun 04 2022
    
  • Maple
    with(combinat):
    T:= (n,k) -> 1/4!*add((-1)^(k-j)*binomial(n+1,k-j)*(j+1)^(n-3)*(j+2)*(j+3)*(j+4),j = 0..k):
    for n from 4 to 12 do
    seq(T(n,k),k = 0..n-4)
    end do;
  • Mathematica
    T[n_, k_] /; 0 < k <= n-4 := T[n, k] = (k+1) T[n-1, k] + (n-k) T[n-1, k-1];
    T[, 0] = 1; T[, _] = 0;
    Table[T[n, k], {n, 4, 12}, {k, 0, n-4}] // Flatten (* Jean-François Alcover, Nov 11 2019 *)
  • SageMath
    m=4 # A144698
    def T(n,k): return (1/m)*sum( (-1)^(k-j)*binomial(n+1,k-j)*binomial(j+m,m-1)*(j+1)^(n-m+1) for j in (0..k) )
    flatten([[T(n,k) for k in (0..n-m)] for n in (m..m+10)]) # G. C. Greubel, Jun 04 2022

Formula

T(n,k) = (1/4!)*Sum_{j = 0..k} (-1)^(k-j)*binomial(n+1,k-j)*(j+1)^(n-3)*(j+2)*(j+3)*(j+4).
T(n,n-k) = (1/4!)*Sum_{j = 4..k} (-1)^(k-j)*binomial(n+1,k-j)*j^(n-3)*(j-1)*(j-2)*(j-3).
Recurrence relation:
T(n,k) = (k + 1)*T(n-1,k) + (n-k)*T(n-1,k-1) with boundary conditions T(n,0) = 1 for n >= 4, T(4,k) = 0 for k >= 1. Special cases: T(n,n-4) = 4^(n-4); T(n,n-5) = 5^(n-3) - 4^(n-3) - (n-3)*4^(n-4).
E.g.f. (with suitable offsets): 1/4*[(1 - x)/(1 - x*exp(t - t*x))]^4 = 1/4 + x*t + (x + 4*x^2)*t^2/2! + (x + 13*x^2 + 16*x^3)*t^3/3! + ... .
The row generating polynomials R_n(x) satisfy the recurrence R_(n+1)(x) = (n*x + 1)*R_n(x) + x*(1 - x)*d/dx(R_n(x)) with R_4(x) = 1. It follows that the polynomials R_n(x) for n >= 5 have only real zeros (apply Corollary 1.2. of [Liu and Wang]).
The (n+3)-th row generating polynomial = (1/4!)*Sum_{k = 1..n} (k+3)!*Stirling2(n,k)*x^(k-1)*(1-x)^(n-k).
For n >= 4,
1/4*(x*d/dx)^(n-3) (1/(1-x)^4) = x/(1-x)^(n+1) * Sum_{k = 0..n-4} T(n,k)*x^k,
1/4*(x*d/dx)^(n-3) (x^4/(1-x)^4) = 1/(1-x)^(n+1) * Sum_{k = 4..n} T(n,n-k)*x^k,
1/(1-x)^(n+1) * Sum {k = 0..n-4} T(n,k)*x^k = (1/4!) * Sum_{m = 0..inf} (m+1)^(n-3)*(m+2)*(m+3)*(m+4)*x^m,
1/(1-x)^(n+1) * Sum {k = 4..n} T(n,n-k)*x^k = (1/4!) * Sum_{m = 4..inf} m^(n-3)*(m-1)*(m-2)*(m-3)*x^m,
Worpitzky-type identities:
Sum_{k = 0..n-4} T(n,k)*binomial(x+k,n) = (1/4!)*x^(n-3)*(x-1)*(x-2)*(x-3).
Sum_{k = 4..n} T(n,n-k)* binomial(x+k,n) = (1/4!)*(x+1)^(n-3)*(x+2)*(x+3)*(x+4).
Relation with Stirling numbers (Frobenius-type identities):
T(n+3,k-1) = (1/4!) * Sum_{j = 0..k} (-1)^(k-j)* (j+3)!* binomial(n-j,k-j)*Stirling2(n,j) for n,k >= 1;
T(n+3,k-1) = 1/4! * Sum_{j = 0..n-k} (-1)^(n-k-j)*(j+3)!* binomial(n-j,k)*S(4;n+4,j+4) for n,k >= 1 and
T(n+4,k) = 1/4! * Sum_{j = 0..n-k} (-1)^(n-k-j)*(j+4)!* binomial(n-j,k)*S(4;n+4,j+4) for n,k >= 0, where S(4;n,k) denotes the 4-Stirling numbers of the second kind A143496(n,k).
For n >=4, the shifted row polynomial t*R(n,t) = (1/4)*D^(n-3)(f(x,t)) evaluated at x = 0, where D is the operator (1-t)*(1+x)*d/dx and f(x,t) = (1+x*t/(t-1))^(-4). - Peter Bala, Apr 22 2012

A152249 Triangle of 4 - restricted Eulerian numbers as polynomials used in exponential data smoothing: m(p,k,x)=((-1)^k*(1 - x)^(p + k)/(k!(p - 1)!))*Sum[(p - 1 + j)!*j^k*x^j/(j!), {j, 0, Infinity}]/x;n=6; t(m,l)=coefficients((-1)^m*m!*M[n, m, x])/n.

Original entry on oeis.org

1, 1, 6, 1, 19, 36, 1, 46, 241, 216, 1, 101, 1091, 2551, 1296, 1, 212, 4182, 18932, 24337, 7776, 1, 435, 14666, 113366, 273141, 217015, 46656, 1, 882, 48783, 600124, 2385999, 3487218, 1845697, 279936, 1, 1777, 156933, 2937109, 17931235, 42397299
Offset: 1

Views

Author

Roger L. Bagula, Nov 30 2008

Keywords

Comments

Row sums are: {1, 7, 56, 504, 5040, 55440, 665280, 8648640, 121080960, 1816214400,...}. The sequences A008292, A144696,A144697,A144698,A144699 and this one, form a matrix of polynomials that are used in data smoothing calculations.

Examples

			{1},
{1, 6},
{1, 19, 36},
{1, 46, 241, 216},
{1, 101, 1091, 2551, 1296},
{1, 212, 4182, 18932, 24337, 7776},
{1, 435, 14666, 113366, 273141, 217015, 46656},
{1, 882, 48783, 600124, 2385999, 3487218, 1845697, 279936},
{1, 1777, 156933, 2937109, 17931235, 42397299, 40817623, 15159367, 1679616},
{1, 3568, 493900, 13631632, 121964374, 433696144, 667299052, 447815920, 121232113,10077696}
		

References

  • Douglas C. Montgomery, Lynwood A, Johnson, Forecasting and Time Series Analysis,McGraw-Hill, New York,1976,page 64.

Crossrefs

Programs

  • Mathematica
    M[p_, k_, x_] = ((-1)^k*(1 - x)^(p + k)/(k!(p - 1)!))*Sum[(p - 1 + j)!*j^k*x^j/(j!), {j, 0, Infinity}]/x;
    Table[Table[CoefficientList[FullSimplify[ExpandAll[(-1)^m*m!*M[n, m, x]]]/n, x], {m, 1, 10}], {n, 1, 10}];
    Table[Flatten[Table[CoefficientList[FullSimplify[ExpandAll[(-1)^m*m!*M[n, m, x]]]/n, x], {m, 1, 10}]], {n, 1, 10}]

Formula

m(p,k,x)=((-1)^k*(1 - x)^(p + k)/(k!(p - 1)!))*Sum[(p - 1 + j)!*j^k*x^j/(j!), {j, 0, Infinity}]/x;n=6;
t(m,l)=coefficients((-1)^m*m!*M[n, m, x])/n
Showing 1-6 of 6 results.