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”).

A191564
Expansion of exp(x*Bernoulli(x)) = 1+sum(n>0, a(n)/(n!)^2*x^n)
0
1, 1, 0, -9, 0, 700, -1800, -188160, 2069760, 114198336, -3503606400, -114527952000, 9096958886400, 105465242488320, -34233337027169280, 553613206386240000, 171717284043841536000, -10454299523104595558400
OFFSET
0,4
FORMULA
a(n)=n!+(n!)^2*(sum(m=1..n-1, sum(k=1..n-m, (k!*stirling1(m+k,m)*stirling2(n-m,k))/(m+k)!)/(n-m)!)), n>0, a(0)=1.
PROG
(Maxima)
a(n):=(n)!+(n!)^2*(sum(sum((k!*stirling1(m+k, m)*stirling2(n-m, k))/(m+k)!, k, 1, n-m)/(n-m)!, m, 1, n-1));
CROSSREFS
Sequence in context: A215484 A177145 A178912 * A067479 A308803 A249489
KEYWORD
sign
AUTHOR
Vladimir Kruchinin, Jun 06 2011
STATUS
approved