OFFSET
1,4
LINKS
Antti Karttunen, Table of n, a(n) for n = 1..20000
FORMULA
a(p) = 1 for p prime.
EXAMPLE
a(30) = Sum_{p|30} p^omega(30/p) = 2^omega(15) + 3^omega(10) + 5^omega(6) = 2^2 + 3^2 + 5^2 = 38.
MATHEMATICA
Table[Sum[k^PrimeNu[n/k] (PrimePi[k] - PrimePi[k - 1]) (1 - Ceiling[n/k] + Floor[n/k]), {k, n}], {n, 100}]
PROG
(PARI) a(n) = my(f=factor(n)); sum(k=1, #f~, f[k, 1]^omega(n/f[k, 1])); \\ Michel Marcus, Jun 16 2021
CROSSREFS
KEYWORD
nonn
AUTHOR
Wesley Ivan Hurt, Jun 15 2021
STATUS
approved