A112486 -id:A112486 - cp's OEIS frontend

A001705 Generalized Stirling numbers: a(n) = n! * Sum_{k=0..n-1} (k+1)/(n-k).

0, 1, 5, 26, 154, 1044, 8028, 69264, 663696, 6999840, 80627040, 1007441280, 13575738240, 196287356160, 3031488633600, 49811492505600, 867718162483200, 15974614352793600, 309920046408806400, 6320046028584960000, 135153868608460800000, 3024476051557847040000

Offset: 0

N. J. A. Sloane

a(n) is also the sum of the positions of the right-to-left minima in all permutations of [n]. Example: a(3)=26 because the positions of the right-to-left minima in the permutations 123,132,213,231,312 and 321 are 123, 13, 23, 3, 23 and 3, respectively and 1 + 2 + 3 + 1 + 3 + 2 + 3 + 3 + 2 + 3 + 3 = 26. - Emeric Deutsch, Sep 22 2008
The asymptotic expansion of the higher order exponential integral E(x,m=2,n=2) ~ exp(-x)/x^2*(1 - 5/x + 26/x^2 - 154/x^3 + 1044/x^4 - 8028/x^5 + 69264/x^6 - ...) leads to the sequence given above. See A163931 and A028421 for more information. - Johannes W. Meijer, Oct 20 2009
a(n) is the total number of cycles (excluding fixed points) in all permutations of [n+1]. - Olivier Gérard, Oct 23 2012; Dec 31 2012
A length n sequence is formed by randomly selecting (one-by-one) n real numbers in (0,1).  a(n)/(n+1)! is the expected value of the sum of the new maximums in such a sequence.  For example for n=3: If we select (in this order): 0.591996, 0.646474, 0.163659  we would add 0.591996 + 0.646474 which would be a bit above the average of a(3)/4! = 26/24. - Geoffrey Critzer, Oct 17 2013

			(1-x)^-2 * (-log(1-x)) = x + 5/2*x^2 + 13/3*x^3 + 77/12*x^4 + ...
Examples: a(6) = 6!*(1/6 + 2/5 + 3/4 + 4/3 + 5/2 + 6/1) = 8028; a(20) = 20!*(1/20 + 2/19 + 3/18 + 4/17 + 5/16 + ... + 16/5 + 17/4 + 18/3 + 19/2 + 20/1) = 135153868608460800000. - _Alexander Adamchuk_, Oct 09 2004
From _Olivier Gérard_, Dec 31 2012: (Start)
The cycle decomposition of all permutations of 4 elements gives the following list: {{{1},{2},{3},{4}}, {{1},{2},{3,4}}, {{1},{2,3},{4}}, {{1},{2,4,3}}, {{1},{2,3,4}}, {{1},{2,4},{3}}, {{1,2},{3},{4}}, {{1,2},{3,4}}, {{1,3,2},{4}},{{1,4,3,2}}, {{1,3,4,2}}, {{1,4,2},{3}}, {{1,2,3},{4}}, {{1,2,4,3}},{{1,3},{2},{4}}, {{1,4,3},{2}}, {{1,3},{2,4}}, {{1,4,2,3}}, {{1,2,3,4}}, {{1,2,4},{3}}, {{1,3,4},{2}}, {{1,4},{2},{3}}, {{1,3,2,4}}, {{1,4},{2,3}}}.
Deleting the fixed points gives the following 26 items: {{3,4}, {2,3}, {2,4,3}, {2,3,4}, {2,4}, {1,2}, {1,2}, {3,4}, {1,3,2}, {1,4,3,2}, {1,3,4,2}, {1,4,2}, {1,2,3}, {1,2,4,3}, {1,3}, {1,4,3}, {1,3}, {2,4}, {1,4,2,3}, {1,2,3,4}, {1,2,4}, {1,3,4}, {1,4}, {1,3,2,4}, {1,4}, {2,3}}. (End)

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).

