結果
問題 | No.1973 Divisor Sequence |
ユーザー |
![]() |
提出日時 | 2022-06-10 22:02:29 |
言語 | PyPy3 (7.3.15) |
結果 |
AC
|
実行時間 | 234 ms / 2,000 ms |
コード長 | 706 bytes |
コンパイル時間 | 211 ms |
コンパイル使用メモリ | 82,108 KB |
実行使用メモリ | 76,532 KB |
最終ジャッジ日時 | 2024-09-21 07:29:48 |
合計ジャッジ時間 | 3,032 ms |
ジャッジサーバーID (参考情報) |
judge2 / judge4 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 3 |
other | AC * 22 |
ソースコード
import collections 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 N,M = map(int, input().split()) ans = 1 mod = 10**9+7 CA = collections.Counter(prime_factorize(M)) CACA = collections.Counter() for k,v in CA.items(): CACA[v]+=1 for k,v in CACA.items(): DP = [1]*(k+1) for i in range(N-1): ACDP = list(DP) for j in range(1,k+1): ACDP[j]=(ACDP[j]+ACDP[j-1])%mod DP = ACDP[::-1] ans = pow(sum(DP),v,mod)*ans%mod print(ans%mod)