login

Year-end appeal: Please make a donation to the OEIS Foundation to support ongoing development and maintenance of the OEIS. We are now in our 61st year, we have over 378,000 sequences, and we’ve reached 11,000 citations (which often say “discovered thanks to the OEIS”).

A001709
Generalized Stirling numbers.
(Formerly M5195 N2259)
4
1, 27, 511, 8624, 140889, 2310945, 38759930, 671189310, 12061579816, 225525484184, 4392554369840, 89142436976320, 1884434077831824, 41471340993035856, 949385215397800224, 22587683825903611680, 557978742043520648256, 14297219701868137003200
OFFSET
0,2
COMMENTS
The asymptotic expansion of the higher order exponential integral E(x,m=6,n=2) ~ exp(-x)/x^6*(1 - 27/x + 511/x^2 - 8624/x^3 + 140889/x^4 - ...) leads to the sequence given above. See A163931 for E(x,m,n) information and A163932 for a Maple procedure for the asymptotic expansion. - Johannes W. Meijer, Oct 20 2009
REFERENCES
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).
LINKS
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]
FORMULA
a(n) = Sum_{k=0..n} (-1)^(n+k)*binomial(k+5, 5)*2^k*Stirling1(n+5, k+5). - Borislav Crstici (bcrstici(AT)etv.utt.ro), Jan 26 2004
E.g.f.: (6-120*log(1-x)+465*log(1-x)^2-580*log(1-x)^3+261*log(1-x)^4-36*log(1-x)^5)/(6*(1-x)^7). - Vladeta Jovovic, Mar 01 2004
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-5) = |f(n,5,2)|, for n>=5. [From Milan Janjic, Dec 21 2008]
MATHEMATICA
nn = 25; t = Range[0, nn]! CoefficientList[Series[-Log[1 - x]^5/(120*(1 - x)^2), {x, 0, nn}], x]; Drop[t, 5] (* T. D. Noe, Aug 09 2012 *)
PROG
(PARI) a(n) = sum(k=0, n, (-1)^(n+k)*binomial(k+5, 5)*2^k*stirling(n+5, k+5, 1)); \\ Michel Marcus, Jan 01 2023
CROSSREFS
Sequence in context: A020568 A021734 A019752 * A016887 A110896 A215039
KEYWORD
nonn
EXTENSIONS
More terms from Borislav Crstici (bcrstici(AT)etv.utt.ro), Jan 26 2004
STATUS
approved