OFFSET
3,1
LINKS
Alois P. Heinz, Table of n, a(n) for n = 3..1000
MAPLE
b:= proc(n, k) option remember; `if`(n=0, 1, add(add(
d*k^d, d=numtheory[divisors](j))*b(n-j, k), j=1..n)/n)
end:
a:= n-> (k-> add(b(n, k-i)*(-1)^i*binomial(k, i), i=0..k))(3):
seq(a(n), n=3..30);
CROSSREFS
KEYWORD
nonn
AUTHOR
Alois P. Heinz, Oct 07 2018
STATUS
approved