hdu5709Claris Loves Painting主席树 奇妙的DFS序

先不考虑层数限制

一棵树上每个点有个颜色，问一棵子树的颜色数

感觉简单多了是吧

考虑每个点的贡献：自己到根的路径上的一个包含自己的连续段

观察最顶端的点的父亲：

它满足有了额外的同色孩子（咦）

这一条链上需要+1（@miaom）

如果差分，就是在这个点+1（@miaom），在最近的有同色孩子的父亲-1，求子树和

最近的有同色孩子的父亲？

根据直觉，它一定是dfs序上和这个点相邻的同色点和这个点的lca

没了

再加上层数限制

相当于一定要在一定的深度以内

所以开主席树保持可持久化（怎么不能离线啊，好烦啊）

 #include <bits/stdc++.h>
 #define LOG 19
 #define mid (l+r>>1)
 using namespace std;
 int NODE,n,m,p,q,E,TIME,x,d;
 int dep[],nex[],fir[],Fir[],Nex[];
 int fa[][],to[];//,ne[200001];
 int tr[],ls[],rs[];
 int c[],root[],start[],en[],pos[];
 void dfs(int now,int fat)
 {
     dep[now]=dep[fat]+;start[now]=++TIME;pos[TIME]=now;
     Nex[now]=Fir[dep[now]];Fir[dep[now]]=now;
     fa[now][]=fat;
     for(int i=;fa[now][i-];i++) fa[now][i]=fa[fa[now][i-]][i-];
     for(int i=fir[now];i;i=nex[i])
     if(to[i]!=fat)
         dfs(to[i],now);
     en[now]=TIME;
 }
 int lca(int x,int y)
 {
     if(dep[x]<dep[y]) swap(x,y);
     for(int delta=dep[x]-dep[y],i=;delta;delta>>=,++i)
         if(delta&) x=fa[x][i];
     for(int i=LOG;i>=;i--)
         if(fa[x][i]!=fa[y][i])
             x=fa[x][i],y=fa[y][i];
     return (x==y)?x:fa[x][];
 }
 int chan(int acc,int l,int r,int x,int y)
 {
     int now=++NODE;
     if(l==r)
     {
         tr[now]=tr[acc]+y;
         return now;
     }
     if(x<=mid) ls[now]=chan(ls[acc],l,mid,x,y),rs[now]=rs[acc];
     else rs[now]=chan(rs[acc],mid+,r,x,y),ls[now]=ls[acc];
     tr[now]=tr[ls[now]]+tr[rs[now]];
     return now;
 }
 void add(int x,int y)
 {
     to[++E]=y;nex[E]=fir[x];fir[x]=E;
 }
 int que(int now,int l,int r,int x,int y)
 {
     if(l==x && r==y)
         return tr[now];
     int ret=;
     if(x<=mid)
         ret+=que(ls[now],l,mid,x,min(y,mid));
     if(y>mid)
         ret+=que(rs[now],mid+,r,max(mid+,x),y);
     return ret;
 }
 set<pair<int,int> > se;
 int main()
 {
 int T;
 for(scanf("%d",&T);T;T--)
 {
     scanf("%d%d",&n,&m);E=;TIME=;NODE=;
     se.clear();
     for(int i=;i<=n;i++)
         scanf("%d",&c[i]);
     for(int i=;i<=n;i++)
         fir[i]=,Fir[i]=;
     for(int i=;i<=n;i++)
         for(int j=;j<=;j++)
             fa[i][j]=;
     for(int i=;i<n;i++)
         scanf("%d",&q),add(q,i+);
     dfs(,);root[]=;tr[]=;ls[]=;rs[]=;
     int MD=;
     for(int i=;i<=n;root[i+]=root[i],i++)
         for(int j=Fir[i];j;j=Nex[j])
         {
             int t=;
             pair <int,int> tem=*se.lower_bound(make_pair(c[j],start[j]));
             int tem1=tem.first;
             int Tem1=tem.second;
             if(tem1==c[j]) t=lca(pos[Tem1],j);
             if(se.lower_bound(make_pair(c[j],start[j]))!=se.begin())
             {
                 tem=*(--se.lower_bound(make_pair(c[j],start[j])));
                 int tem2=tem.first;
                 int Tem2=tem.second;
                 if(tem2==c[j] &&(!t || dep[t]<dep[lca(pos[Tem2],j)])) t=lca(pos[Tem2],j);
             }
             se.insert(make_pair(c[j],start[j]));
             if(t)
                 root[i]=chan(root[i],,n,start[t],-);
             root[i]=chan(root[i],,n,start[j],);
             MD=max(MD,i);
         }
     int lastans=;
     for(int i=;i<=m;i++)
     {
         if(i==)
             int e=;
         scanf("%d%d",&x,&d);
         x^=lastans;d^=lastans;
         int tim=min(dep[x]+d,MD);
         printf("%d\n",lastans=que(root[tim],,n,start[x],en[x]));
     }
 }
     return ;
  }