結果
問題 |
No.1973 Divisor Sequence
|
ユーザー |
![]() |
提出日時 | 2023-04-23 11:33:10 |
言語 | PyPy3 (7.3.15) |
結果 |
AC
|
実行時間 | 513 ms / 2,000 ms |
コード長 | 859 bytes |
コンパイル時間 | 884 ms |
コンパイル使用メモリ | 82,176 KB |
実行使用メモリ | 76,416 KB |
最終ジャッジ日時 | 2024-11-07 19:14:00 |
合計ジャッジ時間 | 5,888 ms |
ジャッジサーバーID (参考情報) |
judge1 / judge5 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 3 |
other | AC * 22 |
ソースコード
def factorization(n): arr = [] temp = n for i in range(2, int(-(-n**0.5//1))+1): if temp%i==0: cnt=0 while temp%i==0: cnt+=1 temp //= i arr.append([i, cnt]) if temp!=1: arr.append([temp, 1]) return arr def solve(K): pre = [0] * (K + 1) pre[0] = 1 Acp = [1] * (K + 2) Acp[0] = 0 for i in range(N): dp = [0] * (K + 1) Ac = [0] * (K + 2) for j in range(K + 1): dp[j] += Acp[K + 1 - j] Ac[j + 1] = Ac[j] + dp[j] dp[j] %= mod Ac[j + 1] %= mod dp, pre = pre, dp Ac, Acp = Acp, Ac return sum(pre) % mod N, M = map(int, input().split()) mod = 10 ** 9 + 7 ans = 1 for p, m in factorization(M): ans *= solve(m) ans %= mod print(ans)