A051523 -id:A051523 - cp's OEIS frontend

A051564 Second unsigned column of triangle A051523.

0, 1, 21, 362, 6026, 101524, 1763100, 31813200, 598482000, 11752855200, 240947474400, 5154170774400, 114942011990400, 2669517204076800, 64496340380102400, 1619153396908185600, 42188624389562112000

Offset: 0

Wolfdieter Lang

The asymptotic expansion of the higher order exponential integral E(x,m=2,n=10) ~ exp(-x)/x^2*(1 - 21/x + 362/x^2 - 6026/x^3 + 101524/x^4 - 1763100/x^5 + 31813200/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 A051523.

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

Cf. A049398 (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[n_] := n!*Sum[(-1)^k*Binomial[-10, k]/(n - k), {k, 0, n - 1}]; Array[f, 17, 0]
Range[0, 16]! CoefficientList[ Series[-Log[(1 - x)]/(1 - x)^10, {x, 0, 16}], x]
(* Or, using elementary symmetric functions: *)
f[k_] := k + 9; 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) = A051523(n, 2)*(-1)^(n-1).
E.g.f.: -log(1-x)/(1-x)^10.
a(n) = n!*Sum_{k=0..n-1}((-1)^k*binomial(-10,k)/(n-k)), for n>=1. - Milan Janjic, Dec 14 2008
a(n) = n!*[9]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

A051565 Third unsigned column of triangle A051523.

Original entry on oeis.org

0, 0, 1, 33, 791, 17100, 358024, 7491484, 159168428, 3463513704, 77559615576, 1792139785920, 42789106278720, 1056302350122240, 26964471256888320, 711643650545422080, 19410244660543737600, 546854985563699289600

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=10) ~ exp(-x)/x^3*(1 - 33/x + 791/x^2 - 17100/x^3 + 358024/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 A051523.

Cf. A049398 (m=0), A051564 (m=1) unsigned columns.

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

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,10)|, for n>=1. - Milan Janjic, Dec 21 2008

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

A087748 Triangle formed by reading triangle of Stirling numbers of the first kind (A048994) mod 2.

Original entry on oeis.org

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

Offset: 0

Philippe Deléham, Oct 02 2003

			Triangle begins:
1,
0, 1,
0, 1, 1,
0, 0, 1, 1,
0, 0, 1, 0, 1,
0, 0, 0, 1, 0, 1,
0, 0, 0, 1, 1, 1, 1,
0, 0, 0, 0, 1, 1, 1, 1,
0, 0, 0, 0, 1, 0, 0, 0, 1,
0, 0, 0, 0, 0, 1, 0, 0, 0, 1,
0, 0, 0, 0, 0, 1, 1, 0, 0, 1, 1,
0, 0, 0, 0, 0, 0, 1, 1, 0, 0, 1, 1,
0, 0, 0, 0, 0, 0, 1, 0, 1, 0, 1, 0, 1,
...

Brand, Neal; Das, Sajal; Jacob, Tom. The number of nonzero entries in recursively defined tables modulo primes. Proceedings of the Twenty-first Southeastern Conference on Combinatorics, Graph Theory, and Computing (Boca Raton, FL, 1990). Congr. Numer. 78 (1990), 47--59. MR1140469 (92h:05004). - From N. J. A. Sloane, Jun 03 2012

Bill Gosper, Colored illustrations of triangle of Stirling numbers of first kind read mod 2, 3, 4, 5, 6, 7

Cf. A008275, A008276, A048994, A087755.

Also parity of triangles A049444, A049459, A051338, A051379, A051523.

T(n, k) = A087755(n, k) = A048994(n, k) mod 2 = A047999([n/2], k-[(n+1)/2]) = T(n-2, k-2) XOR T(n-2, k-1) with T(0, 0) = T(1, 1) = 1 and T(1, 0) = 0; 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

cp's OEIS Frontend

A051564 Second unsigned column of triangle A051523.

Views

Author

Keywords

Comments

References

Links

Crossrefs

Programs

Mathematica

Formula

A051565 Third unsigned column of triangle A051523.

Views

Author

Keywords

Comments

References

Crossrefs

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

A087748 Triangle formed by reading triangle of Stirling numbers of the first kind (A048994) mod 2.

Views

Author

Keywords

Examples

References

Links

Crossrefs

Formula

Extensions

A087745 Numbers A001317 repeated.

Views

Author

Keywords

Comments

Crossrefs

Programs

Python

Extensions