A001717 -id:A001717 - cp's OEIS frontend

A163932 Triangle related to the asymptotic expansion of E(x,m=3,n).

1, 3, 3, 11, 18, 6, 50, 105, 60, 10, 274, 675, 510, 150, 15, 1764, 4872, 4410, 1750, 315, 21, 13068, 39396, 40614, 19600, 4830, 588, 28, 109584, 354372, 403704, 224490, 68040, 11466, 1008, 36, 1026576, 3518100, 4342080, 2693250, 949095, 198450

Offset: 1

Johannes W. Meijer & Nico Baken (n.h.g.baken(AT)tudelft.nl), Aug 13 2009, Oct 22 2009

The higher order exponential integrals E(x,m,n) are defined in A163931. The general formula for the asymptotic expansion E(x,m,n) ~ E(x,m-1,n+1)/x - n*E(x,m-1,n+2)/x^2 + n*(n+1) * E(x,m-1,n+3)/x^3 - n*(n+1)*(n+2)*E(x,m-1,n+4)/x^4 + ...., m >= 1 and n >= 1.
We used this formula and the asymptotic expansion of E(x,m=2,n), see A028421, to determine that E (x,m=3,n) ~ (exp(-x)/x^3)*(1 - (3+3*n)/x + (11+18*n+6*n^2)/x^2 - (50+105*n+ 60*n^2+ 10*n^3)/x^3 + .. ). This formula leads to the triangle coefficients given above.
The asymptotic expansion leads for the values of n from one to ten to known sequences, see the cross-references.
The numerators of the o.g.f.s. of the right hand columns of this triangle lead for z=1 to A001879, see A163938 for more information.
The first Maple program generates the sequence given above and the second program generates the asymptotic expansion of E(x,m=3,n).

			The first few rows of the triangle are:
[1]
[3, 3]
[11, 18, 6]
[50, 105, 60, 10]

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

Cf. A163931 (E(x,m,n)) and A163938.
Cf. A048994 (Stirling1), A000399 (row sums).
A000254, 3*A000399, 6*A000454, 10*A000482, 15*A001233, 21*A001234 equal the first six left hand columns.
A000217, A006011 and A163933 equal the first three right hand columns.
The asymptotic expansion leads to A000399 (n=1), A001706 (n=2), A001712 (n=3), A001717 (n=4), A001722 (n=5), A051525 (n=6), A051546 (n=7), A051561 (n=8), A051563 (n=9) and A051565 (n=10).
Cf. A130534 (m=1), A028421 (m=2) and A163934 (m=4).

nmax:=8; with(combinat): for n1 from 1 to nmax do for m from 1 to n1 do a(n1, m) := (-1)^(n1+m)*binomial(m+1, 2)*stirling1(n1+1, m+1) od: od: seq(seq(a(n1,m), m=1..n1), n1=1..nmax);
# End program 1
with(combinat): imax:=6; EA:=proc(x, m, n) local E, i; E := 0: for i from m-1 to imax+1 do E := E + sum((-1)^(m+k1+1)*binomial(k1, m-1)*n^(k1-m+1)* stirling1(i, k1), k1=m-1..i)/x^(i-m+1) od: E := exp(-x)/x^(m)*E: return(E); end: EA(x, 3, n);
# End program 2

a[n_, m_] /; n >= 1 && 1 <= m <= n = (-1)^(n+m)*Binomial[m+1, 2] * StirlingS1[n+1, m+1]; Flatten[Table[a[n, m], {n, 1, 9}, {m, 1, n}]][[1 ;; 42]] (* Jean-François Alcover, Jun 01 2011, after formula *)

for(n=1,10, for(m=1,n, print1((-1)^(n+m)*binomial(m+1,2) *stirling(n+1,m+1,1), ", "))) \\ G. C. Greubel, Aug 08 2017

a(n,m) = (-1)^(n+m)*binomial(m+1,2)*stirling1(n+1,m+1) for n >= 1 and 1 <= m <= n.

Edited by Johannes W. Meijer, Sep 22 2012

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

A325139 Triangle T(n, k) = [t^n] Gamma(n + k + m + t)/Gamma(k + m + t) for m = 2 and 0 <= k <= n, read by rows.

Original entry on oeis.org

1, 2, 1, 6, 7, 1, 24, 47, 15, 1, 120, 342, 179, 26, 1, 720, 2754, 2070, 485, 40, 1, 5040, 24552, 24574, 8175, 1075, 57, 1, 40320, 241128, 305956, 134449, 24885, 2086, 77, 1, 362880, 2592720, 4028156, 2231012, 541849, 63504, 3682, 100, 1

