[go: up one dir, main page]

login
A366988
The number of prime powers of prime numbers (A053810) that divide n.
4
0, 0, 0, 1, 0, 0, 0, 2, 1, 0, 0, 1, 0, 0, 0, 2, 0, 1, 0, 1, 0, 0, 0, 2, 1, 0, 2, 1, 0, 0, 0, 3, 0, 0, 0, 2, 0, 0, 0, 2, 0, 0, 0, 1, 1, 0, 0, 2, 1, 1, 0, 1, 0, 2, 0, 2, 0, 0, 0, 1, 0, 0, 1, 3, 0, 0, 0, 1, 0, 0, 0, 3, 0, 0, 1, 1, 0, 0, 0, 2, 2, 0, 0, 1, 0, 0, 0
OFFSET
1,8
COMMENTS
a(n) depends only on the prime signature of n.
Every nonnegative number appears in the sequence of record values. k >= 1 first occurs at n = 2^prime(k) (A034785).
LINKS
FORMULA
Additive with a(p^e) = A000720(e).
a(n) = 0 if and only if n is squarefree (A005117).
Asymptotic mean: Limit_{m->oo} (1/m) * Sum_{k=1..m} a(k) = Sum_{p prime} P(p) = 0.67167522222173297323..., where P(s) is the prime zeta function.
MATHEMATICA
f[p_, e_] := PrimePi[e]; a[1] = 0; a[n_] := Plus @@ f @@@ FactorInteger[n]; Array[a, 100]
PROG
(PARI) a(n) = {my(f = factor(n)); sum(i = 1, #f~, primepi(f[i, 2])); }
KEYWORD
nonn,easy
AUTHOR
Amiram Eldar, Oct 31 2023
STATUS
approved