A001716 -id:A001716 - cp's OEIS frontend

A165675 Triangle read by rows. T(n, k) = (n - k + 1)! * H(k, n - k), where H are the hyperharmonic numbers. For 0 <= k <= n.

1, 1, 1, 2, 3, 1, 6, 11, 5, 1, 24, 50, 26, 7, 1, 120, 274, 154, 47, 9, 1, 720, 1764, 1044, 342, 74, 11, 1, 5040, 13068, 8028, 2754, 638, 107, 13, 1, 40320, 109584, 69264, 24552, 5944, 1066, 146, 15, 1, 362880, 1026576, 663696, 241128, 60216, 11274, 1650, 191, 17, 1

Offset: 0

Johannes W. Meijer, Oct 05 2009

Previous name: Extended triangle related to the asymptotic expansions of the E(x, m = 2, n).
For the definition of the hyperharmonic numbers see the formula section.
This triangle is the same as triangle A165674 except for the extra left-hand column T(n, 0) = n!. The T(n) formulas for the right-hand columns generate the coefficients of this extra left-hand column.
Leroy Quet discovered triangle A105954 which is the reversal of our triangle.
In square format, row k gives the (n-1)-st elementary symmetric function of {k, k+1, k+2,..., k+n}, as in the Mathematica section. - Clark Kimberling, Dec 29 2011

			Triangle T(n, k) begins:
  [0]    1;
  [1]    1,     1;
  [2]    2,     3,    1;
  [3]    6,    11,    5,    1;
  [4]   24,    50,   26,    7,   1;
  [5]  120,   274,  154,   47,   9,   1;
  [6]  720,  1764, 1044,  342,  74,  11,  1;
  [7] 5040, 13068, 8028, 2754, 638, 107, 13, 1;
Seen as an array (the triangle arises when read by descending antidiagonals):
  [0] 1,  1,   2,    6,    24,    120,     720,     5040, ...
  [1] 1,  3,  11,   50,   274,   1764,   13068,   109584, ...
  [2] 1,  5,  26,  154,  1044,   8028,   69264,   663696, ...
  [3] 1,  7,  47,  342,  2754,  24552,  241128,  2592720, ...
  [4] 1,  9,  74,  638,  5944,  60216,  662640,  7893840, ...
  [5] 1, 11, 107, 1066, 11274, 127860, 1557660, 20355120, ...
  [6] 1, 13, 146, 1650, 19524, 245004, 3272688, 46536624, ...
  [7] 1, 15, 191, 2414, 31594, 434568, 6314664, 97053936, ...

A105954 is the reversal of this triangle.
A165674, A138771 and A165680 are related triangles.
A080663 equals the third right hand column.
A000142 equals the first left hand column.
A093345 are the row sums.
Columns include A165676, A165677, A165678 and A165679.

nmax := 8; for n from 0 to nmax do a(n, 0) := n! od: for n from 0 to nmax do a(n, n) := 1 od: for n from 2 to nmax do for m from 1 to n-1 do a(n, m) := (n-m+1)*a(n-1, m) + a(n-1, m-1) od: od: seq(seq(a(n, m), m=0..n), n=0..nmax);
# Johannes W. Meijer, revised Nov 27 2012
# Shows the array format, using hyperharmonic numbers.
H := proc(n, k) option remember; if n = 0 then 1/(k+1)
else add(H(n - 1, j), j = 0..k) fi end:
seq(lprint(seq((k + 1)!*H(n, k), k = 0..7)), n = 0..7);
# Shows the array format, using the hypergeometric formula.
A := (n, k) -> (k+1)*((n + k)! / n!)*hypergeom([-k, 1, 1], [2, n + 1], 1):
seq(lprint(seq(simplify(A(n, k)), k = 0..7)), n = 0..7);
# Peter Luschny, Jul 03 2022

