OFFSET
1,4
LINKS
Andrew Howroyd, Table of n, a(n) for n = 1..200
PROG
(PARI) a(n)={n!*polcoeff(sumdiv(n, i, moebius(n/i)*exp(sumdiv(i, j, x^j/j) + O(x*x^n))), n)/eulerphi(n)} \\ Andrew Howroyd, Jul 02 2018
CROSSREFS
KEYWORD
nonn
AUTHOR
Vladeta Jovovic, Sep 30 2002
EXTENSIONS
Terms a(22) and beyond from Andrew Howroyd, Jul 02 2018
STATUS
approved