In this Kata we focus on finding a sum S(n) which is the total number of divisors taken for all natural numbers less or equal to n. More formally, we investigate the sum of n components denoted by d(1) + d(2) + ... + d(n) in which for any i starting from 1 up to n the value of d(i) tells us how many distinct numbers divide i without a remainder. 

Your solution should work for possibly large values of n without a timeout.
Assume n to be greater than zero and not greater than 999 999 999 999 999.
Brute force approaches will not be feasible options in such cases. It is fairly simple to conclude that for every n>1 there holds a recurrence S(n) = S(n-1) + d(n) with initial case S(1) = 1.

For example:
S(1) = 1
S(2) = 3
S(3) = 5
S(4) = 8
S(5) = 10

But is the fact useful anyway? If you find it is rather not, maybe this will help: 

Try to convince yourself that for any natural k, the number S(k) is the same as the number of pairs (m,n) that solve the inequality mn <= k in natural numbers.
Once it becomes clear, we can think of a partition of all the solutions into classes just by saying that a pair (m,n) belongs to the class indexed by n.
The question now arises if it is possible to count solutions of n-th class. If f(n) stands for the number of solutions that belong to n-th class, it means that S(k) = f(1) + f(2) + f(3) + ...

The reasoning presented above leads us to some kind of a formula for S(k), however not necessarily the most efficient one. Can you imagine that all the solutions to inequality mn <= k can be split using sqrt(k) as pivotal item?

Difficulty: interview

Starter code:

def count_divisors(n):
	

Solutions:
def count_divisors(n):
  """Counts the integer points under the parabola xy = n.

  Because the region is symmetric about x = y, it is only necessary to sum up
  to that point (at n^{1/2}), and double it. By this method, a square region is
  counted twice, and thus subtracted off the total.
  """
  r = int(n**(1/2))
  return 2*sum(n // i for i in range(1, r+1)) - r*r

count_divisors=lambda n:2*sum(n//i for i in range(1,int(n**.5)+1))-int(n**.5)**2

def count_divisors(n):
    return 2 * sum(n//k for k in range(int(n ** 0.5), 0, -1)) - int(n ** 0.5) ** 2

# OEIS-A006218
def count_divisors(n):
    x = int(n**0.5)
    return 2*sum(n//y for y in range(1, x+1)) - x**2

def count_divisors(n):
    return 2 * sum(n // i for i in range(1, int(n ** .5) + 1 )) - int(n ** .5)**2

def count_divisors(n):
    sq = int(n**.5)
    return 2 * sum(n // i for i in range(1, sq+1)) - sq**2

def count_divisors(n):
    rt = int(n ** .5)
    su = 0

    for k in range(1, rt + 1): su += n // k

    return 2 * su - rt * rt

import math

def count_divisors(n):
    sum = 0
    h = int(math.sqrt(n))
    for i in range(1, h + 1):
        sum += n // i
    return 2 * sum - h * h

def count_divisors(n): # http://oeis.org/A006218
    sq = int(n ** .5)
    return 2 * sum(n // i for i in range(1, sq + 1)) - sq * sq

# Python 3.8 when
# https://stackoverflow.com/a/53983683
def isqrt(n):
    if n > 0:
        x = 1 << (n.bit_length() + 1 >> 1)
        while True:
            y = (x + n // x) >> 1
            if y >= x:
                return x
            x = y
    elif n == 0:
        return 0
    else:
        raise ValueError("square root not defined for negative numbers")

def count_divisors(n):
    s = isqrt(n)
    return 2*sum(n//i for i in range(1,s+1))-s*s