a[n_] := SymmetricPolynomial[n - 1, t[n]]; z = 10;
t[n_] := Table[k - 1, {k, 1, n}]; t1 = Table[a[n], {n, 1, z}]  (* A000142 *)
t[n_] := Table[k,     {k, 1, n}]; t2 = Table[a[n], {n, 1, z}]  (* A000254 *)
t[n_] := Table[k + 1, {k, 1, n}]; t3 = Table[a[n], {n, 1, z}]  (* A001705 *)
t[n_] := Table[k + 2, {k, 1, n}]; t4 = Table[a[n], {n, 1, z}]  (* A001711 *)
t[n_] := Table[k + 3, {k, 1, n}]; t5 = Table[a[n], {n, 1, z}]  (* A001716 *)
t[n_] := Table[k + 4, {k, 1, n}]; t6 = Table[a[n], {n, 1, z}]  (* A001721 *)
t[n_] := Table[k + 5, {k, 1, n}]; t7 = Table[a[n], {n, 1, z}]  (* A051524 *)
t[n_] := Table[k + 6, {k, 1, n}]; t8 = Table[a[n], {n, 1, z}]  (* A051545 *)
t[n_] := Table[k + 7, {k, 1, n}]; t9 = Table[a[n], {n, 1, z}]  (* A051560 *)
t[n_] := Table[k + 8, {k, 1, n}]; t10 = Table[a[n], {n, 1, z}] (* A051562 *)
t[n_] := Table[k + 9, {k, 1, n}]; t11 = Table[a[n], {n, 1, z}] (* A051564 *)
t[n_] := Table[k + 10, {k, 1, n}];t12 = Table[a[n], {n, 1, z}] (* A203147 *)
t = {t1, t2, t3, t4, t5, t6, t7, t8, t9, t10};
TableForm[t]  (* A165675 in square format *)
m[i_, j_] := t[[i]][[j]];
(* A165675 as a sequence *)
Flatten[Table[m[i, n + 1 - i], {n, 1, 10}, {i, 1, n}]]
(* Clark Kimberling, Dec 29 2011 *)
A[n_, k_] := (k + 1)*((n + k)! / n!)*HypergeometricPFQ[{-k, 1, 1}, {2, n + 1}, 1];
Table[A[n, k], {n, 0, 7}, {k, 0, 7}] // TableForm (* Peter Luschny, Jul 03 2022 *)

from functools import cache
@cache
def Trow(n: int) -> list[int]:
    if n == 0:
        return [1]
    row = Trow(n - 1) + [1]
    for m in range(n - 1, 0, -1):
        row[m] = (n - m + 1) * row[m] + row[m - 1]
    row[0] *= n
    return row
for n in range(9): print(Trow(n))  # Peter Luschny, Feb 27 2025

The hyperharmonic numbers are H(n, k) = Sum_{j=0..k} H(n - 1, j), with base condition H(0, k) = 1/(k + 1).
T(n, k) = (n - k + 1)*T(n - 1, k) + T(n - 1, k - 1), 1 <= k <= n-1, with T(n, 0) = n! and T(n, n) = 1.
From Peter Luschny, Jul 03 2022: (Start)
The rectangular array is given by:
A(n, k) = (k + 1)!*H(n, k).
A(n, k) = (k + 1)*((n + k)! / n!)*hypergeom([-k, 1, 1], [2, n + 1], 1). (End)
From Werner Schulte, Feb 26 2025: (Start)
T(n, k) = n * T(n-1, k) + (n-1)! / (k-1)! for 0 < k < n.
T(n, k) = (Sum_{i=k..n} 1/i) * n! / (k-1)! for 0 < k <= n.
Matrix inverse M = T^(-1) is given by: M(n, n) = 1, M(n, n-1) = 1 - 2 * n for n > 0, M(n, n-2) = (n-1)^2 for n > 1, and M(i, j) = 0 otherwise. (End)

New name from Peter Luschny, Jul 03 2022

A105954 Array read by descending antidiagonals: A(n, k) = (n + 1)! * H(k, n + 1), where H(n, k) is a higher-order harmonic number, H(0, k) = 1/k and H(n, k) = Sum_{j=1..k} H(n-1, j), for 0 <= k <= n.

Original entry on oeis.org

1, 1, 1, 1, 3, 2, 1, 5, 11, 6, 1, 7, 26, 50, 24, 1, 9, 47, 154, 274, 120, 1, 11, 74, 342, 1044, 1764, 720, 1, 13, 107, 638, 2754, 8028, 13068, 5040, 1, 15, 146, 1066, 5944, 24552, 69264, 109584, 40320, 1, 17, 191, 1650, 11274, 60216, 241128, 663696, 1026576, 362880

Offset: 0

Leroy Quet, Jun 26 2005

Antidiagonal sums are A093345 (n! * (1 + Sum_{i=1..n}((1/i)*Sum_{j=0..i-1} 1/j!))). - Gerald McGarvey, Aug 27 2005
A recasting of A093905 and A067176. - R. J. Mathar, Mar 01 2009
The triangular array of this sequence is the reversal of A165675 which is related to the asymptotic expansion of the higher order exponential integral E(x,m=2,n); see also A165674. - Johannes W. Meijer, Oct 16 2009

			A(2, 2) = (1 + (1 + 1/2) + (1 + 1/2 + 1/3))*6 = 26.
Array A(n, k) begins:
  [n\k]  0       1       2        3        4        5          6
  -------------------------------------------------------------------
  [0]    1,      1,      1,       1,       1,       1,         1, ...
  [1]    1,      3,      5,       7,       9,       11,       13, ...
  [2]    2,     11,     26,      47,      74,      107,      146, ...
  [3]    6,     50,    154,     342,     638,     1066,     1650, ...
  [4]   24,    274,   1044,    2754,    5944,    11274,    19524, ...
  [5]  120,   1764,   8028,   24552,   60216,   127860,   245004, ...
  [6]  720,  13068,  69264,  241128,  662640,  1557660,  3272688, ...
  [7] 5040, 109584, 663696, 2592720, 7893840, 20355120, 46536624, ...

