OFFSET
0,3
LINKS
FORMULA
a(n) = [x^n] exp(Sum_{k>=1} ( Sum_{d|k} Sum_{j|d} d*j^n*mu(d/j) ) * x^k/k).
MATHEMATICA
Table[SeriesCoefficient[Product[1/(1 - x^k)^Sum[d^n MoebiusMu[k/d], {d, Divisors[k]}], {k, 1, n}], {x, 0, n}], {n, 0, 16}]
Table[SeriesCoefficient[Exp[Sum[Sum[Sum[d j^n MoebiusMu[d/j], {j, Divisors[d]}], {d, Divisors[k]}] x^k/k, {k, 1, n}]], {x, 0, n}], {n, 0, 16}]
CROSSREFS
KEYWORD
nonn
AUTHOR
Ilya Gutkovskiy, Nov 01 2018
STATUS
approved