OFFSET
1,3
LINKS
Andrew Howroyd, Table of n, a(n) for n = 1..200
FORMULA
a(n) = Sum_{d|n} mu(n/d) * A103293(d+1).
MATHEMATICA
b[n_] := SeriesCoefficient[Exp[(Exp[2*x] - 3)/2 + Exp[x]], {x, 0, n}]*n!;
c[n_] := If[n == 0, 1, (BellB[n - 1] + If[Mod[n, 2] == 1, b[(n - 1)/2], Sum[Binomial[n/2 - 1, k]*b[k], {k, 0, n/2 - 1}]])/2];
a[n_] := DivisorSum[n, MoebiusMu[n/#] c[# + 1]&];
CROSSREFS
KEYWORD
nonn
AUTHOR
Andrew Howroyd, Apr 09 2017
STATUS
approved