結果

問題 No.888 約数の総和
コンテスト
ユーザー yomo3
提出日時 2020-02-13 16:47:09
言語 Python3
(3.13.1 + numpy 2.2.1 + scipy 1.14.1)
結果
AC  
実行時間 81 ms / 2,000 ms
コード長 1,185 bytes
コンパイル時間 116 ms
コンパイル使用メモリ 12,544 KB
実行使用メモリ 10,880 KB
最終ジャッジ日時 2024-10-06 05:12:23
合計ジャッジ時間 2,395 ms
ジャッジサーバーID
（参考情報） 		judge1 / judge4
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 3
other AC * 30
権限があれば一括ダウンロードができます

ソースコード

diff # 

def max_div(n):
    l, r = 1, n + 1
    while r - l > 1:
        m = (l + r) // 2
        if m * m <= n:
            l = m
        else:
            r = m
    return l


def factors(n):
    f = []

    c = 0
    while n % 2 == 0:
        n //= 2
        c += 1
    if c > 0:
        f.append((2, c))

    # lim = int(n**0.5 + 0.1)
    lim = max_div(n)
    p = 3
    while p <= lim:
        c = 0
        while n % p == 0:
            n //= p
            c += 1
        if c > 0:
            f.append((p, c))
        p += 2
    if n != 1:
        f.append((n, 1))

    return f


def factors2(n):
    f = []

    for d in (2, 3):
        c = 0
        while n % d == 0:
            n //= d
            c += 1
        if c > 0:
            f.append((d, c))

    # lim = int(n**0.5 + 0.1)
    lim = max_div(n)
    for k in range(6, lim + 6, 6):
        for d in (k - 1, k + 1):
            c = 0
            while n % d == 0:
                n //= d
                c += 1
            if c > 0:
                f.append((d, c))
    if n != 1:
        f.append((n, 1))

    return f


N = int(input())
ans = 1
for (p, x) in factors(N):
    ans *= sum(p**i for i in range(x + 1))
print(ans)
0