OFFSET
1,2
LINKS
Amiram Eldar, Table of n, a(n) for n = 1..10000
FORMULA
Multiplicative with a(p^e) = (p^(ceiling(e/3)+1) - 1)/(p-1).
Sum_{k=1..n} a(k) ~ c * n^2 / 2, where c = zeta(2) * zeta(5) * Product_{p prime} (1 - 1/p^2 - 1/p^5 + 1/p^6) = 1.01185914843439887476... .
MATHEMATICA
f[p_, e_] := (p^(Ceiling[e/3]+1) - 1)/(p - 1); a[1] = 1; a[n_] := Times @@ f @@@ FactorInteger[n]; Array[a, 100]
PROG
(PARI) a(n) = {my(f = factor(n)); prod(i = 1, #f~, (f[i, 1]^(ceil(f[i, 2]/3)+1) - 1)/(f[i, 1] - 1)); }
CROSSREFS
KEYWORD
nonn,mult,easy
AUTHOR
Amiram Eldar, Nov 14 2025
STATUS
approved
