QTREE

  1. /*
  2.     题目大意：维护一棵树，允许修改边权以及查询链上最大值
  3.     题解：我们将边权转为点权，标记在深度较深的点上，树链剖分后用线段树处理即可
  4. */
  5. #include <cstdio>
  6. #include <cstring>
  7. #include <algorithm>
  8. using namespace std;
  9. const int N=20010;
  10. int top[N],nxt[N],v[N],g[N],d[N],son[N],f[N],size[N],st[N],en[N],dfn,ed;
  11. void add(int x,int y){v[++ed]=y;nxt[ed]=g[x];g[x]=ed;}
  12. void init(){
  13.     memset(g,dfn=ed=0,sizeof(g));
  14.     memset(v,0,sizeof(v));
  15.     memset(nxt,0,sizeof(nxt));
  16.     memset(son,-1,sizeof(son));
  17. }
  18. void dfs(int x){
  19.     size[x]=1;
  20.     for(int i=g[x];i;i=nxt[i])if(v[i]!=f[x]){
  21.         f[v[i]]=x,d[v[i]]=d[x]+1;
  22.         dfs(v[i]),size[x]+=size[v[i]];
  23.         if(size[v[i]]>size[son[x]])son[x]=v[i];
  24.     }
  25. }
  26. void dfs2(int x,int y){
  27.     st[x]=++dfn;top[x]=y;
  28.     if(son[x])dfs2(son[x],y);
  29.     for(int i=g[x];i;i=nxt[i])if(v[i]!=son[x]&&v[i]!=f[x])dfs2(v[i],v[i]);
  30.     en[x]=dfn;
  31. }
  32. struct Node{int l,r;int Max;}Tree[N*3];
  33. void build(int x,int l,int r){
  34.     Tree[x].l=l;Tree[x].r=r;Tree[x].Max=0;
  35.     if(l==r)return;int mid=(l+r)>>1;
  36.     build(x<<1,l,mid);build((x<<1)|1,mid+1,r);
  37. }
  38. void up(int x){Tree[x].Max=max(Tree[x<<1].Max,Tree[(x<<1)|1].Max);}
  39. void update(int x,int k,int val){
  40.     if(Tree[x].l==k&&Tree[x].r==k){Tree[x].Max=val;return;}
  41.     int mid=(Tree[x].l+Tree[x].r)>>1;
  42.     if(k<=mid)update(x<<1,k,val);
  43.     else update((x<<1)|1,k,val);
  44.     up(x);
  45. }
  46. int query(int x,int l,int r){
  47.     if(Tree[x].l==l&&Tree[x].r==r)return Tree[x].Max;
  48.     int mid=(Tree[x].l+Tree[x].r)>>1;
  49.     if(r<=mid)return query(x<<1,l,r);
  50.     else if(l>mid)return query((x<<1)|1,l,r);
  51.     else return max(query(x<<1,l,mid),query((x<<1)|1,mid+1,r));
  52. }
  53. int find(int x,int y){
  54.     int tmp=0;
  55.     for(;top[x]!=top[y];x=f[top[x]]){
  56.         if(d[top[x]]<d[top[y]]){int z=x;x=y;y=z;}
  57.         tmp=max(tmp,query(1,st[top[x]],st[x]));
  58.     }if(x==y)return tmp;
  59.     if(d[x]>d[y]){int z=x;x=y;y=z;}
  60.     return max(tmp,query(1,st[son[x]],st[y]));
  61. }int e[N][3];
  62. int main(int n,int T){
  63.     scanf("%d",&T);
  64.     while(T--){
  65.         init(); scanf("%d",&n);
  66.         for(int i=0;i<n-1;i++){
  67.             scanf("%d%d%d",&e[i][0],&e[i][1],&e[i][2]);
  68.             add(e[i][0],e[i][1]);
  69.             add(e[i][1],e[i][0]);
  70.         }dfs(1);dfs2(1,1);
  71.         build(1,1,dfn);
  72.         for(int i=0;i<n-1;i++){
  73.             if(d[e[i][0]]>d[e[i][1]])swap(e[i][0],e[i][1]);
  74.             update(1,st[e[i][1]],e[i][2]);
  75.         }char op[10]; int u,w;
  76.         while(scanf("%s",op)==1){
  77.             if(op[0]=='D')break;
  78.             scanf("%d%d",&u,&w);
  79.             if(op[0]=='Q')printf("%d\n",find(u,w));
  80.             else update(1,st[e[u-1][1]],w);
  81.         }
  82.     }return 0;
  83. }

