A051380 -id:A051380 - cp's OEIS frontend

A051562 Second unsigned column of triangle A051380.

0, 1, 19, 299, 4578, 71394, 1153956, 19471500, 343976400, 6366517200, 123418922400, 2503748556000, 53091962697600, 1175271048201600, 27123099523027200, 651708291206649600, 16282170039031142400

Offset: 0

Wolfdieter Lang

The asymptotic expansion of the higher order exponential integral E(x,m=2,n=9) ~ exp(-x)/x^2*(1 - 19/x + 299/x^2 - 4578/x^3 + 71394/x^4 - 1153956/x^5 + 19471500/x^6 - ...) leads to the sequence given above. See A163931 and A028421 for more information. - Johannes W. Meijer, Oct 20 2009

Mitrinovic, D. S. and Mitrinovic, R. S. see reference given for triangle A051380.

G. C. Greubel, Table of n, a(n) for n = 0..440

Cf. A049389 (first unsigned column).

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. - Gary Detlefs Jan 04 2011

f[k_] := k + 8; t[n_] := Table[f[k], {k, 1, n}]
a[n_] := SymmetricPolynomial[n - 1, t[n]]
Table[a[n], {n, 1, 16}]
(* Clark Kimberling, Dec 29 2011 *)

a(n) = A051380(n, 2)*(-1)^(n-1).
E.g.f.: -log(1-x)/(1-x)^9.
a(n) = n!*Sum_{k=0..n-1} ((-1)^k*binomial(-9,k)/(n-k)), for n>=1. - Milan Janjic, Dec 14 2008
a(n) = n!*[8]h(n), where [k]h(n) denotes the k-th successive summation of the harmonic numbers from 0 to n. - Gary Detlefs, Jan 04 2011

A051563 Third unsigned column of triangle A051380.

Original entry on oeis.org

0, 0, 1, 30, 659, 13145, 255424, 4985316, 99236556, 2030997852, 42924478536, 938984014584, 21283428847680, 500043968498880, 12176238355176960, 307176581692097280, 8023946251816984320, 216880826334455750400

Offset: 0

Wolfdieter Lang

From Johannes W. Meijer, Oct 20 2009: (Start)
The asymptotic expansion of the higher order exponential integral E(x,m=3,n=9) ~ exp(-x)/x^3*(1 - 30/x + 659/x^2 - 13145/x^3 + 255424/x^4 + ...) leads to the sequence given above. See A163931 and A163932 for more information.
(End)

Mitrinovic, D. S. and Mitrinovic, R. S. see reference given for triangle A051380.

Cf. A049389 (m=0), A051562 (m=1) unsigned columns.

a(n) = A051380(n, 2)*(-1)^n; e.g.f.: ((log(1-x))^2)/(2*(1-x)^9).

If we define f(n,i,a)=sum(binomial(n,k)*stirling1(n-k,i)*product(-a-j,j=0..k-1),k=0..n-i), then a(n) = |f(n,2,9)|, for n>=1. - Milan Janjic, Dec 21 2008

A049389 a(n) = (n+8)!/8!.

Original entry on oeis.org

1, 9, 90, 990, 11880, 154440, 2162160, 32432400, 518918400, 8821612800, 158789030400, 3016991577600, 60339831552000, 1267136462592000, 27877002177024000, 641171050071552000, 15388105201717248000, 384702630042931200000, 10002268381116211200000

Offset: 0

Wolfdieter Lang

The asymptotic expansion of the higher-order exponential integral E(x,m=1,n=9) ~ exp(-x)/x*(1 - 9/x + 90/x^2 - 990/x^3 + 11880/x^4 - 154440/x^5 + ...) leads to the sequence given above. See A163931 and A130534 for more information. - Johannes W. Meijer, Oct 20 2009

Vincenzo Librandi, Table of n, a(n) for n = 0..300

Cf. A000142, A001710, A001715, A001720, A001725, A001730, A049388, A051379, A051380.

Cf. A130534, A163931, A173333, A245334.

a049389 = (flip div 40320) . a000142 . (+ 8)
-- Reinhard Zumkeller, Aug 31 2014

[Factorial(n+8)/40320: n in [0..25]]; // Vincenzo Librandi, Jul 20 2011

a[n_] := (n + 8)!/8!; Array[a, 20, 0] (* Amiram Eldar, Jan 15 2023 *)

PARI
```
a(n) = (n+8)!/8!;
```

