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

P3379 【模板】最近公共祖先（LCA）
提交
129.76k
通过
45.50k
时间限制
1.50s
内存限制
500.00MB
提交答案
加入收藏
题目提供者
HansBug
难度
普及/提高-
历史分数
100
提交记录 查看题解
标签
O2优化
进入讨论版
相关讨论
推荐题目
展开
题目描述
如题，给定一棵有根多叉树，请求出指定两个点直接最近的公共祖先。

输入格式
第一行包含三个正整数 N,M,SN,M,S，分别表示树的结点个数、询问的个数和树根结点的序号。

接下来 N-1N−1 行每行包含两个正整数 x, yx,y，表示 xx 结点和 yy 结点之间有一条直接连接的边（数据保证可以构成树）。

接下来 MM 行每行包含两个正整数 a, ba,b，表示询问 aa 结点和 bb 结点的最近公共祖先。

输出格式
输出包含 MM 行，每行包含一个正整数，依次为每一个询问的结果。

输入输出样例
输入 #1复制
5 5 4
3 1
2 4
5 1
1 4
2 4
3 2
3 5
1 2
4 5
输出 #1复制
4
4
1
4
4
说明/提示
对于 30%30% 的数据，N\leq 10N≤10，M\leq 10M≤10。

对于 70%70% 的数据，N\leq 10000N≤10000，M\leq 10000M≤10000。

对于 100%100% 的数据，N\leq 500000N≤500000，M\leq 500000M≤500000。

样例说明：

该树结构如下：

第一次询问：2, 4 的最近公共祖先，故为 4。

第二次询问：3, 2 的最近公共祖先，故为 4。

第三次询问：3, 5 的最近公共祖先，故为 1。

第四次询问：1, 2的最近公共祖先，故为 4。

第五次询问：4, 5的最近公共祖先，故为 4。

故输出依次为 4, 4, 1, 4, 4

#include <iostream>
#include <string>
#include <cstring>
#include <algorithm>
#include <cstdio>
#include <vector>
#include <queue>
#include <climits>
#include <stack>
#include <map>
#include <cmath>
#include <set>
#include <cstdlib>
#include <unordered_map>
using namespace std;
typedef long long ll;
typedef unsigned long long ull;
#define rep(i,a,b) for(int i=a;i<=b;i++)
#define per(i,a,b) for(int i=a;i>=b;i--)
#define ios ios::sync_with_stdio(false);cin.tie(0);cout.tie(0);
template<typename T>void write(T x)
{
    if (x<0)
    {
        putchar('-');
        x=-x;
    }
    if (x>9)
    {
        write(x/10);
    }
    putchar(x%10+'0');
}

template<typename T> void read(T &x)
{
    x = 0;char ch = getchar();ll f = 1;
    while (!isdigit(ch)) {        
if (ch == '-')f*=-1;ch=getchar();    
}
    while (isdigit(ch)){        
x = x*10+ch-48;ch=getchar();    
}x*=f;
}
const int maxn=5e5+5;
struct node{
    int to,next;
}e[maxn<<1];
int head[maxn],tot;
int dep[maxn],lg[maxn];
int anc[maxn][25];
void add(int x,int y){
   e[++tot]={y,head[x]};
   head[x]=tot;
}
int n,m,s;
void dfs(int u,int fa){
   dep[u]=dep[fa]+1;
   anc[u][0]=fa;
   for(int i=1;(1<<i)<=dep[u];i++){
         anc[u][i]=anc[anc[u][i-1]][i-1];
   }
   for(int i=head[u];i;i=e[i].next){
       int v=e[i].to;
       if(v==fa)continue;
       dfs(v,u);
   }
}
int lca(int x,int y){
  if(dep[x]<dep[y])swap(x,y);
  while(dep[x]>dep[y]){
      x=anc[x][lg[dep[x]-dep[y]]-1];
  }
  if(x==y)return x;
  for(int i=lg[dep[x]]-1;i>=0;i--){
      if(anc[x][i]!=anc[y][i]){
      x=anc[x][i];
      y=anc[y][i];
      }
  }
  return anc[x][0];
}
int main() {
  read(n);
  read(m);
  read(s);
  int a,b;
  rep(i,1,n-1){
     read(a);
     read(b);
     add(a,b);
     add(b,a);
  }
  rep(i,1,n){
      lg[i]=lg[i-1]+((1<<lg[i-1])==i);
  }
  dfs(s,0);
rep(i,1,m){
    read(a);
    read(b);
    write(lca(a,b));
   printf("\n");
}
  return 0;
}