結果
| 問題 |
No.1641 Tree Xor Query
|
| コンテスト | |
| ユーザー |
beet
|
| 提出日時 | 2021-08-06 21:44:16 |
| 言語 | C++23 (gcc 13.3.0 + boost 1.87.0) |
| 結果 |
AC
|
| 実行時間 | 100 ms / 5,000 ms |
| コード長 | 3,971 bytes |
| コンパイル時間 | 3,656 ms |
| コンパイル使用メモリ | 260,672 KB |
| 実行使用メモリ | 17,024 KB |
| 最終ジャッジ日時 | 2024-09-17 01:36:41 |
| 合計ジャッジ時間 | 4,560 ms |
|
ジャッジサーバーID (参考情報) |
judge2 / judge4 |
(要ログイン)
| ファイルパターン | 結果 |
|---|---|
| sample | AC * 3 |
| other | AC * 18 |
ソースコード
#include <bits/stdc++.h>
using namespace std;
using Int = long long;
const char newl = '\n';
template<typename T1,typename T2> inline void chmin(T1 &a,T2 b){if(a>b) a=b;}
template<typename T1,typename T2> inline void chmax(T1 &a,T2 b){if(a<b) a=b;}
template<typename T> void drop(const T &x){cout<<x<<endl;exit(0);}
template<typename T=int>
vector<T> read(size_t n){
vector<T> ts(n);
for(size_t i=0;i<n;i++) cin>>ts[i];
return ts;
}
class EulerTourForVertex{
private:
vector<int> ls,rs;
int pos;
void dfs(int v,int p){
ls[v]=pos++;
for(int u:G[v])
if(u!=p) dfs(u,v);
rs[v]=pos;
}
public:
vector< vector<int> > G;
EulerTourForVertex(int n):ls(n),rs(n),G(n){}
void add_edge(int u,int v){
G[u].emplace_back(v);
G[v].emplace_back(u);
}
void build(int r=0){
pos=0;
dfs(r,-1);
}
int idx(int v){return ls[v];}
template<typename F>
void exec(int v,F f){
f(ls[v],rs[v]);
}
};
template <typename T,typename E>
struct SegmentTree{
using F = function<T(T,T)>;
using G = function<T(T,E)>;
using H = function<E(E,E)>;
int n,height;
F f;
G g;
H h;
T ti;
E ei;
vector<T> dat;
vector<E> laz;
SegmentTree(F f,G g,H h,T ti,E ei):
f(f),g(g),h(h),ti(ti),ei(ei){}
void init(int n_){
n=1;height=0;
while(n<n_) n<<=1,height++;
dat.assign(2*n,ti);
laz.assign(2*n,ei);
}
void build(const vector<T> &v){
int n_=v.size();
init(n_);
for(int i=0;i<n_;i++) dat[n+i]=v[i];
for(int i=n-1;i;i--)
dat[i]=f(dat[(i<<1)|0],dat[(i<<1)|1]);
}
inline T reflect(int k){
return laz[k]==ei?dat[k]:g(dat[k],laz[k]);
}
inline void propagate(int k){
if(laz[k]==ei) return;
laz[(k<<1)|0]=h(laz[(k<<1)|0],laz[k]);
laz[(k<<1)|1]=h(laz[(k<<1)|1],laz[k]);
dat[k]=reflect(k);
laz[k]=ei;
}
inline void thrust(int k){
for(int i=height;i;i--) propagate(k>>i);
}
inline void recalc(int k){
while(k>>=1)
dat[k]=f(reflect((k<<1)|0),reflect((k<<1)|1));
}
void update(int a,int b,E x){
if(a>=b) return;
thrust(a+=n);
thrust(b+=n-1);
for(int l=a,r=b+1;l<r;l>>=1,r>>=1){
if(l&1) laz[l]=h(laz[l],x),l++;
if(r&1) --r,laz[r]=h(laz[r],x);
}
recalc(a);
recalc(b);
}
void set_val(int a,T x){
thrust(a+=n);
dat[a]=x;laz[a]=ei;
recalc(a);
}
T query(int a,int b){
if(a>=b) return ti;
thrust(a+=n);
thrust(b+=n-1);
T vl=ti,vr=ti;
for(int l=a,r=b+1;l<r;l>>=1,r>>=1) {
if(l&1) vl=f(vl,reflect(l++));
if(r&1) vr=f(reflect(--r),vr);
}
return f(vl,vr);
}
template<typename C>
int find(int st,C &check,T &acc,int k,int l,int r){
if(l+1==r){
acc=f(acc,reflect(k));
return check(acc)?k-n:-1;
}
propagate(k);
int m=(l+r)>>1;
if(m<=st) return find(st,check,acc,(k<<1)|1,m,r);
if(st<=l and !check(f(acc,dat[k]))){
acc=f(acc,dat[k]);
return -1;
}
int vl=find(st,check,acc,(k<<1)|0,l,m);
if(~vl) return vl;
return find(st,check,acc,(k<<1)|1,m,r);
}
template<typename C>
int find(int st,C &check){
T acc=ti;
return find(st,check,acc,1,0,n);
}
};
//INSERT ABOVE HERE
signed main(){
cin.tie(0);
ios::sync_with_stdio(0);
int n,q;
cin>>n>>q;
auto cs=read(n);
EulerTourForVertex G(n);
for(int i=1;i<n;i++){
int a,b;
cin>>a>>b;
a--;b--;
G.add_edge(a,b);
}
G.build(0);
using P = pair<int, int>;
auto f=[&](P a,P b){return P(a.first^b.first,a.second+b.second);};
auto g=[&](P a,int b){return P(a.first^(b*(a.second&1)),a.second);};
auto h=[&](int a,int b){return a^b;};
SegmentTree<P, int> seg(f,g,h,P(0,0),0);
vector<P> vp(n);
for(int i=0;i<n;i++) vp[G.idx(i)]=P(cs[i],1);
seg.build(vp);
for(int i=0;i<q;i++){
int t,x,y;
cin>>t>>x>>y;
x--;
if(t==1) seg.set_val(G.idx(x),P(cs[x]^=y,1));
if(t==2) G.exec(x,[&](int l,int r){cout<<seg.query(l,r).first<<newl;});
}
return 0;
}
beet