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)))
 

ZnJvbSBtYXRoIGltcG9ydCBmbG9vciwgc3FydAoKZGVmIEEwMDYyMTgobik6CiAgICByZXR1cm4gMiAqIHN1bShuIC8vIGsgZm9yIGsgaW4gcmFuZ2UoMSwgaW50KHNxcnQobikpICsgMSkpIC0gaW50KHNxcnQobikpICoqIDIKCmRlZiBBMDYxMjAxKG4pOgogICAgbjMgPSBmbG9vcihuICoqICgxIC8gMykpCiAgICBzID0gMAogICAgZm9yIGkgaW4gcmFuZ2UoMSwgbjMgKyAxKToKICAgICAgICBzICs9IDMgKiBBMDA2MjE4KGludChuIC8gaSkpCiAgICAgICAgZm9yIGogaW4gcmFuZ2UoMSwgbjMgKyAxKToKICAgICAgICAgICAgcyAtPSBmbG9vcihuIC8gKGkgKiBqKSkKICAgIHMgKz0gbjMqKjMKICAgIHJldHVybiBzCgpmb3IgbiBpbiByYW5nZSgxLCAxMCk6CiAgICBwcmludCgie30gPSB7fSIuZm9ybWF0KG4sIEEwNjEyMDEobikpKQoKIyBmb3IgbiBpbiByYW5nZSgxLCAxMCk6CiMgICAgIHByaW50KCJ7fSA9IHt9Ii5mb3JtYXQobiwgQTAwNjIxOChuKSkpCg==