QTREE2

  1. /*
  2.     题目大意：给出一棵边权树，要求实现树链求和以及求链上第k个元素
  3.     题解：边权转点权后用LCT维护，注意求Qkth和sum时区别于点权树的讨论
  4. */
  5. #include <cstdio>
  6. #include <algorithm>
  7. #include <vector>
  8. #include <cstring>
  9. using namespace std;
  10. const int N=300010;
  11. namespace Link_Cut_Tree{
  12.     vector<int> v[N],w[N];
  13.     int f[N],son[N][2],tmp[N],size[N],val[N],sum[N];
  14.     void init(){
  15.         memset(val,0,sizeof(val));
  16.         memset(sum,0,sizeof(sum));
  17.         memset(f,0,sizeof(f));
  18.         memset(son,0,sizeof(son));
  19.         for(int i=1;i<N;i++)size[i]=1,v[i].clear(),w[i].clear();
  20.     }
  21.     bool isroot(int x){return !f[x]||son[f[x]][0]!=x&&son[f[x]][1]!=x;}
  22.     void up(int x){
  23.         sum[x]=val[x]+sum[son[x][0]]+sum[son[x][1]];
  24.         size[x]=size[son[x][0]]+size[son[x][1]]+1;
  25.     }
  26.     void rotate(int x){
  27.         int y=f[x],w=son[y][1]==x;
  28.         son[y][w]=son[x][w^1];
  29.         if(son[x][w^1])f[son[x][w^1]]=y;
  30.         if(f[y]){
  31.             int z=f[y];
  32.             if(son[z][0]==y)son[z][0]=x;else if(son[z][1]==y)son[z][1]=x;
  33.         }f[x]=f[y];f[y]=x;son[x][w^1]=y;up(y);
  34.     }
  35.     void splay(int x){
  36.         int s=1,i=x,y;tmp[1]=i;
  37.         while(!isroot(i))tmp[++s]=i=f[i];
  38.         while(!isroot(x)){
  39.             y=f[x];
  40.             if(!isroot(y)){if((son[f[y]][0]==y)^(son[y][0]==x))rotate(x);else rotate(y);}
  41.             rotate(x);
  42.         }up(x);
  43.     }
  44.     void access(int x){for(int y=0;x;y=x,x=f[x])splay(x),son[x][1]=y,up(x);}
  45.     // 求边权树的链和
  46.     int ask(int x,int y){
  47.         access(y),y=0;
  48.         int ans=0;
  49.         while(x){
  50.             splay(x);
  51.             if(!f[x])return sum[son[x][1]]+sum[y];
  52.         	son[x][1]=y; up(x);
  53.             x=f[y=x];
  54.         }
  55.     }
  56.     void dfs(int x,int fx){
  57.         for(int i=0;i<v[x].size();i++){
  58.             int y=v[x][i];
  59.             if(y==fx)continue;
  60.             f[y]=x;
  61.             val[y]=w[x][i];
  62.             dfs(y,x);
  63.         }up(x);
  64.     }
  65.     int kth(int x,int k){
  66.         while(size[son[x][0]]+1!=k){
  67.             if(k<size[son[x][0]]+1)x=son[x][0];
  68.             else{
  69.                 k-=size[son[x][0]]+1;
  70.                 x=son[x][1];
  71.             }
  72.         }return x;
  73.     }
  74.     // 求边权树u->v的第k个点
  75.     int Qkth(int u,int v,int k){
  76.         access(u),u=0;
  77.         while(v){
  78.             splay(v);
  79.             if(!f[v]){
  80.                 if(size[son[v][1]]+1==k)return v;
  81.                 else if(size[son[v][1]]+1>k)return kth(son[v][1],size[son[v][1]]+1-k);
  82.                 else return kth(u,k-(size[son[v][1]]+1));
  83.             }son[v][1]=u;
  84.             up(v); v=f[u=v];
  85.         }
  86.     }
  87. }
  88. using namespace Link_Cut_Tree;
  89. int T,x,y,z,k,n;
  90. char op[10];
  91. int main(){
  92.     scanf("%d",&T);
  93.     while(T--){
  94.         scanf("%d",&n);
  95.         init();
  96.         for(int i=1;i<n;i++){
  97.             scanf("%d%d%d",&x,&y,&z);
  98.               v[x].push_back(y);
  99.               w[x].push_back(z);
  100.               v[y].push_back(x);
  101.               w[y].push_back(z);
  102.         }dfs(1,0);
  103.         while(scanf("%s",op)){
  104.             if(op[1]=='O')break;
  105.             if(op[0]=='D'){
  106.                 scanf("%d%d",&x,&y);
  107.                 printf("%d\n",ask(x,y));
  108.             }
  109.             else{
  110.                 scanf("%d%d%d",&x,&y,&k);
  111.                 printf("%d\n",Qkth(x,y,k));
  112.             }
  113.         }
  114.     }return 0;
  115. }

