login
a(n) = n! * Sum_{k=0..n} k^(3 * (n-k)) / (n-k)!.
3

%I #24 Nov 27 2022 06:44:28

%S 1,1,4,57,1444,61785,4050126,373648513,47101090744,7764843893265,

%T 1630744323319450,426925697290933401,136591846585403311620,

%U 52602030074554601172649,24058544668572618782040022,12916480280574798627072144465

%N a(n) = n! * Sum_{k=0..n} k^(3 * (n-k)) / (n-k)!.

%H Seiichi Manyama, <a href="/A358687/b358687.txt">Table of n, a(n) for n = 0..221</a>

%F E.g.f.: Sum_{k>=0} x^k * exp(k^3 * x).

%F G.f.: Sum_{k>=0} k! * x^k / (1 - k^3 * x)^(k+1).

%F log(a(n)) ~ (6*n*(log(n) - 1) + 3*log(1 + LambertW(n^(2/3))) + 2*n*LambertW(n^(2/3)) * (7*log(n) - 6*log(1 + LambertW(n^(2/3))) + 3*LambertW(n^(2/3)))) / (6*(1 + LambertW(n^(2/3)))). - _Vaclav Kotesovec_, Nov 27 2022

%t Join[{1}, Table[n! * Sum[k^(3*(n-k))/(n-k)!, {k, 0, n}], {n, 1, 20}]] (* _Vaclav Kotesovec_, Nov 27 2022 *)

%o (PARI) a(n) = n!*sum(k=0, n, k^(3*(n-k))/(n-k)!);

%o (PARI) my(N=20, x='x+O('x^N)); Vec(serlaplace(sum(k=0, N, x^k*exp(x)^k^3)))

%o (PARI) my(N=20, x='x+O('x^N)); Vec(sum(k=0, N, k!*x^k/(1-k^3*x)^(k+1)))

%Y Cf. A006153, A193421, A358688.

%Y Cf. A349880, A356673.

%K nonn

%O 0,3

%A _Seiichi Manyama_, Nov 26 2022