OFFSET
1,3
FORMULA
a(n) = Sum_{d|n} mu(n/d) * A084708(d).
MATHEMATICA
u[0, _] = 1; u[k_, j_] := u[k, j] = Sum[Binomial[k - 1, i - 1] Total[u[k - i, j] #^(i - 1) & /@ Divisors[j]], {i, k}];
b[n_] := 1/n*Total[EulerPhi[#] u[Quotient[n, #], #]& /@ Divisors[n] ];
A084708[n_] := b[n]/2 + If[EvenQ[n], u[n/2, 2], Sum[Binomial[n/2 - 1/2, k] u[k, 2], {k, 0, n/2 - 1/2}]]/2;
a[n_] := Sum[MoebiusMu[n/d]*A084708[d], {d, Divisors[n]}];
CROSSREFS
KEYWORD
nonn
AUTHOR
Andrew Howroyd, Apr 09 2017
STATUS
approved