f[i][j]表示从i这个节点出发,向上走2^j步到达的点(超过了最大深度就返回0)
那么显然有。
因为相当于是从i点先向上走2^(j-1)步,再走2^(j-1)步。等价于一共走2^j步。
然后求lca就是:先把深度大的那个点往上跳,使两个点的深度相同。
对应的代码:
for(int i=18;i>=0;i--) if(dep[f[x][i]]>=dep[y]) x=f[x][i];
可以把两点之间的深度差理解为一个二进制数。
举个例子:一个点深度为3,另一个点深度为16。那么深度为16的点向上跳13步它们的深度就相同了。
13转换成二进制为1101。那么这个跳的过程就是先跳(1000)步,再跳(100)步,
然后跳(10)步,发现深度超了,就不会跳。
最后跳(1)步,它们深度就相同了。
之后特判一下这时是否相同。
如果不相同:它们就一起往上跳。
如果发现 f[x][i]==f[y][i] ,那么说明它们跳到了lca或者是超过了lca,就不跳。(当前深度≤lca深度)
但是如果 f[x][i]!=f[y][i] ,那么它们一定还都在lca下面,就继续跳。(当前深度>lca深度)
for(int i=18;i>=0;i--) if(f[x][i]!=f[y][i]) x=f[x][i],y=f[y][i];
这里跳完之后,可以发现,它们一定刚好跳到了lca的下面一个。(可以把深度想成二进制思考。)
那么它们的父亲就是lca。
void dfs(int u,int f){
dep[u]=dep[f]+1,f[u][0]=f;
for(int i=Head[u];i;i=Next[i])
dfs(V[i],u);
}
void pre(){
dfs(1,0);
for(int i=1;(1<<i)<=N;++i)
for(int u=1;u<=N;++u)
f[u][i]=f[f[u][i-1]][i-1];
}
int lca(int x,int y){
if(dep[x]<dep[y]) swap(x,y);
for(int i=18;i>=0;i--) if(dep[f[x][i]]>=dep[y]) x=f[x][i];
if(x==y) return x;
for(int i=18;i>=0;i--) if(f[x][i]!=f[y][i]) x=f[x][i],y=f[y][i];
return f[x][0];
}