bzoj4568 [Scoi2016]幸运数字（树上倍增+线性基）

多次询问路径权值集合的子集的最大异或和。
很容易想到线性基。
我们树上倍增跳lca的同时，暴力合并出我们的线性基即可。
复杂度$O(nlogn60^2+mlogn60^2)$
按道理应该过不去，不过可以卡过去qaq
可以考虑写st表去掉一个logn。

tips：人生一大悲痛，原来我一直写的倍增是有问题的qaq，注意那个预处理时的跳倍增的i，我一直以为写成n还是log[n]是无所谓的，反正会在0退出。但是有可能越界啊啊！！！以后一定要细致一点qaq

#include <cstdio>
#include <cstring>
#include <algorithm>
#include <cmath>
using namespace std;
#define ll long long
#define inf 0x3f3f3f3f
#define N 20010
inline char gc(){
    static char buf[1<<16],*S,*T;
    if(T==S){T=(S=buf)+fread(buf,1,1<<16,stdin);if(S==T) return EOF;}
    return *S++;
}
inline ll read(){
    ll x=0,f=1;char ch=gc();
    while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=gc();}
    while(ch>='0'&&ch<='9') x=x*10+ch-'0',ch=gc();
    return x*f;
}
int n,m,h[N],num=0,fa[N][15],Log[N],dep[N];
ll a[N],p[N][15][65],ans[65];
struct edge{
    int to,next;
}data[N<<1];
inline void ins(ll *p,ll val){
    for(int i=60;i>=0;--i){
        if(!(val>>i&1)) continue;
        if(!p[i]){p[i]=val;break;}
        val^=p[i];
    }
}
inline void merge(ll *p1,ll *p2){
    for(int i=60;i>=0;--i) if(p2[i]) ins(p1,p2[i]);
}
void dfs(int x){
    for(int i=1;i<=Log[n];++i){
        if(!fa[x][i-1]) break;
        fa[x][i]=fa[fa[x][i-1]][i-1];
        memcpy(p[x][i],p[x][i-1],sizeof(p[x][i]));
        merge(p[x][i],p[fa[x][i-1]][i-1]);
    }for(int i=h[x];i;i=data[i].next){
        int y=data[i].to;if(y==fa[x][0]) continue;
        fa[y][0]=x;dep[y]=dep[x]+1;ins(p[y][0],a[x]);dfs(y);
    }
}
inline void lca(int x,int y){
    if(dep[x]<dep[y]) swap(x,y);
    int d=dep[x]-dep[y];ins(ans,a[x]);ins(ans,a[y]);
    for(int i=0;i<=Log[d];++i)
        if(d>>i&1) merge(ans,p[x][i]),x=fa[x][i];
    if(x==y) return;
    for(int i=Log[n];i>=0;--i)
        if(fa[x][i]!=fa[y][i])
            merge(ans,p[x][i]),x=fa[x][i],merge(ans,p[y][i]),y=fa[y][i];
    merge(ans,p[x][0]);
}
int main(){
//  freopen("a.in","r",stdin);
    n=read();m=read();Log[0]=-1;
    for(int i=1;i<=n;++i) a[i]=read(),Log[i]=Log[i>>1]+1;
    for(int i=1;i<n;++i){
        int x=read(),y=read();
        data[++num].to=y;data[num].next=h[x];h[x]=num;
        data[++num].to=x;data[num].next=h[y];h[y]=num;
    }dfs(1);
    while(m--){
        int x=read(),y=read();memset(ans,0,sizeof(ans));
        lca(x,y);ll sum=0;
        for(int i=60;i>=0;--i) if((sum^ans[i])>sum) sum^=ans[i];
        printf("%lld\n",sum);
    }return 0;
}