|
|
|
|
1, 4, 25, 211, 2236, 28471, 422899, 7173580, 136750051, 2893057381, 67241818876, 1702829138209, 46659181547785, 1375237342827076, 43380198327693361, 1458027134026128691, 52014149849253158284, 1962794208713975883415
(list;
graph;
refs;
listen;
history;
text;
internal format)
|
|
|
OFFSET
|
1,2
|
|
LINKS
|
|
|
FORMULA
|
E.g.f. exp(-1+1/sqrt(1-2*x))-1.
Representation of a(n) as n-th moment of a positive function on (0, infinity): a(n)=int(x^n* (x/2)^(-1/2)*exp(-x/2)*(2*hypergeom([], [3/2, 1/2], 1/8*x)/Pi^(1/2)+1/2*sqrt(2)*sqrt(x)*hypergeom([], [2, 3/2], 1/8*x))/(4*exp(1)), x=0..infinity), n=1, 2, ... - Karol A. Penson, Jun 27 2002
Asymptotic expansion for large n: a(n) -> 2^(1/6)*(n^(-1/3) + 2^(-7/3)*n^(-2/3) + O(1/n))*(2*n)^n*exp(-n+(3/2)*(2*n)^(1/3))/(sqrt(3)*exp(1)); (the nature of this approximation of a(n) is the same as that of Stirling approximation of n!). - Karol A. Penson, Sep 02 2002
|
|
MATHEMATICA
|
a[n_, k_] := 2^(n+k)*n!/(4^n*n*k!)*Sum[(j+k)*2^(j)*Binomial[j+k-1, k-1]*Binomial[2*n-j-k-1, n-1], {j, 0, n-k}]; a[n_] := Sum[a[n, k], {k, 1, n}]; Table[a[n], {n, 1, 18}] (* Jean-François Alcover, Jul 05 2013, after Emanuele Munarini *)
Table[Sum[BellY[n, k, (2 Range[n] - 1)!!], {k, n}], {n, 20}] (* Vladimir Reshetnikov, Nov 09 2016 *)
|
|
PROG
|
(Maxima) a(n, k):=2^(n+k)*n!/(4^n*n*k!)*sum((j+k)*2^(j)*binomial(j+k-1, k-1)*binomial(2*n-j-k-1, n-1), j, 0, n-k); makelist(sum(a(n, k), k, 0, n), n, 1, 12); /* Emanuele Munarini, Jun 01 2012 */
|
|
CROSSREFS
|
|
|
KEYWORD
|
easy,nonn
|
|
AUTHOR
|
|
|
STATUS
|
approved
|
|
|
|