G. C. Greubel, Table of n, a(n) for n = 0..445 (terms 0 through 100 were provided by T. D. Noe)
J.-L. Baril and S. Kirgizov, The pure descent statistic on permutations, preprint, 2016.
J.-L. Baril and S. Kirgizov, The pure descent statistic on permutations, Discrete Mathematics, 340(10) (2017), 2550-2558.
Jean-Luc Baril and Sergey Kirgizov, Transformation à la Foata for special kinds of descents and excedances, arXiv:2101.01928 [math.CO], 2021.
Gang Chen, Henrik Johansson, Fei Teng, and Tianheng Wang, Next-to-MHV Yang-Mills kinematic algebra, arXiv:2104.12726 [hep-th], 2021. See p. 46.
Loïc Foissy, Cointeraction on noncrossing partitions and related polynomial invariants, arXiv:2501.18212 [math.CO], 2025. See p. 24.
INRIA Algorithms Project, Encyclopedia of Combinatorial Structures 406.
D. S. Mitrinovic and M. S. Mitrinovic, Tableaux d'une classe de nombres reliés aux nombres de Stirling, Univ. Beograd. Publ. Elektrotehn. Fak. Ser. Mat. Fiz. No. 77 (1962), 1-77.
Robert E. Moritz, On the sum of products of n consecutive integers, Univ. Washington Publications in Math., 1 (No. 3, 1926), 44-49. [Annotated scanned copy]
Motohico Mulase, In Search of a Hidden Curve, arXiv:2501.00716 [math.QA], 2025. See p. 27.
J. Riordan, Letter of 04/11/74.
J. Riordan, Letter, Jul 06 1978.

Cf. A000254 (total number of cycles in permutations, including fixed points).
Cf. A002104 (number of different cycles in permutations, without fixed points).
Cf. A006231 (number of different cycles in permutations, including fixed points).
Cf. A001008, A002805, A006675, A066667, A132159, A143947.
Related to n!*the k-th successive summation of the harmonic numbers:
  (k=0) A000254, (k=1) A001705, (k=2) A001711, (k=3) A001716,
  (k=4) A001721, (k=5) A051524, (k=6) A051545, (k=7) A051560,
  (k=8) A051562, (k=9) A051564.

a := n-> add((n+1)!/k, k=2..n+1): seq(a(n), n=0..21); # Zerinvary Lajos, Jan 22 2008; edited Johannes W. Meijer, Nov 28 2012
a := n -> ((n+1)!*(h(n+1)-1)): h := n-> harmonic(n): seq(a(n), n=0..21); # Gary Detlefs, Dec 18 2009; corrected by Johannes W. Meijer, Nov 28 2012

Table[n!*Sum[Sum[1/k,{k,1,m}], {m,1,n}], {n,0,20}] (* Alexander Adamchuk, Apr 14 2006 *)
a[n_] := (n + 1)! (EulerGamma - 1 + PolyGamma[n + 2]);
Table[a[n], {n, 0, 21}] (* Peter Luschny, Feb 19 2022 *)

a(n):=n!*sum(((-1)^(k+1)*binomial(n+1,k+1))/k,k,1,n); /* Vladimir Kruchinin, Oct 10 2016 */

for(n=0,25, print1(n!*sum(k=0,n-1,(k+1)/(n-k)), ", ")) \\ G. C. Greubel, Jan 20 2017

