OFFSET
1,2
LINKS
Seiichi Manyama, Table of n, a(n) for n = 1..1666
Wikipedia, Möbius function
MATHEMATICA
a[n_] := DivisorSum[n, MoebiusMu[n/#] * (4^# - 3^# - 2^# + 1) &] / n; Array[a, 26] (* Amiram Eldar, Oct 04 2023 *)
PROG
(PARI) {a(n) = sumdiv(n, d, moebius(n/d)*(4^d-3^d-2^d+1))/n}
CROSSREFS
KEYWORD
nonn
AUTHOR
Seiichi Manyama, Nov 23 2017
STATUS
approved