結果
| 問題 | No.2318 Phys Bone Maker |
| コンテスト | |
| ユーザー |
|
| 提出日時 | 2023-05-26 22:45:18 |
| 言語 | PyPy3 (7.3.17) |
| 結果 |
TLE
|
| 実行時間 | - |
| コード長 | 759 bytes |
| 記録 | |
| コンパイル時間 | 151 ms |
| コンパイル使用メモリ | 85,376 KB |
| 実行使用メモリ | 387,328 KB |
| 最終ジャッジ日時 | 2026-05-31 01:02:56 |
| 合計ジャッジ時間 | 5,090 ms |
|
ジャッジサーバーID (参考情報) |
judge1_0 / judge3_1 |
(要ログイン)
| ファイルパターン | 結果 |
|---|---|
| sample | AC * 2 TLE * 1 |
| other | -- * 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)
pos = {d:i for i,d in enumerate(divs)}
divides = [[] for i in range(le)]
for i,d in enumerate(divs):
for j,nd in enumerate(divs):
lcm = d*nd // math.gcd(d,nd)
if lcm > d:
divides[i].append(pos[lcm])
# print(divides,divs)
# print("here")
dp = [0]*le
dp[0] = 1
ans = 0
for i in range(40):
ndp = [0]*le
for j in range(le):
if dp[j] == 0:
continue
for nj in divides[j]:
ndp[nj] += dp[j]
ndp[nj] %= mod
dp = ndp
ans += dp[-1]
ans %= mod
# print(dp)
print(ans)