OFFSET
1,31
LINKS
Amiram Eldar, Table of n, a(n) for n = 1..10000
Sourabhashis Das, Wentang Kuo, and Yu-Ru Liu, On the number of prime factors with a given multiplicity over h-free and h-full numbers, Journal of Number Theory, Vol. 267 (2025), pp. 176-201; arXiv preprint, arXiv:2409.11275 [math.NT], 2024. See Theorem 1.2.
FORMULA
Sum_{A004709(k) <= x} a(k) = c * x + O(sqrt(x)/log(x)), where c = (1/zeta(3)) * Sum_{p prime} ((p-1)/(p^3-1)) = 0.24833233043359932037... (Das et al., 2025).
a(n) = log_2(A382419(n)). - Amiram Eldar, Mar 25 2025
Sum_{k=1..n} a(k) ~ c * n, where c = Sum_{p prime} ((p-1)/(p^3-1)) = 0.29850959207541746... - Vaclav Kotesovec, Mar 25 2025 (according to the above formula)
From Amiram Eldar, Apr 05 2025: (Start)
MATHEMATICA
f[k_] := Module[{e = If[k == 1, {}, FactorInteger[k][[;; , 2]]]}, If[AllTrue[e, # < 3 &], Count[e, 2], Nothing]]; Array[f, 150]
PROG
(PARI) lista(kmax) = {my(e, is); for(k = 1, kmax, e = factor(k)[, 2]; is = 1; for(i = 1, #e, if(e[i] > 2, is = 0; break)); if(is, print1(#select(x -> x == 2, e), ", "))); }
CROSSREFS
KEYWORD
nonn,easy
AUTHOR
Amiram Eldar, Sep 21 2024
STATUS
approved
