結果

問題 No.888 約数の総和
ユーザー yomo3yomo3
提出日時 2020-02-13 16:38:30
言語 Python3
(3.12.2 + numpy 1.26.4 + scipy 1.12.0)
結果
AC  
実行時間 83 ms / 2,000 ms
コード長 1,148 bytes
コンパイル時間 88 ms
コンパイル使用メモリ 13,056 KB
実行使用メモリ 11,008 KB
最終ジャッジ日時 2024-04-15 22:36:04
合計ジャッジ時間 2,302 ms
ジャッジサーバーID
（参考情報） 		judge5 / judge4
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 31 ms
10,880 KB
testcase_01 AC 30 ms
11,008 KB
testcase_02 AC 30 ms
11,008 KB
testcase_03 AC 31 ms
10,880 KB
testcase_04 AC 29 ms
10,880 KB
testcase_05 AC 30 ms
11,008 KB
testcase_06 AC 29 ms
11,008 KB
testcase_07 AC 29 ms
10,880 KB
testcase_08 AC 29 ms
10,880 KB
testcase_09 AC 30 ms
10,880 KB
testcase_10 AC 30 ms
11,008 KB
testcase_11 AC 30 ms
10,880 KB
testcase_12 AC 30 ms
11,008 KB
testcase_13 AC 30 ms
11,008 KB
testcase_14 AC 30 ms
10,880 KB
testcase_15 AC 29 ms
10,880 KB
testcase_16 AC 29 ms
11,008 KB
testcase_17 AC 29 ms
10,880 KB
testcase_18 AC 31 ms
10,880 KB
testcase_19 AC 29 ms
11,008 KB
testcase_20 AC 29 ms
10,880 KB
testcase_21 AC 29 ms
10,880 KB
testcase_22 AC 43 ms
11,008 KB
testcase_23 AC 31 ms
10,880 KB
testcase_24 AC 41 ms
10,880 KB
testcase_25 AC 30 ms
11,008 KB
testcase_26 AC 52 ms
10,880 KB
testcase_27 AC 79 ms
11,008 KB
testcase_28 AC 81 ms
11,008 KB
testcase_29 AC 72 ms
10,880 KB
testcase_30 AC 83 ms
11,008 KB
testcase_31 AC 77 ms
10,880 KB
testcase_32 AC 47 ms
10,880 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)
    for p in range(3, max_div(n) + 1, 2):
        c = 0
        while n % p == 0:
            n //= p
            c += 1
        if c > 0:
            f.append((p, c))
    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)
    for k in range(6, max_div(n) + 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 factors2(N):
    ans *= sum(p**i for i in range(x + 1))
print(ans)
0