a(n)= A051380(n, 0)*(-1)^n (first unsigned column of triangle).
a(n) = (n+8)!/8!.
E.g.f.: 1/(1-x)^9.
a(n) = A173333(n+8,8). - Reinhard Zumkeller, Feb 19 2010
a(n) = A245334(n+8,n) / 9. - Reinhard Zumkeller, Aug 31 2014
From Amiram Eldar, Jan 15 2023: (Start)
Sum_{n>=0} 1/a(n) = 40320*e - 109600.
Sum_{n>=0} (-1)^n/a(n) = 40320/e - 14832. (End)

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

A051523 Generalized Stirling number triangle of first kind.

Original entry on oeis.org

1, -10, 1, 110, -21, 1, -1320, 362, -33, 1, 17160, -6026, 791, -46, 1, -240240, 101524, -17100, 1435, -60, 1, 3603600, -1763100, 358024, -38625, 2335, -75, 1, -57657600, 31813200, -7491484, 976024, -75985, 3535, -91, 1, 980179200, -598482000, 159168428, -24083892, 2267769, -136080, 5082, -108, 1

Offset: 0

Wolfdieter Lang

a(n,m)= ^10P_n^m in the notation of the given reference with a(0,0) := 1. The monic row polynomials s(n,x) := sum(a(n,m)*x^m,m=0..n) which are s(n,x)= product(x-(10+k),k=0..n-1), n >= 1 and s(0,x)=1 satisfy s(n,x+y) = sum(binomial(n,k)*s(k,x)*S1(n-k,y),k=0..n), with the Stirling1 polynomials S1(n,x)=sum(A008275(n,m)*x^m, m=1..n) and S1(0,x)=1. In the umbral calculus (see the S. Roman reference given in A048854) the s(n,x) polynomials are called Sheffer for (exp(10*t),exp(t)-1).

			{1}; {-10,1}; {110,-21,1}; {-1320,362,-331}; ... s(2,x)= 110-21*x+x^2; S1(2,x)= -x+x^2 (Stirling1).

Reinhard Zumkeller, Rows n = 0..125 of triangle, flattened
D. S. Mitrinovic, M. S. Mitrinovic, Tableaux d'une classe de nombres reliés aux nombres de Stirling, Univ. Beograd. Pubi. Elektrotehn. Fak. Ser. Mat. Fiz. 77 (1962).

The first (m=0) unsigned column sequence is A049398. Row sums (signed triangle): A049389(n)*(-1)^n. Row sums (unsigned triangle): A051431(n).
See also A049444 A049458 A049459 A049460 A051338 A051339 A051379 A051380
Cf. A000035 A084938.

a051523 n k = a051523_tabl !! n !! k
a051523_row n = a051523_tabl !! n
a051523_tabl = map fst $ iterate (\(row, i) ->
   (zipWith (-) ([0] ++ row) $ map (* i) (row ++ [0]), i + 1)) ([1], 10)
-- Reinhard Zumkeller, Mar 12 2014

a[n_, m_] := Pochhammer[m + 1, n - m] SeriesCoefficient[Log[1 + x]^m/(1 + x)^10, {x, 0, n}];
Table[a[n, m], {n, 0, 8}, {m, 0, n}] // Flatten (* Jean-François Alcover, Oct 29 2019 *)

a(n, m)= a(n-1, m-1) - (n+9)*a(n-1, m), n >= m >= 0; a(n, m) := 0, n

E.g.f. for m-th column of signed triangle: ((log(1+x))^m)/(m!*(1+x)^10).

Triangle (signed) = [ -10, -1, -11, -2, -12, -3, -13, -14, -4, ...] DELTA A000035; triangle (unsigned) = [10, 1, 11, 2, 12, 3, 13, 4, 14, 5, 15, ...] DELTA A000035; where DELTA is Deléham's operator defined in A084938.

If we define f(n,i,a)=sum(binomial(n,k)*stirling1(n-k,i)*product(-a-j,j=0..k-1),k=0..n-i), then T(n,i) = f(n,i,10), for n=1,2,...;i=0...n. - Milan Janjic, Dec 21 2008

A087755 Triangle read by rows: Stirling numbers of the first kind (A008275) mod 2.

Original entry on oeis.org

