【代码源 Div1 - 105】#451. Dis(倍增求LCA)

举报
小哈里 发表于 2022/05/11 00:20:28 2022/05/11
【摘要】 problem solution 给出 n 个点的一棵树,每个点有各自的点权,m 次询问两个点简单路径所构成点集的异或和。直接在树上求LCA,把每个点权放进去预处理一下即可。 #include&...

problem

在这里插入图片描述

solution

  • 给出 n 个点的一棵树,每个点有各自的点权,m 次询问两个点简单路径所构成点集的异或和。
  • 直接在树上求LCA,把每个点权放进去预处理一下即可。
#include<bits/stdc++.h>
using namespace std;

#define IOS ios::sync_with_stdio(0), cin.tie(0),cout.tie(0)
typedef long long LL;
const LL maxn = 2e5+10;

int n, m;
int a[maxn];
vector<int>G[maxn];

int fa[maxn][20], dep[maxn], val[maxn][20];
void dfs(int u, int f, int d){
    fa[u][0] = f;
    val[u][0] = a[f];
    dep[u] = d;
    for(int x : G[u]){
        if(x==f)continue;
        dfs(x,u,d+1);
    }
}
void init(){
    for(int i = 1; i < 20; i++){
        for(int j = 1; j <= n; j++){
            fa[j][i] = fa[fa[j][i-1]][i-1];
            val[j][i] = val[j][i-1]^val[fa[j][i-1]][i-1];
        }
    }
}
int lca(int x, int y){
    if(dep[x]<dep[y])swap(x,y);
    int ans = 0;
    for(int i = log2(dep[x]-dep[y]); i >= 0; i--){
        if((1<<i)<=dep[x]-dep[y])ans ^= val[x][i], x = fa[x][i];
    }
    if(x == y)return ans^a[y];
    for(int i = log2(dep[x]); i >= 0; i--){
        if(fa[x][i] != fa[y][i]){
            ans ^= val[x][i];
            ans ^= val[y][i];
            x = fa[x][i];
            y = fa[y][i];
        }
    }
    ans ^= val[x][0];
    return ans;
}

int main(){
    IOS;
    cin>>n>>m;
    for(int i = 1; i <= n; i++)cin>>a[i];
    for(int i = 1; i < n; i++){
        int u, v;  cin>>u>>v;
        G[u].push_back(v);
        G[v].push_back(u);
    }
    dfs(1,1,0);
    init();
    while(m--){
        int u, v;  cin>>u>>v;
        int ans = lca(u,v)^a[u]^a[v];
        if(u==v)cout<<a[u]<<"\n";
        else cout<<ans<<'\n';
    }
    return 0;
}


  
 
  • 1
  • 2
  • 3
  • 4
  • 5
  • 6
  • 7
  • 8
  • 9
  • 10
  • 11
  • 12
  • 13
  • 14
  • 15
  • 16
  • 17
  • 18
  • 19
  • 20
  • 21
  • 22
  • 23
  • 24
  • 25
  • 26
  • 27
  • 28
  • 29
  • 30
  • 31
  • 32
  • 33
  • 34
  • 35
  • 36
  • 37
  • 38
  • 39
  • 40
  • 41
  • 42
  • 43
  • 44
  • 45
  • 46
  • 47
  • 48
  • 49
  • 50
  • 51
  • 52
  • 53
  • 54
  • 55
  • 56
  • 57
  • 58
  • 59
  • 60
  • 61
  • 62
  • 63
  • 64
  • 65
  • 66
  • 67
  • 68

文章来源: gwj1314.blog.csdn.net,作者:小哈里,版权归原作者所有,如需转载,请联系作者。

原文链接:gwj1314.blog.csdn.net/article/details/123953490

【版权声明】本文为华为云社区用户转载文章,如果您发现本社区中有涉嫌抄袭的内容,欢迎发送邮件进行举报,并提供相关证据,一经查实,本社区将立刻删除涉嫌侵权内容,举报邮箱: cloudbbs@huaweicloud.com
  • 点赞
  • 收藏
  • 关注作者

评论(0

0/1000
抱歉,系统识别当前为高风险访问,暂不支持该操作

全部回复

上滑加载中

设置昵称

在此一键设置昵称,即可参与社区互动!

*长度不超过10个汉字或20个英文字符,设置后3个月内不可修改。

*长度不超过10个汉字或20个英文字符,设置后3个月内不可修改。