呃,这个常用但是我一直不会
Tarjan
Tarjan 算法基于 dfs ,在 dfs 的过程中,对于每个节点位置的询问做出相应的回答。
dfs 的过程中,当一棵子树被搜索完成之后,就把他和他的父亲合并成同一集合;在搜索当前子树节点的询问时,如果该询问的另一个节点已经被访问过,那么该编号的询问是被标记了的,于是直接输出当前状态下,另一个节点所在的并查集的祖先;如果另一个节点还没有被访问过,那么就做下标记,继续 dfs 。
//tarjian,边建边 边回答问题
#include<cstdio>
#include<cstdlib>
#include<vector>
using namespace std;
int n,m,rt;
const int N=500003,M=500003;
int ans[M];
vector <int> e[N];
struct node
{
int v,id;
node(int vv,int ii)
{ v=vv,id=ii; }
node(){}
};
vector <node> q[N];
int fa[N];
int find(int x)
{ return fa[x]==0?x:fa[x]=find(fa[x]); }
bool vis[N];
void dfs(int x,int f)
{
int sz=e[x].size();
for(int i=0;i<sz;i++)
{
int v=e[x][i];
if(v!=f) dfs(v,x),fa[v]=x;
}
vis[x]=true;
sz=q[x].size();
for(int i=0;i<sz;i++)
{
int v=q[x][i].v ;
if(!ans[q[x][i].id ] && vis[v])
ans[q[x][i].id ]=find(v);
}
}
int main()
{
scanf("%d%d%d",&n,&m,&rt);
int u,v;
for(int i=1;i<n;i++)
scanf("%d%d",&u,&v),e[u].push_back(v),e[v].push_back(u);
for(int i=0;i<m;i++)
scanf("%d%d",&u,&v),q[u].push_back(node(v,i)),q[v].push_back(node(u,i));
dfs(rt,0);
for(int i=0;i<m;i++) printf("%d\n",ans[i]);
return 0;
}
倍增
就从父亲更新到儿子
//tarjian,边建边 边回答问题
#include<cstdio>
#include<cstdlib>
#include<vector>
#include<algorithm>
using namespace std;
int n,m,rt;
const int N=500003,M=500003;
int ans[M];
vector <int> e[N];
int lg[N];
inline int read()
{
int x=0;char c=getchar();
while(c<'0' || c>'9') c=getchar();
while(c>='0' && c<='9') x=(x<<1)+(x<<3) +c-'0',c=getchar();
return x;
}
int dep[N];
int fa[N][19];//2的0-18次方,,,全跳最高能到2^19
void dfs(int x,int f)
{
fa[x][0]=f;
dep[x]=dep[f]+1;
for(int i=1;i<19;i++)
fa[x][i]=fa[fa[x][i-1]][i-1];
int sz=e[x].size();
for(int i=0;i<sz;i++)
{
int v=e[x][i];
if(v!=f) dfs(v,x);
}
}
void get_log()
{
for(int i=1;i<=n;i++) //预先算出log_2(i)+1的值,用的时候直接调用就可以了
lg[i]=lg[i-1]+(1<<lg[i-1]==i); //看不懂的可以手推一下
}//还挺快
int LCA(int u,int v)
{
if(dep[u]<dep[v]) swap(u,v);
int dis=dep[u]-dep[v];
for(int i=0;dis;i++)
if(dis & (1<<i))
{
dis^=(1<<i);
u=fa[u][i];
}
if(u==v) return u;
for(int i=lg[dep[u]]-1;i>=0;i--)
if(fa[u][i]!=fa[v][i])
u=fa[u][i],v=fa[v][i];
return fa[u][0];
}
int main()
{
n=read(),m=read(),rt=read();
int u,v;
for(int i=1;i<n;i++)
u=read(),v=read(),e[u].push_back(v),e[v].push_back(u);
dfs(rt,0);
get_log();
for(int i=0;i<m;i++)
{
u=read(),v=read();
printf("%d\n",LCA(u,v));
}
return 0;
}