[Scoi2016]幸运数字

最新推荐文章于 2019-01-03 19:37:00 发布

quan_tum

最新推荐文章于 2019-01-03 19:37:00 发布

阅读量415

点赞数

CC 4.0 BY-SA版权

分类专栏：线性基

本文链接：https://blog.youkuaiyun.com/quan_tum/article/details/82633131

线性基专栏收录该内容

6 篇文章

订阅专栏

本文介绍了一种解决树上路径查询问题的方法，通过使用线性基合并技术来优化路径上的查询操作。文章详细阐述了如何利用树的倍增、LCA算法以及线性基的数据结构来高效地处理特定类型的查询请求。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

传送门

题意是求树上一条路径上，于是树上倍增lca暴力合并线性基…

听别人说要加上读优否则过不了…

不会求lca的出门左转

不会线性基的出门右转

真的如果都会了就只剩一个暴力合并了…

#include<bits/stdc++.h>
#define il inline
#define ll long long
#define int ll
#define max(a,b) ((a)>(b))?(a):(b)
using namespace std;
const int N=20005;
#define getchar()(p1==p2&&(p2=(p1=buf)+fread(buf,1,1<<21,stdin),p1==p2)?EOF:*p1++)
char buf[1<<21],*p1=buf,*p2=buf;
il ll read(){
    ll x=0,f=1;char c=getchar();
    for(;!isdigit(c);c=getchar()) if(c=='-') f=-1;
    for(;isdigit(c);c=getchar()) x=(x+(x<<2)<<1)+c-48;
    return x*f;
}
char sr[1<<21],z[20];int C=-1,Z;
il void Ot(){fwrite(sr,1,C+1,stdout),C=-1;}
il void print(ll x){
    if(C>1<<20)Ot();if(x<0)sr[++C]=45,x=-x;
    while(z[++Z]=x%10+48,x/=10);
    while(sr[++C]=z[Z],--Z);sr[++C]='\n';
}
struct Edge{int to,nxt;}a[N<<1];
int n,q,k,h[N],f[N][21],w[N],d[N];
ll ans[N],p[N][21][62],sum;
il void Lb(ll *a,ll val){for(int i=61;i>=0;--i)if((val>>i)&1){if(!a[i]){a[i]=val;break;}val^=a[i];}}
il void add(int x,int y){a[++k].to=y,a[k].nxt=h[x],h[x]=k;}
il void dfs(int x,int fa){
    f[x][0]=fa;d[x]=d[fa]+1;
    for(int i=h[x];i;i=a[i].nxt)
    if(a[i].to!=fa) dfs(a[i].to,x);
}
il void mg(ll *a,ll *b){for(int i=61;i>=0;--i) if(b[i]) Lb(a,b[i]);}
il void lca(int x,int y){
    if(d[x]<d[y]) swap(x,y);
    for(int i=20;i>=0;--i)
    if(d[f[x][i]]>=d[y]) mg(ans,p[x][i]),x=f[x][i];
    if(x==y){mg(ans,p[x][0]);return;}
    for(int i=20;i>=0;--i)
    if(f[x][i]!=f[y][i]) mg(ans,p[x][i]),mg(ans,p[y][i]),x=f[x][i],y=f[y][i];
    mg(ans,p[x][0]),mg(ans,p[y][0]),mg(ans,p[f[x][0]][0]);
}
signed main(){
    n=read(),q=read();
    for(int i=1;i<=n;++i) Lb(p[i][0],w[i]=read());
    for(int i=1;i<n;++i){
        int x=read(),y=read();
        add(x,y),add(y,x);
    }
    dfs(1,0);
    for(int j=1;j<20;++j)
    for(int i=1;i<=n;++i){
        f[i][j]=f[f[i][j-1]][j-1];
        memcpy(p[i][j],p[i][j-1],sizeof(p[i][j-1]));
        mg(p[i][j],p[f[i][j-1]][j-1]);
    }
    while(q--){
        memset(ans,0,sizeof(ans));
        int x=read(),y=read();
        lca(x,y);sum=0;
        for(int j=61;j>=0;--j) sum=max(sum,sum^(ll)ans[j]);
        print(sum);
    }Ot();return 0;
}