結果

問題 No.3026 Range LCM (Online Version)
ユーザー t98slider
提出日時 2025-02-14 23:56:21
言語 C++17
(gcc 13.3.0 + boost 1.87.0)
結果
MLE  
実行時間 -
コード長 2,538 bytes
コンパイル時間 4,313 ms
コンパイル使用メモリ 259,952 KB
実行使用メモリ 814,712 KB
最終ジャッジ日時 2025-02-14 23:57:23
合計ジャッジ時間 44,016 ms
ジャッジサーバーID
（参考情報） 		judge7 / judge3
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 1
other AC * 10 MLE * 26
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <bits/stdc++.h>
#include <atcoder/all>
using namespace std;
using namespace atcoder;
using ll = long long;
using mint = atcoder::modint998244353;
#define BL 447
// 59 ~ 447

template <class T, T (*op)(T, T), T (*e)()> struct SparseTable {
    int n, LOGV;
    std::vector<std::vector<T>> table;
    std::vector<int> log_table;
    SparseTable() {}
    SparseTable(const std::vector<T> &v) : n(v.size()), LOGV(std::__lg(std::max(n, 1)) + 1), 
                                           log_table(n + 1), table(LOGV, std::vector<T>(n)) {
        std::copy(v.begin(), v.end(), table[0].begin());
        for(int i = 1; i < LOGV; i++) {
            for(int j = 0; j + (1 << i) <= n; j++) {
                table[i][j] = op(table[i - 1][j], table[i - 1][j + (1 << (i - 1))]);
            }
        }
        for(int i = 2; i <= n; i++) {
            log_table[i] = log_table[i >> 1] + 1;
        }
    }
    T prod(int l, int r){
        assert(0 <= l && l <= r && r <= n);
        if(l == r) return e();
        int b = log_table[r - l];
        return op(table[b][l], table[b][r - (1 << b)]);
    }
};

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<SparseTable<int, op2, e2>> seg3;
    for(int i = 0; i < m; i++) seg3.emplace_back(SparseTable<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';
    }
}
0