SCOI2016幸运数字

题目大意

有一棵有n个节点的树，每个节点有一个权值，有若干组询问，每组询问询问x到y的路径上的点中选出一些点的点权来抑或，求最大的值。

线性基

（属个人理解，如有误，请指出
对于一个集合$S$，它的线性基$S^`$的意义为：集合$S^`$中的若干元素的抑或和得到的集合与$S$中的若干元素的抑或和得到的集合完全相同。

ST表+暴力合并线性基

（如有不懂，详情见标程
由于线性基的特殊性质，我们可以直接从大到小抑或，如果当前答案更优则更新
用$S_{i,j}$表示i节点向上跳长度为$2^j$的数的集合的线性基，对于合并，明显的，我们可以用$S_{i,j-1}$和$S_{fa_{i,j-1},j-1}$来合并，类似高斯消元，我们每次加入一个新的数然后逐个逐个看是否要抑或，如果最后为0则不加入，否则加入当前数。
询问的话，我们可以将一个询问的路径拆成四段，如下图：

然后合并就类似预处理囖。

贴代码哈哈哈哈哈：

#include<cstring>
#include<cstdio>
#include<algorithm>
#include<cmath>
#include<iostream>

#define fo(i,a,b) for(int i=a;i<=b;i++)
#define fd(i,a,b) for(int i=a;i>=b;i--)

using namespace std;

typedef long long LL;
typedef double db;

int get(){
    char ch;
    int s=0;
    bool pd=0;
    while(ch=getchar(),(ch<'0'||ch>'9')&&ch!='-');
    if (ch=='-')pd=1;
    else s=ch-'0';
    while(ch=getchar(),ch>='0'&&ch<='9')s=s*10+ch-'0';
    if (pd)return -s;
    return s;
}

const int lim = 60;
const int N = 20010;

struct linear_basic{
    LL v[lim+2];
}a[N][15],tmp;//线性基
int n,q,fa[N][15],dep[N],tim;
LL g[N];
LL mi[lim+2];
struct edge{
    int x,next;
}e[N*2];
int h[N],tot;

void inse(int x,int y){
    e[++tot].x=y;
    e[tot].next=h[x];
    h[x]=tot;
}

void dfs(int x){
    for(int p=h[x];p;p=e[p].next)
    if (!dep[e[p].x]){
        dep[e[p].x]=dep[x]+1;
        fa[e[p].x][0]=x;
        dfs(e[p].x);
    }
}

void merge(linear_basic a,linear_basic b,linear_basic &c){
    c=a;
    fd(i,lim,0)
    if (b.v[i]){
        LL tp=b.v[i];
        fd(i,lim,0)
        if (tp&mi[i])
            if (!c.v[i]){c.v[i]=tp;break;}
            else tp^=c.v[i];
    }
}//合并线性基

void prepare(){
    fo(i,1,n)
        fd(j,lim,0)
        if (g[i]&mi[j]){
            a[i][0].v[j]=g[i];
            break;
        }
    fo(j,1,tim)
        fo(i,1,n){
            fa[i][j]=fa[fa[i][j-1]][j-1];
            merge(a[i][j-1],a[fa[i][j-1]][j-1],a[i][j]);
        }
}

int getfather(int x,int y){
    if (dep[x]<dep[y])swap(x,y);
    fd(i,tim,0)if (dep[fa[x][i]]>=dep[y])x=fa[x][i];
    if (x==y)return x;
    fd(i,tim,0)
    if (fa[x][i]!=fa[y][i])x=fa[x][i],y=fa[y][i];
    return fa[x][0];
}

int jump(int x,int v){
    fd(i,tim,0)
    if (v>=mi[i]){
        v-=mi[i];
        x=fa[x][i];
    }
    return x;
}

int main(){
    n=get();q=get();
    fo(i,1,n)scanf("%lld",&g[i]);
    fo(i,2,n){
        int x=get(),y=get();
        inse(x,y);
        inse(y,x);
    }
    tim=log(n)/log(2);
    mi[0]=1;
    fo(i,1,lim)mi[i]=mi[i-1]<<1;
    dfs(dep[1]=1);
    prepare();
    fo(i,1,q){
        int x=get(),y=get();
        if (x==y){
            printf("%lld\n",g[x]);
            continue;
        }
//----------------------------我是美丽的分割线^_^--------------------------------
        int t=getfather(x,y);
        fd(v,lim,0)tmp.v[v]=0;
        if (x!=t){
            int v=log(dep[x]-dep[t]+1)/log(2);
            merge(a[x][v],tmp,tmp);
            merge(a[jump(x,dep[x]-dep[t]+1-(1<<v))][v],tmp,tmp);
        }
        if (y!=t){
            int v=log(dep[y]-dep[t]+1)/log(2);
            merge(a[y][v],tmp,tmp);
            merge(a[jump(y,dep[y]-dep[t]+1-(1<<v))][v],tmp,tmp);
        }
//------------------------------以上为合并四段-----------------------------------
        LL ans=0;
        fd(i,lim,0)
        if ((ans^tmp.v[i])>ans)ans^=tmp.v[i];
        printf("%lld\n",ans);
    }
    return 0;
}