結果
問題 |
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)