結果
問題 |
No.1973 Divisor Sequence
|
ユーザー |
![]() |
提出日時 | 2022-07-03 16:01:15 |
言語 | PyPy3 (7.3.15) |
結果 |
AC
|
実行時間 | 352 ms / 2,000 ms |
コード長 | 750 bytes |
コンパイル時間 | 164 ms |
コンパイル使用メモリ | 82,044 KB |
実行使用メモリ | 76,404 KB |
最終ジャッジ日時 | 2024-11-29 09:39:50 |
合計ジャッジ時間 | 3,005 ms |
ジャッジサーバーID (参考情報) |
judge5 / judge2 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 3 |
other | AC * 22 |
ソースコード
from collections import Counter MOD = 10 ** 9 + 7 def prime_factorize(n): a = [] while n % 2 == 0: a.append(2) n //= 2 f = 3 while f * f <= n: if n % f == 0: a.append(f) n //= f else: f += 2 if n != 1: a.append(n) return a def n_fact(N): c = Counter(prime_factorize(N)) return c N, M = map(int, input().split()) X = n_fact(M) x2 = Counter(X.values()) ans = 1 for k, v in x2.items(): dp = [1] * (k + 1) for i in range(N-1): dp2 = [0] * (k+1) sump = 0 for j in range(k+1): sump += dp[j] dp2[k-j] = sump % MOD dp = dp2 ans *= pow(sum(dp),v,MOD) ans %= MOD print(ans)