login
A353895
Expansion of e.g.f. exp( (x * (exp(x) - 1))^3 / 36 ).
3
1, 0, 0, 0, 0, 0, 20, 210, 1400, 7560, 36120, 159390, 850300, 9875580, 170133964, 2688015330, 36706233200, 444802722000, 4939264076016, 52543545234534, 583037908936500, 7645631225897700, 124931080233222340, 2327407301807577066, 44282377224446369800
OFFSET
0,7
FORMULA
a(n) = n! * Sum_{k=0..floor(n/6)} (3*k)! * Stirling2(n-3*k,3*k)/(36^k * k! * (n-3*k)!).
PROG
(PARI) my(N=30, x='x+O('x^N)); Vec(serlaplace(exp((x*(exp(x)-1))^3/36)))
(PARI) a(n) = n!*sum(k=0, n\6, (3*k)!*stirling(n-3*k, 3*k, 2)/(36^k*k!*(n-3*k)!));
CROSSREFS
KEYWORD
nonn
AUTHOR
Seiichi Manyama, May 09 2022
STATUS
approved