A216294 -id:A216294 - cp's OEIS frontend

A048994 Triangle of Stirling numbers of first kind, s(n,k), n >= 0, 0 <= k <= n.

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

N. J. A. Sloane

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

			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)

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.

T. D. Noe, Rows n=0..100 of triangle, flattened
M. Abramowitz and I. A. Stegun, eds., Handbook of Mathematical Functions, National Bureau of Standards, Applied Math. Series 55, Tenth Printing, 1972 [alternative scanned copy].
Paul Barry, Generalized Stirling Numbers, Exponential Riordan Arrays, and Toda Chain Equations, Journal of Integer Sequences, 17 (2014), #14.2.3.
Fatima Zohra Bensaci, Rachid Boumahdi, and Laala Khaldi, Finite Sums Involving Fibonacci and Lucas Numbers, J. Int. Seq. (2024). See p. 9.
R. M. Dickau, Stirling numbers of the first kind. [Illustrates the unsigned Stirling cycle numbers A132393.]
Askar Dzhumadil'daev and Damir Yeliussizov, Walks, partitions, and normal ordering, Electronic Journal of Combinatorics, 22(4) (2015), #P4.10.
Bill Gosper, Colored illustrations of triangle of Stirling numbers of first kind read mod 2, 3, 4, 5, 6, 7.
Gergő Nemes, An Asymptotic Expansion for the Bernoulli Numbers of the Second Kind, J. Int. Seq. 14 (2011), #11.4.8.
A. Hennessy and Paul Barry, Generalized Stirling Numbers, Exponential Riordan Arrays, and Orthogonal Polynomials, J. Int. Seq. 14 (2011), #11.8.2 (A-number typo A048894).
NIST Digital Library of Mathematical Functions, Stirling Numbers
Ken Ono, Larry Rolen, and Florian Sprung, Zeta-Polynomials for modular form periods, p. 6, arXiv:1602.00752 [math.NT], 2016.
Ricardo A. Podestá, New identities for binary Krawtchouk polynomials, binomial coefficients and Catalan numbers, arXiv preprint arXiv:1603.09156 [math.CO], 2016.
Ronald Read, An Introduction to Chromatic Polynomials, Journal of Combinatorial Theory, 4(1968)52-71.
Wikipedia, Stirling numbers and exponential generating functions.

See especially A008275 which is the main entry for this triangle. A132393 is an unsigned version, and A008276 is another version.
Cf. A008277, A039814-A039817, A048993, A084938.
A000142(n) = Sum_{k=0..n} |s(n, k)| for n >= 0.
Row sums give A019590(n+1).
Cf. A002209, A027641/A027642, A216294, A271705.

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

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

Table[StirlingS1[n, m], {n, 0, 9}, {m, 0, n}] (* Peter Luschny, Dec 30 2010 *)

create_list(stirling1(n,k),n,0,12,k,0,n); /* Emanuele Munarini, Mar 11 2011 */

a(n,k) = if(k<0 || k>n,0, if(n==0,1,(n-1)*a(n-1,k)+a(n-1,k-1)))

trg(nn)=for (n=0, nn-1, for (k=0, n, print1(stirling(n,k,1), ", ");); print();); \\ Michel Marcus, Jan 19 2015

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

Offset corrected by R. J. Mathar, Feb 23 2009

Formula corrected by Philippe Deléham, Sep 10 2009

A062147 Row sums of unsigned triangle A062137 (generalized a=3 Laguerre).

Original entry on oeis.org

1, 5, 31, 229, 1961, 19081, 207775, 2501801, 32989969, 472630861, 7307593151, 121247816845, 2148321709561, 40476722545169, 807927483311551, 17028146983530961, 377844723929464865, 8803698102396787861, 214877019857456672479, 5482159931449737760181

Offset: 0

Wolfdieter Lang, Jun 19 2001

Vincenzo Librandi, Table of n, a(n) for n = 0..200
Luis Verde-Star, A Matrix Approach to Generalized Delannoy and Schröder Arrays, J. Int. Seq., Vol. 24 (2021), Article 21.4.1.
Index entries for sequences related to Laguerre polynomials