G. C. Greubel, Table of n, a(n) for the first 27 rows, flattened
Arthur T. Benjamin, David Gaebler and Robert Gaebler, A Combinatorial Approach to Hyperharmonic Numbers, INTEGERS, Electronic Journal of Combinatorial Number Theory, Volum 3, #A15, 2003.

Column 0 = A000142 (factorial numbers).
Column 1 = A000254 (Stirling numbers of first kind s(n, 2)) starting at n=1.
Column 2 = A001705 (Generalized Stirling numbers: a(n) = n!*Sum_{k=0..n-1}(k+1)/(n-k)), starting at n=1.
Column 3 = A001711 (Generalized Stirling numbers: a(n) = Sum_{k=0..n}(-1)^(n+k)*(k+1)*3^k*stirling1(n+1, k+1)).
Column 4 = A001716 (Generalized Stirling numbers: a(n) = Sum_{k=0..n}(-1)^(n+k)*(k+1)*4^k*stirling1(n+1, k+1)).
Column 5 = A001721 (Generalized Stirling numbers: a(n) = Sum_{k=0..n}(-1)^(n+k)*binomial(k+1, 1)*5^k*stirling1(n+1, k+1)).
Column 6 = A051524 (2nd unsigned column of A051338) starting at n=1.
Column 7 = A051545 (2nd unsigned column of A051339) starting at n=1.
Column 8 = A051560 (2nd unsigned column of A051379) starting at n=1.
Column 9 = A051562 (2nd unsigned column of A051380) starting at n=1.
Column 10= A051564 (2nd unsigned column of A051523) starting at n=1.
2nd row is A005408 (2n - 1, starting at n=1).
3rd row is A080663 (3n^2 - 1, starting at n=1).
Main diagonal gives A384024.
Cf. A000254, A165674 and A165675, A028421 and A126671.

H := proc(n, k) option remember; if n = 0 then 1/k else add(H(n - 1, j), j = 1..k) fi end: A := (n, k) -> (n + 1)!*H(k, n + 1):
# Alternative with standard harmonic number:
A := (n, k) -> if k = 0 then n! else (harmonic(n + k) - harmonic(k - 1))*(n + k)! / (k - 1)! fi:
for n from 0 to 7 do seq(A(n, k), k = 0..6) od;
# Alternative with hypergeometric formula:
A := (n, k) -> (n+1)*((n + k)! / k!)*hypergeom([-n, 1, 1], [2, k+1], 1):
seq(print(seq(simplify(A(n, k)), k = 0..6)), n=0..7); # Peter Luschny, Jul 01 2022

H[0, m_] := 1/m; H[n_, m_] := Sum[H[n - 1, k], {k, m}]; a[n_, m_] := m!H[n, m]; Flatten[ Table[ a[i, n - i], {n, 10}, {i, n - 1, 0, -1}]]
Table[ a[n, m], {m, 8}, {n, 0, m + 1}] // TableForm (* to view the table *)
(* Robert G. Wilson v, Jun 27 2005 *)

a(n, k) = polcoef(prod(j=0, n, 1+(j+k)*x), n); \\ Seiichi Manyama, May 19 2025

A(n, k) = (Harmonic(n + k) - Harmonic(k - 1))*(n + k)!/(k - 1)! if k > 0, otherwise n!.
From Gerald McGarvey, Aug 27 2005, edited by Peter Luschny, Jul 02 2022: (Start)
E.g.f. for column k: -log(1 - x)/(x*(1 - x)^k).
Row 3 is r(n) = 4*n^3 + 18*n^2 + 22*n + 6.
Row 4 is r(n) = 5*n^4 + 40*n^3 + 105*n^2 + 100*n + 24.
Row 5 is r(n) = 6*n^5 + 75*n^4 + 340*n^3 + 675*n^2 + 548*n + 120.
Row 6 is r(n) = 7*n^6 + 126*n^5 + 875*n^4 + 2940*n^3 + 4872*n^2 + 3528*n + 720.
Row 7 is r(n) = 8*n^7 + 196*n^6 + 1932*n^5 + 9800*n^4 + 27076*n^3 + 39396*n^2 + 26136*n + 5040.
The sum of the polynomial coefficients for the n-th row is |S1(n, 2)|, which are the unsigned Stirling1 numbers which appear in column 1.
A(m, n) = Sum_{k=1..m} n*A094645(m, n)*(n+1)^(k-1). (A094645 is Generalized Stirling number triangle of first kind, e.g.f.: (1-y)^(1-x).) (End)
In Gerard McGarvey's formulas for the row coefficients we find Wiggen's triangle A028421 and their o.g.f.s lead to Wood's polynomials A126671; see A165674. - Johannes W. Meijer, Oct 16 2009
A(n, k) = (n + 1)*((n + k)! / k!)*hypergeom([-n, 1, 1], [2, k + 1], 1). - Peter Luschny, Jul 01 2022
A(n,k) = [x^n] Product_{j=0..n} (1 + (j+k)*x). - Seiichi Manyama, May 19 2025

