|
|
A071379
|
|
a(n) = (1/e) * Sum_{k>=0} ((k+4)!/k!)^(n-1)/k!.
|
|
8
|
|
|
1, 1, 209, 163121, 326922081, 1346634725665, 9939316337679281, 119802044788535500753, 2205421644124274191535553, 58945667435045762187763602753, 2198513228897522394476415669503377, 110833342180980170285766876408530089329
(list;
graph;
refs;
listen;
history;
text;
internal format)
|
|
|
OFFSET
|
0,3
|
|
COMMENTS
|
This is a Dobinski-type summation formula.
Terms quickly become gigantic: a(15) = 9142140479823239889945170786704021785456107245847570873873. a(n) appears in the process of ordering the n-th power of a product of fourth power of boson creation and fourth power of boson annihilation operators.
Let B_{m}(x) = Sum_{j>=0} (exp(j!/(j-m)!*x-1)/j!) then a(n) = n! [x^n] Taylor(B_{4}(x)), where [x^n] denotes the coefficient of x^k in the Taylor series for B_{4}(x).
a(n) is row 4 of the square array representation of A090210. (End)
|
|
LINKS
|
|
|
FORMULA
|
a(n) = (1/e)*Sum_{k>=4} fallfac(k, 4)^n / k!, n >= 1, with fallfac(n, m) := A008279(n, m) (falling factorials). (From eq.(26) with r=4 of the Schork reference.)
E.g.f. with a(0) := 1: (1/e)*(Sum_{k>=4} e^(fallfac(k, 4)*x)/k! + 8/3). From top of p. 4656 with r=4 of the Schork reference.
|
|
MAPLE
|
if n=0 then 1 else r := [seq(5, i=1..n-1)]; s := [seq(1, i=1..n-1)];
exp(-x)*24^(n-1)*hypergeom(r, s, x); round(evalf(subs(x=1, %), 99)) fi end:
|
|
MATHEMATICA
|
a[n_] := Sum[FactorialPower[k, 4]^n/k!, {k, 4, Infinity}]/E; a[0] = 1; Array[a, 12, 0] (* Jean-François Alcover, Sep 01 2016 *)
|
|
PROG
|
(PARI) default(realprecision, 500); for(n=0, 20, print1(if(n==0, 1, round(exp(-1)*sum(k=0, 500, ((k+4)!/k!)^(n-1)/k!))), ", ")) \\ G. C. Greubel, May 15 2018
|
|
CROSSREFS
|
|
|
KEYWORD
|
nonn
|
|
AUTHOR
|
|
|
EXTENSIONS
|
|
|
STATUS
|
approved
|
|
|
|