[BZOJ4573][ZJOI2016]大森林（LCT建虚点）

题目：

我是超链接

题解：

发现对于询问，每棵树只要把询问的两个点长在该长的位置了，它多长了些什么点，包括不该长的点都无所谓。那么让每棵树先长完了再询问就好。

也就是说问的是有n棵树，每次让所有树长出一个节点，一开始都从1开始长，然后会修改一个区间的树的长的位置。

先让所有节点长出来，给每个更换生长点的操作建一个权值为0的虚点，所有节点就先往它之前的最后一个虚点上长。然后把虚点按前后一个个串起来。

离线所有操作，然后从第左到右处理每棵树。

每个虚点生效的是一段区间l~r，在l的时候虚点开始生效，就把它切下来接到它该接的实点上，在r+1的时候虚点失效，就把它切下来接回它的上一个虚点上，说明这个虚点上面长的点真正该长的地方是它的上一个虚点接的地方，或者上一个点也失效了接到前面去，就是前面某个地方，总会到达一个它该长的地方。

然后求两点距离时不能直接求，需要通过lca，因为存在虚点，若两个点的lca是一个虚点，他们真正的lca其实是虚点的某个祖先，直接算路径就会算少，而求lca可以避免这个问题。

方法是access（x），再access（y），最后一次虚边边实边的x就是lca。

代码：

#include <cstdio>
#include <cstring>
#include <algorithm>
using namespace std;
const int N=300005;
struct hh
{
    int pos,op,x,y; 
    hh (int Pos=0,int Op=0,int X=0,int Y=0){pos=Pos;op=Op;x=X;y=Y;}
}ask[N];
int ch[N][2],f[N],size[N],val[N],delta[N],stack[N],cnt,n,num,m,L[N],id[N],R[N],now,tot,ans[N];
int cmp(hh a,hh b){return a.pos<b.pos || (a.pos==b.pos && a.op<b.op);}
bool isroot(int x){return ch[f[x]][1]!=x && ch[f[x]][0]!=x;}
int get(int x){return ch[f[x]][1]==x;}
void updata(int now)
{
    size[now]=val[now];
    if (ch[now][0]) size[now]+=size[ch[now][0]];
    if (ch[now][1]) size[now]+=size[ch[now][1]];
}
void pushdown(int now)
{
    if (delta[now])
    {
        delta[ch[now][0]]^=1; delta[ch[now][1]]^=1;
        swap(ch[now][0],ch[now][1]);
        delta[now]=0;
    }
}
void rotate(int x)
{
    int old=f[x],oldf=f[old],which=get(x);
    f[x]=oldf; if (!isroot(old)) ch[oldf][get(old)]=x;
    ch[old][which]=ch[x][which^1]; f[ch[x][which^1]]=old;
    ch[x][which^1]=old; f[old]=x;
    updata(old); updata(x);
}
void splay(int x)
{
    int top=0,i;
    for (i=x;!isroot(i);i=f[i]) stack[++top]=i;
    stack[++top]=i;
    for (int i=top;i>=1;i--) pushdown(stack[i]);
    for (;!isroot(x);rotate(x))
      if (!isroot(f[x])) rotate(get(x)==get(f[x])?f[x]:x);
}
int access(int x)
{
    int t=0;
    for (;x;t=x,x=f[x])
    {
        splay(x);
        ch[x][1]=t;
        updata(x);
    }
    return t;
}
void reverse(int x){access(x);splay(x); delta[x]^=1;}
void link(int x,int y){reverse(x);f[x]=y;}
void cut(int x)
{
    access(x); splay(x);
    ch[x][0]=f[ch[x][0]]=0; updata(x);
}
void insert(int v){++cnt;size[cnt]=val[cnt]=v;}
int main()
{
    scanf("%d%d",&n,&m);
    insert(1);num=1;L[num]=id[num]=1; R[num]=n;
    insert(0);now=2;//最近的虚点
    link(2,1); 
    for (int i=1;i<=m;i++)
    {
        int op,k,x,y;
        scanf("%d",&op);
        if (op==0)
        {
            scanf("%d%d",&x,&y); ++num; insert(1);
            L[num]=x; R[num]=y; //多少个实点 
            id[num]=cnt;ask[++tot]=hh(1,i-m,cnt,now);
            //对所有树都长叶子，排在所有2的前面 ,当前点，接到哪个虚点 
        }else if (op==1)
        {
            scanf("%d%d%d",&x,&y,&k);
            x=max(x,L[k]); y=min(y,R[k]);
            if (x<=y)
            {
                insert(0); if (x>1) link(cnt,now);
                ask[++tot]=hh(x,i-m,cnt,id[k]);
                ask[++tot]=hh(y+1,i-m,cnt,now);
                now=cnt;
            }
        }else
        {
            scanf("%d%d%d",&k,&x,&y);
            ask[++tot]=hh(k,i,id[x],id[y]);
        }
    }
    sort(ask+1,ask+tot+1,cmp);
    memset(ans,-1,sizeof(ans));
    int k=1;
    for (int i=1;i<=n;i++)
      for (;k<=tot && ask[k].pos==i;k++)
      {
        if (ask[k].op>0)//2操作 
        {
            access(ask[k].x); splay(ask[k].x); ans[ask[k].op]=size[ask[k].x];
            int t=access(ask[k].y); 
              splay(ask[k].y); 
              ans[ask[k].op]+=size[ask[k].y];
            access(t); splay(t); ans[ask[k].op]-=2*size[t];
        }else cut(ask[k].x),link(ask[k].x,ask[k].y);
      }
    for (int i=1;i<=m;i++) if (ans[i]>=0) printf("%d\n",ans[i]);
}