BZOJ 4568 倍增维护线性基

//By SiriusRen
#include <cstdio>
#include <cstring>
#include <algorithm>
using namespace std;
typedef long long ll;
const int N=20005;
ll Temp[128],ans;
int n,q,first
,next
,v
,tot,xx,yy,deep
;
void add(int x,int y){v[tot]=y,next[tot]=first[x],first[x]=tot++;}
struct Node{int anc,num;ll eli[62];}fa
[16],G
,jy;
int Gauss(Node a,Node b){
int num=a.num+b.num,flag=1;
for(int i=1;i<=a.num;i++)Temp[i]=a.eli[i];
for(int i=1;i<=b.num;i++)Temp[i+a.num]=b.eli[i];
for(int i=60,j;~i;i--){
for(j=flag;j<=num;j++)if(Temp[j]&(1ll<<i))break;
if(j==num+1)continue;
swap(Temp[flag],Temp[j]);
for(int k=1;k<=num;k++)if(k!=flag&&Temp[k]&(1ll<<i))Temp[k]^=Temp[flag];
flag++;
}
return flag-1;
}
void dfs(int x){
for(int i=1;i<=15;i++){
fa[x][i].anc=fa[fa[x][i-1].anc][i-1].anc;
fa[x][i].num=Gauss(fa[x][i-1],fa[fa[x][i-1].anc][i-1]);
for(int j=1;j<=fa[x][i].num;j++)fa[x][i].eli[j]=Temp[j];
}
for(int i=first[x];~i;i=next[i]){
if(v[i]!=fa[x][0].anc){
fa[v[i]][0].anc=x,deep[v[i]]=deep[x]+1;
fa[v[i]][0].num=Gauss(G[v[i]],G[x]);
for(int j=1;j<=fa[v[i]][0].num;j++)fa[v[i]][0].eli[j]=Temp[j];
dfs(v[i]);
}
}
}
void LCA(int x,int y){
if(deep[x]<deep[y])swap(x,y);
for(int i=15;~i;i--)if(deep[x]-(1<<i)>=deep[y]){
jy.num=Gauss(jy,fa[x][i]);
for(int j=1;j<=jy.num;j++)jy.eli[j]=Temp[j];
x=fa[x][i].anc;
}
if(x==y){jy.num=Gauss(jy,G[x]);for(int j=1;j<=jy.num;j++)jy.eli[j]=Temp[j];return;}
for(int i=15;~i;i--){
if(fa[x][i].anc!=fa[y][i].anc){
jy.num=Gauss(jy,fa[x][i]);
for(int j=1;j<=jy.num;j++)jy.eli[j]=Temp[j];
jy.num=Gauss(jy,fa[y][i]);
for(int j=1;j<=jy.num;j++)jy.eli[j]=Temp[j];
x=fa[x][i].anc,y=fa[y][i].anc;
}
}
jy.num=Gauss(jy,fa[x][0]);
for(int j=1;j<=jy.num;j++)jy.eli[j]=Temp[j];
jy.num=Gauss(jy,fa[y][0]);
for(int j=1;j<=jy.num;j++)jy.eli[j]=Temp[j];
}
int main(){
memset(first,-1,sizeof(first));
scanf("%d%d",&n,&q);
for(int i=1;i<=n;i++)scanf("%lld",&G[i].eli[1]),G[i].num=1;
for(int i=1;i<n;i++)scanf("%d%d",&xx,&yy),add(min(xx,yy),max(xx,yy));
deep[1]=1,dfs(1);
for(int i=1;i<=q;i++){
scanf("%d%d",&xx,&yy);
jy.num=ans=0,LCA(xx,yy);
for(int i=1;i<=jy.num;i++)ans^=jy.eli[i];
printf("%lld\n",ans);
}
}