OFFSET
1,1
MAPLE
with(numtheory):
b:= n-> add(mobius(d) *4^(n/d)/n, d=divisors(n)):
a:= n-> add(b(k), k=1..n):
seq(a(n), n=1..30); # Alois P. Heinz, Sep 23 2008
MATHEMATICA
f[n_] := DivisorSum[n, MoebiusMu[#] * 4^(n/#) &] / n; Accumulate[Array[f, 26]] (* Amiram Eldar, Aug 24 2023 *)
PROG
(PARI) a(n)=sum(m=1, n, 1/m* sumdiv(m, d, moebius(d)*4^(m/d) ) ); /* Joerg Arndt, Jul 04 2011 */
CROSSREFS
KEYWORD
nonn
AUTHOR
Gary L Mullen (mullen(AT)math.psu.edu) and Ken Hicks, Jan 06 2006
EXTENSIONS
More terms from Alois P. Heinz, Sep 23 2008
STATUS
approved