結果
| 問題 |
No.2318 Phys Bone Maker
|
| コンテスト | |
| ユーザー |
|
| 提出日時 | 2023-05-26 22:45:18 |
| 言語 | PyPy3 (7.3.15) |
| 結果 |
TLE
|
| 実行時間 | - |
| コード長 | 759 bytes |
| コンパイル時間 | 820 ms |
| コンパイル使用メモリ | 82,176 KB |
| 実行使用メモリ | 334,080 KB |
| 最終ジャッジ日時 | 2024-12-25 09:11:27 |
| 合計ジャッジ時間 | 42,262 ms |
|
ジャッジサーバーID (参考情報) |
judge4 / judge5 |
(要ログイン)
| ファイルパターン | 結果 |
|---|---|
| sample | AC * 2 TLE * 1 |
| other | AC * 37 TLE * 8 |
ソースコード
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)