QTREE3

  1. /*
  2.     题目大意：给出一棵一开始都是白点的树，要求维护操作：
  3.         1.将某个点的颜色取反(黑白互换)
  4.         2.查询1到x的路径上第一个黑色的点
  5.     题解：LCT维护子树中黑点的数量，对于颜色取反，则直接flap，
  6.     对于查询操作通过access(x)将1到x这条链取出，splay上查找即可
  7.     注意是从1到x，所以是makeroot(1)，保证1节点的键值最小
  8. */
  9. #include <cstdio>
  10. #include <algorithm>
  11. #include <cstring>
  12. using namespace std;
  13. const int N=100010;
  14. namespace Link_Cut_Tree{
  15.     int f[N],son[N][2],val[N],sum[N],tmp[N];bool rev[N];
  16.     void init(){
  17.         memset(f,0,sizeof(f));
  18.         memset(son,0,sizeof(son));
  19.         memset(val,0,sizeof(val));
  20.         memset(rev,0,sizeof(rev));
  21.         memset(sum,0,sizeof(sum));
  22.     }
  23.     bool isroot(int x){return !f[x]||son[f[x]][0]!=x&&son[f[x]][1]!=x;}
  24.     void rev1(int x){if(!x)return;swap(son[x][0],son[x][1]);rev[x]^=1;}
  25.     void pb(int x){if(rev[x])rev1(son[x][0]),rev1(son[x][1]),rev[x]=0;}
  26.     void up(int x){sum[x]=sum[son[x][0]]+sum[son[x][1]]+val[x];}
  27.     void rotate(int x){
  28.         int y=f[x],w=son[y][1]==x;
  29.         son[y][w]=son[x][w^1];
  30.         if(son[x][w^1])f[son[x][w^1]]=y;
  31.         if(f[y]){
  32.             int z=f[y];
  33.             if(son[z][0]==y)son[z][0]=x;else if(son[z][1]==y)son[z][1]=x;
  34.         }f[x]=f[y];f[y]=x;son[x][w^1]=y;up(y);
  35.     }
  36.     void splay(int x){
  37.         int s=1,i=x,y;tmp[1]=i;
  38.         while(!isroot(i))tmp[++s]=i=f[i];
  39.         while(s)pb(tmp[s--]);
  40.         while(!isroot(x)){
  41.             y=f[x];
  42.             if(!isroot(y)){if((son[f[y]][0]==y)^(son[y][0]==x))rotate(x);else rotate(y);}
  43.             rotate(x);
  44.         }up(x);
  45.     }
  46.     void access(int x){for(int y=0;x;y=x,x=f[x])splay(x),son[x][1]=y,up(x);}
  47.     int root(int x){access(x);splay(x);while(son[x][0])x=son[x][0];return x;}
  48.     void makeroot(int x){access(x);splay(x);rev1(x);}
  49.     void link(int x,int y){makeroot(x);f[x]=y;access(x);}
  50.     int flip(int x){makeroot(x);val[x]^=1;up(x);}
  51.     // 查询从1->x的第一个黑点
  52.     int ask(int x){
  53.         makeroot(1);
  54.         access(x); splay(x); // splay只保证了x到根节点标记的下传
  55.         if(!sum[x])return -1;
  56.         while(x){
  57.             pb(x); // 因为有makeroot操作，所以需要下传标记
  58.             if(!sum[son[x][0]]&&val[x]==1)return x;
  59.             else if(sum[son[x][0]])x=son[x][0];
  60.             else x=son[x][1];
  61.         }
  62.     }
  63. }
  64. using namespace Link_Cut_Tree;
  65. int n,q,x,y;
  66. int main(){
  67.     while(~scanf("%d%d",&n,&q)){
  68.         init();
  69.         for(int i=1;i<n;i++){
  70.             scanf("%d%d",&x,&y);
  71.             link(x,y);
  72.         }
  73.         while(q--){
  74.             scanf("%d%d",&x,&y);
  75.             if(x)printf("%d\n",ask(y));
  76.             else flip(y);
  77.         }
  78.     }return 0;
  79. }

