DTOJ #2653. 「SCOI2016」幸运数字（luckly）

题目

A 国共有 $n$ 座城市，这些城市由 $n-1$ 条道路相连，使得任意两座城市可以互达，且路径唯一。每座城市都有一个幸运数字，以纪念碑的形式矗立在这座城市的正中心，作为城市的象征。

一些旅行者希望游览 A 国。旅行者计划乘飞机降落在 $x$ 号城市，沿着 $x$ 号城市到 $y$ 号城市之间那条唯一的路径游览，最终从 $y$ 城市起飞离开 A 国。在经过每一座城市时，游览者就会有机会与这座城市的幸运数字拍照，从而将这份幸运保存到自己身上。

然而，幸运是不能简单叠加的，这一点游览者也十分清楚。他们迷信着幸运数字是以异或的方式保留在自己身上的。例如，游览者拍了 $3$ 张照片，幸运值分别是 $5,7,11$，那么最终保留在自己身上的幸运值就是 $9(5xor7xor11)$ 。有些聪明的游览者发现，只要选择性地进行拍照，便能获得更大的幸运值。例如在上述三个幸运值中，只选择 $5$ 和 $11$，可以保留的幸运值为 $14$。

现在，一些游览者找到了聪明的你，希望你帮他们计算出在他们的行程安排中可以保留的最大幸运值是多少。

输入格式

第一行包含 $2$ 个正整数 $n,q$，分别表示城市的数量和旅行者数量。

第二行包含 $n$ 个非负整数，其中第 $i$ 个整数 $G_i$ 表示 $i$ 号城市的幸运值。

随后 $n-1$ 行，每行包含两个正整数 $x,y$，表示 $x$ 号城市和 $y$ 号城市之间有一条道路相连。

随后 $q$ 行，每行包含两个正整数 $x,y$，表示这名旅行者的旅行计划是从 $x$ 号城市到 $y$ 号城市。

输出格式

输出需要包含 $q$ 行，每行包含 $1$ 个非负整数，表示这名旅行者可以保留的最大幸运值。

数据范围

$n\le 2\times 10^4,q\le 2\times 10^5,G_i\le 2^{60}$。

题解

考虑到用线性基来做（最大异或和）。
然后想怎么把 $x-y$ 的简单路径的线性基求出来，因为线性基删除并不好做，所以我们考虑由若干个线性基组合而成。
所以自然而然的想到用倍增：
$a[i][j]$ 表示 $i$ 号点包括自己往上一共 $2^j$ 个点的线性基。
然后用直接倍增求线性基即可。

代码

#include<bits/stdc++.h>
#define ll long long
using namespace std;
inline ll read(){
	char c;ll k=0,f=1;
	while(c<'0'||c>'9'){if(c=='-')f=-1;c=getchar();}
	while(c>='0'&&c<='9')k=k*10+c-'0',c=getchar();
	return k*f;
}
const int N=2e4+5;
ll g[N];
struct A{
	ll s[63];
	A(){memset(s,0,sizeof(s));}
}a[N][20];
int n,q,nxt[N*2],f[N][20],fst[N],d[N],ver[N*2],tot;
inline void add(int x,int y){ver[++tot]=y;nxt[tot]=fst[x];fst[x]=tot;ver[++tot]=x;nxt[tot]=fst[y];fst[y]=tot;}
void AD(A &x,ll k){
	for(int i=62;i>=0;--i){
		if(!((1LL<<i)&k))continue;
		if(!x.s[i]){x.s[i]=k;break;}
		else k^=x.s[i];
	}
}
void merge(A &x,A y){
	for(int j=0;j<62;++j)if(y.s[j])AD(x,y.s[j]);
}
void solve(int x){
	for(int i=1;i<20;++i){
		if(!f[f[x][i-1]][i-1])break;
		f[x][i]=f[f[x][i-1]][i-1];
	}
	AD(a[x][0],g[x]);
	for(int i=1;i<20;++i){
		a[x][i]=a[x][i-1];
		merge(a[x][i],a[f[x][i-1]][i-1]);
	}
}
void dfs(int x,int fa){
	f[x][0]=fa;solve(x);d[x]=d[fa]+1;
	for(int i=fst[x];i;i=nxt[i]){
		if(ver[i]==fa)continue;
		dfs(ver[i],x);
	}
}
int main(){
	n=read();q=read();
	for(int i=1;i<=n;++i)g[i]=read();
	for(int i=1;i<n;++i)add(read(),read());
	dfs(1,0);
	while(q--){
		int x=read(),y=read();if(d[x]<d[y])swap(x,y);
		A k;
		int cha=d[x]-d[y],lca,kx=x,ky=y;
		for(int i=19;i>=0;--i){
			if(cha&(1<<i))kx=f[kx][i];
		}
		for(int i=19;i>=0;--i)if(f[kx][i]!=f[ky][i])kx=f[kx][i],ky=f[ky][i];
		lca=(kx==ky)?kx:f[kx][0];
		int c_x=d[x]-d[lca]+1,c_y=d[y]-d[lca]+1;
		for(int i=19;i>=0;--i){
			if(c_x&(1<<i))merge(k,a[x][i]),x=f[x][i];
			if(c_y&(1<<i))merge(k,a[y][i]),y=f[y][i];
		}
		ll ans=0;
		for(int i=62;i>=0;--i)ans=max(ans,ans^k.s[i]);
		printf("%lld\n",ans);
	}
	return 0;
}