|
%I #20 Oct 13 2022 07:49:22
%S 0,1,1,3,13,55,241,1171,6357,37567,236521,1574331,11068333,82110535,
%T 640794337,5239439011,44723250501,397481121295,3671081354137,
%U 35176098791115,349120380267421,3583273413146647,37975511840454673,415004245048757299,4670891190907818165
%N Stirling transform of A077957 (aerated powers of 2) with 0 prepended [0, 1, 0, 2, 0, 4, 0, 8, ...].
%C a(n) without the leading zero [1, 1, 3, 13, 55, ...] is the binomial transform of A264036.
%H Seiichi Manyama, <a href="/A264037/b264037.txt">Table of n, a(n) for n = 0..567</a>
%H Eric Weisstein's MathWorld, <a href="http://mathworld.wolfram.com/BellPolynomial.html">Bell Polynomial</a>.
%F a(n) = Sum_{k=0..floor(n/2)} 2^k*Stirling2(n,2*k+1).
%F a(n) = (Bell_n(sqrt(2)) - Bell_n(-sqrt(2)))/(2*sqrt(2)), where Bell_n(x) is n-th Bell polynomial.
%F Bell_n(sqrt(2)) = A264036(n) + a(n)*sqrt(2).
%F E.g.f.: sinh(sqrt(2)*(exp(x) - 1))/sqrt(2).
%F a(n) = 0; a(n) = Sum_{k=0..n-1} binomial(n-1, k) * A264036(k). - _Seiichi Manyama_, Oct 12 2022
%e G.f. = x + x^2 + 3*x^3 + 13*x^4 + 55*x^5 + 241*x^7 + 1171*x^8 + 6357*x^9 + ...
%t Table[(BellB[n, Sqrt[2]] - BellB[n, -Sqrt[2]])/(2 Sqrt[2]), {n, 0, 24}]
%o (PARI) vector(100, n, n--; sum(k=0, n\2, 2^k*stirling(n, 2*k+1, 2))) \\ _Altug Alkan_, Nov 01 2015
%Y Cf. A077957, A264036.
%Y Cf. A024429, A357572, A357598.
%K nonn
%O 0,4
%A _Vladimir Reshetnikov_, Nov 01 2015
|
