結果

問題 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)
הההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההה
XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX
0