結果

問題 No.1973 Divisor Sequence
ユーザー shobonvip
提出日時 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
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 3
other AC * 22
権限があれば一括ダウンロードができます

ソースコード

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