#line 1 "library/Template/template.hpp" #include using namespace std; #define rep(i,a,b) for(int i=(int)(a);i<(int)(b);i++) #define ALL(v) (v).begin(),(v).end() using ll=long long int; const int inf = 0x3fffffff; const ll INF = 0x1fffffffffffffff; templateinline bool chmax(T& a,T b){if(ainline bool chmin(T& a,T b){if(a>b){a=b;return 1;}return 0;} #line 2 "library/Utility/fastio.hpp" #include class FastIO{ static constexpr int L=1<<16; char rdbuf[L]; int rdLeft=0,rdRight=0; inline void reload(){ int len=rdRight-rdLeft; memmove(rdbuf,rdbuf+rdLeft,len); rdLeft=0,rdRight=len; rdRight+=fread(rdbuf+len,1,L-len,stdin); } inline bool skip(){ for(;;){ while(rdLeft!=rdRight and rdbuf[rdLeft]<=' ')rdLeft++; if(rdLeft==rdRight){ reload(); if(rdLeft==rdRight)return false; } else break; } return true; } template::value,int> =0>inline bool _read(T& x){ if(!skip())return false; if(rdLeft+20>=rdRight)reload(); bool neg=false; if(rdbuf[rdLeft]=='-'){ neg=true; rdLeft++; } x=0; while(rdbuf[rdLeft]>='0' and rdLeft::value,int> =0>inline bool _read(T& x){ if(!skip())return false; if(rdLeft+20>=rdRight)reload(); bool neg=false; if(rdbuf[rdLeft]=='-'){ neg=true; rdLeft++; } x=0; while(rdbuf[rdLeft]>='0' and rdbuf[rdLeft]<='9' and rdLeft='0' and rdbuf[rdLeft]<='9' and rdLeft=rdRight)reload(); x=rdbuf[rdLeft++]; return true; } inline bool _read(string& x){ if(!skip())return false; for(;;){ int pos=rdLeft; while(pos' ')pos++; x.append(rdbuf+rdLeft,pos-rdLeft); if(rdLeft==pos)break; rdLeft=pos; if(rdLeft==rdRight)reload(); else break; } return true; } templateinline bool _read(vector& v){ for(auto& x:v){ if(!_read(x))return false; } return true; } char wtbuf[L],tmp[50]; int wtRight=0; inline void flush(){ fwrite(wtbuf,1,wtRight,stdout); wtRight=0; } inline void _write(const char& x){ if(wtRight>L-32)flush(); wtbuf[wtRight++]=x; } inline void _write(const string& x){ for(auto& c:x)_write(c); } template::value,int> =0>inline void _write(T x){ if(wtRight>L-32)flush(); if(x==0){ _write('0'); return; } else if(x<0){ _write('-'); if (__builtin_expect(x == std::numeric_limits::min(), 0)) { switch (sizeof(x)) { case 2: _write("32768"); return; case 4: _write("2147483648"); return; case 8: _write("9223372036854775808"); return; } } x=-x; } int pos=0; while(x!=0){ tmp[pos++]=char((x%10)|48); x/=10; } rep(i,0,pos)wtbuf[wtRight+i]=tmp[pos-1-i]; wtRight+=pos; } templateinline void _write(const vector& v){ rep(i,0,v.size()){ if(i)_write(' '); _write(v[i]); } } public: FastIO(){} ~FastIO(){flush();} inline void read(){} template inline void read(Head& head,Tail&... tail){ assert(_read(head)); read(tail...); } templateinline void write(){if(ln)_write('\n');} template inline void write(const Head& head,const Tail&... tail){ if(space)_write(' '); _write(head); write(tail...); } }; /** * @brief Fast IO */ #line 3 "sol.cpp" #line 2 "library/Graph/hld.hpp" struct HLD{ using P=pair; vector> g; vector sz,in,out,rev,hs,par,dist; void dfs(int v,int p){ par[v]=p; sz[v]=1; if(p!=-1)dist[v]=dist[p]+1; if(!g[v].empty() and g[v][0]==p)swap(g[v][0],g[v].back()); for(auto& to:g[v])if(to!=p){ dfs(to,v); sz[v]+=sz[to]; if(sz[g[v][0]]in[v])swap(u,v); if(hs[u]==hs[v])return u; } } vector

