luogu题解 P3379 【模板】最近公共祖先（LCA）

看了看题解的dalao们,都是什么RMQ啊,在线的倍增啊(然后我都不会)

嗯,就发一篇tarjan的离线写法把(刚好最近写tarjan)

(第一篇题解)

嗯,大致的思路如下:

先选择一个节点u为根节点，从根节点开始搜索。（标记u已访问过）

遍历该点u的所有儿子节点v，并标记v已访问过。

若v还有儿子节点，对v重复ii操作，否则进入下一操作。

把v合并到u上（并查集）。

把当前的点设为u，遍历与u有询问关系的节点v。

如果v在之前已经被访问过，那么u和v的最近公共祖先就是v通过并查集合并后的父亲节点（注意是合并后），即当前的find（v）。

#include <cstdio>
#include <cstring>
#include <cmath>
#include <climits>
#include <iostream>
#include <algorithm>
#define ri register int 

using namespace std;

struct node{                                        //前向星 
    int next,to;
}edge[1000005];
int head[500005],cnt;

void add(int u,int v){
    edge[++cnt].to = v;
    edge[cnt].next = head[u];
    head[u] = cnt;
}

struct node2{
    int id;                                         //编号,方便记录答案 
    int to;
    int next;
}que[1000005];                                      //记录询问的结构体 
int ans[500005];
int headq[500005],cnt2;

void add_q(int u,int v,int i){
    que[++cnt2].to = v;
    que[cnt2].id = i;
    que[cnt2].next = headq[u];
    headq[u] = cnt2;
}

int n,m,s;
int fa[500005];
bool vis[500005];
int lca[500005];

int find(int x){
    if(x == fa[x])
        return x;
    return fa[x] = find(fa[x]);
}

void tarjan(int u){
    vis[u] = 1;
    for(ri i = head[u];i;i = edge[i].next){
        int v = edge[i].to;
        if(!vis[v]){                                        //如果没有被访问过,就tagjan这个点 
            tarjan(v);
            fa[v] = find(u);
        }
    }
    for(ri i = headq[u];i;i = que[i].next){
        int v = que[i].to;
        if(vis[v])                                          //找答案时,已经访问过的节点就可以记录lca了 
            lca[que[i].id] = find(v);
    }
}

int main(){
    ios::sync_with_stdio(false);
    cin >> n >> m >> s;
    int a,b;
    for(ri i = 1;i <= n-1;i ++){
        cin >> a >> b;
        add(a,b);
        add(b,a);
    }
    for(ri i = 1;i <= m;i ++){
        cin >> a >> b;
        add_q(a,b,i);
        add_q(b,a,i);
    }
    for(ri i = 1;i <= n;i ++)
        fa[i] = i;
    tarjan(s);
    for(ri i = 1;i <= m;i ++)
        cout << lca[i] << endl;
    return 0;
}