from math import factorial
def A001705(n):
    f = factorial(n)
    return sum(f*(k+1)//(n-k) for k in range(n)) # Chai Wah Wu, Jun 23 2022

Partial sum of first n harmonic numbers multiplied by n!.
a(n) = n!*Sum_{m=1..n} Sum_{k=1..m} 1/k = n!*Sum_{m=1..n} H(m), where H(m) = Sum_{k=1..m} 1/k = A001008(m)/A002805(m) is m-th Harmonic number.
E.g.f.: - log (1 - x) / (1 - x)^2.
a(n) = (n+1)! * H(n) - n*n!, H(n) = Sum_{k=1..n} (1/k).
a(n) = A112486(n, 1).
a(n) = a(n-1)*(n+1) + n! = A000254(n+1) - A000142(n+1) = A067176(n+1, 1). - Henry Bottomley, Jan 09 2002
a(n) = Sum_{k=0..n-1} ((-1)^(n-1+k) * (k+1) * 2^k * Stirling1(n, k+1)). - Borislav Crstici (bcrstici(AT)etv.utt.ro), Jan 26 2004
With alternating signs: Ramanujan polynomials psi_2(n, x) evaluated at 0. - Ralf Stephan, Apr 16 2004
a(n) = Sum_{k=1..n} (k*StirlingCycle(n+1,k+1)). - David Callan, Sep 25 2006
a(n) = Sum_{k=n..n*(n+1)/2} k*A143947(n,k). - Emeric Deutsch, Sep 22 2008
For n >= 1, a(n) = Sum_{j=0..n-1} ((-1)^(n-j-1) * 2^j * (j+1) * Stirling1(n,j+1)). - Milan Janjic, Dec 14 2008
a(n) = (2*n+1)*a(n-1) - n^2*a(n-2). - Gary Detlefs, Nov 27 2009
a(n) = (n+1)!*(H(n+1) - 1) where H(n) is the n-th harmonic number. - Gary Detlefs, Dec 18 2009
a(n) = n!*Sum_{k=1..n} (-1)^(k+1)*binomial(n+1,k+1)/k. - Vladimir Kruchinin, Oct 10 2016
a(n) = (n+1)!*Sum_{k = 1..n} (-1)^(k+1)*binomial(n+1,k+1)*k/(k+1). - Peter Bala, Feb 15 2022
a(n) = Gamma(n + 2) * (Digamma(n + 2) + EulerGamma - 1). - Peter Luschny, Feb 19 2022
From Mélika Tebni, Jun 22 2022: (Start)
a(n) = -Sum_{k=0..n} k!*A066667(n, k+1).
a(n) = Sum_{k=0..n} k!*A132159(n, k+1). (End)
a(n) = n*(n + 1)!*hypergeom([1, 1, 1 - n], [2, 3], 1)/2. - Peter Luschny, Jun 22 2022

More terms from Sascha Kurz, Mar 22 2002

A008276 Triangle of Stirling numbers of first kind, s(n, n-k+1), n >= 1, 1 <= k <= n. Also triangle T(n,k) giving coefficients in expansion of n!*binomial(x,n)/x in powers of x.

Original entry on oeis.org

1, 1, -1, 1, -3, 2, 1, -6, 11, -6, 1, -10, 35, -50, 24, 1, -15, 85, -225, 274, -120, 1, -21, 175, -735, 1624, -1764, 720, 1, -28, 322, -1960, 6769, -13132, 13068, -5040, 1, -36, 546, -4536, 22449, -67284, 118124, -109584, 40320, 1, -45

Offset: 1

N. J. A. Sloane

n-th row of the triangle = charpoly of an (n-1) X (n-1) matrix with (1,2,3,...) in the diagonal and the rest zeros. - Gary W. Adamson, Mar 19 2009
From Daniel Forgues, Jan 16 2016: (Start)
For n >= 1, 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)
The moment generating function of the probability density function p(x, m=q, n=1, mu=q) = q^q*x^(q-1)*E(x, q, 1)/(q-1)!, with q >= 1, is M(a, m=q, n=1, mu=q) = Sum_{k=0..q}(A000312(q) / A000142(q-1)) * A008276(q, k) * polylog(k, a) / a^q , see A163931 and A274181. - Johannes W. Meijer, Jun 17 2016
Triangle of coefficients of the polynomial x(x-1)(x-2)...(x-n+1), also denoted as falling factorial (x)n, expanded into decreasing powers of x. - _Ralf Stephan, Dec 11 2016

			3!*binomial(x,3) = x*(x-1)*(x-2) = x^3 - 3*x^2 + 2*x.
Triangle begins
  1;
  1,  -1;
  1,  -3,   2;
  1,  -6,  11,   -6;
  1, -10,  35,  -50,  24;
  1, -15,  85, -225, 274, -120;
...

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.
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, 2nd ed. (Addison-Wesley, 1994), p. 257.

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].
T. Copeland, Generators, Inversion, and Matrix, Binomial, and Integral Transforms
Yahia Djemmada, Abdelghani Mehdaoui, László Németh, and László Szalay, The Fibonacci-Fubini and Lucas-Fubini numbers, arXiv:2407.04409 [math.CO], 2024. See pp. 10, 12.
Bill Gosper, Colored illustrations of triangle of Stirling numbers of first kind read mod 2, 3, 4, 5, 6, 7
Eric Weisstein's World of Mathematics, Stirling Number of the First Kind
Wikipedia, Stirling numbers and exponential generating functions

See A008275 and A048994, which are the main entries for this triangle of numbers.
See A008277 triangle of Stirling numbers of the second kind, S2(n,k).
Cf. A054654, A054655, A084938, A145324, A094216, A003422, A000166, A000110, A000204, A000045, A000108.

a008276 n k = a008276_tabl !! (n-1) !! (k-1)
a008276_row n = a008276_tabl !! (n-1)
a008276_tabl = map init $ tail a054654_tabl
-- Reinhard Zumkeller, Mar 18 2014

seq(seq(coeff(expand(n!*binomial(x,n)),x,j),j=n..1,-1),n=1..15); # Robert Israel, Jan 24 2016
A008276 := proc(n, k): combinat[stirling1](n, n-k+1) end: seq(seq(A008276(n, k), k=1..n), n=1..9); # Johannes W. Meijer, Jun 17 2016

