OFFSET
0,3
LINKS
Alois P. Heinz, Table of n, a(n) for n = 0..450
FORMULA
a(n) = n! * [x^n] exp(Sum_{k divides n} x^k/k).
a(n) = Sum_{d|n} A057731(n,d) for n >= 1. - Alois P. Heinz, Jul 05 2021
MAPLE
A:= proc(n, k) option remember; `if`(n<0, 0, `if`(n=0, 1,
add(mul(n-i, i=1..j-1)*A(n-j, k), j=numtheory[divisors](k))))
end:
a:= n-> A(n, n):
seq(a(n), n=0..25); # Alois P. Heinz, Feb 14 2013
MATHEMATICA
Table[a = Sum[x^i/i, {i, Divisors[n]}]; Part[Range[0, 20]! CoefficientList[Series[Exp[a], {x, 0, 20}], x], n + 1], {n, 0, 20}] (* Geoffrey Critzer, Dec 04 2011 *)
CROSSREFS
KEYWORD
easy,nonn
AUTHOR
Vladeta Jovovic, Sep 28 2002
STATUS
approved