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][20];
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]+1;
        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=19;i>=0;i--)
    {
        if(deep[y]-deep[x]>=(1<<i))
        {
            y=s[y][i];
        }
    }
    if(x==y)return x;
    for(i=19;i>=0;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=0;
    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=0;
    if (l>=r) return;
    m=(l+r)>>1;
    for (i=l;i<=m;i++)
    if (b[i].typ==2)
    {
        w[++tot].id=i;
        w[tot].v=m-i;
    }
    for (i=m+1;i<=r;i++)
    if (b[i].typ==1)
    {
        w[++tot].id=i;
        w[tot].v=b[i].w-(i-m);
    }
    sort(w+1,w+1+tot,cmp);
    
    for (i=1;i<=tot;i++)
    if (w[i].id<=m)
    {
        cc(L[b[w[i].id].w],1);
        cc(R[b[w[i].id].w],-1);
    }
    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]-1);
    }
    for (i=1;i<=tot;i++)
    if (w[i].id<=m)
    {
        cc(L[b[w[i].id].w],-1);
        cc(R[b[w[i].id].w],1);
    }

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

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