OFFSET
1,3
FORMULA
MAPLE
with(numtheory):
a:= n-> add(mobius(n/d)*tau(n^d-1), d=divisors(n)):
seq(a(n), n=1..30);
MATHEMATICA
a[n_] := DivisorSum[n, MoebiusMu[n/#]*DivisorSigma[0, n^#-1]&]; a[1] = 0;
Table[an = a[n]; Print["a(", n, ") = ", an]; an, {n, 1, 46}] (* Jean-François Alcover, Mar 25 2017, translated from Maple *)
PROG
(PARI) a(n) = if (n==1, 1, sumdiv(n, d, moebius(n/d)*numdiv(n^d-1))); \\ Michel Marcus, Mar 25 2017
CROSSREFS
KEYWORD
nonn
AUTHOR
Alois P. Heinz, Dec 21 2014
STATUS
approved