結果

問題 No.2341 Triple Tree Query (Medium)
ユーザー SSRSSSRS
提出日時 2023-05-25 23:27:29
言語 C++14
(gcc 12.3.0 + boost 1.83.0)
結果
WA  
実行時間 -
コード長 4,055 bytes
コンパイル時間 2,461 ms
コンパイル使用メモリ 179,812 KB
実行使用メモリ 30,060 KB
最終ジャッジ日時 2023-08-28 00:45:24
合計ジャッジ時間 18,540 ms
ジャッジサーバーID
（参考情報） 		judge14 / judge13
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 WA -
testcase_01 WA -
testcase_02 RE -
testcase_03 RE -
testcase_04 RE -
testcase_05 RE -
testcase_06 RE -
testcase_07 RE -
testcase_08 RE -
testcase_09 RE -
testcase_10 RE -
testcase_11 RE -
testcase_12 RE -
testcase_13 RE -
testcase_14 RE -
testcase_15 RE -
testcase_16 RE -
testcase_17 RE -
testcase_18 RE -
testcase_19 RE -
testcase_20 RE -
testcase_21 RE -
testcase_22 RE -
testcase_23 RE -
testcase_24 RE -
testcase_25 RE -
testcase_26 RE -
testcase_27 RE -
testcase_28 RE -
testcase_29 RE -
testcase_30 RE -
testcase_31 RE -
testcase_32 RE -
testcase_33 RE -
testcase_34 RE -
testcase_35 RE -
testcase_36 RE -
testcase_37 RE -
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <bits/stdc++.h>
using namespace std;
const long long MOD = 998244353;
struct affine{
  long long a, b;
  affine(){
    a = 1;
    b = 0;
  }
  affine(int a, int b): a(a), b(b){
  }
};
affine composite(affine A, affine B){
  return affine(A.a * B.a % MOD, (A.b * B.a + B.b) % MOD);
}
int value(affine A, int x){
  return (A.a * x + A.b) % MOD;
}
template <typename T>
struct dual_segment_tree{
  int N;
  vector<T> ST;
  function<T(T, T)> f;
  T E;
  dual_segment_tree(int n, function<T(T, T)> f, T E): f(f), E(E){
    N = 1;
    while (N < n){
      N *= 2;
    }
    ST = vector<T>(N * 2 - 1, E);
  }
  void push(int i){
    if (i < N - 1){
      ST[i * 2 + 1] = f(ST[i * 2 + 1], ST[i]);
      ST[i * 2 + 2] = f(ST[i * 2 + 2], ST[i]);
      ST[i] = E;
    }
  }
  T operator [](int k){
    int v = 0;
    for (int i = N / 2; i >= 1; i >>= 1){
      push(v);
      if ((k & i) == 0){
        v = v * 2 + 1;
      } else {
        v = v * 2 + 2;
      }
    }
    return ST[v];
  }
  void range_apply(int L, int R, T x, int i, int l, int r){
    if (r <= L || R <= l){
    } else if (L <= l && r <= R){
      ST[i] = f(ST[i], x);
    } else {
      push(i);
      int m = (l + r) / 2;
      range_apply(L, R, x, i * 2 + 1, l, m);
      range_apply(L, R, x, i * 2 + 2, m, r);
    }
  }
  void range_apply(int L, int R, T x){
    range_apply(L, R, x, 0, 0, N);
  }
};
int main(){
  int N, Q;
  cin >> N >> Q;
  vector<vector<int>> E(N);
  for (int i = 0; i < N - 1; i++){
    int A, B;
    cin >> A >> B;
    A--;
    B--;
    E[A].push_back(B);
    E[B].push_back(A);
  }
  vector<int> X(N);
  for (int i = 0; i < N; i++){
    cin >> X[i];
  }
  vector<int> p(N, -1);
  vector<vector<int>> c(N);
  vector<int> b;
  queue<int> q;
  q.push(0);
  while (!q.empty()){
    int v = q.front();
    q.pop();
    b.push_back(v);
    for (int w : E[v]){
      if (w != p[v]){
        p[w] = v;
        c[v].push_back(w);
        q.push(w);
      }
    }
  }
  reverse(b.begin(), b.end());
  vector<int> sz(N, 1);
  for (int v : b){
    for (int &w : c[v]){
      sz[v] += sz[w];
      if (sz[w] > sz[c[v][0]]){
        swap(w, c[v][0]);
      }
    }
  }
  vector<int> in(N);
  in[0] = 0;
  int t = 1;
  vector<int> next(N);
  for (int i = 0; i < N; i++){
    next[i] = i;
  }
  auto dfs = [&](auto dfs, int v = 0) -> void {
    if (v != 0){
      if (in[p[v]] == in[v] - 1){
        next[v] = next[p[v]];
      }
    }
    if (!c[v].empty()){
      in[c[v][0]] = t;
      t++;
      dfs(dfs, c[v][0]);
      int cnt = c[v].size();
      for (int i = 1; i < cnt; i++){
        in[c[v][i]] = t;
        t++;
      }
      for (int i = 1; i < cnt; i++){
        dfs(dfs, c[v][i]);
      }
    }
  };
  dfs(dfs);
  dual_segment_tree<affine> ST(N, composite, affine());
  for (int i = 0; i < Q; i++){
    int T;
    cin >> T;
    if (T == 1){
      int V;
      cin >> V;
      V--;
      cout << value(ST[in[V]], X[V]) << endl;
    }
    if (T == 2){
      int V, K;
      long long B, C;
      cin >> V >> K >> B >> C;
      V--;
      if (V != 0){
        ST.range_apply(in[p[V]], in[p[V]] + 1, affine(B, C));
      }
      ST.range_apply(in[V], in[V] + 1, affine(B, C));
      if (!c[V].empty()){
        ST.range_apply(in[c[V][0]], in[c[V][0]] + 1, affine(B, C));
      }
      if (c[V].size() >= 2){
        ST.range_apply(in[c[V][1]], in[c[V].back()] + 1, affine(B, C));
      }
    }
    if (T == 3){
      int U, V;
      long long B, C;
      cin >> U >> V >> B >> C;
      U--;
      V--;
      while (true){
        if (in[U] > in[V]){
          swap(U, V);
        }
        if (next[U] == next[V]){
          ST.range_apply(in[U], in[V] + 1, affine(B, C));
          break;
        }
        ST.range_apply(in[next[V]], in[V] + 1, affine(B, C));
        V = p[next[V]];
      }
    }
    if (T == 4){
      int V;
      long long B, C;
      cin >> V >> B >> C;
      V--;
      ST.range_apply(in[V], in[V] + 1, affine(B, C));
      if (!c[V].empty()){
        ST.range_apply(in[c[V][0]], in[c[V][0]] + sz[V] - 1, affine(B, C));
      }
    }
  }
}
0