OFFSET
1,2
FORMULA
a(n) = n! * [x^n] Sum_{k>=1} (exp(n*x^k) - 1).
MATHEMATICA
Table[n! Sum[n^d/d!, {d, Divisors[n]}], {n, 1, 18}]
Table[n! SeriesCoefficient[Sum[(Exp[n x^k] - 1), {k, 1, n}], {x, 0, n}], {n, 1, 18}]
PROG
(PARI) a(n) = n! * sumdiv(n, d, n^d/d!); \\ Michel Marcus, Aug 12 2020
CROSSREFS
KEYWORD
nonn
AUTHOR
Ilya Gutkovskiy, Aug 10 2020
STATUS
approved