結果
問題 | No.2318 Phys Bone Maker |
ユーザー |
|
提出日時 | 2023-05-26 21:43:42 |
言語 | C++23 (gcc 13.3.0 + boost 1.87.0) |
結果 |
AC
|
実行時間 | 926 ms / 3,000 ms |
コード長 | 1,293 bytes |
コンパイル時間 | 3,769 ms |
コンパイル使用メモリ | 252,220 KB |
実行使用メモリ | 5,248 KB |
最終ジャッジ日時 | 2024-12-25 05:56:16 |
合計ジャッジ時間 | 12,103 ms |
ジャッジサーバーID (参考情報) |
judge5 / judge1 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 3 |
other | AC * 45 |
ソースコード
#include <bits/stdc++.h> using namespace std; using i64 = int64_t; constexpr i64 mod = 998244353; int main() { ios::sync_with_stdio(false); cin.tie(nullptr); i64 n; cin >> n; vector<i64> p, d; for (i64 i = 1; i * i <= n; i += 1) { if (n % i == 0) { d.push_back(i); if (i * i != n) { d.push_back(n / i); } } } sort(d.begin(), d.end()); for (i64 i = 2; i * i <= n; i += 1) { if (n % i == 0) { p.push_back(i); while (n % i == 0) { n /= i; } } } if (n > 1) { p.push_back(n); } vector<i64> f(d.size()); f[0] = 1; for (int i = 0; i < d.size(); i += 1) { for (int j = 0; j < i; j += 1) { if (d[i] % d[j] == 0) { i64 x = d[i] / d[j], y = 1; for (i64 pi : p) { if (x % pi != 0) { int z = 1; for (i64 di = d[i]; di % pi == 0; di /= pi) { z += 1; } y *= z; } } f[i] = (y % mod * f[j] + f[i]) % mod; } } } cout << f.back(); }