len = 47; m = Ceiling[Sqrt[2*len]]; t[n_, k_] = StirlingS1[n, n-k+1]; Flatten[Table[t[n, k], {n, 1, m}, {k, 1, n}]][[1 ;; len]] (* Jean-François Alcover, May 31 2011 *)
Flatten@Table[CoefficientList[Product[1-k x, {k, 1, n}], x], {n, 0, 8}] (* Oliver Seipel, Jun 14 2024 *)
Flatten@Table[Coefficient[Product[x-k, {k, 0, n-1}], x, Reverse@Range[n]], {n, Range[9]}] (* Oliver Seipel, Jun 14 2024, after  Ralf Stephan *)

T(n,k)=if(n<1,0,n!*polcoeff(binomial(x,n),n-k+1))

T(n,k)=if(n<1,0,n!*polcoeff(polcoeff(y*(1+y*x+x*O(x^n))^(1/y),n),k))

def T(n,k): return falling_factorial(x,n).expand().coefficient(x,n-k+1) # Ralf Stephan, Dec 11 2016

n!*binomial(x, n) = Sum_{k=1..n-1} T(n, k)*x^(n-k).
|A008276(n, k)| = T(n-1, k-1) where T(n, k) is the triangle, read by rows, given by [1, 0, 1, 0, 1, 0, 1, 0, 1, ...] DELTA [1, 1, 2, 2, 3, 3, 4, 4, 5, 5, ...]; A008276(n, k) = T(n-1, k-1) where T(n, k) is the triangle, read by rows, given by [1, 0, 1, 0, 1, 0, 1, 0, 1, ...] DELTA [ -1, -1, -2, -2, -3, -3, -4, -4, -5, -5, ...]. Here DELTA is the operator defined in A084938. - Philippe Deléham, Dec 30 2003
|T(n, k)| = Sum_{m=0..n} A008517(k, m+1)*binomial(n+m, 2*(k-1)), n >= k >= 1. A008517 is the second-order Eulerian triangle. See the Graham et al. reference p. 257, eq. (6.44).
A094638 formula for unsigned T(n, k).
|T(n, k)| = Sum_{m=0..min(k-1, n-k)} A112486(k-1, m)*binomial(n-1, k-1+m) if n >= k >= 1, else 0. - Wolfdieter Lang, Sep 12 2005, see A112486.
|T(n, k)| = (f(n-1, k-1)/(2*(k-1))!)* Sum_{m=0..min(k-1, n-k)} A112486(k-1, m)*f(2*(k-1), k-1-m)*f(n-k, m) if n >= k >= 1, else 0, where f(n, k) stands for the falling factorial n*(n-1)*...*(n-(k-1)) and f(n, 0):=1. - Wolfdieter Lang, Sep 12 2005, see A112486.
With P(n,t) = Sum_{k=0..n-1} T(n,k+1) * t^k = (1-t)*(1-2*t)*...*(1-(n-1)t) and P(0,t) = 1, exp(P(.,t)*x) = (1+t*x)^(1/t) . Compare A094638. T(n,k+1) = (1/k!) (D_t)^k (D_x)^n ( (1+t*x)^(1/t) - 1 ) evaluated at t=x=0 . - Tom Copeland, Dec 09 2007
Product_{i=1..n} (x-i) = Sum_{k=0..n} T(n,k)*x^k. - Reinhard Zumkeller, Dec 29 2007
E.g.f.: Sum_{n>=0} (Sum_{k=0..n} T(n,n-k)*t^k)/n!) = Sum_{n>=0} (x)n * t^k/n! = exp(x * log(1+t)), with (x)_n the n-th falling factorial polynomial. - _Ralf Stephan, Dec 11 2016
Sum_{j=0..m} T(m, m-j)*s2(j+k+1, m) =  m^k, where s2(j, m) are Stirling numbers of the second kind. - Tony Foster III, Jul 25 2019
For n>=2, Sum_{k=1..n} k*T(n,k) = (-1)^(n-1)*(n-2)!. - Zizheng Fang, Dec 27 2020

A001303 Stirling numbers of first kind, s(n+3, n), negated.

Original entry on oeis.org

6, 50, 225, 735, 1960, 4536, 9450, 18150, 32670, 55770, 91091, 143325, 218400, 323680, 468180, 662796, 920550, 1256850, 1689765, 2240315, 2932776, 3795000, 4858750, 6160050, 7739550, 9642906, 11921175, 14631225, 17836160, 21605760, 26016936, 31154200

Offset: 1

N. J. A. Sloane