More terms from Robert G. Wilson v, Jun 27 2005

Edited by Peter Luschny, Jul 02 2022

A052774 a(n) = (4*n+1)^(n-1).

Original entry on oeis.org

1, 1, 9, 169, 4913, 194481, 9765625, 594823321, 42618442977, 3512479453921, 327381934393961, 34050628916015625, 3909821048582988049, 491258904256726154641, 67046038752496061076057, 9876832533361318095112441, 1562069488955406402587890625

Offset: 0

encyclopedia(AT)pommard.inria.fr, Jan 25 2000

G. C. Greubel, Table of n, a(n) for n = 0..320
INRIA Algorithms Project, Encyclopedia of Combinatorial Structures 731
J.-C. Novelli, J.-Y. Thibon, Hopf Algebras of m-permutations,(m+1)-ary trees, and m-parking functions, arXiv:1403.5962 [math.CO], 2014.

Cf. A000169, A052750, A052752, A052782, A000272, A001716.

[(4*n+1)^(n-1): n in [0..30]]; // G. C. Greubel, Nov 14 2017

spec := [S,{B=Prod(Z,S,S,S,S),S=Set(B)},labeled]: seq(combstruct[count](spec,size=n), n=0..20);
A052774:=n->(4*n+1)^(n-1); seq(A052774(n), n=0..20); # Wesley Ivan Hurt, Dec 03 2013

Table[(4n+1)^(n-1), {n,0,20}] (* Wesley Ivan Hurt, Dec 03 2013 *)
With[{nmax = 50}, CoefficientList[Series[Exp[-LambertW[-4*x]/4], {x, 0, nmax}], x]*Range[0, nmax]!] (* G. C. Greubel, Nov 14 2017 *)

for(n=0, 30, print1((4*n+1)^(n-1), ", ")) \\ G. C. Greubel, Nov 14 2017

x='x+O('x^50); Vec(serlaplace(exp(-lambertw(-4*x)/4))) \\ G. C. Greubel, Nov 14 2017

E.g.f.: exp(-1/4*LambertW(-4*x)).
a(n) = A016813(n)^A023443(n). - Wesley Ivan Hurt, Dec 03 2013
From Peter Bala, Dec 19 2013: (Start)
The e.g.f. A(x) = 1 + x + 9*x^2/2! + 169*x^3/3! + 4913*x^4/4! + ... satisfies:
1) A(x*exp(-4*x)) = exp(x) = 1/A(-x*exp(4*x));
2) A^4(x) = 1/x*series reversion(x*exp(-4*x));
3) A(x^4) = 1/x*series reversion(x*exp(-x^4));
4) A(x) = exp(x*A(x)^4);
5) A(x) = 1/A(-x*A(x)^8). (End)
E.g.f.: (-LambertW(-4*x)/(4*x))^(1/4). - Vaclav Kotesovec, Dec 07 2014
Related to A001716 by Sum_{n >= 1} a(n)*x^n/n! = series reversion( 1/(1 + x)^4*log(1 + x) ) = series reversion(x - 9*x^2/2! + 74*x^3/3! - 638*x^4/4! + ...). Cf. A000272, A052750. - Peter Bala, Jun 15 2016

Better description from Vladeta Jovovic, Sep 02 2003

A067176 A triangle of generalized Stirling numbers: sum of consecutive terms in the harmonic sequence multiplied by the product of their denominators.

Original entry on oeis.org

0, 1, 0, 3, 1, 0, 11, 5, 1, 0, 50, 26, 7, 1, 0, 274, 154, 47, 9, 1, 0, 1764, 1044, 342, 74, 11, 1, 0, 13068, 8028, 2754, 638, 107, 13, 1, 0, 109584, 69264, 24552, 5944, 1066, 146, 15, 1, 0, 1026576, 663696, 241128, 60216, 11274, 1650, 191, 17, 1, 0, 10628640

Offset: 0

Henry Bottomley, Jan 09 2002

In the Coupon Collector's Problem with n types of coupon, the expected number of coupons required until there are only k types of coupon uncollected is a(n,k)*k!/(n-1)!.
If n+k is even, then a(n,k) is divisible by (n+k+1). For n>=k and k>= 0, a(n,k) = (n-k)!*H(k+1,n-k), where H(m,n) is a generalized harmonic number, i.e., H(0,n) = 1/n and H(m,n) = Sum_{j=1..n} H(m-1,j). - Leroy Quet, Dec 01 2006
This triangle is the same as triangle A165674, which is generated by the asymptotic expansion of the higher order exponential integral E(x,m=2,n), minus the first right hand column. - Johannes W. Meijer, Oct 16 2009

			Rows start 0; 1,0; 3,1,0; 11,5,1,0; 50,26,7,1,0; 274,154,47,9,1,0 etc. a(5,2) = 3*4*5*(1/3 + 1/4 + 1/5) = 4*5 + 3*5 + 3*4 = 20 + 15 + 12 = 47.

