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

A008275 Triangle read by rows of Stirling numbers of first kind, s(n,k), n >= 1, 1 <= k <= n.

Original entry on oeis.org

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

Views

Author

Keywords

Comments

The unsigned numbers are also called Stirling cycle numbers: |s(n,k)| = number of permutations of n objects with exactly k cycles.
The unsigned numbers (read from right to left) also give the number of permutations of 1..n with complexity k, where the complexity of a permutation is defined to be the sum of the lengths of the cycles minus the number of cycles. In other words, the complexity equals the sum of (length of cycle)-1 over all cycles. For n=5, the numbers of permutations with complexity 0,1,2,3,4 are 1, 10, 35, 50, 24. - N. J. A. Sloane, Feb 08 2019
The unsigned numbers are also the number of permutations of 1..n with k left to right maxima (see Khovanova and Lewis, Smith).
With P(n) = the number of integer partitions of n, T(i,n) = the number of parts of the i-th partition of n, D(i,n) = the number of different parts of the i-th partition of n, p(j,i,n) = the j-th part of the i-th partition of n, m(j,i,n) = multiplicity of the j-th part of the i-th partition of n, Sum_[T(i,n)=k]{i=1}^{P(n)} = sum running from i=1 to i=p(n) but taking only partitions with T(i,n)=k parts into account, Product{j=1..T(i,n)} = product running from j=1 to j=T(i,n), Product_{j=1..D(i,n)} = product running from j=1 to j=D(i,n) one has S1(n,k) = Sum_[T(i,n)=k]{i=1}^{P(n)} (n!/Product{j=1..T(i,n)} p(j,i,n))* (1/Product_{j=1..D(i,n)} m(j,i,n)!). For example, S1(6,3) = 225 because n=6 has the following partitions with k=3 parts: (114), (123), (222). Their complexions are: (114): (6!/1*1*4)*(1/2!*1!) = 90, (123): (6!/1*2*3)*(1/1!*1!*1!) = 120, (222): (6!/2*2*2)*(1/3!) = 15. The sum of the complexions is 90+120+15 = 225 = S1(6,3). - Thomas Wieder, Aug 04 2005
Row sums equal 0. - Jon Perry, Nov 14 2005
|s(n,k)| enumerates unordered n-vertex forests composed of k increasing non-plane (unordered) trees. Proof from the e.g.f. of the first column and the F. Bergeron et al. reference, especially Table 1, last row (non-plane "recursive"), given in A049029. - Wolfdieter Lang, Oct 12 2007
|s(n,k)| enumerates unordered increasing n-vertex k-forests composed of k unary trees (out-degree r from {0,1}) whose vertices of depth (distance from the root) j >= 0 come in j+1 colors (j=0 for the k roots). - Wolfdieter Lang, Oct 12 2007, Feb 22 2008
A refinement of the unsigned array is A036039. For an association to forests of "naturally grown" rooted non-planar trees, dispositions of flags on flagpoles, and colorings of the vertices of the complete graphs K_n, see A130534. - Tom Copeland, Mar 30 and Apr 05 2014
The Stirling numbers of the first kind were related to the falling factorial and the convolved, or generalized, Bernoulli numbers B_n by Norlund in 1924 through Sum_{k=1..n+1} T(n+1,k) * x^(k-1) = (x-1)!/(x-1-n)! = (x + B.(0))^n = B_n(x), umbrally evaluated with (B.(0))^k = B_k(0) and the associated Appell polynomial B_n(x) defined by the e.g.f. (t/(exp(t) - 1))^(n+1) * exp(x*t) = exp(B.(x)t). - Tom Copeland, Sep 29 2015
With x = e^z, D_x = d/dx, D_z = d/dz, and p_n(x) the row polynomials of this entry, x^n (D_x)^n = p_n(D_z) = (D_z)! / (D_z - n)! = (xD_x)! / (xD_x - n)!. - Tom Copeland, Nov 27 2015
From the operator relation z + Psi(1) + sum_{n > 0} (-1)^n (-1/n) binomial(D,n) = z + Psi(1+D) with D = d/dz and Psi the digamma function, Zeta(n+1) = Sum_{k > n-1} (1/k) |S(k,n)| / k! for n > 0 and Zeta the Riemann zeta function. - Tom Copeland, Aug 12 2016
Let X_1,...,X_n be i.i.d. random variables with exponential distribution having mean = 1. Let Y = max{X_1,...,X_n}. Then (-1)^n*n!/( Sum_{k=1..n+1} a(n+1,k) t^(k-1) ) is the moment generating function of Y. The expectation of Y is the n-th harmonic number. - Geoffrey Critzer, Dec 25 2018
In the Ewens sampling theory describing the multivariate probability distribution of the sizes of the allelic classes in a sample of size n under the Infinite Alleles Model, |s(n,k)| gives the coefficient in the formula for the probability that a sample of n alleles has exactly k distinct types. - Noah A Rosenberg, Feb 10 2019
Named by Nielson (1906) after the Scottish mathematician James Stirling (1692-1770). - Amiram Eldar, Jun 11 2021 and Oct 02 2023
The first few row polynomials along with a recursion formula are found in a manuscript by Newton written in 1664 or 1665 (p. 169 of Turnbull) giving a geometric presentation of the binomial theorem for rational powers. - Tom Copeland, Dec 10 2022

