[SCOI2016] bzoj 4568 幸运数字 - 线性基 - 拟阵？

题目大意：多组询问点带权的树上链中选出一些点使得异或和最大。
题解：线性基具有拟阵性质，dfs预处理即可，询问暴力合并。
没了？
当然不是，首先把询问分成从LCA出发的两段，分别处理出来然后合并起来即可。
考虑类似保留区间的边问联通块个数那个题的维护方式，维护从每一个点向根的方向走的时候线性基以及每个数字的最早加入时间；但直接做是saygoodbye的，我们考虑维护从根开始dfs，然后每次有矛盾（即两个数字该位都是1）的时候，保留深度更大的（因为二者任一皆可，而询问的时候显然深度越大越有可能在线性基里面）。然后询问的时候查询线性基里面深度>=LCA的深度的所有数字即可，复杂度瓶颈在于合并线性基。当然也可以点分治，维护每个重心到其控制子树的每个点的线性基，然后询问的时候求一下两个点在点分树上的LCA，用LCA到这两个点的线性基合并起来询问也可以。
代码（第一种方法）：

#include<iostream>
#include<cstring>
#include<cstdio>
#include<algorithm>
#define LB LinearBase
#define N 20010
#define lint long long
#define gc getchar()
#define debug(x) cerr<<#x<<"="<<x
#define sp <<" "
#define ln <<endl
using namespace std;
lint v[N];int Log[N],d[N],up[N][20],LOG;
inline int inn()
{
    int x,ch;while((ch=gc)<'0'||ch>'9');
    x=ch^'0';while((ch=gc)>='0'&&ch<='9')
        x=(x<<1)+(x<<3)+(ch^'0');return x;
}
struct LinearBase{
    lint p[60+5];int d[60+5];
    LinearBase()
    {   memset(p,0,sizeof(p)),memset(d,0,sizeof(d));    }
    inline int insert(lint x,int dep)
    {
        if(!x) return 0;
        for(int i=LOG;i>=0;i--)
            if((x>>i)&1)
            {
                if(dep>d[i]) (p[i]?swap(p[i],x),0:p[i]=x),swap(d[i],dep);
                x^=p[i];if(!x) return 0;
            }
        return 0;
    }
    inline LB operator=(const LB &b)
    {
        memcpy(p,b.p,sizeof(lint)*(LOG+1)),
        memcpy(d,b.d,sizeof(int)*(LOG+1));
        return *this;
    }
    inline int show()
    {
        for(int i=1;i>=0;i--) debug(i)sp,debug(p[i])sp,debug(d[i])ln;
        cerr ln ln;return 0;
    }
}b[N];
inline LB merge_LB(LB a,const LB &b)
{
    for(int i=0;i<=LOG;i++)
        if(b.p[i]) a.insert(b.p[i],b.d[i]);
    return a;
}
inline lint query(const LB &b,int d,lint ans=0ll)
{
    for(int i=LOG;i>=0;i--)
        if(b.d[i]>=d&&(!((ans>>i)&1))) ans^=b.p[i];
    return ans;
}
struct edges{
    int to,pre;
}e[N<<1];int h[N],etop;
inline int add_edge(int u,int v)
{   return e[++etop].to=v,e[etop].pre=h[u],h[u]=etop;   }
void dfs(int x,int fa)
{
    b[x]=b[up[x][0]=fa],b[x].insert(v[x],d[x]=d[fa]+1);
    for(int i=1;i<=Log[d[x]];i++) up[x][i]=up[up[x][i-1]][i-1];
    for(int i=h[x],y;i;i=e[i].pre) if((y=e[i].to)^fa) dfs(y,x);
}
inline int getLCA(int x,int y)
{
    if(d[x]<d[y]) swap(x,y);
    for(int i=Log[d[x]];i>=0;i--)
        if(d[up[x][i]]>=d[y]) x=up[x][i];
    if(x==y) return x;
    for(int i=Log[d[x]];i>=0;i--)
        if(up[x][i]^up[y][i]) x=up[x][i],y=up[y][i];
    return up[x][0];
}
int main()
{
    int n=inn(),q=inn();
    for(int i=2;i<=n;i++) Log[i]=Log[i>>1]+1;
    for(int i=1;i<=n;i++) scanf("%lld",&v[i]);
    for(int i=1;i<=n;i++)
        for(int j=60;j>=0;j--)
            if(v[i]&(1ll<<j)) { LOG=max(LOG,j);break; }
    for(int i=1,u,v;i<n;i++) u=inn(),v=inn(),add_edge(u,v),add_edge(v,u);
    dfs(1,0);
//  for(int i=1;i<=n;i++) debug(i)sp,debug(d[i])sp,debug(up[i][0])sp,debug(up[i][1])ln;
    for(int i=1,x,y,c;i<=q;i++) c=getLCA(x=inn(),y=inn()),
//      debug(x)ln,b[x].show(),debug(y)ln,b[y].show(),merge_LB(b[x],b[y]).show(),
        printf("%lld\n",query(merge_LB(b[x],b[y]),d[c]));
    return 0;
}