結果
| 問題 |
No.1973 Divisor Sequence
|
| コンテスト | |
| ユーザー |
H20
|
| 提出日時 | 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)
H20