結果

問題 No.1973 Divisor Sequence
ユーザー pitPpitP
提出日時 2022-06-11 01:25:05
言語 PyPy3
(7.3.15)
結果
AC  
実行時間 549 ms / 2,000 ms
コード長 731 bytes
コンパイル時間 245 ms
コンパイル使用メモリ 81,696 KB
実行使用メモリ 231,592 KB
最終ジャッジ日時 2023-10-21 07:34:47
合計ジャッジ時間 5,378 ms
ジャッジサーバーID
（参考情報） 		judge11 / judge9
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 40 ms
53,544 KB
testcase_01 AC 40 ms
53,544 KB
testcase_02 AC 210 ms
106,208 KB
testcase_03 AC 71 ms
74,196 KB
testcase_04 AC 150 ms
93,504 KB
testcase_05 AC 87 ms
79,360 KB
testcase_06 AC 153 ms
94,632 KB
testcase_07 AC 85 ms
78,716 KB
testcase_08 AC 114 ms
85,632 KB
testcase_09 AC 140 ms
91,200 KB
testcase_10 AC 167 ms
97,560 KB
testcase_11 AC 79 ms
77,440 KB
testcase_12 AC 138 ms
90,840 KB
testcase_13 AC 88 ms
80,596 KB
testcase_14 AC 142 ms
93,876 KB
testcase_15 AC 193 ms
102,184 KB
testcase_16 AC 186 ms
101,696 KB
testcase_17 AC 149 ms
93,888 KB
testcase_18 AC 68 ms
72,636 KB
testcase_19 AC 176 ms
99,792 KB
testcase_20 AC 85 ms
79,040 KB
testcase_21 AC 173 ms
98,952 KB
testcase_22 AC 166 ms
97,164 KB
testcase_23 AC 549 ms
231,592 KB
testcase_24 AC 459 ms
173,348 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

def fact(num):
	i=2
	ans=dict()
	while i*i <=num:
		while num%i==0:
			num=num//i
			if i in ans:
				ans[i]+=1
			else:
				ans[i]=1
		i+=1
	if num!=1:
		ans[num]=1
	return ans

mod = 10**9+7
n,m = map(int,input().split())

primes = fact(m)
ans = 1

for _,value in primes.items():
    dp = [[0 for _ in range(n+1)] for _ in range(value + 1)]
    dp[0][0] = 1

    for i in range(n):
        s = 0
        for j in range(value+1):
            s += dp[j][i]
            s %= mod
        for j in range(value+1):
            dp[j][i+1] += s
            dp[j][i+1] %= mod
            s -= dp[value-j][i]
    
    cnt = 0
    for j in range(value+1):
        cnt += dp[j][n]
        cnt += mod
    ans *= cnt
    ans %= mod

print(ans)
0