結果
| 問題 |
No.2318 Phys Bone Maker
|
| コンテスト | |
| ユーザー |
|
| 提出日時 | 2023-05-27 10:20:54 |
| 言語 | PyPy3 (7.3.15) |
| 結果 |
AC
|
| 実行時間 | 921 ms / 3,000 ms |
| コード長 | 974 bytes |
| コンパイル時間 | 221 ms |
| コンパイル使用メモリ | 82,304 KB |
| 実行使用メモリ | 78,208 KB |
| 最終ジャッジ日時 | 2024-12-25 20:18:09 |
| 合計ジャッジ時間 | 11,516 ms |
|
ジャッジサーバーID (参考情報) |
judge1 / judge2 |
(要ログイン)
| ファイルパターン | 結果 |
|---|---|
| sample | AC * 3 |
| other | AC * 45 |
ソースコード
import math
n = int(input())
mod = 998244353
divs = set()
for i in range(1,int(n**0.5)+1):
if n%i:
continue
divs.add(i)
divs.add(n//i)
divs = sorted(divs)
le = len(divs)
primes = []
for i in range(2,int(n**0.5)+1):
if n%i:
continue
primes.append(i)
while n % i == 0:
n //= i
if n != 1:
primes.append(n)
div_count = [[0]*len(primes) for i in range(le)]
for i,div in enumerate(divs):
now = div
for j,prime in enumerate(primes):
while now%prime == 0:
now //= prime
div_count[i][j] += 1
dp = [0]*le
dp[0] = 1
lp = len(primes)
for i in range(le-1):
divi = div_count[i]
for j in range(i+1,le):
if divs[j]%divs[i]:
continue
divj = div_count[j]
pat = 1
for x in range(lp):
if divi[x] == divj[x]:
pat *= divi[x]+1
pat %= mod
dp[j] += dp[i]*pat
dp[j] %= mod
print(dp[-1])