Offset: 0

Peter Luschny, Apr 15 2019

			0:        1;
1:        2,        1;
2:        6,        7,        1;
3:       24,       47,       15,        1;
4:      120,      342,      179,       26,        1;
5:      720,     2754,     2070,      485,       40,       1;
6:     5040,    24552,    24574,     8175,     1075,      57,      1;
7:    40320,   241128,   305956,   134449,    24885,    2086,     77,    1;
8:   362880,  2592720,  4028156,  2231012,   541849,   63504,   3682,  100,   1;
9:  3628800, 30334320, 56231712, 37972304, 11563650, 1768809, 142632, 6054, 126, 1;
A:  A000142,  A001711,  A001717,  A001723, ...

Row sums are A325140.
Columns are: A000142,  A001711,  A001717,  A001723.
Family: A307419 (m=0), A325137 (m=1), this sequence (m=2).

T := (n, k) -> add(binomial(j+k, k)*(k+2)^j*abs(Stirling1(n, j+k)), j=0..n-k):
seq(seq(T(n, k), k=0..n), n=0..8);
# Note that for n > 16 Maple fails (at least in some versions) to compute the
# terms properly. Inserting 'simplify' or numerical evaluation might help.
A325139Row := proc(n) local ogf, ser; ogf := (n, k) -> GAMMA(n+k+2+x)/GAMMA(k+2+x);
ser := (n, k) -> series(ogf(n,k), x, k+2); seq(coeff(ser(n,k), x, k), k=0..n) end:
seq(A325139Row(n), n=0..9);

T(n, k) = Sum_{j=0..n-k} binomial(j+k, k)*abs(Stirling1(n, j+k))*(k+2)^j.

A376634 Triangle read by rows: T(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 = 2.

Original entry on oeis.org

1, 9, 1, 71, 12, 1, 580, 119, 15, 1, 5104, 1175, 179, 18, 1, 48860, 12154, 2070, 251, 21, 1, 509004, 133938, 24574, 3325, 335, 24, 1, 5753736, 1580508, 305956, 44524, 5000, 431, 27, 1, 70290936, 19978308, 4028156, 617624, 74524, 7155, 539, 30, 1, 924118272, 270074016, 56231712, 8969148, 1139292, 117454, 9850, 659, 33, 1, 13020978816, 3894932448, 832391136, 136954044, 18083484, 1961470, 176554, 13145, 791, 36, 1

Offset: 0

Igor Victorovich Statsenko, Sep 30 2024

The columns of the triangle T(m,n,k) represent the coefficients of the asymptotic expansion of the higher order exponential integral E(x,m+1,k+2), for m=2, k>=0. For reference see. A163931.

			Triangle starts:
 [0]          1;
 [1]          9,          1;
 [2]         71,         12,          1;
 [3]        580,        119,         15,        1;
 [4]       5104,       1175,        179,       18,        1;
 [5]      48860,      12154,       2070,      251,       21,       1;
 [6]     509004,     133938,      24574,     3325,      335,      24,     1;
 [7]    5753736,    1580508,     305956,    44524,     5000,     431,    27,     1;

Igor Victorovich Statsenko, Relationships of P-generalized Stirling numbers of the first kind with other generalized Stirling numbers, Innovation science No 10-1, State Ufa, Aeterna Publishing House, 2024, pp. 19-22. In Russian.

Column k: A001706 (k=0), A001712 (k=1), A001717 (k=2), A001722 (k=3), A051525 (k=4), A051546 (k=5), A051561 (k=6).

Cf. A094587 and A173333 for m=0, A376582 for m=1.

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:=2:seq(seq(T(m,n,k), k=0..n),n=0..10);

cp's OEIS Frontend

A163932 Triangle related to the asymptotic expansion of E(x,m=3,n).

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Maple

Mathematica

PARI

Formula

Extensions

A001718 Generalized Stirling numbers.

Views

Author

Keywords

Comments

References

Links

Crossrefs

Programs

Mathematica

PARI

Formula

Extensions

A325139 Triangle T(n, k) = [t^n] Gamma(n + k + m + t)/Gamma(k + m + t) for m = 2 and 0 <= k <= n, read by rows.

Views

Author

Keywords

Examples

Crossrefs

Programs

Maple

Formula

A376634 Triangle read by rows: T(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 = 2.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Maple

A376634 Triangle read by rows: T(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 = 2.