|
|
A068939
|
|
a(n) = Bell(n^2), where Bell(n) are the Bell numbers, cf. A000110.
|
|
2
|
|
|
|
OFFSET
|
0,3
|
|
LINKS
|
|
|
FORMULA
|
Integral representation as n-th moment of a positive function on a positive half-axis, in Maple notation: a(n)=int(x^n*(sum(exp(-ln(x)^2/ (4*ln(k)))/(k!*sqrt(ln(k))), k=2..infinity)/ (2*exp(1)*sqrt(Pi)*x) +Dirac(1-x)/exp(1)), x=0..infinity), n=0, 1, ...
|
|
MATHEMATICA
|
|
|
PROG
|
(PARI) for(n=0, 50, print1(round(sum(i=0, 1000, i^(n^2)/(i)!)/exp(1)), ", "))
(Python)
from sympy import bell
|
|
CROSSREFS
|
|
|
KEYWORD
|
nonn
|
|
AUTHOR
|
|
|
STATUS
|
approved
|
|
|
|