QTREE4

  1. /*
  2.     题目大意：给出一棵边权树，一开始所有点都是白点，
  3.     要求维护两种操作：
  4.         1.颜色取反（黑白互变）
  5.         2.查询树上最远的两个白点之间的距离
  6.     题解：我们维护重心树，对于每个重心维护分治层所有白点到其距离，
  7.     将其保存在对应方向子重心的优先队列中，
  8.     重心处另外维护一个队列保存该分治层每个子重心所产生的最大答案，即最大值加次大值，
  9.     那么每个重心处产生的最大答案的最值就是答案，
  10.     考虑修改问题，等价于在处理优先队列的删除问题，
  11.     对于每个需要删除操作的优先队列，我们额外创建一个优先队列将删除元素加入其中，
  12.     当两个队列top元素相同时同时删去即可。
  13. */
  14. #include <cstdio>
  15. #include <algorithm>
  16. #include <cstring>
  17. #include <queue>
  18. using namespace std;
  19. const int N=200010;
  20. const int INF=0x3f3f3f3f;
  21. int tot,g[N],nxt[N<<1],v[N<<1],w[N<<1];
  22. int TOT,G[N],NXT[N*20],V[N*20],W[N*20];
  23. void init(){
  24.     memset(g,-1,sizeof(g));tot=0;
  25.     memset(G,-1,sizeof(G));TOT=0;
  26. }
  27. void add_edge(int x,int y,int z){v[tot]=y,w[tot]=z,nxt[tot]=g[x],g[x]=tot++;}
  28. void ADD_EDGE(int x,int y,int z){V[TOT]=y,W[TOT]=z,NXT[TOT]=G[x],G[x]=TOT++;}
  29. int sum,root,size[N],dp[N],vis[N];
  30. void getroot(int x,int fx){
  31.     size[x]=1; dp[x]=0;
  32.     for(int i=g[x];~i;i=nxt[i]){
  33.         if(!vis[v[i]]&&v[i]!=fx){
  34.             getroot(v[i],x);
  35.             size[x]+=size[v[i]];
  36.             dp[x]=max(dp[x],size[v[i]]);
  37.         }
  38.     }dp[x]=max(dp[x],sum-size[x]);
  39.     if(dp[x]<dp[root])root=x;
  40. }
  41. struct Que{
  42.     priority_queue<int> Q,D;
  43.     void clear(){
  44.         while(!Q.empty())Q.pop();
  45.         while(!D.empty())D.pop();
  46.     }
  47.     void add(int x){Q.push(x);}
  48.     void del(int x){D.push(x);}
  49.     int top(){
  50.         for(;;){
  51.             if(Q.empty())return -INF;
  52.             else if(!D.empty()&&Q.top()==D.top())Q.pop(),D.pop();
  53.             else return Q.top();
  54.         }
  55.     }
  56.     int toptwo(){
  57.         int x=top();
  58.         del(x);
  59.         int y=top();
  60.         add(x);
  61.         if(y==-INF)return x==-INF?x:0;
  62.         return max(x+y,0);
  63.     }
  64. }P[N],Q[N],ans;
  65. int dis[N],col[N];
  66. void getdis(int rt,int x,int fx){
  67.     Q[rt].add(dis[x]);
  68.     for(int i=g[x];~i;i=nxt[i]){
  69.         if(v[i]==fx||vis[v[i]])continue;
  70.         dis[v[i]]=dis[x]+w[i];
  71.         getdis(rt,v[i],x);
  72.     }
  73. }
  74. void getdeep(int rt,int x,int fx){
  75.     ADD_EDGE(x,rt,dis[x]);
  76.     for(int i=g[x];~i;i=nxt[i]){
  77.         if(v[i]==fx||vis[v[i]])continue;
  78.         dis[v[i]]=dis[x]+w[i];
  79.         getdeep(rt,v[i],x);
  80.     }
  81. }
  82. int belong[N];
  83. void build(int x,int fx,int from,int d){
  84.     belong[x]=fx;
  85.     dis[x]=0; getdeep(x,x,fx); P[x].add(0);
  86.     if(fx){
  87.         dis[from]=d;
  88.         getdis(x,from,fx);
  89.         P[fx].add(Q[x].top());
  90.     }
  91.     vis[x]=1;
  92.     for(int i=g[x];~i;i=nxt[i])if(!vis[v[i]]){
  93.         root=0; sum=size[v[i]];
  94.         getroot(v[i],x);
  95.         build(root,x,v[i],w[i]);
  96.     }
  97.     ans.add(P[x].toptwo());
  98. }
  99. void change(int rt){
  100.     int a=P[rt].toptwo();
  101.     if(col[rt])P[rt].add(0);
  102.     else P[rt].del(0);
  103.     int b=P[rt].toptwo();
  104.     if(a!=b)ans.del(a),ans.add(b);
  105.     int x=rt;
  106.     for(int i=G[rt];~NXT[i];i=NXT[i]){
  107.         int y=V[NXT[i]],len=W[NXT[i]];
  108.         x=V[i];
  109.         a=Q[x].top();
  110.         if(col[rt])Q[x].add(len);
  111.         else Q[x].del(len);
  112.         b=Q[x].top();
  113.         if(a!=b){
  114.             int c=P[y].toptwo();
  115.             P[y].del(a),P[y].add(b);
  116.             int d=P[y].toptwo();
  117.             if(c!=d)ans.del(c),ans.add(d);
  118.         }
  119.     }
  120. }
  121. int n,m,a,b,c,x;
  122. char op[10];
  123. int main(){
  124.     scanf("%d",&n); init();
  125.     for(int i=1;i<=n;i++)col[i]=1;
  126.     for(int i=1;i<n;i++){
  127.         scanf("%d%d%d",&a,&b,&c);
  128.         add_edge(a,b,c);
  129.         add_edge(b,a,c);
  130.     }
  131.     memset(vis,0,sizeof(vis));
  132.     dp[root=0]=sum=n;
  133.     getroot(1,0);
  134.     build(root,0,0,0);
  135.     scanf("%d",&m);
  136.     int cnt=n;
  137.     while(m--){
  138.         scanf("%s",op);
  139.         if(op[0]=='A'){
  140.             if(cnt)printf("%d\n",ans.top());
  141.             else puts("They have disappeared.");
  142.         }else{
  143.             scanf("%d",&x);
  144.             col[x]^=1;
  145.             if(col[x])cnt++;
  146.             else cnt--;
  147.             change(x);
  148.         }
  149.     }return 0;
  150. }

