結果

問題 No.888 約数の総和
ユーザー yomo3yomo3
提出日時 2020-02-13 16:47:09
言語 Python3
(3.12.2 + numpy 1.26.4 + scipy 1.12.0)
結果
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
このコードへのチャレンジ
(要ログイン)

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 25 ms
10,752 KB
testcase_01 AC 26 ms
10,752 KB
testcase_02 AC 26 ms
10,624 KB
testcase_03 AC 25 ms
10,880 KB
testcase_04 AC 25 ms
10,880 KB
testcase_05 AC 26 ms
10,880 KB
testcase_06 AC 27 ms
10,880 KB
testcase_07 AC 26 ms
10,880 KB
testcase_08 AC 26 ms
10,752 KB
testcase_09 AC 27 ms
10,880 KB
testcase_10 AC 27 ms
10,752 KB
testcase_11 AC 26 ms
10,880 KB
testcase_12 AC 26 ms
10,880 KB
testcase_13 AC 26 ms
10,880 KB
testcase_14 AC 27 ms
10,752 KB
testcase_15 AC 26 ms
10,880 KB
testcase_16 AC 26 ms
10,752 KB
testcase_17 AC 26 ms
10,880 KB
testcase_18 AC 27 ms
10,880 KB
testcase_19 AC 26 ms
10,880 KB
testcase_20 AC 26 ms
10,880 KB
testcase_21 AC 29 ms
10,880 KB
testcase_22 AC 44 ms
10,752 KB
testcase_23 AC 27 ms
10,880 KB
testcase_24 AC 38 ms
10,752 KB
testcase_25 AC 28 ms
10,880 KB
testcase_26 AC 46 ms
10,880 KB
testcase_27 AC 73 ms
10,880 KB
testcase_28 AC 77 ms
10,880 KB
testcase_29 AC 64 ms
10,880 KB
testcase_30 AC 81 ms
10,880 KB
testcase_31 AC 69 ms
10,880 KB
testcase_32 AC 72 ms
10,752 KB
権限があれば一括ダウンロードができます

ソースコード

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