A009599 Expansion of e.g.f. sinh(sinh(x)*exp(x)).

0, 1, 2, 5, 20, 117, 782, 5441, 39496, 306921, 2602682, 24116413, 241121564, 2561633245, 28613237382, 334511450617, 4089814554384, 52302564139985, 699179303859698, 9751200460426357, 141494250613386916

Offset

0,3

Links

Vincenzo Librandi, Table of n, a(n) for n = 0..200

Formula

a(n) = Sum_{k=0..n} Stirling2(n, k)*(1-(-1)^k)/2*2^(n-k). - Vladeta Jovovic, Sep 26 2003

G.f.: Sum_{k>=0} x^(2*k+1)/Product_{i=0..2*k+1} (1 - 2*i*x). - Sergei N. Gladkovskii, Jan 06 2013

G.f.: x/( G(0)-x^2 ) where G(k) = x^2 + (4*x*k-1)*(4*x*k+2*x-1) - x^2*(4*x*k-1)*(4*x*k+2*x-1)/G(k+1); (recursively defined continued fraction). - Sergei N. Gladkovskii, Jan 06 2013

Mathematica

Table[Sum[StirlingS2[n, k]*(1-(-1)^k)/2*2^(n-k), {k, 0, n}], {n, 0, 20}] (* Vaclav Kotesovec, Aug 06 2014 after Vladeta Jovovic *)
Table[(BellB[n, 1/2] - BellB[n, -1/2]) 2^(n-1), {n, 0, 20}] (* Vladimir Reshetnikov, Nov 01 2015 *)
With[{nn=20}, CoefficientList[Series[Sinh[Sinh[x]Exp[x]], {x, 0, nn}], x] Range[ 0, nn]!] (* Harvey P. Dale, Jun 02 2017 *)

Prog

(PARI) x='x+O('x^30); concat([0], Vec(serlaplace(sinh(sinh(x)*exp(x))))) \\ G. C. Greubel, Jan 22 2018
(Magma) [(&+[2^(n-k)*StirlingSecond(n, k)*(1 - (-1)^k)/2: k in [0..n]]): n in [0..30]]; // G. C. Greubel, Jan 22 2018

Crossrefs

Cf. A065143.

Sequence in context: A052850 A000130 A288841 * A112833 A144503 A012321

Adjacent sequences: A009596 A009597 A009598 * A009600 A009601 A009602

Keyword

nonn,easy

Author

R. H. Hardin

Extensions

Extended and signs tested by Olivier Gérard, Mar 15 1997

Status

approved