Cf. A062137, A216294.

[Factorial(n)*(&+[Binomial(n+3,n-m)/Factorial(m): m in [0..n]]): n in [0..30]]; // G. C. Greubel, Feb 06 2018

A062147 := n -> n!*simplify(LaguerreL(n,3,-1), 'LaguerreL');
seq(A062147(n), n = 0 .. 30); # G. C. Greubel, Mar 10 2021

Table[Sum[n!*Binomial[n+3,n-k]/k!,{k,0,n}],{n,0,20}]
(* or *)
Table[n!*SeriesCoefficient[E^(x/(1-x))/(1-x)^4,{x,0,n}],{n,0,20}] (* Vaclav Kotesovec, Oct 11 2012 *)

my(x='x+O('x^66)); Vec(serlaplace(exp(x/(1-x))/(1-x)^4)) \\ Joerg Arndt, May 06 2013

a(n) = vecsum(apply(abs,Vec(n!*pollaguerre(n, 3)))); \\ Michel Marcus, Feb 06 2021

[factorial(n)*gen_laguerre(n, 3, -1) for n in (0..30)] # G. C. Greubel, Mar 10 2021

E.g.f.: exp(x/(1-x))/(1-x)^4.
a(n) = Sum_{m=0..n} n!*binomial(n+3, n-m)/m!.
a(n) = (2*n+3)*a(n-1) - (n-1)*(n+2)*a(n-2). - Vaclav Kotesovec, Oct 11 2012
a(n) ~ exp(2*sqrt(n)-n-1/2)*n^(n+7/4)/sqrt(2). - Vaclav Kotesovec, Oct 11 2012
a(n) = n!*LaguerreL(n, 3, -1). - G. C. Greubel, Mar 10 2021

A062192 Row sums of unsigned triangle A062138 (generalized a=5 Laguerre).

Original entry on oeis.org

1, 7, 57, 529, 5509, 63591, 805597, 11109337, 165625929, 2654025319, 45481765921, 829903882017, 16062421776397, 328634683136839, 7086337847838789, 160604998959958441, 3816483607166825617

Offset: 0

Wolfdieter Lang, Jun 19 2001

Seiichi Manyama, Table of n, a(n) for n = 0..441 (terms 0..100 from Harry J. Smith)
Index entries for sequences related to Laguerre polynomials

Cf. A062138, A216294.

[Factorial(n)*(&+[Binomial(n+5,n-m)/Factorial(m): m in [0..n]]): n in [0..30]]; // G. C. Greubel, Feb 06 2018

A062192:= n -> n!*simplify(LaguerreL(n,5,-1), 'LaguerreL');
seq(A062192(n), n = 0 .. 30); # G. C. Greubel, Mar 10 2021

Table[Sum[n! Binomial[n+5,n-m]/m!,{m,0,n}],{n,0,20}] (* Harvey P. Dale, Apr 11 2012 *)

{ f=1; for (n=0, 100, if (n>1, f*=n); a=f*binomial(n+5, n); g=1; a+=sum(m=1, n, f*binomial(n+5, n-m)/g*=m); write("b062192.txt", n, " ", a) ) } \\ Harry J. Smith, Aug 02 2009

my(x='x+O('x^30)); Vec(serlaplace(exp(x/(1-x))/(1-x)^6)) \\ G. C. Greubel, Feb 06 2018

a(n) = vecsum(apply(abs, Vec(n!*pollaguerre(n, 5)))); \\ Michel Marcus, Feb 06 2021

[factorial(n)*gen_laguerre(n, 5, -1) for n in (0..30)] # G. C. Greubel, Mar 10 2021

E.g.f.: exp(x/(1-x))/(1-x)^6.
a(n) = Sum_{m=0..n} n!*binomial(n+5, n-m)/m!.
a(n) = (2*n+5)*a(n-1) - (n-1)*(n+4)*a(n-2). - Vaclav Kotesovec, Oct 11 2012
a(n) ~ exp(2*sqrt(n)-n-1/2)*n^(n+11/4)/sqrt(2). - Vaclav Kotesovec, Oct 11 2012
a(n) = n!*LaguerreL(n, 5, -1). - G. C. Greubel, Mar 10 2021

