OFFSET
0,3
PROG
(PARI) a(n) = n!*3^(n*(n-1)/2)*polcoef(exp(sum(k=1, n, x^k/(k!*3^(k*(k-1)/2)))+x*O(x^n)), n);
(PARI) T(n, k) = if(k==1, 1, sum(j=1, n-1, 3^(j*(n-j))*binomial(n-1, j)*T(j, k-1)));
a(n) = if(n==0, 1, sum(k=1, n, T(n, k)));
CROSSREFS
KEYWORD
nonn
AUTHOR
Seiichi Manyama, Jun 18 2022
STATUS
approved