a(n) is equal to the sum of the products of each distinct grouping of 3 members of the set {1, 2, 3, ..., n + 2} (a(1) = 1*2*3, a(2) = 1*2*3 + 1*2*4 + 1*3*4 + 2*3*4, a(3) = 1*2*3 + 1*2*4 + 1*2*5 + 1*3*4 + 1*3*5 + 1*4*5 + 2*3*4 + 2*3*5 + 2*4*5 + 3*4*5). - Jeffreylee R. Snow, Sep 23 2013

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.
Louis Comtet, Advanced Combinatorics, Reidel, 1974, p. 227, #16.
F. N. David, M. G. Kendall and D. E. Barton, Symmetric Function and Allied Tables, Cambridge, 1966, p. 226.
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).

T. D. Noe, Table of n, a(n) for n = 1..1000
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].
Karl Dienger, Beiträge zur Lehre von den arithmetischen und geometrischen Reihen höherer Ordnung, Jahres-Bericht Ludwig-Wilhelm-Gymnasium Rastatt, Rastatt, 1910. [Annotated scanned copy]
G. C. Greubel, A Note on Jain basis functions, arXiv:1612.09385 [math.CA], 2016.
Robert E. Moritz, On the sum of products of n consecutive integers, Univ. Washington Publications in Math., Vol. 1, No. 3 (1926), pp. 44-49. [Annotated scanned copy]
Simon Plouffe, Approximations de séries génératrices et quelques conjectures, Dissertation, Université du Québec à Montréal, 1992, arXiv:0911.4975 [math.NT], 2009.
Simon Plouffe, 1031 Generating Functions, Appendix to Thesis, Montreal, 1992.
Index entries for linear recurrences with constant coefficients, signature (7,-21,35,-35,21,-7,1).

Cf. A000217, A001297, A008275.

seq(numbperm (n,2)*numbperm (n,4)/48, n=4..33); # Zerinvary Lajos, Apr 26 2007
seq(15*binomial(n+2,6)-10*binomial(n+1,5)+binomial(n,4),n=4..30); # Miklos Kristof, Nov 04 2007
A001303 := proc(n)
    -combinat[stirling1](n+3,n) ;
end proc: # R. J. Mathar, May 19 2016

Table[-StirlingS1[n + 3, n], {n, 100}] (* T. D. Noe, Jun 27 2012 *)
a[ n_] := n (n + 1) (n + 2)^2 (n + 3)^2 / 48; (* Michael Somos, Sep 04 2017 *)

a(n) = n*(n+1)*(n+2)^2*(n+3)^2/48; \\ Altug Alkan, Aug 29 2017

[stirling_number1(n,n-3) for n in range(4, 34)] # Zerinvary Lajos, May 16 2009

a(n) = binomial(n+3, 4)*binomial(n+3, 2).
G.f.: x*(6 + 8*x + x^2)/(1 - x)^7. - Simon Plouffe in his 1992 dissertation
E.g.f. with offset 3: exp(x)*(6*(x^3)/3! + 26*(x^4)/4! + 35*(x^5)/5! + 15*(x^6)/6!). See row k=3 of A112486 for the coefficients [6, 26, 35, 15].
a(n) = (f(n+2, 3)/6!)*Sum_{m=0..min(3, n)} A112486(3,m)*f(6, 3-m)*f(n-1, m), with the falling factorials notation f(n, m):=n*(n-1)*...*(n-(m-1)).
From Jason Lang, Oct 03 2006: (Start)
a(n) = A000217(n) * n! / ( 4! * (n-4)! ) [for n > 4 and A000217 = the triangular numbers];
a(n) = ((n+4)! / n! ) ^2 / ( (n+2) * (n+1) * 2*4!);
a(n) = (n-0)^2 * (n-1)^2 * (n-2) * (n-3) / (2*4!). (End)
From Miklos Kristof, Nov 04 2007: (Start)
a(n) = 15*binomial(n+5,6) - 10*binomial(n+4,5) + binomial(n+3,4).
E.g.f. with offset 4: exp(x)*((1/4)*x^4 + (1/6)*x^5 + (1/48)*x^6). (End)
a(n) = n*(n+1)(n+2)^2*(n+3)^2/48. - Jeremy Galvagni, Mar 03 2009
From Gary Detlefs, Jun 06 2010: (Start)
a(n) = (n+3)^2/(n^2-1)*a(n-1), n > 1;
a(n) = 6*Product_{k=2..n} (k+3)^2/(k^2 - 1). (End)
a(n) = A001297(-3-n) for all n in Z. - Michael Somos, Sep 04 2017
From Amiram Eldar, Jan 10 2022: (Start)
Sum_{n>=1} 1/a(n) = 16*Pi^2/3 - 472/9.
Sum_{n>=1} (-1)^(n+1)/a(n) = 4*Pi^2/3 + 16/9 - 64*log(2)/3. (End)