G. C. Greubel, Table of n, a(n) for the first 50 rows, flattened

Columns are A000254, A001705, A001711, A001716, A001721, A051524, A051545, A051560, A051562, A051564, etc.

Cf. A093905 and A165674. - Johannes W. Meijer, Oct 16 2009

T[0, k_] := 1; T[n_, k_] := T[n, k] = Sum[ i*k^(i - 1)*Abs[StirlingS1[n - k, i]], {i, 1, n - k}]; Table[T[n,k], {n,1,10}, {k,1,n}] (* G. C. Greubel, Jan 21 2017 *)

a(n, k) = (n!/k!)*Sum_{j=k+1..n} 1/j = (A000254(n) - A000254(k)*A008279(n, n-k))/A000142(k) = a(n-1, k)*n + (n-1)!/k! = (a(n, k-1)-n!/k!)/k.

a(n, k) = Sum_{i=1..n-k} i*k^(i-1)*abs(stirling1(n-k, i)). - Vladeta Jovovic, Feb 02 2003

A093905 Triangle read by rows: for 0 <= k < n, a(n, k) is the sum of the products of all subsets of {n-k, n-k+1, ..., n} with k members.

Original entry on oeis.org

1, 1, 3, 1, 5, 11, 1, 7, 26, 50, 1, 9, 47, 154, 274, 1, 11, 74, 342, 1044, 1764, 1, 13, 107, 638, 2754, 8028, 13068, 1, 15, 146, 1066, 5944, 24552, 69264, 109584, 1, 17, 191, 1650, 11274, 60216, 241128, 663696, 1026576, 1, 19, 242, 2414, 19524, 127860

Offset: 1

Amarnath Murthy, Apr 24 2004

Triangle A165674, which is the reversal of this triangle, is generated by the asymptotic expansion of the higher order exponential integral E(x,m=2,n). - Johannes W. Meijer, Oct 16 2009

			Triangle begins:
1
1 3
1 5 11
1 7 26 50
1 9 47 154 274
...
a(5, 3) = 4*3*2+5*3*2+5*4*2+5*4*3 = 154.

G. C. Greubel, Table of n, a(n) for the first 50 rows, flattened

The leading diagonal is given by A000254, Stirling numbers of first kind. The next nine diagonals are A001705, A001711, A001716, A001721, A051524, A051545, A051560, A051562 and A051564, generalized Stirling numbers.
Cf. A001705, A001711, A067176, A093344, A105954.
A165674 is the reversal of this triangle. - Johannes W. Meijer, Oct 16 2009

T[n_, 0] := 1; T[n_, k_]:= Product[i, {i, n - k, n}]*Sum[1/i, {i, n - k, n}]; Table[T[n, k], {n, 1, 10}, {k, 0, n - 1}] (* G. C. Greubel, Jan 21 2017 *)

a(n, k) = prod(i=n-k, n, i)*sum(i=n-k,n,1/i);
tabl(nn) = for (n=1, nn, for (k=0, n-1, print1(a(n,k), ", ")); print()); \\ Michel Marcus, Jan 21 2017

a(n, k) = (Product_{i=n-k..n} i)*(Sum_{i=n-k..n} 1/i), where a(n, 0) = 1.

a(n, k) = A067176(n, n-k-1) = A105954(k+1, n-k). Row sums are given by A093344.

Edited and extended by David Wasserman, Apr 24 2007

A001717 Generalized Stirling numbers.

Original entry on oeis.org

1, 15, 179, 2070, 24574, 305956, 4028156, 56231712, 832391136, 13051234944, 216374987520, 3785626465920, 69751622298240, 1350747863435520, 27437426560500480, 583506719443584000, 12969079056388224000, 300749419818102528000, 7265204785551331584000

Offset: 0

N. J. A. Sloane

nonn

The asymptotic expansion of the higher order exponential integral E(x,m=3,n=4) ~ exp(-x)/x^3*(1 - 15/x + 179/x^2 - 2070/x^3 + 24574/x^4 - 305956/x^5 + ...) leads to the sequence given above. See A163931 and A163932 for more information. - Johannes W. Meijer, Oct 20 2009
From Petros Hadjicostas, Jun 25 2020: (Start)
For nonnegative integers n, m and complex numbers a, b (with b <> 0), the numbers R_n^m(a,b) were introduced by Mitrinovic (1961) and Mitrinovic and Mitrinovic (1962) using slightly different notation.
These numbers are defined via the g.f. Product_{r=0..n-1} (x - (a + b*r)) = Sum_{m=0..n} R_n^m(a,b)*x^m for n >= 0.
As a result, R_n^m(a,b) = R_{n-1}^{m-1}(a,b) - (a + b*(n-1))*R_{n-1}^m(a,b) for n >= m >= 1 with R_0^0(a,b) = 1, R_1^0(a,b) = a, R_1^1(a,b) = 1, and R_n^m(a,b) = 0 for n < m.
We have R_n^m(a,b) = Sum_{k=0}^{n-m} (-1)^k * a^k * b^(n-m-k) * binomial(m+k, k) * S1(n, m+k) for n >= m >= 0.
For the current sequence, a(n) = R_{n+2}^2(a=-4, b=-1) for n >= 0. (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).