get(int u,int p,bool es=0){ assert(in[p]<=in[u] and out[u]<=out[p]); vector

res; while(hs[u]!=hs[p]){ res.push_back({in[hs[u]],in[u]+1}); u=par[hs[u]]; } res.push_back({in[p]+es,in[u]+1}); return res; } int jump(int u,int v,int k){ if(k<0)return -1; int g=lca(u,v); int d0=dist[u]+dist[v]-dist[g]*2; if(d0=in[to])return rev[in[st]-k]; k-=in[st]-in[to]+1; st=par[to]; } } }; /** * @brief Heavy Light Decomposition */ #line 2 "library/DataStructure/lazysegtree.hpp" template class LazySegmentTree{ int sz,height; vector data; vector lazy; void update(int k){data[k]=f(data[k*2],data[k*2+1]);} void apply(int k,N x){ data[k]=g(data[k],x); if(k(n,m1())){} LazySegmentTree(const vector& a){ sz=1,height=0; while(sz<(int)a.size())sz<<=1,height++; data.assign(2*sz,m1()); lazy.assign(sz,n1()); rep(i,0,a.size())data[sz+i]=a[i]; for(int i=sz-1;i;i--)update(i); } void set(int k,M x){ k+=sz; for(int i=height;i;i--)down(k>>i); data[k]=x; for(int i=1;i<=height;i++)update(k>>i); } void update(int L,int R,N x){ if(L>=R)return; L+=sz,R+=sz; for(int i=height;i;i--){ if(((L>>i)<>i); if(((R>>i)<>i); } int lb=L,rb=R; while(L>=1; R>>=1; } L=lb,R=rb; for(int i=1;i<=height;i++){ if(((L>>i)<>i); if(((R>>i)<>i); } } M query(int L,int R){ if(L>=R)return m1(); L+=sz,R+=sz; for(int i=height;i;i--){ if(((L>>i)<>i); if(((R>>i)<>i); } M lb=m1(),rb=m1(); while(L>=1; R>>=1; } return f(lb,rb); } }; /** * @brief Lazy Segment Tree */ #line 2 "library/DataStructure/bit.hpp" templatestruct BIT{ int n; T all=0; vector val; BIT(int _n=0):n(_n),val(_n+10){} void clear(){val.assign(n+10,0); all=T();} void add(int i,T x){ for(i++;i<=n;i+=(i&-i))val[i]=val[i]+x; all+=x; } T sum(int i){ T res=0; for(;i;i-=(i&-i))res+=val[i]; return res; } T sum(int L,int R){return sum(R)-sum(L);} // [L,R) int lower_bound(T x){ int ret=0,len=1; while(2*len<=n)len<<=1; for(;len>=1;len>>=1){ if(ret+len<=n and val[ret+len] cnt(n),sum(n); vector tmp(n); tmp[0]=M(1); rep(i,1,n){ tmp[hld.in[i]]=M((i+1)-(hld.par[i]+1)); } LazySegmentTree seg(tmp); vector S(n); auto change=[&](int x)->void{ if(S[x]==0){ S[x]=1; cnt.add(hld.in[x],1); sum.add(hld.in[x],x+1); auto ps=hld.get(x,0); for(auto& [L,R]:ps)seg.update(L,R,1); } else{ S[x]=0; cnt.add(hld.in[x],-1); sum.add(hld.in[x],-(x+1)); auto ps=hld.get(x,0); for(auto& [L,R]:ps)seg.update(L,R,-1); } }; auto fixroot=[&](int x)->ll{ ll ret=seg.query(hld.in[x],hld.out[x]).a2bsum; ret-=sum.sum(hld.in[x],hld.out[x]); if(x!=0){ ll sz=cnt.sum(hld.in[x],hld.out[x]); ret+=(hld.par[x]+1)*sz*sz; } assert(ret%2==0); return ret/2; }; auto same=[&](int x)->ll{ ll ret=fixroot(0); auto ps=hld.get(x,0,1); M done=M(); for(auto& [L,R]:ps){ done=f(done,seg.query(L,R)); } ret+=done.absum*cnt.all; ret-=done.a2bsum; return ret; }; while(q--){ int u,r,v; io.read(u,r,v); u--; r--; v--; change(u); ll ret; if(r==v)ret=same(r); else if(hld.in[r]