A062266 Row sums of unsigned triangle A062140 (generalized a=4 Laguerre).

Original entry on oeis.org

1, 6, 43, 358, 3393, 36046, 424051, 5470158, 76751233, 1163391958, 18941512731, 329604456886, 6103575192193, 119823200043678, 2485452283923043, 54309931242376606, 1246803623807490561, 29999359707124127398, 754865494585690965643, 19824604328577866107398

Offset: 0

Wolfdieter Lang, Jun 19 2001

Vincenzo Librandi, Table of n, a(n) for n = 0..200
Luis Verde-Star, A Matrix Approach to Generalized Delannoy and Schröder Arrays, J. Int. Seq., Vol. 24 (2021), Article 21.4.1.
Index entries for sequences related to Laguerre polynomials

Cf. A062140, A216294.

[(&+[Factorial(n)*Binomial(n+4,n-m)/Factorial(m): m in [0..n]]): n in [0..20]]; // G. C. Greubel, Feb 06 2018

A062266:= n -> n!*simplify(LaguerreL(n,4,-1), 'LaguerreL');
seq(A062266(n), n = 0 .. 30); # G. C. Greubel, Mar 10 2021

Table[n!*SeriesCoefficient[E^(x/(1-x))/(1-x)^5,{x,0,n}],{n,0,20}] (* Vaclav Kotesovec, Oct 11 2012 *)
With[{nn=20},CoefficientList[Series[Exp[x/(1-x)]/(1-x)^5,{x,0,nn}],x] Range[0,nn]!] (* Harvey P. Dale, Jul 11 2025 *)

my(x='x+O('x^66)); Vec(serlaplace(exp(x/(1-x))/(1-x)^5)) \\ Joerg Arndt, May 06 2013

a(n) = vecsum(apply(abs,Vec(n!*pollaguerre(n, 4)))); \\ Michel Marcus, Feb 06 2021

[factorial(n)*gen_laguerre(n, 4, -1) for n in (0..30)] # G. C. Greubel, Mar 10 2021

E.g.f.: exp(x/(1-x))/(1-x)^5.
a(n) = Sum_{m=0..n} n!*binomial(n+4, n-m)/m!.
a(n) = 2*(n+2)*a(n-1) - (n-1)*(n+3)*a(n-2). - Vaclav Kotesovec, Oct 11 2012
a(n) ~ exp(2*sqrt(n)-n-1/2)*n^(n+9/4)/sqrt(2). - Vaclav Kotesovec, Oct 11 2012
a(n) = n!*LaguerreL(n, 4, -1). - G. C. Greubel, Mar 10 2021

cp's OEIS Frontend

A048994 Triangle of Stirling numbers of first kind, s(n,k), n >= 0, 0 <= k <= n.

Views

Author

Keywords

Comments

Examples

References

Links

Crossrefs

Programs

Haskell

Maple

Mathematica

Maxima

PARI

PARI

Formula

Extensions

A062147 Row sums of unsigned triangle A062137 (generalized a=3 Laguerre).

Views

Author

Keywords

Links

Crossrefs

Programs

Magma

Maple

Mathematica

PARI

PARI

Sage

Formula

A062192 Row sums of unsigned triangle A062138 (generalized a=5 Laguerre).

Views

Author

Keywords

Links

Crossrefs

Programs

Magma

Maple

Mathematica

PARI

PARI

PARI

Sage

Formula

A062266 Row sums of unsigned triangle A062140 (generalized a=4 Laguerre).

Views

Author

Keywords

Links

Crossrefs

Programs

Magma

Maple

Mathematica

PARI

PARI

Sage

Formula

A216313 Total number of cycles in all partial permutations of {1,2,...,n}.

Views

Author

Keywords

Crossrefs

Programs

Mathematica

Formula