結果

問題 No.1641 Tree Xor Query
ユーザー vwxyz
提出日時 2024-12-06 11:59:10
言語 C++23
(gcc 13.3.0 + boost 1.87.0)
結果
AC  
実行時間 400 ms / 5,000 ms
コード長 7,730 bytes
コンパイル時間 4,163 ms
コンパイル使用メモリ 274,544 KB
実行使用メモリ 39,936 KB
最終ジャッジ日時 2024-12-06 11:59:17
合計ジャッジ時間 6,363 ms
ジャッジサーバーID
（参考情報） 		judge2 / judge1
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 3
other AC * 18
権限があれば一括ダウンロードができます

ソースコード

diff #
プレゼンテーションモードにする

#include <bits/stdc++.h>
using namespace std;
using uint = unsigned int;
using ll = long long;
using ull = unsigned long long;
template<class T> using V = vector<T>;
template<class T> using VV = V<V<T>>;
constexpr ll TEN(int n) { return (n == 0) ? 1 : 10 * TEN(n-1); }
#define FOR(i, a, b) for(int i=(int)(a);i<(int)(b);i++)
#define rep(i,N) for(int i=0;i<(int)(N);i++)
#define rep1(i,N) for(int i=1;i<=(int)(N);i++)
#define fs first
#define sc second
#define eb emplace_back
#define pb eb
#define all(x) x.begin(),x.end()
template<class T, class U> void chmin(T& t, const U& u) { if (t > u) t = u; }
template<class T, class U> void chmax(T& t, const U& u) { if (t < u) t = u; }
// bit op
int popcnt(uint x) { return __builtin_popcount(x); }
int popcnt(ull x) { return __builtin_popcountll(x); }
int bsr(uint x) { return 31 - __builtin_clz(x); }
int bsr(ull x) { return 63 - __builtin_clzll(x); }
int bsf(uint x) { return __builtin_ctz(x); }
int bsf(ull x) { return __builtin_ctzll(x); }
template <class S, S (*op)(S, S), S (*e)(), class F,
          S (*mapping)(F, S), F (*composition)(F, F), F (*id)()>
