OFFSET
1,4
COMMENTS
LINKS
Amiram Eldar, Table of n, a(n) for n = 1..10000
FORMULA
Asymptotic mean: Limit_{m->oo} (1/m) * Sum_{k=1..m} a(k) = zeta(4) * Product_{p prime} (1 + 1/p^2 + 1/p^3 - 3/p^4) = 1.57226906210272200398... .
In general, the asymptotic mean of the product of exponents in the prime factorization of the k-free numbers (numbers that are not divisible by a k-th power other than 1), for k >= 2, is zeta(k) * Product_{p prime} (1 + 1/p^2 + 1/p^3 + ... + 1/p^(k-1) - (k-1)/p^k). For k = 2 (squarefree numbers), the mean is 1 since the sequence contains only 1's. The limit when k->oo is zeta(2)*zeta(3)/zeta(6) (A082695).
MATHEMATICA
s[n_] := Times @@ FactorInteger[n][[;; , 2]]; biqFreeQ[n_] := Max[FactorInteger[n][[;; , 2]]] < 4; s /@ Select[Range[100], biqFreeQ]
PROG
(PARI) list(kmax) = {my(e); print1(1, ", "); for(k = 2, kmax, e = factor(k)[, 2]; if(vecmax(e) < 4, print1(vecprod(e), ", "))); }
CROSSREFS
KEYWORD
nonn,easy
AUTHOR
Amiram Eldar, Mar 25 2025
STATUS
approved