1, 1, 1, 0, 1, 1, 0, 1, 0, 1, 0, 0, 1, 0, 1, 0, 0, 1, 1, 1, 1, 0, 0, 0, 1, 1, 1, 1, 0, 0, 0, 1, 0, 0, 0, 1, 0, 0, 0, 0, 1, 0, 0, 0, 1, 0, 0, 0, 0, 1, 1, 0, 0, 1, 1, 0, 0, 0, 0, 0, 1, 1, 0, 0, 1, 1, 0, 0, 0, 0, 0, 1, 0, 1, 0, 1, 0, 1, 0, 0, 0, 0, 0, 0, 1, 0, 1, 0, 1, 0, 1, 0, 0, 0, 0, 0, 0, 1, 1, 1, 1, 1, 1, 1, 1

Offset: 1

Philippe Deléham, Oct 02 2003

Essentially also parity of Mitrinovic's triangles A049458, A049460, A051339, A051380.

			Triangle begins:
1
1 1
0 1 1
0 1 0 1
0 0 1 0 1
0 0 1 1 1 1
0 0 0 1 1 1 1
0 0 0 1 0 0 0 1
0 0 0 0 1 0 0 0 1
0 0 0 0 1 1 0 0 1 1
0 0 0 0 0 1 1 0 0 1 1
0 0 0 0 0 1 0 1 0 1 0 1
0 0 0 0 0 0 1 0 1 0 1 0 1
0 0 0 0 0 0 1 1 1 1 1 1 1 1

Das, Sajal K., Joydeep Ghosh, and Narsingh Deo. "Stirling networks: a versatile combinatorial topology for multiprocessor systems." Discrete applied mathematics 37 (1992): 119-146. See p. 122. - N. J. A. Sloane, Nov 20 2014

p = 2; s=14; S1T = matrix(s,s,n,k, if(k==1,(-1)^(n-1)*(n-1)!)); for(n=2,s,for(k=2,n, S1T[n,k]=-(n-1)*S1T[n-1,k]+S1T[n-1,k-1]));
S1TMP = matrix(s,s,n,k, S1T[n,k]%p);
for(n=1,s,for(k=1,n,print1(S1TMP[n,k]," "));print()) /* Gerald McGarvey, Oct 17 2009 */

T(n, k) = A087748(n, k) = A008275(n, k) mod 2 = A047999([n/2], k-[(n+1)/ 2]) = T(n-2, k-2) XOR T(n-2, k-1) with T(1, 1) = T(2, 1) = T(2, 2) = 1; T(2n, k) = T(2n-1, k-1) XOR T(2n-1, k); T(2n+1, k) = T(2n, k-1). - Henry Bottomley, Dec 01 2003

Edited and extended by Henry Bottomley, Dec 01 2003

A087756 a(n) = A087745(n+1).

Original entry on oeis.org

1, 3, 3, 5, 5, 15, 15, 17, 17, 51, 51, 85, 85, 255, 255, 257, 257, 771, 771, 1285, 1285, 3855, 3855, 4369, 4369, 13107, 13107, 21845, 21845, 65535, 65535, 65537, 65537, 196611, 196611, 327685, 327685, 983055, 983055, 1114129, 1114129

Offset: 0

Philippe Deléham, Oct 02 2003

Essentially a duplicate of A087745.
Mitrinovic's triangles A049458, A049460, A051339, A051380 and triangle of Stirling numbers of first kind (A008275) mod 2 converted to decimal.
See also A001317 = [1, 3, 5, 15, 17, ...].

Cf. A001317, A008275, A049458, A049460, A051339, A051380.

Edited by Omar E. Pol and N. J. A. Sloane, Dec 26 2008

Showing 1-7 of 7 results.

cp's OEIS Frontend

A051562 Second unsigned column of triangle A051380.

Views

Author

Keywords

Comments

References

Links

Crossrefs

Programs

Mathematica

Formula

A051563 Third unsigned column of triangle A051380.

Views

Author

Keywords

Comments

References

Crossrefs

Formula

A049389 a(n) = (n+8)!/8!.

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

Haskell

Magma

Mathematica

PARI

Formula

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

A051523 Generalized Stirling number triangle of first kind.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Haskell

Mathematica

Formula

A087755 Triangle read by rows: Stirling numbers of the first kind (A008275) mod 2.

Views

Author

Keywords

Comments

Examples

References

Programs

PARI

Formula

Extensions

A087756 a(n) = A087745(n+1).

Views

Author

Keywords

Comments

Crossrefs

Extensions