[BZOJ4573][UOJ#195][Zjoi2016][LCT][离线]大森林

http://blog.csdn.net/lych_cys/article/details/53515748#

神想法啊…

至于求两点间距离为什么不能直接把两点的链拉出来

可以考虑因为树上虚点，两个点的LCA可以是一个虚点，虚点的父节点才是它们真正的LCA，所以直接拉一条链会导致少算一个点，而算LCA的话不会有这个问题（可以脑补一下）

#include <cstdio>
#include <iostream>
#include <algorithm>
#include <string>
#include <cstring>
#include <vector>
#define N 100010

using namespace std;

int n,m,cnt,top,tot,g,cnt0;
int Ans
,stc
,L
,R
;
struct lef{
int Size,rev,key;
lef *ch[2],*fa;
int isr(){ return (!fa)||(fa->ch[0]!=this&&fa->ch[1]!=this); }
int isl(){ return fa&&fa->ch[1]==this; }
void Up(){ Size=key+(ch[0]?ch[0]->Size:0)+(ch[1]?ch[1]->Size:0); }
void reverse(){ access(); splay(); rev^=1; }
void Push(){
if(!rev) return ;
swap(ch[0],ch[1]);
if(ch[0]) ch[0]->rev^=1;
if(ch[1]) ch[1]->rev^=1;
rev=0;
}
lef *access(){
lef *t,*x,*p;
for(t=0,x=this;x;x=x->fa)
p=x,x->splay(),x->ch[1]=t,t=x,x->Up();
return p;
}
void splay(){
Pushtop();
for(lef *x=this;!x->isr();x->rot())
if(!x->fa->isr()) ((x->isl())^(x->fa->isl()))?x->rot():x->fa->rot();
}
void rot(){
lef *x=this,*y=x->fa,*z=y->fa; int wh=x->isl();
if(!y->isr()) z->ch[z->ch[1]==y]=x; x->fa=z;
if(y->ch[wh]=x->ch[wh^1]) y->ch[wh]->fa=y;
x->ch[wh^1]=y; y->fa=x;
y->Up(); x->Up();
}
void Pushtop(){ if(!isr()) this->fa->Pushtop(); Push(); }
}A[N<<2];
struct stp{int op,x,y,z,g;};
vector<stp> work
;

inline void link(int x,int y){ A[x].reverse(); A[x].fa=&A[y]; A[x].access(); }
inline void cut(int x,int y){ A[x].reverse(); A[y].access(); A[y].splay(); A[y].ch[0]=A[x].fa=0; A[y].Up(); }
inline int query(int x,int y){
int ret=0; lef *lca; A[1].reverse();
A[x].access(); A[x].splay(); ret=A[x].Size;
lca=A[y].access(); A[y].splay(); ret+=A[y].Size;
lca->access(); lca->splay(); ret-=2*(lca->Size);
return ret;
}

inline char nc(){
static char buf[100000],*p1=buf,*p2=buf;
return p1==p2&&(p2=(p1=buf)+fread(buf,1,100000,stdin),p2==p1)?EOF:*p1++;
}

inline void rea(int &x){
char c=nc(); x=0;
for(;c>'9'||c<'0';c=nc());for(;c>='0'&&c<='9';x=x*10+c-'0',c=nc());
}

inline lef *newd(int x,int y){
A[x].Size=y; A[x].ch[0]=A[x].ch[1]=A[x].fa=0; A[x].key=y;
return &A[x];
}

int main(){
rea(n); rea(m);
newd(cnt=1,1); L[cnt]=1; R[cnt]=n;
newd(tot=m+1,0); link(tot,cnt);
for(int i=1;i<=m;i++){
int op; rea(op);
if(op==0){
++cnt; rea(L[cnt]); rea(R[cnt]);
newd(cnt,1); link(cnt,tot);
}
else if(op==1){
int l,r,x;
rea(l); rea(r); rea(x);
l=max(l,L[x]); r=min(r,R[x]);
if(l<=r){
newd(++tot,0); link(tot,tot-1);
work[l].push_back((stp){1,tot,tot-1,x});
work[r+1].push_back((stp){1,tot,x,tot-1});
}
}
else{
int x,y,z; rea(x); rea(y); rea(z);
work[x].push_back((stp){2,x,y,z,++cnt0});
}
}
for(int i=1;i<=n;i++){
stp now;
for(int j=0;j<work[i].size();j++){
now=work[i][j];
if(now.op==1){
cut(now.x,now.y); link(now.x,now.z);
}
else
Ans[now.g]=query(now.y,now.z);
}
}
for(int i=1;i<=cnt0;i++) printf("%d\n",Ans[i]);
return 0;
}