OFFSET
0,3
FORMULA
a(n) = Sum_{k=1..n} k^(k*n) * binomial(n-1,k-1) for n > 0.
PROG
(PARI) my(N=10, x='x+O('x^N)); Vec(sum(k=0, N, (k^k*x/(1-k^k*x))^k))
(PARI) a(n) = if(n==0, 1, sum(k=1, n, k^(k*n)*binomial(n-1, k-1)));
CROSSREFS
KEYWORD
nonn
AUTHOR
Seiichi Manyama, Jul 03 2022
STATUS
approved