QTREE5

  1. /*
  2.     题目大意：给出一棵树，边权均为1，一开始所有点都是白点，
  3.     要求维护两种操作：
  4.         1.颜色取反（黑白互变）
  5.         2.查询树上某点距离最近的白点与其距离
  6.     题解：我们维护重心树，对于每个重心维护分治层所有白点到其距离，
  7.     将其保存在对应方向子重心的优先队列中，
  8.     重心处另外维护一个队列表示每个子重心的最小答案，
  9.     考虑修改问题，等价于在处理优先队列的删除问题，
  10.     对于每个需要删除操作的优先队列，我们额外创建一个优先队列将删除元素加入其中，
  11.     当两个队列top元素相同时同时删去即可。
  12.     对于查询操作，我们沿着重心链用经过每个重心的最优值更新答案
  13. */
  14. #include <cstdio>
  15. #include <algorithm>
  16. #include <cstring>
  17. #include <queue>
  18. using namespace std;
  19. const int N=200010;
  20. const int INF=0x3f3f3f3f;
  21. int tot,g[N],nxt[N<<1],v[N<<1],w[N<<1];
  22. int TOT,G[N],NXT[N*20],V[N*20],W[N*20];
  23. void init(){
  24.     memset(g,-1,sizeof(g));tot=0;
  25.     memset(G,-1,sizeof(G));TOT=0;
  26. }
  27. void add_edge(int x,int y,int z){v[tot]=y,w[tot]=z,nxt[tot]=g[x],g[x]=tot++;}
  28. void ADD_EDGE(int x,int y,int z){V[TOT]=y,W[TOT]=z,NXT[TOT]=G[x],G[x]=TOT++;}
  29. int sum,root,size[N],dp[N],vis[N];
  30. void getroot(int x,int fx){
  31.     size[x]=1; dp[x]=0;
  32.     for(int i=g[x];~i;i=nxt[i]){
  33.         if(!vis[v[i]]&&v[i]!=fx){
  34.             getroot(v[i],x);
  35.             size[x]+=size[v[i]];
  36.             dp[x]=max(dp[x],size[v[i]]);
  37.         }
  38.     }dp[x]=max(dp[x],sum-size[x]);
  39.     if(dp[x]<dp[root])root=x;
  40. }
  41. struct Que{
  42.     priority_queue<int,vector<int>,greater<int> > Q,D;
  43.     void clear(){
  44.         while(!Q.empty())Q.pop();
  45.         while(!D.empty())D.pop();
  46.     }
  47.     void add(int x){Q.push(x);}
  48.     void del(int x){D.push(x);}
  49.     int top(){
  50.         for(;;){
  51.             if(Q.empty())return INF;
  52.             else if(!D.empty()&&Q.top()==D.top())Q.pop(),D.pop();
  53.             else return Q.top();
  54.         }
  55.     }
  56. }P[N],Q[N];
  57. int dis[N],col[N];
  58. void getdis(int rt,int x,int fx){
  59.     for(int i=g[x];~i;i=nxt[i]){
  60.         if(v[i]==fx||vis[v[i]])continue;
  61.         dis[v[i]]=dis[x]+w[i];
  62.         getdis(rt,v[i],x);
  63.     }
  64. }
  65. void getdeep(int rt,int x,int fx){
  66.     ADD_EDGE(x,rt,dis[x]);
  67.     for(int i=g[x];~i;i=nxt[i]){
  68.         if(v[i]==fx||vis[v[i]])continue;
  69.         dis[v[i]]=dis[x]+w[i];
  70.         getdeep(rt,v[i],x);
  71.     }
  72. }
  73. int belong[N];
  74. void build(int x,int fx,int from,int d){
  75.     belong[x]=fx;
  76.     dis[x]=0; getdeep(x,x,fx);
  77.     if(fx){
  78.         dis[from]=d;
  79.         getdis(x,from,fx);
  80.     }
  81.     vis[x]=1;
  82.     for(int i=g[x];~i;i=nxt[i])if(!vis[v[i]]){
  83.         root=0; sum=size[v[i]];
  84.         getroot(v[i],x);
  85.         build(root,x,v[i],w[i]);
  86.     }
  87. }
  88. void change(int rt){
  89.     if(col[rt])P[rt].add(0);
  90.     else P[rt].del(0);
  91.     int x=rt;
  92.     for(int i=G[rt];~NXT[i];i=NXT[i]){
  93.         int y=V[NXT[i]],len=W[NXT[i]];
  94.         x=V[i];
  95.         int a=Q[x].top();
  96.         if(col[rt])Q[x].add(len);
  97.         else Q[x].del(len);
  98.         int b=Q[x].top();
  99.         if(a!=b)P[y].del(a),P[y].add(b);
  100.     }
  101. }
  102. int ask(int rt){
  103.     int ans=INF;
  104.     for(int i=G[rt];~i;i=NXT[i]){
  105.         int x=V[i],len=W[i];
  106.         int a=P[x].top();
  107.         if(a+len<ans)ans=a+len;
  108.     }
  109.     return ans;
  110. }
  111. int n,m,a,b,c,x,op;
  112. int main(){
  113.     scanf("%d",&n); init();
  114.     for(int i=1;i<=n;i++)col[i]=0;
  115.     for(int i=1;i<n;i++){
  116.         scanf("%d%d",&a,&b);
  117.         add_edge(a,b,1);
  118.         add_edge(b,a,1);
  119.     }
  120.     memset(vis,0,sizeof(vis));
  121.     dp[root=0]=sum=n;
  122.     getroot(1,0);
  123.     build(root,0,0,0);
  124.     scanf("%d",&m);
  125.     int cnt=0;
  126.     while(m--){
  127.         scanf("%d%d",&op,&x);
  128.         if(op==1){
  129.             if(cnt)printf("%d\n",ask(x));
  130.             else puts("-1");
  131.         }else{
  132.             col[x]^=1;
  133.             if(col[x])cnt++;
  134.             else cnt--;
  135.             change(x);
  136.         }
  137.     }return 0;
  138. }

