OFFSET
1,9
COMMENTS
Sum of indices of non-unitary prime factors of n (prime factors for which the exponent exceeds 1).
LINKS
FORMULA
Additive with a(p^e) = primepi(p) = A000720(p) if e >= 2, and 0 otherwise. - Amiram Eldar, Jul 24 2024
EXAMPLE
a(450) = a(2 * 3^2 * 5^2) = a(prime(1) * prime(2)^2 * prime(3)^2) = 2 + 3 = 5.
MATHEMATICA
nmax = 104; CoefficientList[Series[Sum[k x^(Prime[k]^2)/(1 - x^(Prime[k]^2)), {k, 1, nmax}], {x, 0, nmax}], x] // Rest
f[p_, e_] := If[e == 1, 0, PrimePi[p]]; a[1] = 0; a[n_] := Plus @@ f @@@ FactorInteger[n]; Array[a, 100] (* Amiram Eldar, Jul 24 2024 *)
PROG
(PARI) A333842(n) = { my(f=factor(n)); sum(k=1, #f~, if(1==f[k, 2], 0, 1)*primepi(f[k, 1])); }; \\ Antti Karttunen, Jun 12 2020
CROSSREFS
KEYWORD
nonn
AUTHOR
Ilya Gutkovskiy, Apr 07 2020
STATUS
approved