【代码源 Div1 - 105】#451. Dis(倍增求LCA)
【摘要】
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)