QTREE6

  1. /*
  2.     题目大意：给出一棵树，一开始树上点均为黑色，要求维护操作
  3.         1.翻转某点的颜色(黑白互换)，
  4.         2.查询某个点所在连通块(颜色相同则连通)大小
  5.     题解：我们分别维护黑点LCT和白点LCT，当一个点从黑变白时，
  6.     将其从黑点LCT中与父节点断开，然后在白点LCT中与父节点相连，
  7.     这样我们就保证了每个连通块在LCT中只有父节点是不同色的而其余一定是连通的。
  8.     考虑用LCT维护子树信息，根据在LCT上的连接情况，我们将点的儿子分为实儿子和虚儿子
  9.     实儿子是原本树结构上相连的点，实儿子的子树为实子树，虚儿子的子树为虚子树
  10.     x的LCT子树的信息和等于x的实儿子的LCT子树信息和加上x的虚子树的信息和加上x自己
  11.     在进行access操作时，我们会有更换一个点的x右儿子的操作，
  12.     这时我们要把x原来的右儿子的LCT子树信息加入x的虚子树信息，
  13.     把x的新的右儿子的LCT子树信息从x的虚子树信息中减去
  14.     同时在link x到y的时候，我们需要对y进行access再splay
  15.     这样就只会对y的虚子树信息和LCT子树信息产生影响，而不会影响到y的祖先节点
  16. */
  17. #include <cstdio>
  18. #include <algorithm>
  19. #include <set>
  20. const int N=100010;
  21. using namespace std;
  22. int n,m,i,k,x,y,c[N],fa[N],g[N],v[N<<1],nxt[N<<1],ed;
  23. struct LCT{
  24.     int son[N][2],f[N],sum[N],s[N];
  25.     bool isroot(int x){return !f[x]||son[f[x]][0]!=x&&son[f[x]][1]!=x;}
  26.     void up(int x){sum[x]=1+sum[son[x][0]]+sum[son[x][1]]+s[x];}
  27.     void rotate(int x){
  28.         int y=f[x],w=son[y][1]==x;
  29.         son[y][w]=son[x][w^1];
  30.         if(son[x][w^1])f[son[x][w^1]]=y;
  31.         if(f[y]){
  32.             int z=f[y];
  33.             if(son[z][0]==y)son[z][0]=x;
  34.             if(son[z][1]==y)son[z][1]=x;
  35.         }f[x]=f[y];f[y]=x;son[x][w^1]=y;up(y);
  36.     }
  37.     void splay(int x){
  38.         while(!isroot(x)){
  39.             int y=f[x];
  40.             if(!isroot(y)){
  41.                 if((son[f[y]][0]==y)^(son[y][0]==x))rotate(x);
  42.                 else rotate(y);
  43.             }rotate(x);
  44.         }up(x);
  45.     }
  46.     void access(int x){
  47.         for(int y=0;x;y=x,x=f[x]){
  48.             splay(x);
  49.             if(son[x][1])s[x]+=sum[son[x][1]];
  50.             if(son[x][1]=y)s[x]-=sum[y];
  51.             up(x);
  52.         }
  53.     }
  54.     int root(int x){access(x);splay(x);while(son[x][0])x=son[x][0];return x;}
  55.     void link(int x){access(fa[x]);splay(fa[x]);splay(x);son[fa[x]][1]=x;up(f[x]=fa[x]);}
  56.     void cut(int x){access(x);splay(x);f[son[x][0]]=0;son[x][0]=0;up(x);}
  57.     int ask(int x){splay(x=root(x));return sum[son[x][1]];}
  58. }T[2];
  59. void add(int x,int y){v[++ed]=y;nxt[ed]=g[x];g[x]=ed;}
  60. void dfs(int x){
  61.     for(int i=g[x];i;i=nxt[i])if(v[i]!=fa[x]){
  62.         fa[v[i]]=T[c[v[i]]].f[v[i]]=x;
  63.         dfs(v[i]);
  64.         T[c[v[i]]].s[x]+=T[c[v[i]]].sum[v[i]];
  65.     }T[0].up(x),T[1].up(x);
  66. }
  67. void read(int&a){
  68.     char c;
  69.     while(!(((c=getchar())>='0')&&(c<='9')));a=c-'0';
  70.     while(((c=getchar())>='0')&&(c<='9'))(a*=10)+=c-'0';
  71. }
  72. int main(){
  73.      read(n);
  74.      for(i=1;i<n;i++)read(x),read(y),add(x,y),add(y,x);
  75.      for(i=1;i<=n;i++)c[i]=0;
  76.      fa[1]=T[c[1]].f[1]=n+1;dfs(1);
  77.      T[c[1]].s[n+1]+=T[c[1]].sum[1];read(m);
  78.      while(m--){
  79.          read(k);read(x);
  80.          if(!k)printf("%d\n",T[c[x]].ask(x));
  81.          else T[c[x]].cut(x),T[c[x]^=1].link(x);
  82.      }return 0;
  83. }

