洛谷P3292 [SCOI2016]幸运数字（倍增+线性基）

传送门

 

 

不知道线性基是什么东西的可以看看蒟蒻的总结

第一眼：这不会是个倍增LCA暴力合并线性基吧……

打了一发……A了？

所以这真的是个暴力倍增LCA合并线性基么……

ps：据某大佬说其实可以离线之后用点分做，那样的话因为每次只要合并两个线性基，复杂度可以减一个$log$

//minamoto
#include<iostream>
#include<cstring>
#include<cstdio>
#define ll long long
#define max(a,b) ((a)>(b)?(a):(b))
#define getc() (p1==p2&&(p2=(p1=buf)+fread(buf,1,1<<21,stdin),p1==p2)?EOF:*p1++)
char buf[1<<21],*p1=buf,*p2=buf;
template<class T>inline bool cmax(T&a,const T&b){return a<b?a=b,1:0;}
using namespace std;
inline ll read(){
    #define num ch-'0'
    char ch;bool flag=0;ll res;
    while((ch=getc())>'9'||ch<'0')
    (ch=='-')&&(flag=true);
    for(res=num;(ch=getc())<='9'&&ch>='0';res=res*10+num);
    (flag)&&(res=-res);
    #undef num
    return res;
}
char sr[1<<21],z[20];int C=-1,Z;
inline void Ot(){fwrite(sr,1,C+1,stdout),C=-1;}
inline 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';
}
const int N=20005;
int n,q,tot,head[N],Next[N<<1],ver[N<<1],dep[N];
ll fa[N][21],b[N][21][62],sum,ans[62],val[N];
inline void add(int u,int v){
    ver[++tot]=v,Next[tot]=head[u],head[u]=tot;
}
inline void get(ll *b,ll x){
    for(int i=61;i>=0;--i)
    if(x>>i&1){
        if(!b[i]) return (void)(b[i]=x);
        x^=b[i];
    }
}
inline void merge(ll *b,ll *x){
    for(int i=61;i>=0;--i)
    if(x[i]) get(b,x[i]);
}
inline void init(int i){
    for(int j=1;j<20;++j){
        fa[i][j]=fa[fa[i][j-1]][j-1];
        memcpy(b[i][j],b[i][j-1],sizeof(b[i][j-1]));
        merge(b[i][j],b[fa[i][j-1]][j-1]);
    }
}
void dfs(int u,int f){
    fa[u][0]=f,dep[u]=dep[f]+1,init(u);
    for(int i=head[u];i;i=Next[i])
    if(ver[i]!=f) dfs(ver[i],u);
}
void LCA(int u,int v){
    if(dep[u]<dep[v]) swap(u,v);
    for(int i=20;i>=0;--i)
    if(dep[fa[u][i]]>=dep[v])
    merge(ans,b[u][i]),u=fa[u][i];
    if(u==v) return (void)(merge(ans,b[u][0]));
    for(int i=20;i>=0;--i)
    if(fa[u][i]!=fa[v][i]){
        merge(ans,b[u][i]),merge(ans,b[v][i]);
        u=fa[u][i],v=fa[v][i];
    }
    merge(ans,b[u][0]),merge(ans,b[v][0]),merge(ans,b[fa[u][0]][0]);
}
int main(){
//    freopen("testdata.in","r",stdin);
    n=read(),q=read();
    for(int i=1;i<=n;++i)
    get(b[i][0],val[i]=read());
    for(int i=1,u,v;i<n;++i)
    u=read(),v=read(),add(u,v),add(v,u);
    dfs(1,0);
    while(q--){
        memset(ans,0,sizeof(ans));
        int u=read(),v=read();
        LCA(u,v);
        sum=0;
        for(int i=61;i>=0;--i)
        cmax(sum,sum^ans[i]);
        print(sum);
    }
    Ot();
    return 0;
}