T. D. Noe, Table of n, a(n) for n = 0..100
D. S. Mitrinovic, Sur une classe de nombres reliés aux nombres de Stirling, Comptes rendus de l'Académie des sciences de Paris, t. 252 (1961), 2354-2356. [The numbers R_n^m(a,b) are introduced.]
D. S. Mitrinovic and R. 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 [jstor stable version].
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.

nn = 20; t = Range[0, nn]! CoefficientList[Series[(1 - 9*Log[1 - x] + 10*Log[1 - x]^2)/(1 - x)^6, {x, 0, nn}], x] (* T. D. Noe, Aug 09 2012 *)

a(n) = sum(k=0, n, (-1)^(n+k)*binomial(k+2, 2)*4^k*stirling(n+2, k+2, 1)); \\ Michel Marcus, Jan 20 2016

a(n) = Sum_{k=0..n} (-1)^(n+k) * binomial(k+2, 2) * 4^k * Stirling1(n+2, k+2). - Borislav Crstici (bcrstici(AT)etv.utt.ro), Jan 26 2004
E.g.f.: (1 - 9*log(1 - x) + 10*log(1 - x)^2)/(1 - x)^6. - Vladeta Jovovic, Mar 01 2004
If we define f(n,i,a) = Sum_{k=0..n-i} binomial(n,k) * Stirling1(n-k,i) * Product_{j=0..k-1} (-a-j), then a(n-2) = |f(n,2,4)| for n>=2. - Milan Janjic, Dec 21 2008
From Petros Hadjicostas, Jun 26 2020: (Start)
a(n) = [x^2] Product_{r=0..n+1} (x + 4 + r) = (Product_{r=0..n+1} (4 + r)) * Sum_{0 <= i < j <= n+1} 1/((4 + i)*(4 + j)).
Since a(n) = R_{n+2}^2(a=-4, b=-1) and R_n^m(a,b) = R_{n-1}^{m-1}(a,b) - (a + b*(n-1))*R_{n-1}^m(a,b), we conclude that:
(i) a(n) = A001716(n) + (n+5)*a(n-1) for n >= 1;
(ii) a(n) = (n+3)!/6 + (2*n+9)*a(n-1) - (n+4)^2*a(n-2) for n >= 2.
(iii) a(n) = 3*(n+4)*a(n-1) - (3*n^2+21*n+37)*a(n-2) + (n+3)^3*a(n-3) for n >= 3. (End)

More terms from Borislav Crstici (bcrstici(AT)etv.utt.ro), Jan 26 2004

A001718 Generalized Stirling numbers.

Original entry on oeis.org

1, 22, 355, 5265, 77224, 1155420, 17893196, 288843260, 4876196776, 86194186584, 1595481972864, 30908820004608, 626110382381184, 13246845128678016, 292374329134060800, 6723367631258860800, 160883166944083161600, 4001062259532015244800

Offset: 0

N. J. A. Sloane

nonn

The asymptotic expansion of the higher order exponential integral E(x,m=4,n=4) ~ exp(-x)/x^4*(1 - 22/x + 355/x^2 - 5265/x^3 + 77224/x^4 - 1155420/x^5 + ...) leads to the sequence given above. See A163931 and A163934 for more information. - Johannes W. Meijer, Oct 20 2009
From Petros Hadjicostas, Jun 26 2020: (Start)
For nonnegative integers n, m and complex numbers a, b (with b <> 0), the numbers R_n^m(a,b) were introduced by Mitrinovic (1961) and Mitrinovic and Mitrinovic (1962) using slightly different notation.
These numbers are defined via the g.f. Product_{r=0..n-1} (x - (a + b*r)) = Sum_{m=0..n} R_n^m(a,b)*x^m for n >= 0.
As a result, R_n^m(a,b) = R_{n-1}^{m-1}(a,b) - (a + b*(n-1))*R_{n-1}^m(a,b) for n >= m >= 1 with R_0^0(a,b) = 1, R_1^0(a,b) = a, R_1^1(a,b) = 1, and R_n^m(a,b) = 0 for n < m.
We have R_n^m(a,b) = Sum_{k=0}^{n-m} (-1)^k * a^k * b^(n-m-k) * binomial(m+k, k) * S1(n, m+k) for n >= m >= 0.
For the current sequence, a(n) = R_{n+3}^3(a=-4, b=-1) for n >= 0. (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).

