前往小程序,Get更优阅读体验!
立即前往
发布
社区首页 >专栏 >RMQ求LCA

RMQ求LCA

作者头像
attack
发布2018-09-30 09:59:59
发布2018-09-30 09:59:59
58100
代码可运行
举报
运行总次数:0
代码可运行

题目链接

rmq求LCA,interesting。

一直没有学这玩意儿是因为CTSC的Day1T2,当时我打的树剖LCA 65分,gxb打的rmq LCA 45分。。。

不过rmq理论复杂度还是小一点的,就学一下把。

RMQ求LCA

我们要用到三个数组

$dfn[i]$:第$i$个节点位置的时间戳

$id[i][j]$:在欧拉序中$i$到$i + 2^j - 1$这段区间内深度最小的节点编号

$dep[i]$:第$i$个节点的深度

实际上用到了一个性质:

对于任意两点的$LCA$,一定是它们欧拉序中两点之间的最小值

欧拉序是什么?就是把dfs中遍历到每一个一个节点(包括回溯时遍历到)加到一个序列里,最终得到的就是欧拉序

时空复杂度

设$T = 2 * n - 1$

时间复杂度:

预处理:$O(TlogT)$

查询:$O(1)$

空间复杂度:

考虑欧拉序中有多少个点,首先每个点被访问到的时候会做出$1$的贡献

其次在遍历每条边时会多出$1$的共贡献

因此总空间复杂度为:$O(T)$

代码语言:javascript
代码运行次数:0
复制
// luogu-judger-enable-o2
#include<bits/stdc++.h>
const int MAXN = 1e6 + 10;
using namespace std;
inline int read() {
    char c = getchar(); int x = 0, f = 1;
    while(c < '0' || c > '9') {if(c == '-') f = -1; c = getchar();}
    while(c >= '0' && c <= '9') x = x * 10 + c - '0', c = getchar();
    return x * f;
}
int N, Q, S, tot, dfn[MAXN], rev[MAXN], dep[MAXN], id[MAXN][21], lg2[MAXN], rd[MAXN];
vector<int> v[MAXN];
void dfs(int x, int fa) {
    dfn[x] = ++tot; dep[x] = dep[fa] + 1; id[tot][0] = x; 
    for(int i = 0, to; i < v[x].size(); i++) {
        if((to = v[x][i]) == fa) continue;
        dfs(to, x);
        id[++tot][0] = x;
    }
}
void RMQ() {
    for(int i = 2; i <= tot; i++) lg2[i] = lg2[i >> 1] + 1;
    for(int j = 1; j <= 20; j++) {
        for(int i = 1; (i + (1 << j) - 1) <= tot; i++) {
            int r = i + (1 << (j - 1));
            id[i][j] = dep[id[i][j - 1]] < dep[id[r][j - 1]] ? id[i][j - 1] : id[r][j - 1];
        }
    }
}
int Query(int l, int r) {
    if(l > r) swap(l, r);
    int k = lg2[r - l + 1];
    return dep[id[l][k]] < dep[id[r - (1 << k) + 1][k]] ? id[l][k] : id[r - (1 << k) + 1][k];
}
int main() {
    freopen("a.in", "r", stdin);
    N = read(); Q = read(); S = read();
    for(int i = 1; i <= N - 1; i++) {
        int x = read(), y = read();
        v[x].push_back(y); v[y].push_back(x);
    }
    dfs(S, 0);
    RMQ();
    while(Q--) {
        int x = read(), y = read();
        printf("%d\n", Query(dfn[x], dfn[y]));
    }
    return 0;
}
本文参与 腾讯云自媒体同步曝光计划,分享自作者个人站点/博客。
原始发表:2018-09-28 ,如有侵权请联系 cloudcommunity@tencent.com 删除

本文分享自 作者个人站点/博客 前往查看

如有侵权,请联系 cloudcommunity@tencent.com 删除。

本文参与 腾讯云自媒体同步曝光计划  ,欢迎热爱写作的你一起参与!

评论
登录后参与评论
0 条评论
热度
最新
推荐阅读
目录
  • RMQ求LCA
  • 时空复杂度
领券
问题归档专栏文章快讯文章归档关键词归档开发者手册归档开发者手册 Section 归档