from math import floor, sqrt

def A006218(n):
    return 2 * sum(n // k for k in range(1, int(sqrt(n)) + 1)) - int(sqrt(n)) ** 2

def A061201(n):
    n3 = floor(n ** (1 / 3))
    s = 0
    for i in range(1, n3 + 1):
        s += 3 * A006218(int(n / i))
        for j in range(1, n3 + 1):
            s -= floor(n / (i * j))
    s += n3**3
    return s

for n in range(1, 10):
    print("{} = {}".format(n, A061201(n)))

# for n in range(1, 10):
#     print("{} = {}".format(n, A006218(n)))