struct lazy_segtree {
  int _n, size, log;
  vector<S> d;
  vector<F> lz;
  void update(int k) {
    d[k] = op(d[2 * k], d[2 * k + 1]);
  }
  void all_apply(int k, F f) {
    d[k] = mapping(f, d[k]);
    if (k < size) lz[k] = composition(f, lz[k]);
  }
  void push(int k) {
    all_apply(2 * k, lz[k]);
    all_apply(2 * k + 1, lz[k]);
    lz[k] = id();
  }
  lazy_segtree() : lazy_segtree(0) {
  }
  explicit lazy_segtree(int n)
      : lazy_segtree(vector<S>(n, e())) {
  }
  explicit lazy_segtree(const vector<S>& v)
      : _n(int(v.size())) {
    size = bit_ceil((unsigned int)(_n));
    log = countr_zero((unsigned int)size);
    d = vector<S>(2 * size, e());
    lz = vector<F>(size, id());
    for (int i = 0; i < _n; i++) d[size + i] = v[i];
    for (int i = size - 1; i >= 1; i--) {
      update(i);
    }
  }
  void set(int p, S x) {
    p += size;
    for (int i = log; i >= 1; i--) push(p >> i);
    d[p] = x;
    for (int i = 1; i <= log; i++) update(p >> i);
  }
  S get(int p) {
    p += size;
    for (int i = log; i >= 1; i--) push(p >> i);
    return d[p];
  }
  S prod(int l, int r) {
    if (l == r) return e();
    l += size;
    r += size;
    for (int i = log; i >= 1; i--) {
      if (((l >> i) << i) != l) push(l >> i);
      if (((r >> i) << i) != r) push((r - 1) >> i);
    }
    S sml = e(), smr = e();
    while (l < r) {
      if (l & 1) sml = op(sml, d[l++]);
      if (r & 1) smr = op(d[--r], smr);
      l >>= 1;
      r >>= 1;
    }
    return op(sml, smr);
  }
  S all_prod() {
    return d[1];
  }
  void apply(int p, F f) {
    p += size;
    for (int i = log; i >= 1; i--) push(p >> i);
    d[p] = mapping(f, d[p]);
    for (int i = 1; i <= log; i++) update(p >> i);
  }
  void apply(int l, int r, F f) {
    if (l == r) return;
    l += size;
    r += size;
    for (int i = log; i >= 1; i--) {
      if (((l >> i) << i) != l) push(l >> i);
      if (((r >> i) << i) != r) push((r - 1) >> i);
    }
    {
      int l2 = l, r2 = r;
      while (l < r) {
        if (l & 1) all_apply(l++, f);
        if (r & 1) all_apply(--r, f);
        l >>= 1;
        r >>= 1;
      }
      l = l2;
      r = r2;
    }
    for (int i = 1; i <= log; i++) {
      if (((l >> i) << i) != l) update(l >> i);
      if (((r >> i) << i) != r) update((r - 1) >> i);
    }
  }
  template <bool (*g)(S)>
  int max_right(int l) {
    return max_right(l, [](S x) { return g(x); });
  }
  template <class G>
  int max_right(int l, G g) {
    if (l == _n) return _n;
    l += size;
    for (int i = log; i >= 1; i--) push(l >> i);
    S sm = e();
    do {
      while (l % 2 == 0) l >>= 1;
      if (!g(op(sm, d[l]))) {
        while (l < size) {
          push(l);
          l = (2 * l);
          if (g(op(sm, d[l]))) {
            sm = op(sm, d[l]);
            l++;
          }
        }
        return l - size;
      }
      sm = op(sm, d[l]);
      l++;
    } while ((l & -l) != l);
    return _n;
  }
  template <bool (*g)(S)>
  int min_left(int r) {
    return min_left(r, [](S x) { return g(x); });
  }
  template <class G>
  int min_left(int r, G g) {
    if (r == 0) return 0;
    r += size;
    for (int i = log; i >= 1; i--) push((r - 1) >> i);
    S sm = e();
    do {
      r--;
      while (r > 1 && (r % 2)) r >>= 1;
      if (!g(op(d[r], sm))) {
        while (r < size) {
          push(r);
          r = (2 * r + 1);
          if (g(op(d[r], sm))) {
            sm = op(d[r], sm);
            r--;
          }
        }
        return r + 1 - size;
      }
      sm = op(d[r], sm);
    } while ((r & -r) != r);
    return 0;
  }
};
//頂点vに関する情報はgrp[v]番目のセグ木のpos[v]番目に入れる
//g番目のセグ木の長さはgsz[g]
struct HL_Decomp{
  int N,gord;
  vector<vector<int>> tr;
  vector<int> par,dep,sub,grp,pos,gsz,gpar,gdep;
  HL_Decomp(int N,vector<vector<int>> graph):gord(),par(N),dep(N),sub(N),grp(N),pos(N),gsz(N),gpar(N),gdep(N){
    tr=graph;
    dfs1(0,-1,0);
    dfs2(0,-1,gord++,0);
  }
  int dfs1(int v,int p,int d){
    par[v]=p;
    sub[v]=1;
    dep[v]=d;
    for(auto to:tr[v])if(to!=p)sub[v]+=dfs1(to,v,d+1);
    return sub[v];
  }
  void dfs2(int v,int p,int g,int d){
    grp[v]=g;
    if(gsz[g]==0){
      gpar[g]=p;
      gdep[g]=d;
    }
    pos[v]=gsz[g]++;
    pair<int,int> si(0,-1);
    for(auto to:tr[v])if(to!=p)chmax(si,pair<int,int>(sub[to],to));
    for(auto to:tr[v])if(to!=p){
      if(to==si.second)dfs2(to,v,g,d);
      else dfs2(to,v,gord++,d+1);
    }
  }
  int LCA(int a,int b){
    if(gdep[grp[a]]>gdep[grp[b]])swap(a,b);
    while(gdep[grp[a]]<gdep[grp[b]])b=gpar[grp[b]];
    while(grp[a]!=grp[b]){
      a=gpar[grp[a]];
      b=gpar[grp[b]];
    }
    if(dep[a]>dep[b])swap(a,b);
    return a;
  }
  //{g,l,r}はg番目のセグ木のl~rという意味
  //l<rなら[l:r)、l>rなら[r:l)の逆順
  //edge=trueの場合はlcaを含まない区間を返す
  vector<tuple<int,int,int>> query(int a,int b,bool edge=false){
    vector<tuple<int,int,int>> l,r;
    while(gdep[grp[a]]>gdep[grp[b]]){
      l.emplace_back(grp[a],pos[a]+1,0);
      a=gpar[grp[a]];
    }
    while(gdep[grp[a]]<gdep[grp[b]]){
      r.emplace_back(grp[b],0,pos[b]+1);
      b=gpar[grp[b]];
    }
    while(grp[a]!=grp[b]){
      l.emplace_back(grp[a],pos[a]+1,0);
      r.emplace_back(grp[b],0,pos[b]+1);
      a=gpar[grp[a]];
      b=gpar[grp[b]];
    }
    int e=edge;
    if(dep[a]>dep[b])l.emplace_back(grp[a],pos[a]+1,pos[b]+e);
    else if(dep[a]<dep[b])r.emplace_back(grp[b],pos[a]+e,pos[b]+1);
    else if(!edge)l.emplace_back(grp[a],pos[a]+1,pos[a]);
    reverse(r.begin(),r.end());
    l.insert(l.end(),r.begin(),r.end());
    return l;
  }
};
int op(int x,int y){return x^y;}
int e(){return 0;}
int mapping(int x,int a){return x^a;}
int composition (int a,int b){return a^b;}
int id(){return 0;}
int main(){
  int N,Q;cin>>N>>Q;
  vector<int> C(N);
  for(int i=0;i<N;i++){
    cin>>C[i];
  }
  vector<vector<int>> graph(N);
  for(int i=0;i<N-1;i++){
    int a,b;cin>>a>>b;
    a-=1;b-=1;
    graph[a].push_back(b);
    graph[b].push_back(a);
  }
  HL_Decomp HL(N,graph);
  vector<lazy_segtree<int,op,e,int,mapping,composition,id>> LST(N);
  for(int g=0;g<N;g++){
    LST[g]=lazy_segtree<int,op,e,int,mapping,composition,id>(HL.gsz[g]);
  }
  for(int x=0;x<N;x++){
    for(auto [g,l,r]:HL.query(x,0)){
      LST[g].apply(min(l,r),max(l,r),C[x]);
    }
  }
  while(Q--){
    int T,x,y;cin>>T>>x>>y;
    x-=1;
    if(T==1){
      for(auto [g,l,r]:HL.query(x,0)){
        LST[g].apply(min(l,r),max(l,r),y);
      }
    }
    else{
      int ans=LST[HL.grp[x]].prod(HL.pos[x],HL.pos[x]+1);
      cout<<ans<<"\n";
    }
  }
}
הההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההה
XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX
0