#include using namespace std; using ll=long long; using ull=unsigned long long; using P=pair; templateusing minque=priority_queue,greater>; templatebool chmax(T &a,const T &b){return (abool chmin(T &a,const T &b){return (a>b?(a=b,true):false);} templateistream &operator>>(istream &is,pair&p){is>>p.first>>p.second;return is;} templateistream &operator>>(istream &is,tuple&a){is>>std::get<0>(a)>>std::get<1>(a)>>std::get<2>(a);return is;} templateistream &operator>>(istream &is,array&a){for(auto&i:a)is>>i;return is;} templateistream &operator>>(istream &is,vector &a){for(auto &i:a)is>>i;return is;} templatevoid operator++(pair&a,int n){a.first++,a.second++;} templatevoid operator--(pair&a,int n){a.first--,a.second--;} templatevoid operator++(vector&a,int n){for(auto &i:a)i++;} templatevoid operator--(vector&a,int n){for(auto &i:a)i--;} #define overload3(_1,_2,_3,name,...) name #define rep1(i,n) for(int i=0;i<(int)(n);i++) #define rep2(i,l,r) for(int i=(int)(l);i<(int)(r);i++) #define rep(...) overload3(__VA_ARGS__,rep2,rep1)(__VA_ARGS__) #define reps(i,l,r) rep2(i,l,r) #define all(x) x.begin(),x.end() #define pcnt(x) __builtin_popcountll(x) #define fin(x) return cout<<(x)<<'\n',static_cast(0) #define yn(x) cout<<((x)?"Yes\n":"No\n") #define uniq(x) sort(all(x)),x.erase(unique(all(x)),x.end()) template inline int fkey(vector&z,T key){return lower_bound(z.begin(),z.end(),key)-z.begin();} ll myceil(ll a,ll b){return (a+b-1)/b;} template auto vec(const int (&d)[n],const T &init=T()){ if constexpr (id(d,init)); else return init; } #ifdef LOCAL #include #define SWITCH(a,b) (a) #else #define debug(...) static_cast(0) #define debugg(...) static_cast(0) #define SWITCH(a,b) (b) templateostream &operator<<(ostream &os,const pair&p){os<>testcase; for(int i=0;i struct Edge{ int from,to; T weight; int index; Edge(int from_,int to_,T weight_=T(),int index_=-1):from(from_),to(to_),weight(weight_),index(index_){} Edge():from(-1),to(-1),weight(),index(-1){} friend std::ostream &operator<<(std::ostream &os,const Edge&e){ os<<'['; os<<"from:"< struct Tree{ int n,r; std::vector>edge; std::vector>g; std::vectorptr; struct tree_range{ using iterator=typename std::vector>::iterator; iterator l,r; iterator begin()const{return l;} iterator end()const{return r;} int size()const{return r-l;} bool empty()const{return !size();} Edge &operator[](int i)const{return l[i];} }; struct const_tree_range{ using iterator=typename std::vector>::const_iterator; iterator l,r; iterator begin()const{return l;} iterator end()const{return r;} int size()const{return r-l;} bool empty()const{return !size();} const Edge &operator[](int i)const{return l[i];} }; explicit Tree(int n_):n(n_){ edge.reserve(n-1); } Tree():n(0){} Tree(int n_,const std::vector>&e,bool dir=false):n(n_),r(-1),edge(e){ if(!dir)build(); else{ std::vectorseen(n,false); ptr.resize(n+1); for(const Edge&i:edge)ptr[i.from]++,ptr[i.to]++,seen[i.to]=true; for(int i=1;i<=n;i++)ptr[i]+=ptr[i-1]; r=std::find(seen.begin(),seen.end(),false)-seen.begin(); assert(ptr[n]==n*2-2); g.resize(ptr[n]); for(const Edge&i:edge)g[--ptr[i.to]]=Edge(i.to,i.from,i.weight,i.index); for(const Edge&i:edge)g[--ptr[i.from]]=i; } } template void read(){ for(int i=0;i>u>>v; if constexpr(index)u--,v--; if constexpr(weighted)std::cin>>w; else w=1; edge.emplace_back(u,v,w,i); } build(); } template void readp(){ ptr.resize(n+1); for(int i=1;i>p; if constexpr(index)p--; edge.emplace_back(p,i,1,i-1); ptr[p]++; ptr[i]++; } for(int i=1;i<=n;i++)ptr[i]+=ptr[i-1]; g.resize(n*2-2); for(auto&&[u,v,w,i]:edge)g[--ptr[v]]=Edge(v,u,w,i); for(int i=0;ipar(n,-1); par[root]=-1; std::queueque; que.push(root); while(!que.empty()){ int x=que.front(); que.pop(); for(int i=ptr[x];i&e=g[i]; if(e.to!=par[x]){ par[e.to]=x; assert(e.indexbfs_order()const{ assert(is_directed()); std::vectorbfs(n); int p=0,q=0; bfs[q++]=root(); while(p&e:(*this)[x])bfs[q++]=e.to; } return bfs; } std::vectordfs_order()const{ assert(is_directed()); std::vectorres; res.reserve(n); std::vectorst(n); int p=0; st[p++]=root(); while(p){ int x=st[--p]; res.push_back(x); p+=(*this)[x].size(); for(const Edge&e:(*this)[x])st[--p]=e.to; p+=(*this)[x].size(); } return res; } std::vectorrbfs_order()const{ std::vectorbfs=bfs_order(); std::reverse(bfs.begin(),bfs.end()); return bfs; } void hld(){ assert(is_directed()); std::vectorsub(n); for(int x:rbfs_order()){ sub[x]=1; int mx=-1; for(Edge&e:(*this)[x]){ sub[x]+=sub[e.to]; if(mx,std::vector>in_out_order(){ assert(is_directed()); std::vectorin(n),out(n); int p=0; auto dfs=[&](auto self,int x)->void { in[x]=p++; for(const Edge&e:(*this)[x]){ self(self,e.to); } out[x]=p; }; dfs(dfs,root()); return std::make_pair(in,out); } std::pair>diameter()const{ assert(!is_directed()); static constexpr T inf=std::numeric_limits::max(); std::vectordst(n,inf); dst[0]=0; std::vectorque(n); int p=0,q=1; que[0]=0; while(p&e:(*this)[x])if(dst[e.to]==inf){ dst[e.to]=dst[x]+e.weight; que[q++]=e.to; } } int u=std::max_element(dst.begin(),dst.end())-dst.begin(); std::fill(dst.begin(),dst.end(),inf); dst[u]=0; p=0,q=1; que[0]=u; while(p&e:(*this)[x])if(dst[e.to]==inf){ dst[e.to]=dst[x]+e.weight; que[q++]=e.to; } } int v=std::max_element(dst.begin(),dst.end())-dst.begin(); T weight=dst[v]; std::vectorres; while(u!=v){ res.push_back(v); for(const Edge&e:(*this)[v])if(dst[e.to]& get_edge(int i)const{return edge[i];} inline int parent(int i)const{return i==r?-1:g[ptr[i+1]-1].to;} const Edge& parent_edge(int i)const{ assert(r!=-1&&i!=r); return g[ptr[i+1]-1]; } inline int root()const{return r;} typename std::vector>::iterator begin(){return edge.begin();} typename std::vector>::iterator end(){return edge.end();} typename std::vector>::const_iterator begin()const{return edge.begin();} typename std::vector>::const_iterator end()const{return edge.end();} }; struct StaticTopTree{ private: template void build(Treet){ int n=t.size(); left.reserve(n*2-1),right.reserve(n*2-1),par.reserve(n*2-1),A.reserve(n*2-1),B.reserve(n*2-1); left.resize(n,-1),right.resize(n,-1),par.resize(n,-1),A.resize(n),B.resize(n); for(int i=0;istd::pair { std::vector>vs{{0,x}}; while(t[x].size()>=1){ std::priority_queue,std::vector>,std::greater>>que; int heavy=t[x][0].to; que.emplace(0,heavy); for(int i=1;i=2){ auto [d1,v1]=que.top();que.pop(); auto [d2,v2]=que.top();que.pop(); if(B[v2]==heavy)std::swap(d1,d2),std::swap(v1,v2); int nv=left.size(); left.push_back(v1),right.push_back(v2),par.push_back(-1),A.push_back(x),B.push_back(B[v1]); par[v1]=par[v2]=nv; que.emplace(std::max(d1,d2)+1,nv); } vs.push_back(que.top()); while(true){ int sz=vs.size(); if(sz>=3&&(vs[sz-3].first==vs[sz-2].first||vs[sz-3].first<=vs[sz-1].first)){ int nv=left.size(); left.push_back(vs[sz-3].second),right.push_back(vs[sz-2].second),par.push_back(-1),A.push_back(A[vs[sz-3].second]),B.push_back(B[vs[sz-2].second]); par[vs[sz-3].second]=par[vs[sz-2].second]=nv; vs[sz-3].first=std::max(vs[sz-3].first,vs[sz-2].first)+1; vs[sz-3].second=nv; vs[sz-2]=vs[sz-1]; vs.pop_back(); } else if(sz>=2&&vs[sz-2].first<=vs[sz-1].first){ int nv=left.size(); left.push_back(vs[sz-2].second),right.push_back(vs[sz-1].second),par.push_back(-1),A.push_back(A[vs[sz-2].second]),B.push_back(B[vs[sz-1].second]); par[vs[sz-2].second]=par[vs[sz-1].second]=nv; vs[sz-2].first=std::max(vs[sz-2].first,vs[sz-1].first)+1; vs[sz-2].second=nv; vs.pop_back(); } else break; } x=heavy; } while((int)vs.size()>=2){ int sz=vs.size(); int nv=left.size(); left.push_back(vs[sz-2].second),right.push_back(vs[sz-1].second),par.push_back(-1),A.push_back(A[vs[sz-2].second]),B.push_back(B[vs[sz-1].second]); par[vs[sz-2].second]=par[vs[sz-1].second]=nv; vs[sz-2].first=std::max(vs[sz-2].first,vs[sz-1].first)+1; vs[sz-2].second=nv; vs.pop_back(); } return vs[0]; }; dfs(dfs,t.root()); } public: std::vectorleft,right,par,A,B; StaticTopTree(){} template explicit StaticTopTree(Treet,int){ build(std::move(t)); } template explicit StaticTopTree(Treet){ assert(t.is_directed()); t.hld(); build(std::move(t)); } }; template struct DynamicTreeDP{ private: using Data=typename DP::Data; std::vectordata; StaticTopTree stt; inline void update(int x){ int l=stt.left[x],r=stt.right[x]; if(stt.A[l]==stt.A[r])data[x]=DP::rake(data[l],data[r]); else data[x]=DP::compress(data[l],data[r]); } public: DynamicTreeDP(){} template explicit DynamicTreeDP(Treet,std::vectorinit){ assert(t.size()==(int)init.size()); int n=t.size(); stt=StaticTopTree(std::move(t)); data=std::move(init); data.resize(n*2-1); for(int i=n;i>n>>r>>b; Tree t(n); double phi=(sqrt(5)+1)/2; rep(i,1,n+1){ int j=i*phi; if(iinit(n); vectorcol(n); rep(i,n){ int c=0; if(i!=0)c=t.parent_edge(i).weight; init[i]=M::Data(0,c); } DynamicTreeDPdp(t,init); int q; cin>>q; while(q--){ int type,x; cin>>type>>x; if(type==1){ x--; col[x]^=1; int c=0; if(x!=0)c=t.parent_edge(x).weight; init[x]=M::Data(col[x],c); dp.set(x,init[x]); } else if(type==2)r=x; else b=x; debug(dp.get().red.sum,dp.get().blue.sum); cout<