Examples

			|s(3,2)| = 3 for the three unordered 2-forest with 3 vertices and two increasing (nonplane) trees: ((1),(2,3)), ((2),(1,3)), ((3),(1,2)).
Triangle begins:
                                      1
                                 -1,      1
                               2,    -3,      1
                          -6,    11,     -6,     1
                      24,    -50,    35,    -10,    1
                -120,    274,  -225,     85,   -15,    1
             720,  -1764,   1624,  -735,    175,  -21,   1
       -5040,  13068, -13132,  6769,  -1960,   322,  -28,  1
  40320, -109584, 118124, -67284, 22449,  -4536,  546, -36,  1
Another version of the same triangle, from _Joerg Arndt_, Oct 05 2009: (Start)
s(n,k) := number of permutations of n elements with exactly k cycles ("Stirling cycle numbers")
  n|  total   m=1      2      3     4     5    6   7  8 9
  -+-----------------------------------------------------
  1|      1     1
  2|      2     1      1
  3|      6     2      3      1
  4|     24     6     11      6     1
  5|    120    24     50     35    10     1
  6|    720   120    274    225    85    15    1
  7|   5040   720   1764   1624   735   175   21   1
  8|  40320  5040  13068  13132  6769  1960  322  28  1
  9| 362880 40320 109584 118124 67284 22449 4536 546 36 1
(End)
|s(4,2)| = 11 for the eleven unordered 2-forest with 4 vertices, composed of two increasing (nonplane) trees: ((1),((23)(24))), ((2),((13)(14))), ((3),((12)(14))), ((4),((12)(13))); ((1),(2,3,4)),((2),(1,2,3)), ((3), (1,2,4)), ((4),(1,2,3)); ((1,2),(3,4)), ((1,3),(2,4)), ((1,4),(2,3)). - _Wolfdieter Lang_, Feb 22 2008
		

References

  • M. Abramowitz and I. A. Stegun, eds., Handbook of Mathematical Functions, National Bureau of Standards Applied Math. Series 55, 1964 (and various reprintings), p. 833.
  • Arthur T. Benjamin and Jennifer Quinn, Proofs that really count: the art of combinatorial proof, M.A.A. 2003, p. 93ff.
  • Boris A. Bondarenko, Generalized Pascal Triangles and Pyramids (in Russian), FAN, Tashkent, 1990, ISBN 5-648-00738-8.
  • George Boole, Finite Differences, 5th ed. New York, NY: Chelsea, 1970.
  • Louis Comtet, Advanced Combinatorics, Reidel, 1974; Chapter V, also p. 310.
  • John H. Conway and Richard K. Guy, The Book of Numbers, Copernicus Press, NY, 1996, p. 93.
  • Florence Nightingale David, Maurice George Kendall and David Elliot Barton, Symmetric Function and Allied Tables, Cambridge, 1966, p. 226.
  • Saber N. Elaydi, An Introduction to Difference Equations, 3rd ed. Springer, 2005.
  • Herman H. Goldstine, A History of Numerical Analysis, Springer-Verlag, 1977; Section 2.7.
  • Ronald L. Graham, Donald E. Knuth and Oren Patashnik, Concrete Mathematics. Addison-Wesley, Reading, MA, 1990, p. 245. In the second edition, see Chapter 6, especially p. 259.
  • M. Miyata and J. W. Son, On the complexity of permutations and the metric space of bijections, Tensor, 60 (1998), No. 1, 109-116 (MR1768839).
  • Isaac Newton, A Method whereby to find ye areas of Those Lines wch can be squared, pp. 168-171 of Turnbull below.
  • John Riordan, An Introduction to Combinatorial Analysis, p. 48.
  • Robert Sedgewick and Phillipe Flajolet, An Introduction to the Analysis of Algorithms, Addison-Wesley, Reading, MA, 1996.
  • H. Turnbull (editor), The Correspondence of Isaac Newton Vol. II 1676-1687, Cambridge Univ. Press, 1960.

