結果

問題 No.1973 Divisor Sequence
ユーザー shobonvipshobonvip
提出日時 2022-06-11 18:42:22
言語 PyPy3
(7.3.15)
結果
AC  
実行時間 1,181 ms / 2,000 ms
コード長 717 bytes
コンパイル時間 155 ms
コンパイル使用メモリ 82,328 KB
実行使用メモリ 388,840 KB
最終ジャッジ日時 2024-09-22 05:37:56
合計ジャッジ時間 7,981 ms
ジャッジサーバーID
（参考情報） 		judge4 / judge2
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 38 ms
53,396 KB
testcase_01 AC 37 ms
52,820 KB
testcase_02 AC 418 ms
247,780 KB
testcase_03 AC 67 ms
75,464 KB
testcase_04 AC 268 ms
173,540 KB
testcase_05 AC 103 ms
95,616 KB
testcase_06 AC 282 ms
183,796 KB
testcase_07 AC 96 ms
88,780 KB
testcase_08 AC 187 ms
132,904 KB
testcase_09 AC 235 ms
153,776 KB
testcase_10 AC 325 ms
185,024 KB
testcase_11 AC 87 ms
84,112 KB
testcase_12 AC 243 ms
158,620 KB
testcase_13 AC 113 ms
95,456 KB
testcase_14 AC 235 ms
139,896 KB
testcase_15 AC 368 ms
223,760 KB
testcase_16 AC 366 ms
222,816 KB
testcase_17 AC 287 ms
169,344 KB
testcase_18 AC 89 ms
80,800 KB
testcase_19 AC 359 ms
199,464 KB
testcase_20 AC 100 ms
87,308 KB
testcase_21 AC 343 ms
207,568 KB
testcase_22 AC 305 ms
190,208 KB
testcase_23 AC 680 ms
283,484 KB
testcase_24 AC 1,181 ms
388,840 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

def pfact(m):
	pf = {}
	for i in range(2,int(m**0.5)+1):
		while m%i == 0:
			pf[i] = pf.get(i,0) + 1
			m //= i
	if m>1 : pf[m]=1
	return pf


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

#pairwiseな素因数ごとに独立
# p^nについて考えよう

# 1   -> 1, p, p^2, p^3
# p   -> 1, p, p^2
# p^2 -> 1, p
# p^3 -> 1

ans = 1
pf = pfact(m)
for num, cnt in pf.items():
	dp = [[0] * (cnt+1) for i in range(n)]
	for j in range(cnt+1):
		dp[0][j] = 1
	for i in range(n-1):
		for j in range(cnt+1):
			dp[i+1][cnt-j] = dp[i][j]
		for j in range(cnt-1,-1,-1):
			dp[i+1][j] += dp[i+1][j+1]
			dp[i+1][j] %= mod
	tmp = 0
	for j in range(cnt+1):
		tmp += dp[n-1][j]
		tmp %= mod
	ans *= tmp
	ans %= mod

print(ans)
0