%I #11 Oct 21 2021 01:25:34
%S 0,1,8,51,312,1965,13248,97839,800208,7260921,72806040,801515979,
%T 9620317512,125071036389,1751016829968,26265324194055,420245416687392,
%U 7144172815479921,128595113003161512,2443307154421058019,48866143111666389720,1026189005418216656541,22576158119430894214368
%N a(n) = n! * Sum_{k=0..n-1} 3^k / k!.
%F E.g.f.: x * exp(3*x) / (1 - x).
%F a(0) = 0; a(n) = n * (a(n-1) + 3^(n-1)).
%F a(n) ~ exp(3)*n!. - _Stefano Spezia_, Oct 11 2021
%t Table[n! Sum[3^k/k!, {k, 0, n - 1}], {n, 0, 22}]
%t nmax = 22; CoefficientList[Series[x Exp[3 x]/(1 - x), {x, 0, nmax}], x] Range[0, nmax]!
%o (PARI) a(n) = n!*sum(k=0, n-1, 3^k/k!); \\ _Michel Marcus_, Oct 11 2021
%Y Cf. A007526, A053486, A066534, A348314.
%K nonn
%O 0,3
%A _Ilya Gutkovskiy_, Oct 11 2021
|