QTREE7

  1. /*
  2.     题目大意：给出一棵树，树上每个点为黑色或者白色，每个点有个点权，
  3.     我们认为树上颜色相同并通过边连接的块为连通块
  4.     要求维护操作：
  5.         0.查询x所在连通块中权值最大的点
  6.         1.翻转某点的颜色(黑白互换)
  7.         2.改变某点的权值
  8.     题解：我们分别维护黑点LCT和白点LCT，当一个点从黑变白时，
  9.     将其从黑点LCT中与父节点断开，然后在白点LCT中与父节点相连，
  10.     这样我们就保证了每个连通块在LCT中只有父节点是不同色的而其余一定是连通的。
  11.     考虑用LCT维护子树信息，根据在LCT上的连接情况，我们将点的儿子分为实儿子和虚儿子
  12.     实儿子是原本树结构上相连的点，实儿子的子树为实子树，虚儿子的子树为虚子树
  13.     x的LCT子树的信息和等于x的实儿子的LCT子树信息和加上x的虚子树的信息和加上x自己
  14.     在进行access操作时，我们会有更换一个点的x右儿子的操作，
  15.     这时我们要把x原来的右儿子的LCT子树信息加入x的虚子树信息，
  16.     把x的新的右儿子的LCT子树信息从x的虚子树信息中减去
  17.     同时在link x到y的时候，我们需要对y进行access再splay
  18.     这样就只会对y的虚子树信息和LCT子树信息产生影响，而不会影响到y的祖先节点
  19.     因为要维护的信息是极值，因此我们要用multiset来维护虚子树信息
  20. */
  21. #include <cstdio>
  22. #include <algorithm>
  23. #include <ctype.h>
  24. #include <set>
  25. const int N=100010;
  26. using namespace std;
  27. int n,m,k,x,y,c[N],fa[N],g[N],v[N<<1],nxt[N<<1],ed;
  28. struct LCT{
  29.     int son[N][2],f[N],val[N],mx[N];
  30.     multiset<int> s[N];
  31.     bool isroot(int x){return !f[x]||son[f[x]][0]!=x&&son[f[x]][1]!=x;}
  32.     void up(int x){
  33.         mx[x]=val[x];
  34.         if(s[x].size())mx[x]=max(mx[x],*s[x].rbegin());
  35.         if(son[x][0])mx[x]=max(mx[x],mx[son[x][0]]);
  36.         if(son[x][1])mx[x]=max(mx[x],mx[son[x][1]]);
  37.     }
  38.     void rotate(int x){
  39.         int y=f[x],w=son[y][1]==x;
  40.         son[y][w]=son[x][w^1];
  41.         if(son[x][w^1])f[son[x][w^1]]=y;
  42.         if(f[y]){
  43.             int z=f[y];
  44.             if(son[z][0]==y)son[z][0]=x;
  45.             if(son[z][1]==y)son[z][1]=x;
  46.         }f[x]=f[y];f[y]=x;son[x][w^1]=y;up(y);
  47.     }
  48.     void splay(int x){
  49.         while(!isroot(x)){
  50.             int y=f[x];
  51.             if(!isroot(y)){
  52.                 if((son[f[y]][0]==y)^(son[y][0]==x))rotate(x);
  53.                 else rotate(y);
  54.             }rotate(x);
  55.         }up(x);
  56.     }
  57.     void access(int x){
  58.         for(int y=0;x;y=x,x=f[x]){
  59.             splay(x);
  60.             if(son[x][1])s[x].insert(mx[son[x][1]]);
  61.             if(son[x][1]=y)s[x].erase(s[x].find(mx[y]));
  62.             up(x);
  63.         }
  64.     }
  65.     int root(int x){access(x);splay(x);while(son[x][0])x=son[x][0];return x;}
  66.     void link(int x){access(fa[x]);splay(fa[x]);splay(x);son[fa[x]][1]=x;up(f[x]=fa[x]);}
  67.     void cut(int x){access(x);splay(x);f[son[x][0]]=0;son[x][0]=0;up(x);}
  68.     int ask(int x){
  69. 		int t=c[x]; splay(x=root(x));
  70. 		return t==c[x]?mx[x]:mx[son[x][1]];
  71. 	}
  72.     void modify(int x,int v){access(x),splay(x),val[x]=v,up(x);}
  73. }T[2];
  74. void add(int x,int y){v[++ed]=y;nxt[ed]=g[x];g[x]=ed;}
  75. void dfs(int x){
  76.     for(int i=g[x];i;i=nxt[i])if(v[i]!=fa[x]){
  77.         fa[v[i]]=T[c[v[i]]].f[v[i]]=x;
  78.         dfs(v[i]);
  79.         T[c[v[i]]].s[x].insert(T[c[v[i]]].mx[v[i]]);
  80.     }T[0].up(x),T[1].up(x);
  81. }
  82. void read(int &x){
  83.     char ch;
  84.     for(ch=getchar();!isdigit(ch)&&ch!='-';ch=getchar());
  85.     int t;if(ch=='-')t=-1,ch=getchar();else t=1;
  86.     for(x=0;isdigit(ch);ch=getchar())x=x*10+ch-'0';
  87.     x=x*t;
  88. }
  89. int main(){
  90.      read(n);
  91.      for(int i=1;i<n;i++)read(x),read(y),add(x,y),add(y,x);
  92.      for(int i=1;i<=n;i++)read(c[i]);
  93.      for(int i=1;i<=n;i++)read(T[0].val[i]),T[1].val[i]=T[0].val[i];
  94.      dfs(1); read(m);
  95.      while(m--){
  96.          read(k);read(x);
  97.          if(k==0)printf("%d\n",T[c[x]].ask(x));
  98.          else if(k==1)T[c[x]].cut(x),T[c[x]^=1].link(x);
  99.          else if(k==2)read(y),T[0].modify(x,y),T[1].modify(x,y);
  100.      }return 0;
  101. }

