OFFSET
1,1
LINKS
Alois P. Heinz, Table of n, a(n) for n = 1..1000
FORMULA
a(n) ~ 5^(n+1) / (4*n). - Vaclav Kotesovec, Sep 05 2014
MAPLE
with(numtheory):
b:= n-> add(mobius(d) *5^(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[#] * 5^(n/#) &] / n; Accumulate[Array[f, 30]] (* Amiram Eldar, Aug 24 2023 *)
PROG
(PARI) a(n)=sum(m=1, n, 1/m* sumdiv(m, d, moebius(d)*5^(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