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-4 of 4 results.

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

Original entry on oeis.org

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

Views

Author

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

Keywords

Comments

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

Examples

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

References

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

Crossrefs

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

Programs

  • Haskell
    a132393 n k = a132393_tabl !! n !! k
    a132393_row n = a132393_tabl !! n
    a132393_tabl = map (map abs) a048994_tabl
    -- Reinhard Zumkeller, Nov 06 2013
    
  • Maple
    a132393_row := proc(n) local k; seq(coeff(expand(pochhammer (x,n)),x,k),k=0..n) end: # Peter Luschny, Nov 28 2010
  • Mathematica
    p[t_] = 1/(1 - t)^x; Table[ ExpandAll[(n!)SeriesCoefficient[ Series[p[t], {t, 0, 30}], n]], {n, 0, 10}]; a = Table[(n!)* CoefficientList[SeriesCoefficient[ Series[p[t], {t, 0, 30}], n], x], {n, 0, 10}]; Flatten[a] (* Roger L. Bagula, Apr 18 2008 *)
    Flatten[Table[Abs[StirlingS1[n,i]],{n,0,10},{i,0,n}]] (* Harvey P. Dale, Feb 04 2014 *)
  • Maxima
    create_list(abs(stirling1(n,k)),n,0,12,k,0,n); /* Emanuele Munarini, Mar 11 2011 */
    
  • PARI
    column(n,k) = my(v1, v2); v1 = vector(n-1, i, 0); v2 = vector(n, i, 0); v2[1] = 1; for(i=1, n-1, v1[i] = (i+k)*(i+k-1)/2*v2[i]; for(j=1, i-1, v1[j] *= (i-j)*(i+k)/(i-j+2)); v2[i+1] = vecsum(v1)/i); v2 \\ generates n first elements of the k-th column starting from the first nonzero element. - Mikhail Kurkov, Mar 05 2025

Formula

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

A000369 Triangle of numbers related to triangle A049213; generalization of Stirling numbers of second kind A008277, Bessel triangle A001497.

Original entry on oeis.org

1, 3, 1, 21, 9, 1, 231, 111, 18, 1, 3465, 1785, 345, 30, 1, 65835, 35595, 7650, 825, 45, 1, 1514205, 848925, 196245, 24150, 1680, 63, 1, 40883535, 23586255, 5755050, 775845, 62790, 3066, 84, 1, 1267389585, 748471185, 190482705, 27478710
Offset: 1

Views

Author

Keywords

Comments

a(n,m) := S2p(-3; n,m), a member of a sequence of triangles including S2p(-1; n,m) := A001497(n-1,m-1) (Bessel triangle) and ((-1)^(n-m))*S2p(1; n,m) := A008277(n,m) (Stirling 2nd kind). a(n,1)= A008545(n-1).
a(n,m), n>=m>=1, enumerates unordered n-vertex m-forests composed of m increasing plane (aka ordered) trees, with one vertex of out-degree r=0 (leafs or a root) and each vertex with out-degree r>=1 comes in r+2 types (like for an (r+2)-ary vertex). Proof from the e.g.f. of the first column Y(z):=1-(1-4*x)^(1/4) and the F. Bergeron et al. reference given in A001498, eq. (8), Y'(z)= phi(Y(z)), Y(0)=0, with out-degree o.g.f. phi(w)=1/(1-w)^3. - Wolfdieter Lang, Oct 12 2007
Also the Bell transform of the quadruple factorial numbers Product_{k=0..n-1} (4*k+3) (A008545) adding 1,0,0,0,... as column 0. For the definition of the Bell transform see A264428 and for cross-references A265606. - Peter Luschny, Dec 31 2015

Examples

			Triangle begins:
  1;
  3, 1;
  21, 9, 1;
  231, 111, 18, 1;
  3465, 1785, 345, 30, 1; ...
Tree combinatorics for a(3,2)=9: there are three m=2 forests each with one tree a root (with out-degree r=0) and the other tree a root and a leaf coming in three versions (like for a 3-ary vertex). Each such forest can be labeled increasingly in three ways (like (1,(23)), (2,(13)) and (3,(12))) yielding 9 such forests. - _Wolfdieter Lang_, Oct 12 2007
		

Crossrefs

Row sums give A016036. Cf. A004747.
Columns include A008545.
Alternating row sums A132163.

Programs

Formula

a(n, m) = n!*A049213(n, m)/(m!*4^(n-m)); a(n+1, m) = (4*n-m)*a(n, m) + a(n, m-1), n >= m >= 1; a(n, m) := 0, n
E.g.f. of m-th column: ((1-(1-4*x)^(1/4))^m)/m!.
From Peter Bala, Jun 08 2016: (Start)
With offset 0, the e.g.f. is 1/(1 - 4*x)^(3/4)*exp(t*(1 - (1 - 4*x)^(1/4))) = 1 + (3 + t)*x + (21 + 9*t + t^2)*x^2/2! + ....
Thus with row and column numbering starting at 0, this triangle is the exponential Riordan array [d/dx(F(x)), F(x)], belonging to the Derivative subgroup of the exponential Riordan group, where F(x) = 1 - (1 - 4*x)^(1/4).
Row polynomial recurrence: R(n+1,t) = t*Sum_{k = 0..n} binomial(n,k)*A008545(k)*R(n-k,t) with R(0,t) = 1. (End)

A051142 Generalized Stirling number triangle of first kind.

Original entry on oeis.org

1, -4, 1, 32, -12, 1, -384, 176, -24, 1, 6144, -3200, 560, -40, 1, -122880, 70144, -14400, 1360, -60, 1, 2949120, -1806336, 415744, -47040, 2800, -84, 1, -82575360, 53526528, -13447168, 1732864, -125440, 5152, -112, 1, 2642411520, -1795424256, 483835904
Offset: 1