SPOJ QTREE的更多相关文章

  1. QTREE 树链剖分---模板 spoj QTREE

    <树链剖分及其应用> 一文讲得非常清楚,我一早上就把他学会了并且A了这题的入门题. spoj QTREE 题目: 给出一棵树,有两种操作: 1.修改一条边的边权. 2.询问节点a到b的最大 ...

  2. SPOJ QTREE 系列解题报告

    题目一 : SPOJ 375 Query On a Tree http://www.spoj.com/problems/QTREE/ 给一个树,求a,b路径上最大边权,或者修改a,b边权为t. #in ...

  3. SPOJ QTREE Query on a tree 树链剖分+线段树

    题目链接:http://www.spoj.com/problems/QTREE/en/ QTREE - Query on a tree #tree You are given a tree (an a ...

  4. spoj QTREE - Query on a tree(树链剖分+线段树单点更新，区间查询)

    传送门:Problem QTREE https://www.cnblogs.com/violet-acmer/p/9711441.html 题解: 树链剖分的模板题,看代码比看文字解析理解来的快~~~ ...

  5. SPOJ QTREE Query on a tree ——树链剖分 线段树

    [题目分析] 垃圾vjudge又挂了. 树链剖分裸题. 垃圾spoj,交了好几次,基本没改动却过了. [代码](自带常数,是别人的2倍左右) #include <cstdio> #incl ...

  6. 【学术篇】SPOJ QTREE 树链剖分

    发现链剖这东西好久不写想一遍写对是有难度的.. 果然是熟能生巧吧.. WC的dalao们都回来了 然后就用WC的毒瘤题荼毒了我们一波, 本来想打个T1 44分暴力 然后好像是特判写挂了还是怎么的就只能 ...

  7. SPOJ QTree【树链剖分】

    一 题目 QTREE 二 分析 第一道树链剖分的题,写的好艰难啊. 题意还是比较好理解的,就是在树上操作. 对于修改,题中要求的是单点修改,就算是直接树上操作也是非常简单的. 对于查询,查询的时候,是 ...

  8. SPOJ QTREE Query on a tree

    题意:给一颗n个点的树,有两种操作CHANGE i ti : 把第i条边的权变为tiQUERY a b : 问点a 到 点b 之间的边的最大权 思路:树剖处理边权.由于是边,所以只需要把边权处理到子节 ...

  9. SPOJ QTREE Query on a tree --树链剖分

    题意:给一棵树,每次更新某条边或者查询u->v路径上的边权最大值. 解法:做过上一题,这题就没太大问题了,以终点的标号作为边的标号,因为dfs只能给点分配位置,而一棵树每条树边的终点只有一个. ...

随机推荐

  1. 2016.6.17——Remove Duplicates from Sorted Array

    Remove Duplicates from Sorted Array 本题收获: 1.“删除”数组中元素 2.数组输出 题目: Given a sorted array, remove the du ...

  2. Dream------scala--类的属性和对象私有字段实战详解

    Scala类的属性和对象私有字段实战详解 一.类的属性 scala类的属性跟java有比较大的不同,需要注意的是对象的私有(private)字段 1.私有字段:字段必须初始化(当然即使不是私有字段也要 ...

  3. Flask--wtforms快速使用和表单验证(附示例)

    一.Form类 表单提供WTForms中最高级别的API.它们包含您的字段定义,委托验证,获取输入,聚合错误,并且通常用作将所有内容组合在一起的粘合剂. class wtforms.form.Form ...

  4. Django 内置模板标签和过滤器

    一.内置模板标签 语法:{%  %} autoescape : 是否转义,on或off作为参数,并确定自动转义是否在块内有效.该块以endautoescape结束 {% autoescape on % ...

  5. python网络编程--线程锁(互斥锁Mutex)

    一:为什么需要线程锁 一个进程下可以启动多个线程,多个线程共享父进程的内存空间,也就意味着每个线程可以访问同一份数据,此时,如果2个线程同时要修改同一份数据,会出现什么状况? 很简单,假设你有A,B两 ...

  6. Python 正则表达式提高

    re模块的高级用法 search re.search(pattern, string[, flags]) ​ 若string中包含pattern子串,则返回Match对象,否则返回None,注意,如果 ...

  7. 如何阻止点击scrollviewer里面的单位内容时，自动滚动

    <Style TargetType="{x:Type ListBoxItem}"> <Setter Property="FocusVisualStyle ...

  8. JMeter接口&性能测试

    JMeter接口测试 目前最新版本发展到5.0版本,需要Java7以上版本环境,下载解压目录后,进入\apache-jmeter-5.0\bin\,双击ApacheJMeter.jar文件启动JMem ...

  9. Unix IPC之Posix消息队列(2)

    /* Query status and attributes of message queue MQDES. */ extern int mq_getattr (mqd_t __mqdes, stru ...

  10. Linux命令之cp命令

    cp命令:用来将一个或多个源文件或者目录复制到指定的目的文件或目录.它可以将单个源文件复制成一个指定文件名的具体的文件或一个已经存在的目录下.cp命令还支持同时复制多个文件,当一次复制多个文件时,目标 ...