結果

問題 No.1973 Divisor Sequence
ユーザー shobonvipshobonvip
提出日時 2022-06-11 18:42:22
言語 PyPy3
(7.3.15)
結果
AC  
実行時間 1,085 ms / 2,000 ms
コード長 717 bytes
コンパイル時間 158 ms
コンパイル使用メモリ 81,688 KB
実行使用メモリ 387,912 KB
最終ジャッジ日時 2023-10-22 04:16:55
合計ジャッジ時間 7,703 ms
ジャッジサーバーID
（参考情報） 		judge15 / judge12
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 38 ms
53,648 KB
testcase_01 AC 35 ms
53,444 KB
testcase_02 AC 387 ms
247,252 KB
testcase_03 AC 61 ms
74,820 KB
testcase_04 AC 249 ms
172,996 KB
testcase_05 AC 98 ms
94,928 KB
testcase_06 AC 258 ms
183,088 KB
testcase_07 AC 89 ms
88,068 KB
testcase_08 AC 170 ms
132,336 KB
testcase_09 AC 217 ms
153,208 KB
testcase_10 AC 291 ms
184,500 KB
testcase_11 AC 79 ms
83,608 KB
testcase_12 AC 224 ms
158,060 KB
testcase_13 AC 103 ms
94,784 KB
testcase_14 AC 217 ms
139,240 KB
testcase_15 AC 341 ms
223,328 KB
testcase_16 AC 336 ms
222,236 KB
testcase_17 AC 269 ms
168,932 KB
testcase_18 AC 79 ms
78,976 KB
testcase_19 AC 343 ms
198,788 KB
testcase_20 AC 90 ms
86,652 KB
testcase_21 AC 314 ms
206,952 KB
testcase_22 AC 280 ms
189,596 KB
testcase_23 AC 622 ms
282,960 KB
testcase_24 AC 1,085 ms
387,912 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