|
%I #10 Jan 23 2024 16:19:31
%S 0,1,1,0,1,2,1,0,0,2,1,1,1,2,2,1,1,1,1,1,2,2,1,1,0,2,0,1,1,3,1,0,2,2,
%T 2,0,1,2,2,1,1,3,1,1,1,2,1,2,0,1,2,1,1,1,2,1,2,2,1,2,1,2,1,0,2,3,1,1,
%U 2,3,1,0,1,2,1,1,2,3,1,2,1,2,1,2,2,2,2
%N The number of exponents in the prime factorization of n that are squares.
%C First differs from A125070 at n = 64.
%H Amiram Eldar, <a href="/A369428/b369428.txt">Table of n, a(n) for n = 1..10000</a>
%F Additive with a(p^e) = A010052(e).
%F Sum_{k=1..n} a(k) ~ n * (log(log(n)) + B - C), where B is Mertens's constant (A077761) and C = P(2) + Sum_{k>=2} (P(k^2+1) - P(k^2)) = 0.40999077396414387641..., and P(s) is the prime zeta function.
%p a:= n-> add(`if`(issqr(i[2]), 1, 0), i=ifactors(n)[2]):
%p seq(a(n), n=1..100); # _Alois P. Heinz_, Jan 23 2024
%t f[p_, e_] := If[IntegerQ[Sqrt[e]], 1, 0]; a[1] = 0; a[n_] := Plus @@ f @@@ FactorInteger[n]; Array[a, 100]
%o (PARI) a(n) = vecsum(apply(x -> if(issquare(x), 1, 0), factor(n)[, 2]));
%Y Cf. A000290, A010052, A077761, A197680.
%Y Similar sequences: A125070, A293439, A366074, A367512, A369426, A369427.
%K nonn,easy
%O 1,6
%A _Amiram Eldar_, Jan 23 2024
|