T. D. Noe, Table of n, a(n) for n = 0..100
D. S. Mitrinovic, Sur une classe de nombres reliés aux nombres de Stirling, Comptes rendus de l'Académie des sciences de Paris, t. 252 (1961), 2354-2356. [The numbers R_n^m(a,b) are introduced.]
D. S. Mitrinovic and R. 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 [jstor stable version].
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.

Cf. A001716, A001717, A163931, A163934.

nn = 20; t = Range[0, nn]! CoefficientList[Series[(1 - 15*Log[1 - x] + 37*Log[1 - x]^2 - 20*Log[1 - x]^3)/(1 - x)^7, {x, 0, nn}], x] (* T. D. Noe, Aug 09 2012 *)

a(n) = sum(k=0, n, (-1)^(n+k)*binomial(k+3, 3)*4^k*stirling(n+3, k+3, 1)); \\ Michel Marcus, Jan 20 2016

a(n) = Sum_{k=0..n} (-1)^(n+k) * binomial(k+3, 3) * 4^k * Stirling1(n+3, k+3). - Borislav Crstici (bcrstici(AT)etv.utt.ro), Jan 26 2004
E.g.f.: (1 - 15*log(1 - x) + 37*log(1 - x)^2 - 20*log(1 - x)^3)/(1 - x)^7. - Vladeta Jovovic, Mar 01 2004
If we define f(n,i,a) = Sum_{k=0..n-i} binomial(n,k) * Stirling1(n-k,i) * Product_{j=0..k-1} (-a-j), then a(n-3) = |f(n,3,4)| for n >= 3. - Milan Janjic, Dec 21 2008
From Petros Hadjicostas, Jun 26 2020: (Start)
a(n) = [x^3] Product_{r=0..n+2} (x + 4 + r) = (Product_{r=0..n+2} (4 + r)) * Sum_{0 <= i < j < k <= n+2} 1/((4 + i)*(4 + j)*(4 + k)).
E.g.f.: Sum_{n >= 0} a(n)/(n+3)!*x^(n+3) = -(log(1 - x))^3/(6*(1 - x)^4).
Since a(n) = R_{n+3}^3(a=-4, b=-1) and R_n^m(a,b) = R_{n-1}^{m-1}(a,b) - (a + b*(n-1))*R_{n-1}^m(a,b), we conclude that:
(i) a(n) = A001717(n) + (n+6)*a(n-1) for n >= 1;
(ii) a(n) = A001716(n) + (2*n+11)*a(n-1) - (n+5)^2*a(n-2) for n >= 2.
(iii) a(n) = (n+3)!/6 + 3*(n+5)*a(n-1) - (3*n^2+27*n+61)*a(n-2) + (n+4)^3*a(n-3) for n >= 3.
(iv) a(n) = 2*(2*n+9)*a(n-1) - (6*n^2+48*n+97)*a(n-2) + (2*n+7)*(2*n^2+14*n+25)*a(n-3) - (n+3)^4*a(n-4) = 0 for n >= 4. (End)

More terms from Borislav Crstici (bcrstici(AT)etv.utt.ro), Jan 26 2004

A138771 Triangle read by rows: T(n,k) is the number of permutations of {1,2,...,n} whose 2nd cycle has k entries; each cycle is written with the smallest element first and cycles are arranged in increasing order of their first elements (n>=1; 0<=k<=n-1). For example, 1432=(1)(24)(3) has 2 entries in the 2nd cycle; 3421=(1324) has 0 entries in the 2nd cycle.

Original entry on oeis.org

1, 1, 1, 2, 3, 1, 6, 11, 5, 2, 24, 50, 26, 14, 6, 120, 274, 154, 94, 54, 24, 720, 1764, 1044, 684, 444, 264, 120, 5040, 13068, 8028, 5508, 3828, 2568, 1560, 720, 40320, 109584, 69264, 49104, 35664, 25584, 17520, 10800, 5040

Offset: 1

Emeric Deutsch, Apr 10 2008

T(n,0)=(n-1)!=A000142(n-1).
T(n,1)=A000254(n-1).
T(n,2)=A001705(n-2).
T(n,3)=2*A001711(n-4).
T(n,4)=6*A001716(n-5).
T(n,n-1)=(n-2)! (n>=2).
Sum(kT(n,k),k=0..n-1)=(n-1)!(n-1)(n+2)/4=A138772(n).

			T(4,2)=5 because we have (1)(23)(4), (1)(24)(3), (13)(24), (12)(34) and (14)(23).
Triangle starts;
1;
1,1;
2,3,1;
6,11,5,2;
24,50,26,14,6;
120,274,154,94,54,24;

Cf. A000142, A000254, A001705, A001711, A001716, A138772.
From Johannes W. Meijer, Oct 16 2009: (Start)
A000142 equals for n=>1 the row sums.
a(n) = A165680(n) * A165675(n-1).
(End)