More terms from Klaus Strassburger (strass(AT)ddfi.uni-duesseldorf.de), Jan 17 2000

Notation of the polynomial formula edited by R. J. Mathar, Sep 15 2009

A000915 Stirling numbers of first kind s(n+4, n).

Original entry on oeis.org

24, 274, 1624, 6769, 22449, 63273, 157773, 357423, 749463, 1474473, 2749747, 4899622, 8394022, 13896582, 22323822, 34916946, 53327946, 79721796, 116896626, 168423871, 238810495, 333685495, 460012995, 626334345, 843041745, 1122686019, 1480321269, 1933889244

Offset: 1

N. J. A. Sloane

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, p. 227, #16.
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, 2nd. ed., 1994, p. 259.
J. Riordan, An Introduction to Combinatorial Analysis, Wiley, 1958, p. 48.
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).

T. D. Noe, Table of n, a(n) for n = 1..1000
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].
Index entries for linear recurrences with constant coefficients, signature (9,-36,84,-126,126,-84,36,-9,1).

Cf. A008275, A094216, A001303 for s(n+3,n), A053567 for s(n+5,n).

Cf. A001298.

A000915 := proc(n)
    combinat[stirling1](n+4,n) ;
end proc:
seq(A000915(n),n=1..10) ; # R. J. Mathar, May 19 2016

Table[Binomial[n + 4, 5]*(15*n^3 + 150*n^2 + 485*n + 502)/48, {n, 50}] (* T. D. Noe, Jun 20 2012 *)
a[ n_] := n (n + 1) (n + 2) (n + 3) (n + 4) (15 n^3 + 150 n^2 + 485 n + 502) / 5760; (* Michael Somos, Sep 04 2017 *)

{a(n) = n * (n+1) * (n+2) * (n+3) * (n+4) * (15*n^3+ 150*n^2 + 485*n + 502) / 5760}; /* Michael Somos, Sep 04 2017 */

[stirling_number1(n,n-4) for n in range(5, 30)] # Zerinvary Lajos, May 16 2009

a(n) = binomial(n+4, 5)*(15*n^3 + 150*n^2 + 485*n + 502)/48. - André F. Labossière, Sep 30 2004
Stirling1(n+1, n-3) = Sum_{L=1..n} (Sum_{k=L+1..n} (Sum_{j=k+1..n} (Sum_{i=j+1..n} i*j*k*L))), cf. A001298. - Vladeta Jovovic, Jan 31 2005
E.g.f. with offset 4: exp(x)*(Sum_{m=0..4} A112486(4,m)*(x^(4+m))/(4+m)!).
a(n) = (f(n+3, 4)/8!)*Sum_{m=0..min(4, n-1)} A112486(4,m)*f(8, 4-m)*f(n-1, m), with the falling factorials f(n, m):=n*(n-1)*...*(n-(m-1)).
G.f.: x*(24 + 58*x + 22*x^2 + x^3)/(1 - x)^9, see the k=3 row of triangle A112007 for [24, 58, 22, 1].
a(n) = A001298(-4-n) for all n in Z. - Michael Somos, Sep 04 2017

More terms from Klaus Strassburger (strass(AT)ddfi.uni-duesseldorf.de), Jan 17 2000

A269939 Triangle read by rows, Ward numbers T(n, k) = Sum_{m=0..k} (-1)^(m + k) * binomial(n + k, n + m) * Stirling2(n + m, m), for n >= 0, 0 <= k <= n.

Original entry on oeis.org

1, 0, 1, 0, 1, 3, 0, 1, 10, 15, 0, 1, 25, 105, 105, 0, 1, 56, 490, 1260, 945, 0, 1, 119, 1918, 9450, 17325, 10395, 0, 1, 246, 6825, 56980, 190575, 270270, 135135, 0, 1, 501, 22935, 302995, 1636635, 4099095, 4729725, 2027025

Offset: 0

Peter Luschny, Mar 26 2016

We propose to call this sequence the 'Ward set numbers' and sequence A269940 the 'Ward cycle numbers'. - Peter Luschny, Nov 25 2022

			Triangle starts:
  1;
  0, 1;
  0, 1,   3;
  0, 1,  10,   15;
  0, 1,  25,  105,   105;
  0, 1,  56,  490,  1260,    945;
  0, 1, 119, 1918,  9450,  17325,  10395;
  0, 1, 246, 6825, 56980, 190575, 270270, 135135;

