OFFSET
1,4
PROG
(PARI) {a(n) = my(u='u); n!*polcoef(polcoef(prod(k=1, n, (1-x^k+x*O(x^n))^(-u/k)), n), 3)}
CROSSREFS
KEYWORD
nonn
AUTHOR
Seiichi Manyama, Nov 10 2020
STATUS
approved