T:=proc (n,k) if k = 0 then factorial(n-1) elif n <= k then 0 else (n-1)*T(n-1, k)+factorial(n-2) end if end proc: for n to 9 do seq(T(n, k), k=0..n-1) end do;

T(n,k)=(n-1)T(n-1,k)+(n-2)! (1<=k<=n-1). The row generating polynomials P[n](t) satisfy: P[n+1](t)=nP[n](t)+(n-1)!(t+t^2+...+t^n).

A346846 E.g.f.: log(1 + x) / (1 - x)^4.

Original entry on oeis.org

1, 7, 50, 386, 3304, 31176, 323280, 3656880, 44890560, 594463680, 8453128320, 128473430400, 2079045964800, 35692494566400, 648044312832000, 12406994498304000, 249834635947008000, 5278539223415808000, 116768100285720576000, 2699047267616544768000, 65071515565786447872000

Offset: 1

Ilya Gutkovskiy, Aug 05 2021

nonn

Cf. A000292, A001716, A024167, A109792, A346845, A346847.

nmax = 21; CoefficientList[Series[Log[1 + x]/(1 - x)^4, {x, 0, nmax}], x] Range[0, nmax]! // Rest
Table[n! Sum[(-1)^(k + 1) Binomial[n - k + 3, 3]/k , {k, 1, n}], {n, 1, 21}]

my(x='x+O('x^25)); Vec(serlaplace(log(1+x)/(1-x)^4)) \\ Michel Marcus, Aug 06 2021

a(n) = n! * Sum_{k=1..n} (-1)^(k+1) * binomial(n-k+3,3) / k.

a(n) ~ log(2) * n^3 * n! / 6. - Vaclav Kotesovec, Aug 06 2021

A376582 Triangle of generalized Stirling numbers.

Original entry on oeis.org

1, 5, 1, 26, 7, 1, 154, 47, 9, 1, 1044, 342, 74, 11, 1, 8028, 2754, 638, 107, 13, 1, 69264, 24552, 5944, 1066, 146, 15, 1, 663696, 241128, 60216, 11274, 1650, 191, 17, 1, 6999840, 2592720, 662640, 127860, 19524, 2414, 242, 19, 1, 80627040, 30334320, 7893840, 1557660, 245004, 31594, 3382, 299, 21, 1

Offset: 0

Igor Victorovich Statsenko, Sep 29 2024

			Triangle starts:
[0]       1;
[1]       5,       1;
[2]      26,       7,       1;
[3]     154,      47,       9,        1;
[4]    1044,     342,      74,       11,       1;
[5]    8028,    2754,     638,      107,      13,     1;
[6]   69264,   24552,    5944,     1066,     146,    15,    1;
[7]  663696,  241128,   60216,    11274,    1650,   191,   17,    1;

Column k: A001705 (k=0), A001711 (k=1), A001716 (k=2), A001721 (k=3), A051524 (k=4), A051545 (k=5), A051560 (k=6).

Cf. A094587 and A173333 for m=0.

T:=(m,n,k)->add(Stirling1(i+m,m)*binomial(n+m+1,n-k-i)*(n+m-k)!/(i+m)!,i=0..n-k): m:=1: seq(seq(T(m,n,k), k=0..n), n=0..10);

T(m,n,k) = Sum_{i=0..n-k} Stirling1(i+m,m)*binomial(n+m+1,n-k-i)*(n+m-k)!/(i+m)!, for m=1.

cp's OEIS Frontend

A165675 Triangle read by rows. T(n, k) = (n - k + 1)! * H(k, n - k), where H are the hyperharmonic numbers. For 0 <= k <= n.

Views

Author

Keywords

Comments

Examples

Crossrefs

Programs

Maple

Mathematica

Python

Formula

Extensions

A105954 Array read by descending antidiagonals: A(n, k) = (n + 1)! * H(k, n + 1), where H(n, k) is a higher-order harmonic number, H(0, k) = 1/k and H(n, k) = Sum_{j=1..k} H(n-1, j), for 0 <= k <= n.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Maple

Mathematica

PARI

Formula

Extensions

A052774 a(n) = (4*n+1)^(n-1).

Views

Author

Keywords

Links

Crossrefs

Programs

Magma

Maple

Mathematica

PARI

PARI

Formula

Extensions

A067176 A triangle of generalized Stirling numbers: sum of consecutive terms in the harmonic sequence multiplied by the product of their denominators.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Mathematica

Formula

A093905 Triangle read by rows: for 0 <= k < n, a(n, k) is the sum of the products of all subsets of {n-k, n-k+1, ..., n} with k members.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Mathematica

PARI

Formula

Extensions

A001717 Generalized Stirling numbers.

Views

Author

Keywords

Comments

References

Links

Programs

Mathematica

PARI

Formula

Extensions

A001718 Generalized Stirling numbers.

Views