bzoj4448 [Scoi2015]情报传递

　　第一问不解释，对于第二问的处理，可以使用cdq分治，假设分治的询问区间是[L,R],那么我们对于标号在[L,mid]的修改操作赋予一个权值，因为在当前[L,R]中[L,mid]的修改操作只会对[mid+1，R]的询问操作，所以第i修改操作至少经过m-i的时间，因此赋予的权值是m-i,而对于[mid+1,R]区间中的询问操作，也赋予一个权值w-m,这里w为询问的数值，我们可以预处理出树的dfs序并维护一个树状数组，这样我们就可以把这些操作按权值从大到小插入或者询问，时间复杂度O(nlognlogn)。

　　代码，运行时间差不多垫底。。(看了其他博主的题解貌似我想的复杂了。。。)

 #include<cstdio>
 #include<algorithm>
 #define N 500010
 #define lb(x) (x&-x)
 using namespace std;
 int dp,p[N],pre[N],tt[N],n,a,i,m,o,c[N],deep[N],fa[N];
 int L[N],R[N],ans[N];
 int s[N][];
 struct g{
     int typ,l,r,w;
 }b[N];
 struct gg{
     int id,v;
 }w[N];
 void cc(int x,int w)
 {
     while (x<=o)
     {
         c[x]+=w;
         x+=lb(x);
     }
 }
 void dfs(int x)
 {
     int i=p[x];
     L[x]=++o;
     while (i)
     {
         fa[tt[i]]=x;
         deep[tt[i]]=deep[x]+;
         dfs(tt[i]);
         i=pre[i];
     }
     R[x]=++o;
 }
 int lca(int x,int y)
 {
     if(deep[x]>deep[y])x^=y^=x^=y;
     int i;
     for(i=;i>=;i--)
     {
         if(deep[y]-deep[x]>=(<<i))
         {
             y=s[y][i];
         }
     }
     if(x==y)return x;
     for(i=;i>=;i--)
     {
         if(s[x][i]!=s[y][i])
         {
             x=s[x][i];
             y=s[y][i];
         }
     }
     return fa[x];
 }
 int sum(int x)
 {
     int ans=;
     while (x)
     {
         ans+=c[x];
         x-=lb(x);
     }
     return ans;
 }
 void link(int x,int y)
 {
     dp++;pre[dp]=p[x];p[x]=dp;tt[dp]=y;
 }
 bool cmp(gg a,gg b)
 {
     if (a.v==b.v)
     return a.id>b.id;
     return a.v>b.v;
 }
 void solve(int l,int r)
 {
     int m,i,tot=;
     if (l>=r) return;
     m=(l+r)>>;
     for (i=l;i<=m;i++)
     if (b[i].typ==)
     {
         w[++tot].id=i;
         w[tot].v=m-i;
     }
     for (i=m+;i<=r;i++)
     if (b[i].typ==)
     {
         w[++tot].id=i;
         w[tot].v=b[i].w-(i-m);
     }
     sort(w+,w++tot,cmp);

     for (i=;i<=tot;i++)
     if (w[i].id<=m)
     {
         cc(L[b[w[i].id].w],);
         cc(R[b[w[i].id].w],-);
     }
     else
     {
         int u=b[w[i].id].l;
         int v=b[w[i].id].r;
         int LCA=lca(u,v);
         ans[w[i].id]+=sum(L[u])+sum(L[v])-sum(L[LCA])-sum(L[LCA]-);
     }
     for (i=;i<=tot;i++)
     if (w[i].id<=m)
     {
         cc(L[b[w[i].id].w],-);
         cc(R[b[w[i].id].w],);
     }

     solve(l,m);solve(m+,r);
 }

 int main()
 {
     scanf("%d",&n);
     for (i=;i<=n;i++)
     {
         scanf("%d",&a);
         if (a) link(a,i);
     }
     scanf("%d",&m);
     for (i=;i<=m;i++)
     {
         scanf("%d",&b[i].typ);
         if (b[i].typ==)
         scanf("%d%d%d",&b[i].l,&b[i].r,&b[i].w);
         else scanf("%d",&b[i].w);
     }
     dfs();
     for(i=;i<=n;i++)s[i][]=fa[i];
     for(int h=;h<;h++)
     {
         for(i=;i<=n;i++)
         {
             s[i][h]=s[s[i][h-]][h-];
         }
     }
     solve(,m);
     for (i=;i<=m;i++)
     if (b[i].typ==)
     {
         int Ans=lca(b[i].l,b[i].r);
         Ans=deep[b[i].l]+deep[b[i].r]-*deep[Ans]+;
         printf("%d %d\n",Ans,ans[i]);
     }
 }