Bishal Deb and Alan D. Sokal, Higher-order Stirling cycle and subset triangles: Total positivity, continued fractions and real-rootedness, arXiv:2507.18959 [math.CO], 2025. See p. 6.
Peter Luschny, The P-transform.
Andrew Elvey Price and Alan D. Sokal, Phylogenetic trees, augmented perfect matchings, and a Thron-type continued fraction (T-fraction) for the Ward polynomials, arXiv:2001.01468 [math.CO], 2020.
Marko Riedel, Math Stackexchange, Upper Stirling numbers and Ward numbers.
Aleks Žigon Tankosič, Recurrence Relations for Some Integer Sequences Related to Ward Numbers, arXiv:2508.04754 [math.CO], 2025. See p. 2.
Elena L. Wang and Guoce Xin, On Ward Numbers and Increasing Schröder Trees, arXiv:2507.15654 [math.CO], 2025. See pp. 1, 12.

Variants: A134991 (main entry for this triangle), A181996.
Row sums are A000311.
Alternating row sums are signed factorials A133942.
Cf. A269940 (Stirling1 counterpart), A268437.

# first version
A269939 := (n,k) -> add((-1)^(m+k)*binomial(n+k,n+m)*Stirling2(n+m, m), m=0..k):
seq(seq(A269939(n,k), k=0..n), n=0..8);
# Alternatively:
T := proc(n,k) option remember;
    `if`(k=0 and n=0, 1,
    `if`(k<=0 or k>n, 0,
    k*T(n-1,k)+(n+k-1)*T(n-1,k-1))) end:
for n from 0 to 6 do seq(T(n,k),k=0..n) od;
# simple, third version
T := (n,k)->  (n+k)!*coeftayl((exp(z)-z-1)^k/k!, z=0, n+k); # Marko Riedel, Apr 14 2016

Table[Sum[(-1)^(m + k) Binomial[n + k, n + m] StirlingS2[n + m, m], {m, 0, k}], {n, 0, 8}, {k, 0, n}] // Flatten (* Michael De Vlieger, Apr 15 2016 *)

