login
A253769
Sum of number of divisors of all positive integers <= prime(n).
1
3, 5, 10, 16, 29, 37, 52, 60, 76, 103, 113, 142, 160, 170, 188, 219, 249, 263, 294, 314, 328, 358, 379, 413, 461, 484, 494, 516, 530, 554, 637, 659, 697, 707, 768, 782, 822, 858, 878, 919, 953, 973, 1033, 1049, 1072, 1086, 1168, 1240, 1267, 1281, 1307, 1343, 1365, 1423, 1468, 1504, 1544, 1562, 1604, 1632, 1642, 1709
OFFSET
1,1
COMMENTS
a(n) is the index of the first position of prime(n) in A027750, the sequence that lists the divisors of all integers. - Michel Marcus, Oct 17 2015
FORMULA
a(n) = A006218(A000040(n)).
EXAMPLE
For n = 3 the third prime number is 5 and the sum of the number of divisors of the first five positive integers is 1 + 2 + 2 + 3 + 2 = 10, so a(3) = 10.
PROG
(PARI) a(n) = sum(i=1, prime(n), numdiv(i)); \\ Michel Marcus, Jan 15 2015
(Python)
from math import isqrt
from sympy import prime
def A253769(n): return (lambda m, p: 2*sum(p//k for k in range(1, m+1))-m*m)(isqrt(prime(n)), prime(n)) # Chai Wah Wu, Oct 09 2021
CROSSREFS
Partial sums of A139140.
Sequence in context: A070558 A233758 A301653 * A070559 A320788 A000990
KEYWORD
nonn
AUTHOR
Omar E. Pol, Jan 14 2015
STATUS
approved