close
login
A398528
a(n) = Sum_{k=1..n} omega(floor(n/k)), where omega(n) = A001221(n) is the number of distinct primes dividing n.
2
0, 1, 1, 2, 2, 4, 3, 4, 4, 6, 5, 8, 7, 8, 8, 8, 8, 11, 10, 13, 12, 12, 11, 15, 14, 15, 14, 16, 15, 20, 18, 18, 18, 19, 18, 22, 21, 22, 21, 24, 23, 27, 25, 26, 26, 26, 25, 29, 27, 29, 29, 30, 29, 32, 31, 33, 32, 32, 31, 40, 38, 38, 37, 36, 36, 39, 37, 39, 38, 41
OFFSET
1,4
LINKS
Olivier Bordellès, On certain sums of number theory, International Journal of Number Theory, Vol. 18, No. 9 (2022), pp. 2053-2074; arXiv preprint, arXiv:2009.05751 [math.NT], 2020.
Mihoub Bouderbala and Meselem Karras, On a sum involving the number of distinct prime factors function related to the integer part function, Notes on Number Theory and Discrete Mathematics, Vol. 26, No. 4 (2020), pp. 52-56.
Kui Liu, Jie Wu, and Zhishan Yang, On some sums involving the integral part function, International Journal of Number Theory, Vol. 20, No. 3 (2024), pp. 831-847; arXiv preprint, arXiv:2109.01382 [math.NT], 2021.
FORMULA
a(n) = c * n + O(n^(1/2)*log(n)), where c = 0.591845... = A398529 (Bouderbala and Karras, 2020).
The error term was further improved to O(n^(455/914 + eps)), for any eps > 0 (Bordellès, 2022), and O(n^(53/110 + eps)), for any eps > 0 (Liu et al., 2024).
MATHEMATICA
a[n_] := Sum[PrimeNu[Floor[n/k]], {k, 1, n}]; Array[a, 70]
PROG
(PARI) a(n) = sum(k = 1, n, omega(n\k));
CROSSREFS
Sequence in context: A383213 A094950 A087874 * A363341 A166267 A117484
KEYWORD
nonn,easy
AUTHOR
Amiram Eldar, Aug 01 2026
STATUS
approved