ST表求LCA
好多天前写的……
dfs一遍,得到特殊的长度为 2*n-1 dfs序,然后维护一个在dfs序上的区间深度最小值,而拥有这个最小深度值的节点就是我们要求的LCA。
简单易懂。
#include<bits/stdc++.h>
#define RG register
using namespace std;
const int N=500000+5,P=20;
int to[N<<1],head[N<<1],nxt[N<<1],cn;
int logbx[N<<1],stu[N<<1][P+5];
int seq[N<<1],dep[N<<1],pos[N<<1];
int n,m,root,x,y,idc;
void create(int u,int v){
cn++;
to[cn]=v;
nxt[cn]=head[u];
head[u]=cn;
}
void init(){
logbx[0]=-1;
for(RG int i=1;i<N;i++)
logbx[i]=logbx[i>>1]+1;
for(RG int p=0;p<=P;p++)
for(RG int i=1;i+(1<<p)-1<=idc;i++)//i+(1<<p)-1<=idc惨痛教训!
if(p==0)
stu[i][p]=seq[i];//stu[i][p]=seq[i]惨痛教训!
else{
int l=stu[i][p-1],r=stu[i+(1<<(p-1))][p-1];//位运算一定要打上括号!吃屎吧优先级!
stu[i][p]=dep[l]<dep[r]?l:r;
}
}
void dfs(int u,int f){
int v;
for(RG int i=head[u];i;i=nxt[i]){
v=to[i];
if(v==f) continue;
seq[++idc]=v;//!!
pos[v]=idc;//!!
dep[v]=dep[u]+1;
dfs(v,u);
seq[++idc]=u;
pos[u]=idc;
}
}
int LCA(int u,int v){
u=pos[u],v=pos[v];
if(u>v) u^=v^=u^=v;
int d=logbx[v-u+1];
int x=stu[u][d],y=stu[v-(1<<d)+1][d];
return dep[x]<dep[y]?x:y;
}
int main(){
scanf("%d%d%d",&n,&m,&root);
for(RG int i=1;i<=n-1;i++){
scanf("%d%d",&x,&y);
create(x,y);
create(y,x);
}
seq[++idc]=root;
dep[root]=1;
dfs(root,root);
init();
for(RG int i=1;i<=m;i++){
scanf("%d%d",&x,&y);
printf("%d\n",LCA(x,y));
}
return 0;
}
ok,我还要更好多csdn啊。idy002都讲数论了我还在更原来的博客。
慢工细活
努力嘻嘻