#1193862 (C++23(gcc16)) No.2571 列辞書順列

• 問題: No.2571 列辞書順列
• ユーザー: detteiuu
• 提出日時: 2026-10-20 20:27:51
• 言語: C++23(gcc16) (gcc 16.2.0 + boost 1.92.0 + ACL)
• 結果: AC
• 実行時間: 437 ms / 3,000 ms
• コード長: 3,673 bytes
• コンパイル時間: 6,144 ms
• コンパイル使用メモリ: 396,440 KB
• 実行使用メモリ: 10,912 KB
• 最終ジャッジ日時: 2026-10-20 20:28:15
• 合計ジャッジ時間: 18,817 ms
• ジャッジサーバーID （参考情報）: judge6_0 / judge2_0

ファイルパターン	結果
sample	AC * 2
other	AC * 29

ソースコード

#ifndef ONLINE_JUDGE
#define _GLIBCXX_DEBUG
#endif
#include <bits/stdc++.h>
using namespace std;

#define pass (void)0
#define INF (1<<30)-1
#define INFLL (1LL<<60)-1
#define rep(i, n) for (int i = 0; i < (int)(n); i++)
#define repr(i, n) for (int i = (int)(n) - 1; i >= 0; i--)
#define rep2(i, a, b) for (int i = (int)(a); i < (int)(b); i++)
#define repr2(i, a, b) for (int i = (int)(b) - 1; i >= (int)(a); i--)
#define all(x) (x).begin(), (x).end()
#define rall(x) (x).rbegin(), (x).rend()
#define sz(x) ((int)(x).size())
#define YesNo(cond) cout << ((cond) ? "Yes\n" : "No\n")
#define YESNO(cond) cout << ((cond) ? "YES\n" : "NO\n")

using ll = long long;
using pii = pair<int,int>;
using pll = pair<ll,ll>;
using vi = vector<int>;
using vl = vector<ll>;
using vvi = vector<vi>;
using vvl = vector<vl>;

template <typename T> void print(const T& value) { cout << value << "\n"; }
template <typename T> void print(const vector<T>& vec) { for (auto& v : vec) cout << v << " "; cout << "\n"; }
template <typename T> void input(vector<T>& vec) { for (auto& v : vec) cin >> v; };
template <typename T> bool chmin(T& a, const T& b) { if (a > b) { a = b; return true; } return false; }
template <typename T> bool chmax(T& a, const T& b) { if (a < b) { a = b; return true; } return false; }

#include <atcoder/all>
using namespace atcoder;
using mint = modint998244353;
using vm = vector<mint>;

vm C;
vm D;

struct Node {
    mint val;
    mint val2;
    ll size;
};

Node op(Node left, Node right) {
    if (left.size == 0) return right;
    if (right.size == 0) return left;
    return {left.val+left.val2*(D[right.size]-1)+right.val,
            left.val2*C[right.size]+right.val2,
            left.size+right.size};
}

Node e() { return {1, 0, 0}; }

mint pow2(mint x, long long n) {
    mint ret = 1;
    while (n > 0) {
        if (n & 1) ret = ret * x;  // n の最下位bitが 1 ならば x^(2^i) をかける
        x = x * x;
        n >>= 1;  // n を1bit 左にずらす
    }
    return ret;
}

vector<vector<mint>> matrix_mul(const vector<vector<mint>> &a,
                                const vector<vector<mint>> &b) {
    int H1 = a.size();
    int W1 = a[0].size();
    int W2 = b[0].size();

    vector<vector<mint>> ans(H1, vector<mint>(W2, 0));

    for (int i = 0; i < H1; i++) {
        for (int k = 0; k < W1; k++) {
            mint aik = a[i][k];
            for (int j = 0; j < W2; j++) {
                ans[i][j] += aik * b[k][j];
            }
        }
    }
    return ans;
}

int main() {
    ios::sync_with_stdio(false);
    cin.tie(nullptr);
    cout << fixed << setprecision(10);

    ll N, M, K, Q;
    cin >> N >> M >> K >> Q;
    vl A(N); input(A);

    C.push_back(1);
    rep (_, N) C.push_back(C.back()*M);
    D.push_back(1);
    rep2 (i, 1, N+1) D.push_back(D.back()+C[i]);

    vector<Node> def(N);
    rep (i, N) def[i] = {A[i], A[i]-1, 1};
    segtree<Node, op, e> seg(def);
    while (Q--) {
        ll t; cin >> t;
        if (t == 1) {
            ll l, r; cin >> l >> r;
            l--; r--;
            auto res = seg.prod(l, r+1);
            mint ans = res.val;
            if (r-l+1 < K) {
                vector<vm> mat = {{M, 0}, {1, 1}};
                vector<vm> ansM = {{M}, {0}};
                rep (i, 30) {
                    if ((1<<i) & (K-(r-l+1))) {
                        ansM = matrix_mul(mat, ansM);
                    }
                    mat = matrix_mul(mat, mat);
                }
                ans += res.val2*ansM[1][0];
            }
            print(ans.val());
        } else {
            ll x, y; cin >> x >> y;
            x--;
            seg.set(x, {y, y-1, 1});
        }
    }
}