好久不更了。但是学新算法,就算做作业时间再紧张,也要记一下。
初学lct。几个要点提一下
1.同一个联通块是同一棵。不是一个splay或者一条链;
2.根节点有很多。刚开始打还想着用变量rt我是傻了吗;
3.那个access还算好理解。关键是代码里面那个rever是把x变为x所在树(不是splay)的根节点,理解了好久;
差不多就这么多了吧。
代码是抄hzwer黄学长的。orz~:
#include<iostream>
#include<cstdio>
#include<cstring>
#define N 30025
using namespace std;
int n,m,c[N][2],fa[N],h[N];
bool rev[N];
inline int read(){
int x=0,f=1; char ch=getchar();
while (ch==' ' || ch=='\r' || ch=='\n') ch=getchar();
if (ch=='+'){ f=1; ch=getchar(); }
if (ch=='-'){ f=-1; ch=getchar(); }
while (ch>='0' && ch<='9'){ x=x*10+ch-'0'; ch=getchar(); }
return x*f;
}
bool isrt(int x){
return c[fa[x]][0]!=x && c[fa[x]][1]!=x;
}
void pushdown(int x){
int l=c[x][0],r=c[x][1];
if (rev[x]){
rev[x]=0; rev[l]^=1; rev[r]^=1;
swap(c[x][0],c[x][1]);
}
}
void rotate(int x){
int y=fa[x],z=fa[y],l,r;
l=(c[y][0]==x)?0:1; r=l^1;
if (!isrt(y)){
if (c[z][0]==y) c[z][0]=x; else c[z][1]=x;
}
fa[x]=z; fa[y]=x; fa[c[x][r]]=y;
c[y][l]=c[x][r]; c[x][r]=y;
}
void splay(int x){
int top=1,p,i,y,z; h[1]=x;
for (p=x; !isrt(p); p=fa[p]) h[++top]=fa[p];
for (i=top; i; i--) pushdown(h[i]);
while (!isrt(x)){
y=fa[x]; z=fa[y];
if (!isrt(y)){
if ((c[y][0]==x)^(c[z][0]==y))
rotate(x); else rotate(y);
}
rotate(x);
}
}
void access(int x){
int y=0;
while (x){
splay(x); c[x][1]=y;
y=x; x=fa[x];
}
}
void rever(int x){
access(x); splay(x); rev[x]^=1;
}
void link(int x,int y){
rever(x); fa[x]=y; splay(x);
}
void cut(int x,int y){
rever(x); access(y); splay(y); fa[x]=c[y][0]=0;
}
int find(int x){
access(x); splay(x);
int y=x;
while (c[y][0]) y=c[y][0];
return y;
}
int main(){
n=read(); m=read();
int i,x,y; char ch[10];
for (i=1; i<=m; i++){
scanf("%s",ch);
x=read(); y=read();
if (ch[0]=='C') link(x,y);
else if (ch[0]=='D') cut(x,y);
else puts((find(x)==find(y))?"Yes":"No");
}
return 0;
}
2015.8.2
by lych