結果
問題 |
No.1973 Divisor Sequence
|
ユーザー |
![]() |
提出日時 | 2022-06-10 21:39:59 |
言語 | PyPy3 (7.3.15) |
結果 |
TLE
(最新)
AC
(最初)
|
実行時間 | - |
コード長 | 1,050 bytes |
コンパイル時間 | 168 ms |
コンパイル使用メモリ | 81,908 KB |
実行使用メモリ | 76,128 KB |
最終ジャッジ日時 | 2024-09-21 07:27:48 |
合計ジャッジ時間 | 6,021 ms |
ジャッジサーバーID (参考情報) |
judge5 / judge3 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 3 |
other | AC * 21 TLE * 1 |
ソースコード
mod = 1000000007 def main(): import sys input = sys.stdin.readline def PrimeDecomposition(N): ret = {} n = int(N ** 0.5) for d in range(2, n + 1): while N % d == 0: if d not in ret: ret[d] = 1 else: ret[d] += 1 N //= d if N == 1: break if N != 1: ret[N] = 1 return ret N, M = map(int, input().split()) P = PrimeDecomposition(M) ans = 1 for p in P: K = P[p] dp = [0] * (K+1) dp[0] = 1 for _ in range(N): dp_new = [0] * (K+1) for k in range(K+1): for k2 in range(K+1): if k + k2 > K: break dp_new[k2] = (dp_new[k2] + dp[k]) % mod dp = dp_new S = 0 for d in dp: S = (S + d)%mod ans = (ans * S)%mod print(ans) if __name__ == '__main__': main()