%I #10 Sep 02 2021 15:39:10
%S 1,1,4,20,123,902,7656,73509,785154,9213324,117624569,1621028312,
%T 23959376436,377730250003,6322478398476,111904530008040,
%U 2087093471665987,40891426070289970,839329531471890724,18004595602417946685,402747680140030433886,9376084240910510840672,226760664399026618376569
%N E.g.f.: exp( (exp(x) - 1) * exp(exp(x) - 1) ).
%C Exponential transform of A138378.
%C Stirling transform of A000248.
%H Alois P. Heinz, <a href="/A347339/b347339.txt">Table of n, a(n) for n = 0..465</a>
%F a(0) = 1; a(n) = Sum_{k=1..n} binomial(n-1,k-1) * (Bell(k+1) - Bell(k)) * a(n-k).
%F a(n) = Sum_{k=0..n} Stirling2(n,k) * A000248(k).
%p g:= proc(n) option remember; `if`(n=0, 1,
%p add(g(n-j)*j*binomial(n-1, j-1), j=1..n))
%p end:
%p b:= proc(n, m) option remember; `if`(n=0,
%p g(m), m*b(n-1, m)+b(n-1, m+1))
%p end:
%p a:= n-> b(n, 0):
%p seq(a(n), n=0..22); # _Alois P. Heinz_, Aug 27 2021
%t nmax = 22; CoefficientList[Series[Exp[(Exp[x] - 1) Exp[Exp[x] - 1]], {x, 0, nmax}], x] Range[0, nmax]!
%t a[0] = 1; a[n_] := a[n] = Sum[Binomial[n - 1, k - 1] (BellB[k + 1] - BellB[k]) a[n - k], {k, 1, n}]; Table[a[n], {n, 0, 22}]
%o (PARI) my(x='x+O('x^25)); Vec(serlaplace(exp((exp(x)-1)*exp(exp(x)-1)))) \\ _Michel Marcus_, Aug 27 2021
%Y Cf. A000110, A000248, A000258, A005493, A005727, A087761, A138378.
%K nonn
%O 0,3
%A _Ilya Gutkovskiy_, Aug 27 2021
|