OFFSET
1,3
LINKS
Amiram Eldar, Table of n, a(n) for n = 1..10000
FORMULA
Multiplicative with a(p^e) = p^ceiling(e/3) - 1.
Dirichlet g.f.: zeta(s) * zeta(3*s-1) * Product_{p prime} (1 + 1/p^(s-1) - 2/p^s - 1/p^(3*s-1) + 1/p^(4*s-1)).
Sum_{k=1..n} a(k) ~ c * n^2 / 2, where c = zeta(2) * zeta(5) * Product_{p prime} (1 - 3/p^2 + 2/p^3 - 1/p^5 + 1/p^6 + 1/p^7 - 1/p^8) = 0.47590270384210734089... .
MATHEMATICA
f[p_, e_] := p^Ceiling[e/3] - 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); }
CROSSREFS
KEYWORD
nonn,mult,easy
AUTHOR
Amiram Eldar, Dec 30 2025
STATUS
approved
