|
|
A295659
|
|
Number of exponents larger than 2 in the prime factorization of n.
|
|
13
|
|
|
0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 1, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 1, 0, 1, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 1, 1, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 1, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 1
(list;
graph;
refs;
listen;
history;
text;
internal format)
|
|
|
OFFSET
|
1,216
|
|
LINKS
|
|
|
FORMULA
|
Additive with a(p^e) = 1 if e>2, 0 otherwise.
Asymptotic mean: lim_{n->oo} (1/n) * Sum_{k=1..n} a(k) = Sum_{p prime} 1/p^3 = 0.174762... (A085541). - Amiram Eldar, Nov 01 2020
|
|
EXAMPLE
|
For n = 120 = 2^3 * 3^1 * 5^1 there is only one exponent larger than 2, thus a(120) = 1.
For n = 216 = 2^3 * 3^3 there are two exponents larger than 2, thus a(216) = 2.
|
|
MATHEMATICA
|
Array[Count[FactorInteger[#][[All, -1]], _?(# > 2 &)] &, 105] (* Michael De Vlieger, Nov 28 2017 *)
|
|
PROG
|
(PARI) a(n) = { my(v = factor(n)[, 2], i=0); for(x=1, length(v), if(v[x]>2, i++)); i; } \\ Iain Fox, Nov 29 2017
|
|
CROSSREFS
|
Cf. A001221, A003557, A053150, A056169, A056170, A085541, A162642, A295657, A295662, A295883, A295884.
|
|
KEYWORD
|
nonn
|
|
AUTHOR
|
|
|
STATUS
|
approved
|
|
|
|