Keywords

Comments

a(n,m) = R_n^m(a=0, b=4) in the notation of the given 1961 and 1962 references.
a(n,m) is a Jabotinsky matrix, i.e., the monic row polynomials E(n,x) := Sum_{m=1..n} a(n,m)*x^m = Product_{j=0..n-1} (x - 4*j), n >= 1, and E(0,x) := 1 are exponential convolution polynomials (see A039692 for the definition and a Knuth reference).
This is the signed Stirling1 triangle with diagonal d >= 0 (main diagonal d = 0) scaled with 4^d.
Also the Bell transform of the quadruple factorial numbers Product_{k=0..n-1} (4*k+4) (A047053) giving unsigned values and adding 1, 0, 0, 0, ... as column 0. For the definition of the Bell transform, see A264428 and for cross-references A265606. - Peter Luschny, Dec 31 2015

Examples

			Triangle a(n,m) (with rows n >= 1 and columns m = 1..n) begins:
        1;
       -4,     1;
       32,   -12,      1;
     -384,   176,    -24,    1;
     6144, -3200,    560,  -40,   1,
  -122880, 70144, -14400, 1360, -60, 1;
  ...
3rd row o.g.f.: E(3,x) = 32*x - 12*x^2 + x^3.
		

Crossrefs

First (m=1) column sequence is: A047053(n-1).
Row sums (signed triangle): A008545(n-1)*(-1)^(n-1).
Row sums (unsigned triangle): A007696(n).
Cf. A008275 (Stirling1 triangle, b=1), A039683 (b=2), A051141 (b=3).

Programs

  • Mathematica
    Table[StirlingS1[n, m] 4^(n - m), {n, 9}, {m, n}] // Flatten (* Michael De Vlieger, Dec 31 2015 *)
  • Sage
    # uses[bell_transform from A264428]
    # Unsigned values and an additional first column (1,0,0,0, ...).
    def A051142_row(n):
        multifact_4_4 = lambda n: prod(4*k + 4 for k in (0..n-1))
        mfact = [multifact_4_4(k) for k in (0..n)]
        return bell_transform(n, mfact)
    [A051142_row(n) for n in (0..9)] # Peter Luschny, Dec 31 2015

Formula

a(n, m) = a(n-1, m-1) - 4*(n-1)*a(n-1, m) for n >= m >= 1; a(n, m) := 0 for n < m; a(n, 0) := 0 for n >= 1; a(0, 0) = 1.
E.g.f. for the m-th column of the signed triangle: (log(1 + 4*x)/4)^m/m!.
a(n, m) = S1(n, m)*4^(n-m), with S1(n, m) := A008275(n, m) (signed Stirling1 triangle).

A119274 Triangle of coefficients of numerators in Padé approximation to exp(x).

Original entry on oeis.org

1, 2, 1, 12, 6, 1, 120, 60, 12, 1, 1680, 840, 180, 20, 1, 30240, 15120, 3360, 420, 30, 1, 665280, 332640, 75600, 10080, 840, 42, 1, 17297280, 8648640, 1995840, 277200, 25200, 1512, 56, 1, 518918400, 259459200, 60540480, 8648640, 831600, 55440, 2520
Offset: 0

Author

Paul Barry, May 12 2006

Keywords

Comments

n-th numerator of Padé approximation is (1/n!)*sum{j=0..n, C(n,j)(2n-j)!x^j}. Reversal of A113025. Row sums are A001517. First column is A001813. Inverse is A119275.
Also the Bell transform of the quadruple factorial numbers Product_{k=0..n-1} (4*k+2) (A001813) adding 1,0,0,0,... as column 0. For the definition of the Bell transform see A264428 and for cross-references A265606. - Peter Luschny, Dec 31 2015
Dividing each diagonal by its initial element generates A054142. - Tom Copeland, Oct 10 2016

Examples

			Triangle begins
1,
2, 1,
12, 6, 1,
120, 60, 12, 1,
1680, 840, 180, 20, 1,
30240, 15120, 3360, 420, 30, 1
		

Crossrefs

Programs

  • Maple
    # The function BellMatrix is defined in A264428.
    # Adds (1,0,0,0, ..) as column 0.
    BellMatrix(n -> (2*n)!/n!, 9); # Peter Luschny, Jan 27 2016
  • Mathematica
    BellMatrix[f_Function, len_] := With[{t = Array[f, len, 0]}, Table[BellY[n, k, t], {n, 0, len - 1}, {k, 0, len - 1}]];
    rows = 12;
    M = BellMatrix[(2#)!/#!&, rows];
    Table[M[[n, k]], {n, 2, rows}, {k, 2, n}] // Flatten (* Jean-François Alcover, Jun 24 2018, after Peter Luschny *)
  • Sage
    # uses[bell_transform from A264428]
    # Adds a column 1,0,0,0,... at the left side of the triangle.
    def A119274_row(n):
        multifact_4_2 = lambda n: prod(4*k + 2 for k in (0..n-1))
        mfact = [multifact_4_2(k) for k in (0..n)]
        return bell_transform(n, mfact)
    [A119274_row(n) for n in (0..9)] # Peter Luschny, Dec 31 2015

Formula

Number triangle T(n,k) = C(n,k)(2n-k)!/n!.
After adding a leading column (1,0,0,0,...), the triangle gives the coefficients of the Sheffer associated sequence (binomial-type polynomials) for the delta (lowering) operator D(1-D) with e.g.f. exp[ x * (1 - sqrt(1-4t)) / 2 ] . See Mathworld on Sheffer sequences. See A134685 for relation to Catalan numbers. - Tom Copeland, Feb 09 2008
Showing 1-4 of 4 results.