T(n) = {[Vecrev(Pol(p)) | p<-Vec(serlaplace(1/((1+y)*(1 + lambertw(-y/(1+y)*exp((x-y)/(1+y) + O(x*x^n)))))))]}
{ my(A=T(8)); for(n=1, #A, print(A[n])) } \\ Andrew Howroyd, Jan 14 2022

T = lambda n,k: sum((-1)^(m+k)*binomial(n+k,n+m)*stirling_number2(n+m,m) for m in (0..k))
for n in (0..6): print([T(n,k) for k in (0..n)])

# uses[PtransMatrix from A269941]
PtransMatrix(8, lambda n: 1/(n+1), lambda n, k: (-1)^k*falling_factorial(n+k,n))

T(n,k) = (-1)^k*FF(n+k,n)*P[n,k](1/(n+1)) where P is the P-transform and FF the falling factorial function. For the definition of the P-transform see the link.
T(n,k) = A268437(n,k)*FF(n+k,n)/(2*n)!.
T(n,k) = (n+k)! [z^{n+k}] (exp(z)-z-1)^k/k!. - Marko Riedel, Apr 14 2016
From Fabián Pereyra, Jan 12 2022: (Start)
T(n,k) = k*T(n-1,k) + (n+k-1)*T(n-1,k-1) for n > 0, T(0,0) = 1, T(n,0) = 0 for n > 0. (See the second Maple program.)
E.g.f.: A(x,t) = 1/((1+t)*(1 + W(-t/(1+t)*exp((x-t)/(1+t))))), where W(x) is the Lambert W-function.
T(n,k) = Sum_{j=0..k} E2(n,j)*binomial(n-j,k-j), where E2(n,k) are the second-order Eulerian numbers A340556.
T(n,k) = Sum_{j=k..n} (-1)^(n-j)*A112486(n,j)*binomial(j,k). (End)

A053567 Stirling numbers of first kind, s(n+5, n).

Original entry on oeis.org

-120, 1764, -13132, 67284, -269325, 902055, -2637558, 6926634, -16669653, 37312275, -78558480, 156952432, -299650806, 549789282, -973941900, 1672280820, -2792167686, 4546047198, -7234669596, 11276842500, -17247104875, 25922927745, -38343278610, 55880640270

Offset: 1

Klaus Strassburger (strass(AT)ddfi.uni-duesseldorf.de), Jan 17 2000

a(n) is equal to (-1)^n times the sum of the products of each distinct grouping of 5 members of the set {1, 2, 3, ..., n + 4}. So, a(1) = (-1)*1*2*3*4*5 = -120, and a(2) = 1*2*3*4*5 + 1*2*3*4*6 + 1*2*3*5*6 + 1*2*4*5*6 + 1*3*4*5*6 + 2*3*4*5*6 = 1764. See comment at A001303. - Greg Dresden, Aug 26 2019

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.

Vincenzo Librandi, Table of n, a(n) for n = 1..200
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].
G. C. Greubel, A Note on Jain basis functions, arXiv:1612.09385 [math.CA], 2016.
Index entries for linear recurrences with constant coefficients, signature (-11,-55,-165,-330,-462,-462,-330,-165,-55,-11,-1).

Next |Stirling1| diagonal A112002, 5th diagonal of A130534.

[(-1)^n*Binomial(n+5, 6)*Binomial(n+5, 2)*(3*n^2+23*n+38)/8: n in [1..30]]; // Vincenzo Librandi, Jun 09 2011

A053567 := proc(n) (-1)^(n+1)*combinat[stirling1](n+5,n) ; end proc: # R. J. Mathar, Jun 08 2011

Table[StirlingS1[n+5,n](-1)^(n-1),{n,30}] (* Harvey P. Dale, Sep 21 2011 *)
(* or *)
CoefficientList[Series[-x*(120 - 444*x + 328*x^2 - 52*x^3 + x^4)/(1+x)^11, {x, 0, 27}], x] (* Georg Fischer, May 19 2019 *)

a(n) = (-1)^(n-1)*stirling(n+5, n, 1); \\ Michel Marcus, Aug 29 2017

[stirling_number1(n,n-5)*(-1)^(n+1) for n in range(6, 26)] # Zerinvary Lajos, May 16 2009

a(n) = (-1)^n*binomial(n+5, 6)*binomial(n+5, 2)*(3*n^2 + 23*n + 38)/8.
G.f.: -x*(120 - 444*x + 328*x^2 - 52*x^3 + x^4)/(1+x)^11. See row k=4 of triangle A112007 for the coefficients. [G.f. corrected by Georg Fischer, May 19 2019]
E.g.f. with offset 5: exp(x)*(Sum_{m=0..5} A112486(5, m)*(x^(5+m)/(5+m)!).
a(n) = (f(n+4, 5)/10!)*Sum_{m=0..min(5, n-1)} A112486(5, m)*f(10, 5-m)*f(n-1, m)), with the falling factorials f(n, m):=n*(n-1)*, ..., *(n-(m-1)). From the e.g.f.

Definition edited by Eric M. Schmidt, Aug 29 2017

Incorrect formula removed by Greg Dresden, Aug 26 2019

cp's OEIS Frontend

A112488 Third column of triangle A112486 used for e.g.f.s of |Stirling1| = |A008275| diagonals.

Views

Author

Keywords

Examples

Links

Programs

Mathematica

Formula

A112489 Fourth column of triangle A112486 used for e.g.f.s of |Stirling1|=|A008275| diagonals.

Views

Author

Keywords

Examples

Formula

A112490 Fifth column of triangle A112486 used for e.g.f.s of |Stirling1|=|A008275| diagonals.

Views

Author

Keywords

Examples

Formula

A112491 Sixth column of triangle A112486 used for e.g.f.s of |Stirling1|=|A008275| diagonals.

Views

Author

Keywords

Examples

Formula

A001705 Generalized Stirling numbers: a(n) = n! * Sum_{k=0..n-1} (k+1)/(n-k).

Views

Author

Keywords

Comments

Examples

References

Links

Crossrefs

Programs

Maple

Mathematica

Maxima

PARI

Python

Formula

Extensions

A008276 Triangle of Stirling numbers of first kind, s(n, n-k+1), n >= 1, 1 <= k <= n. Also triangle T(n,k) giving coefficients in expansion of n!*binomial(x,n)/x in powers of x.

Views

Author

Keywords

Comments

Examples

References

Links

Crossrefs

Programs

Haskell

Maple

Mathematica

PARI

PARI

Sage

Formula

A001303 Stirling numbers of first kind, s(n+3, n), negated.

Views

Author

Keywords

Comments

References

Links

Crossrefs

Programs

Maple

Mathematica

PARI

Sage

Formula

Extensions

A000915 Stirling numbers of first kind s(n+4, n).

Views

Author