結果

問題 No.2265 Xor Range Substring Sum Query
ユーザー SSRSSSRS
提出日時 2023-04-04 20:24:12
言語 C++14
(gcc 12.3.0 + boost 1.83.0)
結果
TLE  
実行時間 -
コード長 2,917 bytes
コンパイル時間 1,659 ms
コンパイル使用メモリ 181,264 KB
実行使用メモリ 305,776 KB
最終ジャッジ日時 2024-10-01 10:52:23
合計ジャッジ時間 8,634 ms
ジャッジサーバーID
（参考情報） 		judge4 / judge3
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 1 ms
5,248 KB
testcase_01 AC 1 ms
5,248 KB
testcase_02 AC 2 ms
5,248 KB
testcase_03 AC 3 ms
5,248 KB
testcase_04 TLE -
testcase_05 -- -
testcase_06 -- -
testcase_07 -- -
testcase_08 -- -
testcase_09 -- -
testcase_10 -- -
testcase_11 -- -
testcase_12 -- -
testcase_13 -- -
testcase_14 -- -
testcase_15 -- -
testcase_16 -- -
testcase_17 -- -
testcase_18 -- -
testcase_19 -- -
testcase_20 -- -
testcase_21 -- -
testcase_22 -- -
testcase_23 -- -
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <bits/stdc++.h>
using namespace std;
const long long MOD = 998244353;
struct node{
  long long x, p2, p11;
  node(): x(0), p2(1), p11(1){
  }
  node(int x): x(x), p2(2), p11(11){
  }
};
node op(node L, node R){
  node ans;
  ans.x = (L.x * R.p11 + L.p2 * R.x) % MOD;
  ans.p2 = L.p2 * R.p2 % MOD;
  ans.p11 = L.p11 * R.p11 % MOD;
  return ans;
}
node e(){
  return node();
}
template <typename T, T (*op)(T, T), T (*e)()>
struct xor_segment_tree{
  int N;
  vector<vector<T>> ST;
  xor_segment_tree(vector<T> &A){
    N = A.size();
    ST = vector<vector<T>>(N * 2 - 1);
    for (int i = 0; i < N; i++){
      ST[N - 1 + i].push_back(A[i]);
    }
    for (int i = N - 2; i >= 0; i--){
      int cnt = ST[i * 2 + 1].size();
      for (int j = 0; j < cnt; j++){
        ST[i].push_back(op(ST[i * 2 + 1][j], ST[i * 2 + 2][j]));
      }
      for (int j = 0; j < cnt; j++){
        ST[i].push_back(op(ST[i * 2 + 2][j], ST[i * 2 + 1][j]));
      }
    }
  }
  T range_fold(int L, int R, int x, int i, int l, int r){
    if (r <= L || R <= l){
      return e();
    } else if (L <= l && r <= R){
      return ST[i][x];
    } else {
      int p = (r - l) / 2;
      int m = (l + r) / 2;
      if ((x & p) == 0){
        T resL = range_fold(L, R, x, i * 2 + 1, l, m);
        T resR = range_fold(L, R, x, i * 2 + 2, m, r);
        return op(resL, resR);
      } else {
        T resL = e();
        if (R >= m){
          resL = range_fold(max(L, m) - p, R - p, x ^ p, i * 2 + 1, l, m);
        }
        T resR = e();
        if (L < m){
          resR = range_fold(L + p, min(R, m) + p, x ^ p, i * 2 + 2, m, r);
        }
        return op(resR, resL);
      }
    }
  }
  T range_fold(int L, int R, int x){
    return range_fold(L, R, x, 0, 0, N);
  }
};
int main(){
  ios::sync_with_stdio(false);
  cin.tie(nullptr);
  int n;
  cin >> n;
  string S;
  cin >> S;
  vector<node> A(1 << n);
  for (int i = 0; i < (1 << n); i++){
    A[i] = node(S[i] - '0');
  }
  int Q;
  cin >> Q;
  vector<int> pos;
  xor_segment_tree<node, op, e> ST(A);
  for (int i = 0; i < Q; i++){
    int t;
    cin >> t;
    if (t == 1){
      int x, y;
      cin >> x >> y;
      A[x] = node(y);
      pos.push_back(x);
    }
    if (t == 2){
      int L, R, X;
      cin >> L >> R >> X;
      R++;
      vector<int> pos2 = {L - 1, R};
      for (int j : pos){
        if (L <= (j ^ X) && (j ^ X) < R){
          pos2.push_back(j ^ X);
        }
      }
      sort(pos2.begin(), pos2.end());
      pos2.erase(unique(pos2.begin(), pos2.end()), pos2.end());
      int cnt2 = pos2.size() - 1;
      node ans;
      for (int j = 0; j < cnt2; j++){
        ans = op(ans, ST.range_fold(pos2[j] + 1, pos2[j + 1], X));
        if (j < cnt2 - 1){
          ans = op(ans, A[pos2[j + 1] ^ X]);
        }
      }
      cout << ans.x << '\n';
    }
    if (pos.size() == 50){
      ST = xor_segment_tree<node, op, e>(A);
      pos.clear();
    }
  }
}
0