結果

問題 No.1973 Divisor Sequence
ユーザー gr1msl3ygr1msl3y
提出日時 2022-06-11 10:05:08
言語 PyPy3
(7.3.15)
結果
AC  
実行時間 239 ms / 2,000 ms
コード長 556 bytes
コンパイル時間 261 ms
コンパイル使用メモリ 82,048 KB
実行使用メモリ 76,160 KB
最終ジャッジ日時 2024-09-21 18:44:40
合計ジャッジ時間 3,684 ms
ジャッジサーバーID
（参考情報） 		judge3 / judge2
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 37 ms
51,968 KB
testcase_01 AC 38 ms
51,840 KB
testcase_02 AC 121 ms
75,904 KB
testcase_03 AC 59 ms
65,024 KB
testcase_04 AC 105 ms
76,032 KB
testcase_05 AC 72 ms
73,088 KB
testcase_06 AC 100 ms
76,160 KB
testcase_07 AC 71 ms
71,552 KB
testcase_08 AC 87 ms
76,160 KB
testcase_09 AC 99 ms
75,904 KB
testcase_10 AC 111 ms
76,160 KB
testcase_11 AC 64 ms
69,376 KB
testcase_12 AC 98 ms
76,032 KB
testcase_13 AC 76 ms
75,904 KB
testcase_14 AC 106 ms
76,160 KB
testcase_15 AC 118 ms
75,904 KB
testcase_16 AC 115 ms
75,904 KB
testcase_17 AC 107 ms
76,032 KB
testcase_18 AC 70 ms
72,064 KB
testcase_19 AC 126 ms
76,160 KB
testcase_20 AC 74 ms
72,704 KB
testcase_21 AC 113 ms
75,776 KB
testcase_22 AC 107 ms
75,904 KB
testcase_23 AC 169 ms
76,032 KB
testcase_24 AC 239 ms
76,160 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

N, M = map(int, input().split())
if M == 1:
    print(1)
    exit()
MOD = 10**9+7

fact = []
v = M
for i in range(2, int(M**0.5)+1):
    if v % i:
        continue
    ct = 0
    while v % i == 0:
        v //= i
        ct += 1
    fact.append([i, ct])
if v > 1:
    fact.append([v, 1])


def solve(p, K):
    dp = [1]*(K+1)
    for _ in range(N-1):
        for i in range(K):
            dp[i+1] += dp[i]
            dp[i] %= MOD
        dp = dp[::-1]
    return sum(dp) % MOD


ans = 1
for p, c in fact:
    ans *= solve(p, c)
    ans %= MOD
print(ans)
0