結果
問題 |
No.3026 Range LCM (Online Version)
|
ユーザー |
|
提出日時 | 2025-02-14 23:33:14 |
言語 | C++17 (gcc 13.3.0 + boost 1.87.0) |
結果 |
TLE
|
実行時間 | - |
コード長 | 1,585 bytes |
コンパイル時間 | 4,742 ms |
コンパイル使用メモリ | 259,436 KB |
実行使用メモリ | 265,556 KB |
最終ジャッジ日時 | 2025-02-14 23:36:02 |
合計ジャッジ時間 | 114,850 ms |
ジャッジサーバーID (参考情報) |
judge5 / judge1 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 1 |
other | AC * 10 TLE * 26 |
ソースコード
#include <bits/stdc++.h> #include <atcoder/all> using namespace std; using namespace atcoder; using ll = long long; using mint = atcoder::modint998244353; #define BL 447 mint op(mint lhs, mint rhs){ return lhs * rhs; } mint e(){return 1;} int op2(int lhs, int rhs){return max(lhs, rhs);} constexpr int e2(){return 1;} int main(){ ios::sync_with_stdio(false); cin.tie(0); int n; cin >> n; vector<mint> tmp(n, 1); constexpr int r = 200000; vector<int> id(r + 1), tb(r + 1); tb[0] = tb[1] = -1; vector<int> p; for(int i = 2; i <= r; i++){ if(tb[i] != 0) continue; if(i < BL){ id[i] = p.size(); p.emplace_back(i); } for(int j = i; j <= r; j += i){ if(tb[j] == 0) tb[j] = i; } } const int m = p.size(); vector<vector<int>> seg2(m, vector<int>(n, 1)); for(int i = 0; i < n; i++){ int v; cin >> v; while(v != 1){ if(tb[v] <= BL) seg2[id[tb[v]]][i] *= tb[v]; else tmp[i] *= tb[v]; v /= tb[v]; } } segtree<mint, op, e> seg(tmp); vector<segtree<int, op2, e2>> seg3(m); for(int i = 0; i < m; i++) seg3[i] = segtree<int, op2, e2>(seg2[i]); int Q; cin >> Q; mint ans = 1; while(Q--){ int a, b; cin >> a >> b; int l = (a * ans).val() % n, r = (b * ans).val() % n; if(l > r) swap(l, r); r++; ans = seg.prod(l, r); for(int j = 0; j < m; j++) ans *= seg3[j].prod(l, r); cout << ans.val() << '\n'; } }