原篇：（ACM算法）tarjan算法求LCA - 知乎 (zhihu.com)

顾名思义，就是求两个节点最近的共同祖先，就好比下图，2和3的共同祖先为3，2和4的共同祖先为1。

关于LCA求解有3种算法。

1.标记回溯法（也称暴力枚举，从一个点开始向上标记他的父节点，直接标记到根节点为止，然后另另一个点也开始向上回溯，当这个点被标记过，则找到这两个点的共同祖先），当数据过大时，很明显这个算法会超时，所以本文不过多解释此方法。

2.ST倍增法（本文也不重点介绍次方法，读者可看其他博客阅读此方法）

3.targan算法（本文重点介绍）

tarjan算法就这几个步骤：
1.标记当前节点x已访问，但是还没有回溯

2.枚举所有没有访问的过的子节点y，继续遍历子节点

3.合并y到x上

4.访问所有与当前节点有询问并且回溯过的y x,y的lca就是find(y)

5.标记当前节点回溯

首先targan是基于并查集实现的，没学过的读者可点击此传送门：何为并查集？-CSDN博客

具体过程看以下代码：

#include<bits/stdc++.h>
#define int long long 
#define endl "\n"
#define fi first
#define se second
#define PII pair<int,int> 
using namespace std;
const int N=5e5+5;
int fa[N],vis[N];
vector<int> v[N];
vector<PII> query[N];
int ans[N]; 
int find(int x){//并查集模板 if(fa[x]==x)  return x;else return fa[x]=find(fa[x]);
}
void tarjan(int x){vis[x]++;//第一次走过的路径标记为1 for(auto i :v[x]){if(!vis[i]){//若子节点没走过 继续遍历子节点 tarjan(i);fa[i]=x;//合并子节点到当前节点 }		}for(auto t : query[x]){int a=t.fi,b=t.se;if(vis[a]==2&&!ans[b]) ans[b]=find(a);// 若已经回溯并且b下标为空 }vis[x]++;//回溯的路径标记为2 
}
void solve(){	int n,m,q;cin >> n >> m >> q;for(int i=1;i<n;++i){int x,y;cin >> x >> y;v[x].push_back(y);//将x，y存入数组中 v[y].push_back(x);}for(int i=1;i<=n;++i) fa[i]=i;//预处理所有节点都是自己的父节点 for(int i=1;i<=m;++i){int x,y;cin >> x >> y;if(x==y){//特判 ans[i]=x;continue;}query[x].push_back({y,i});//存入查询的数字 query[y].push_back({x,i});}tarjan(q);for(int i=1;i<=m;++i){cout << ans[i] << endl;}return ;
}
signed main(){ios::sync_with_stdio(false);cin.tie(0);cout.tie(0);int t=1;//cin >> t;while(t--) solve();return 0;
}

经典例题：P3884 [JLOI2009] 二叉树问题 - 洛谷 | 计算机科学教育新生态 (luogu.com.cn)

解题思路：用层序遍历的思想求出树的深度和宽度，同时标记每个节点所在的层数，然后用tarjan算法求出a和b的共同祖先，最后套公式输出距离即可。

实现代码：

#include<bits/stdc++.h>
#define int long long 
#define endl "\n"
#define fi first
#define se second
#define PII pair<int,int> 
using namespace std;
const int N=105;
vector<int> v[N];
vector<int> s[N];
queue<PII> q;
int a,b;
int vis[N],fa[N],c[N];
vector<PII> query[N];
int ans1,ans2;
int k;
PII pa;
void check(int x){q.push({1,1});c[1]=1;int sum=0;int t=1;while(!q.empty()){pa=q.front();q.pop();if(pa.se==t) sum++;else{ans2=max(ans2,sum);t++;sum=1;}c[pa.fi]=pa.se;for(int i=0;i<s[pa.fi].size();++i){q.push({s[pa.fi][i],pa.se+1});}} ans1=t;ans2=max(ans2,sum);
}
int find(int x){if(fa[x]==x)  return x;else return fa[x]=find(fa[x]);
}
void tarjan(int x){vis[x]++; for(auto i :v[x]){if(!vis[i]){tarjan(i);fa[i]=x;}		}for(auto t : query[x]){int a=t.fi;if(vis[a]==2){k=find(a);break;} }vis[x]++;
}
void solve(){int n;cin >> n;for(int i=1;i<n;++i){int x,y;cin >> x >> y;v[x].push_back(y);v[y].push_back(x);s[x].push_back(y);}	for(int i=1;i<=n;++i) fa[i]=i; cin >> a >> b;query[a].push_back({b,1});query[b].push_back({a,1});check(1);cout << ans1 << endl << ans2 << endl;tarjan(1);//for(int i=1;i<=n;++i) cout << c[i] << endl; cout << (c[a]-c[k])*2+(c[b]-c[k]);return ;
}
signed main(){ios::sync_with_stdio(false);cin.tie(0);cout.tie(0);int t=1;//cin >> t;while(t--) solve();return 0;
}