Crossrefs

Diagonals: A000217, A000914, A001303, A000915, A053567, etc.
Cf. A048994, A008277 (Stirling numbers of second kind), A039814, A039815, A039816, A039817, A048993, A087748.
Cf. A084938, A094216, A008276 (row reversed), A008277, A008278, A094262, A121632, A130534 (unsigned version), A087755 (triangle mod 2), A000142 (row sums of absolute values).

Programs

  • Haskell
    a008275 n k = a008275_tabl !! (n-1) !! (k-1)
    a008275_row n = a008275_tabl !! (n-1)
    a008275_tabl = map tail $ tail a048994_tabl
    -- Reinhard Zumkeller, Mar 18 2013
  • Maple
    with (combinat):seq(seq(stirling1(n, k), k=1..n), n=1..10); # Zerinvary Lajos, Jun 03 2007
    for i from 0 to 9 do seq(stirling1(i, j), j = 1 .. i) od; # Zerinvary Lajos, Nov 29 2007
  • Mathematica
    Flatten[Table[StirlingS1[n, k], {n, 1, 10}, {k, 1, n}]] (* Jean-François Alcover, May 18 2011 *)
    Flatten@Table[Coefficient[Product[x-k, {k, 0, n-1}], x, Range[n]], {n, Range[10]}] (* Oliver Seipel, Jun 11 2024 *)
    a[n_, n_] := 1; a[n_, 0] := 0; a[0, k_] := 0;
    a[n_, k_] := a[n, k] = a[n-1, k-1] + (n-1) a[n-1, k];
    Flatten@Table[(-1)^(n-k) a[n, k], {n, 1, 10}, {k, 1, n}] (* Oliver Seipel, Jun 11 2024 *)
  • Maxima
    create_list(stirling1(n+1,k+1),n,0,30,k,0,n); /* Emanuele Munarini, Jun 01 2012 */
    
  • PARI
    T(n,k)=if(n<1,0,n!*polcoeff(binomial(x,n),k))
    
  • PARI
    T(n,k)=if(n<1,0,n!*polcoeff(polcoeff((1+x+x*O(x^n))^y,n),k))
    
  • PARI
    vecstirling(n)=Vec(factorback(vector(n-1,i,1-i*'x))) /* (A function that returns all the s(n,k) as a vector) */ \\ Bill Allombert (Bill.Allombert(AT)math.u-bordeaux1.fr), Mar 16 2009
    

Formula

s(n, k) = s(n-1, k-1) - (n-1)*s(n-1, k), n, k >= 1; s(n, 0) = s(0, k) = 0; s(0, 0) = 1.
The unsigned numbers a(n, k)=|s(n, k)| satisfy a(n, k) = a(n-1, k-1) + (n-1)*a(n-1, k), n, k >= 1; a(n, 0) = a(0, k) = 0; a(0, 0) = 1.
E.g.f.: for m-th column (unsigned): ((-log(1-x))^m)/m!.
s(n, k) = T(n-1, k-1), n>1 and k>1, where T(n, k) is the triangle [ -1, -1, -2, -2, -3, -3, -4, -4, -5, -5, -6, -6, ...] DELTA [1, 0, 1, 0, 1, 0, 1, 0, 1, ...] and DELTA is Deléham's operator defined in A084938. The unsigned numbers are also |s(n, k)| = T(n-1, k-1), for n>0 and k>0, where T(n, k) = [1, 1, 2, 2, 3, 3, 4, 4, 5, 5, ...] DELTA [1, 0, 1, 0, 1, 0, 1, 0, ...].
Sum_{i=0..n} (-1)^(n-i) * StirlingS1(n, i) * binomial(i, k) = (-1)^(n-k) * StirlingS1(n+1, k+1). - Carlo Wood (carlo(AT)alinoe.com), Feb 13 2007
G.f. for row n: Product_{j=1..n} (x-j) (e.g., (x-1)*(x-2)*(x-3) = x^3 - 6*x^2 + 11*x - 6). - Jon Perry, Nov 14 2005
s(n,k) = A048994(n,k), for k=1..n. - Reinhard Zumkeller, Mar 18 2013 (Corrected by N. J. A. Sloane, May 07 2025 at the suggestion of Manfred Boergens, May 07 2025)
As lower triangular matrices A008277*A008275 = I, the identity matrix. - Tom Copeland, Apr 25 2014
a(n,k) = s(n,k) = lim_{y -> 0} Sum_{j=0..k} (-1)^j*binomial(k,j)*((-j*y)!/(-j*y-n)!)*y^(-k)/k! = Sum_{j=0..k} (-1)^(n-j)*binomial(k,j)*((j*y - 1 + n)!/(j*y-1)!)*y^(-k)/k!. - Tom Copeland, Aug 28 2015
From Daniel Forgues Jan 16 2016: (Start)
Let x_(0) := 1 (empty product), and for n >= 1:
x_(n) := Product_{k=0..n-1} (x-k), called a factorial term (Boole, 1970) or a factorial polynomial (Elaydi, 2005: p. 60), and also x_(-n) := 1 / [Product_{k=0..n-1} (x+k)].
Then, for n >= 1: x_(n) = Sum_{k=1..n} T(n,k) * x^k, 1 / [x_(-n)] = Sum_{k=1..n} |T(n,k)| * x^k, x^n = Sum_{k=1..n} A008277(n,k) * x_(k), where A008277(n,k) are Stirling numbers of the second kind.
The row sums (of either signed or absolute values) are Sum_{k=1..n} T(n,k) = 0^(n-1), Sum_{k=1..n} |T(n,k)| = T(n+1,1) = n!. (End)
s(n,m) = ((-1)^(n-m)/n)*Sum_{i=0..m-1} C(2*n-m-i, m-i-1)*A008517(n-m+1,n-m-i+1). - Vladimir Kruchinin, Feb 14 2018
Orthogonal relation: Sum_{i=0..n} i^p*Sum_{j=k..n} (-1)^(i+j) * binomial(j,i) * Stirling1(j,k)/j! = delta(p,k), i,k,p <= n, n >= 1. - Leonid Bedratyuk, Jul 27 2020
From Zizheng Fang, Dec 28 2020: (Start)
Sum_{k=1..n} (-1)^k * k * T(n, k) = -T(n+1, 2).
Sum_{k=1..n} k * T(n, k) = (-1)^n * (n-2)! = T(n-1, 1) for n>=2. (End)
n-th row polynomial = n!*Sum_{k = 0..2*n} (-1)^(n+k)*binomial(x, k)*binomial(x-1, 2*n-k) = n!*Sum_{k = 0..2*n+1} (-1)^(n+k+1)*binomial(x, k)*binomial(x-1, 2*n+1-k). - Peter Bala, Mar 29 2024

A048994 Triangle of Stirling numbers of first kind, s(n,k), n >= 0, 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

Keywords

Comments

The unsigned numbers are also called Stirling cycle numbers: |s(n,k)| = number of permutations of n objects with exactly k cycles.
Mirror image of the triangle A054654. - Philippe Deléham, Dec 30 2006
Also the triangle gives coefficients T(n,k) of x^k in the expansion of C(x,n) = (a(k)*x^k)/n!. - Mokhtar Mohamed, Dec 04 2012
From Wolfdieter Lang, Nov 14 2018: (Start)
This is the Sheffer triangle of Jabotinsky type (1, log(1 + x)). See the e.g.f. of the triangle below.
This is the inverse Sheffer triangle of the Stirling2 Sheffer triangle A008275.
The a-sequence of this Sheffer triangle (see a W. Lang link in A006232)
is from the e.g.f. A(x) = x/(exp(x) -1) a(n) = Bernoulli(n) = A027641(n)/A027642(n), for n >= 0. The z-sequence vanishes.
The Boas-Buck sequence for the recurrences of columns has o.g.f. B(x) = Sum_{n>=0} b(n)*x^n = 1/((1 + x)*log(1 + x)) - 1/x. b(n) = (-1)^(n+1)*A002208(n+1)/A002209(n+1), b = {-1/2, 5/12, -3/8, 251/720, -95/288, 19087/60480,...}. For the Boas-Buck recurrence of Riordan and Sheffer triangles see the Aug 10 2017 remark in A046521, adapted to the Sheffer case, also for two references. See the recurrence and example below. (End)
Let G(n,m,k) be the number of simple labeled graphs on [n] with m edges and k components. Then T(n,k) = Sum (-1)^m*G(n,m,k). See the Read link below. Equivalently, T(n,k) = Sum mu(0,p) where the sum is over all set partitions p of [n] containing k blocks and mu is the Moebius function in the incidence algebra associated to the set partition lattice on [n]. - Geoffrey Critzer, May 11 2024

Examples

			Triangle begins:
  n\k 0     1       2       3      4      5      6    7    8   9 ...
  0   1
  1   0     1
  2   0    -1       1
  3   0     2      -3       1
  4   0    -6      11      -6      1
  5   0    24     -50      35    -10      1
  6   0  -120     274    -225     85    -15      1
  7   0   720   -1764    1624   -735    175    -21    1
  8   0 -5040   13068  -13132   6769  -1960    322  -28    1
  9   0 40320 -109584  118124 -67284  22449  -4536  546  -36   1
  ... - _Wolfdieter Lang_, Aug 22 2012
------------------------------------------------------------------
From _Wolfdieter Lang_, Nov 14 2018: (Start)
Recurrence: s(5,2)= s(4, 1) - 4*s(4, 2) = -6 - 4*11 = -50.
Recurrence from the a- and z-sequences: s(6, 3) = 2*(1*1*(-50) + 3*(-1/2)*35 + 6*(1/6)*(-10) + 10*0*1) = -225.
Boas-Buck recurrence for column k = 3, with b = {-1/2, 5/12, -3/8, ...}:
s(6, 3) = 6!*((-3/8)*1/3! + (5/12)*(-6)/4! + (-1/2)*35/5!) = -225. (End)
		

References

  • M. Abramowitz and I. A. Stegun, eds., Handbook of Mathematical Functions, National Bureau of Standards Applied Math. Series 55, 1964 (and various reprintings), p. 833.
  • L. Comtet, Advanced Combinatorics, Reidel, 1974; Chapter V, also p. 310.
  • J. H. Conway and R. K. Guy, The Book of Numbers, Copernicus Press, NY, 1996, p. 93.
  • F. N. David, M. G. Kendall and D. E. Barton, Symmetric Function and Allied Tables, Cambridge, 1966, p. 226.
  • R. L. Graham, D. E. Knuth and O. Patashnik, Concrete Mathematics. Addison-Wesley, Reading, MA, 1990, p. 245.
  • J. Riordan, An Introduction to Combinatorial Analysis, p. 48.

Crossrefs

See especially A008275 which is the main entry for this triangle. A132393 is an unsigned version, and A008276 is another version.
A000142(n) = Sum_{k=0..n} |s(n, k)| for n >= 0.
Row sums give A019590(n+1).

Programs

  • Haskell
    a048994 n k = a048994_tabl !! n !! k
    a048994_row n = a048994_tabl !! n
    a048994_tabl = map fst $ iterate (\(row, i) ->
    (zipWith (-) ([0] ++ row) $ map (* i) (row ++ [0]), i + 1)) ([1], 0)
    -- Reinhard Zumkeller, Mar 18 2013
  • Maple
    A048994:= proc(n,k) combinat[stirling1](n,k) end: # R. J. Mathar, Feb 23 2009
    seq(print(seq(coeff(expand(k!*binomial(x,k)),x,i),i=0..k)),k=0..9); # Peter Luschny, Jul 13 2009
    A048994_row := proc(n) local k; seq(coeff(expand(pochhammer(x-n+1,n)), x,k), k=0..n) end: # Peter Luschny, Dec 30 2010
  • Mathematica
    Table[StirlingS1[n, m], {n, 0, 9}, {m, 0, n}] (* Peter Luschny, Dec 30 2010 *)
  • Maxima
    create_list(stirling1(n,k),n,0,12,k,0,n); /* Emanuele Munarini, Mar 11 2011 */
    
  • PARI
    a(n,k) = if(k<0 || k>n,0, if(n==0,1,(n-1)*a(n-1,k)+a(n-1,k-1)))
    
  • PARI
    trg(nn)=for (n=0, nn-1, for (k=0, n, print1(stirling(n,k,1), ", ");); print();); \\ Michel Marcus, Jan 19 2015
    

Formula

s(n, k) = A008275(n,k) for n >= 1, k = 1..n; column k = 0 is {1, repeat(0)}.
s(n, k) = s(n-1, k-1) - (n-1)*s(n-1, k), n, k >= 1; s(n, 0) = s(0, k) = 0; s(0, 0) = 1.
The unsigned numbers a(n, k)=|s(n, k)| satisfy a(n, k)=a(n-1, k-1)+(n-1)*a(n-1, k), n, k >= 1; a(n, 0) = a(0, k) = 0; a(0, 0) = 1.
Triangle (signed) = [0, -1, -1, -2, -2, -3, -3, -4, -4, ...] DELTA [1, 0, 1, 0, 1, 0, 1, 0, ...]; Triangle(unsigned) = [0, 1, 1, 2, 2, 3, 3, 4, 4, ...] DELTA [1, 0, 1, 0, 1, 0, 1, 0, ...]; where DELTA is Deléham's operator defined in A084938.
Sum_{k=0..n} (-m)^(n-k)*s(n, k) = A000142(n), A001147(n), A007559(n), A007696(n), ... for m = 1, 2, 3, 4, ... .- Philippe Deléham, Oct 29 2005
A008275*A007318 as infinite lower triangular matrices. - Gerald McGarvey, Aug 20 2009
T(n,k) = n!*[x^k]([t^n]exp(x*log(1+t))). - Peter Luschny, Dec 30 2010, updated Jun 07 2020
From Wolfdieter Lang, Nov 14 2018: (Start)
Recurrence from the Sheffer a-sequence (see a comment above): s(n, k) = (n/k)*Sum_{j=0..n-k} binomial(k-1+j, j)*Bernoulli(j)*s(n-1, k-1+j), for n >= 1 and k >= 1, with s(n, 0) = 0 if n >= 1, and s(0,0) = 1.
Boas-Buck type recurrence for column k: s(n, k) = (n!*k/(n - k))*Sum_{j=k..n-1} b(n-1-j)*s(j, k)/j!, for n >= 1 and k = 0..n-1, with input s(n, n) = 1. For sequence b see the Boas-Buck comment above. (End)
T(n,k) = Sum_{j=k..n} (-1)^(n-j)*A271705(n,j)*A216294(j,k). - Mélika Tebni, Feb 23 2023

Extensions

Offset corrected by R. J. Mathar, Feb 23 2009
Formula corrected by Philippe Deléham, Sep 10 2009

A039814 Matrix square of Stirling-1 triangle A008275.

Original entry on oeis.org

1, -2, 1, 7, -6, 1, -35, 40, -12, 1, 228, -315, 130, -20, 1, -1834, 2908, -1485, 320, -30, 1, 17582, -30989, 18508, -5005, 665, -42, 1, -195866, 375611, -253400, 81088, -13650, 1232, -56, 1, 2487832, -5112570, 3805723, -1389612, 279048, -32130, 2100, -72, 1
Offset: 1

Views

Author

Christian G. Bower, Feb 15 1999

Keywords

Comments

Exponential Riordan array [1/((1 + x)*(1 + log(1 + x))), log(1 + log(1 + x))]. The row sums of the unsigned array give A007840 (apart from the initial term). - Peter Bala, Jul 22 2014
Also the Bell transform of (-1)^n*A003713(n+1). For the definition of the Bell transform see A264428. - Peter Luschny, Jan 28 2016

Examples

			Triangle begins:
      1;
     -2,    1;
      7,   -6,     1;
    -35,   40,   -12,   1;
    228, -315,   130, -20,   1;
  -1834, 2908, -1485, 320, -30, 1;
...
		

Crossrefs

Column k=1..3 give (-1)^(n-1) * A003713(n), (-1)^n * A341587(n), (-1)^(n-1) * A341588(n).
Cf. A007840.

Programs

  • Maple
    # The function BellMatrix is defined in A264428.
    # Adds (1,0,0,0, ..) as column 0.
    BellMatrix(n -> (-1)^n*add(k!*abs(Stirling1(n+1,k+1)), k=0..n), 10); # Peter Luschny, Jan 28 2016
  • Mathematica
    max = 9; t = Table[StirlingS1[n, k], {n, 1, max}, {k, 1, max}]; t2 = t.t; Table[t2[[n, k]], {n, 1, max}, {k, 1, n}] // Flatten (* Jean-François Alcover, Feb 01 2013 *)
    rows = 9;
    t = Table[(-1)^n*Sum[k!*Abs[StirlingS1[n+1, k+1]], {k,0,n}], {n, 0, rows}];
    T[n_, k_] := BellY[n, k, t];
    Table[T[n, k], {n, 1, rows}, {k, 1, n}] // Flatten (* Jean-François Alcover, Jun 22 2018, after Peter Luschny *)
  • PARI
    T(n, k) = sum(j=0, n, stirling(n, j, 1)*stirling(j, k, 1)); \\ Seiichi Manyama, Feb 13 2022

Formula

E.g.f. of k-th column: ((log(1+log(1+x)))^k)/k!.
E.g.f.: 1/(1 + t)*( 1 + log(1 + t) )^(x-1) = 1 + (-2 + x)*t + (7 - 6*x + x^2)*t^2/2! + .... - Peter Bala, Jul 22 2014
T(n,k) = Sum_{j=0..n} Stirling1(n,j) * Stirling1(j,k). - Seiichi Manyama, Feb 13 2022

A000359 Coefficients of iterated exponentials.

Original entry on oeis.org

1, 5, 40, 440, 6170, 105315, 2120610, 49242470, 1296133195, 38152216495, 1242274374380, 44345089721923, 1722416374173854, 72330102999829054, 3265871028909088036, 157797437377747327987, 8124524883679977475839, 444098724261935142753430
Offset: 1

Views

Author

Keywords

References

  • J. Ginsburg, Iterated exponentials, Scripta Math., 11 (1945), 340-353.
  • N. J. A. Sloane, A Handbook of Integer Sequences, Academic Press, 1973 (includes this sequence).
  • N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).

Crossrefs

a(n) = |A039817(n, 1)| (first column of triangle). Cf. A003713, A000268, A000310, A000406, A001765.

Programs

  • Mathematica
    max = 20; CoefficientList[-Log[1 + Log[1 + Log[1 + Log[1 + Log[1 - x]]]]]/x + O[x]^max, x]*Range[max]! (* Jean-François Alcover, Feb 08 2016 *)
  • PARI
    T(n, k) = if(k==1, (n-1)!, sum(j=1, n, abs(stirling(n, j, 1))*T(j, k-1)));
    a(n) = T(n, 5); \\ Seiichi Manyama, Feb 11 2022
    
  • PARI
    my(N=20, x='x+O('x^N)); Vec(serlaplace(-log(1+log(1+log(1+log(1+log(1-x))))))) \\ Seiichi Manyama, Feb 11 2022

Formula

E.g.f.: -log(1+log(1+log(1+log(1+log(1-x))))).

A039815 Triangle read by rows: matrix cube of the Stirling-1 triangle A008275.

Original entry on oeis.org

1, -3, 1, 15, -9, 1, -105, 87, -18, 1, 947, -975, 285, -30, 1, -10472, 12657, -4680, 705, -45, 1, 137337, -188090, 82887, -15960, 1470, -63, 1, -2085605, 3159699, -1598954, 370237, -43890, 2730, -84, 1, 36017472, -59326371, 33613353, -9009294, 1292067, -103950, 4662, -108, 1
Offset: 1

Views

Author

Christian G. Bower, Feb 15 1999

Keywords

Examples

			Triangle begins:
       1;
      -3,     1;
      15,    -9,     1;
    -105,    87,   -18,   1;
     947,  -975,   285, -30,   1;
  -10472, 12657, -4680, 705, -45, 1;
  ...
		

Crossrefs

Cf. A000268 (first column), A008275.

Programs

  • Maple
    T:= Matrix(10,10,(i,j) -> `if`(i>= j, combinat:-stirling1(i,j),0)):
    M:= T^3:
    seq(seq(M[i,j],j=1..i),i=1..10); # Robert Israel, Sep 12 2022
  • Mathematica
    Flatten[Table[SeriesCoefficient[(Log[1+Log[1+Log[1+x]]])^k, {x,0,n}] n!/k!, {n,9}, {k,n}]] (* Stefano Spezia, Sep 12 2022 *)

Formula

E.g.f. of k-th column: ((log(1+log(1+log(1+x))))^k)/k!.

A039816 Triangle read by rows: matrix 4th power of the Stirling-1 triangle A008275.

Original entry on oeis.org

1, -4, 1, 26, -12, 1, -234, 152, -24, 1, 2696, -2210, 500, -40, 1, -37919, 36976, -10710, 1240, -60, 1, 630521, -704837, 245896, -36750, 2590, -84, 1, -12111114, 15132932, -6120324, 1109696, -101500, 4816, -112, 1, 264051201, -362099010, 165387680, -34990620, 3901296, -241164, 8232, -144, 1
Offset: 1

Views

Author

Christian G. Bower, Feb 15 1999

Keywords

Examples

			Triangle begins:
       1;
      -4,     1;
      26,   -12,      1;
    -234,   152,    -24,    1;
    2696, -2210,    500,  -40,   1;
  -37919, 36976, -10710, 1240, -60, 1;
  ...
		

Crossrefs

Cf. A000310 (first column), A008275.

Programs

  • Maple
    T:= Matrix(10,10,(i,j) -> `if`(i>= j, combinat:-stirling1(i,j),0)):
    M:= T^4:
    seq(seq(M[i,j],j=1..i),i=1..10); # Robert Israel, Sep 12 2022
  • Mathematica
    Flatten[Table[SeriesCoefficient[(Log[1+Log[1+Log[1+Log[1+x]]]])^k,{x,0,n}] n!/k!, {n,9}, {k,n}]] (* Stefano Spezia, Sep 12 2022 *)

Formula

E.g.f. of k-th column: ((log(1+log(1+log(1+log(1+x)))))^k)/k!.

A351527 Expansion of e.g.f. (log(1 + log(1 + log(1 + log(1 + log(1+ x))))))^2 / 2.

Original entry on oeis.org

1, -15, 235, -4200, 86020, -2001055, 52305780, -1520815230, 48747603100, -1709228504170, 65115320810260, -2679459929923699, 118482699493123571, -5604477255138004835, 282449438671531808676, -15111729578894643263239
Offset: 2

Views

Author

Seiichi Manyama, Feb 13 2022

Keywords

Crossrefs

Programs

  • PARI
    my(N=20, x='x+O('x^N)); Vec(serlaplace(log(1+log(1+log(1+log(1+log(1+x)))))^2/2))
    
  • PARI
    T(n, k) = if(k==0, n==1, sum(j=0, n, abs(stirling(n, j, 1))*T(j, k-1)));
    a(n) = (-1)^n*sum(k=1, n-1, binomial(n-1, k)*T(k, 5)*T(n-k, 5));

Formula

a(n) = (-1)^n * Sum_{k=1..n-1} binomial(n-1,k) * A000359(k) * A000359(n-k).
Showing 1-7 of 7 results.