YBTOJ&洛谷P3292：幸运数字（线性基、点分治/倍增）

解析

虽然使用三个log的倍增算法艹过去了
但是我们还是来聊聊正解吧

考虑点分治
对于当前的根，dfs求出联通块内每个点到当前根的线性基
一条路径的答案应该在路径出现上第一个成为根的点时统计到
具体来说，就是路径的两端点在同一个solve函数的不同子树内
通过奇怪的时间戳标记可以实现
对于询问，每个结点开个vector暴力遍历即可
由于每个结点最多被dfslog遍，所以复杂度是对的
复杂度瓶颈在dfs求线性基上，时间复杂度$nlo{g}^{2}n$
然后就ok啦
但是懒得再写一遍了，所以还是贴的3log代码

代码

#include<bits/stdc++.h>
using namespace std;
#define ll long long
#define inf (n+1)
//#define debug(...) fprintf(stderr,__VA_ARGS__)
const int N=2e4+100;
const int M=2e5+10500;
const double eps=1e-5;
inline ll read(){ll x(0),f(1);char c=getchar();
	while(!isdigit(c)){if(c=='-')f=-1;c=getchar();}
	while(isdigit(c)){x=(x<<1)+(x<<3)+c-'0';c=getchar();}
	return x*f;
}

int n,m;

struct node{
	int to,nxt;
}p[N<<1];
int fi[N],cnt;
inline void addline(int x,int y){
	p[++cnt]=(node){y,fi[x]};fi[x]=cnt;
	return;
}

ll mi[65];

struct basic{
	ll v[62];
	basic(){memset(v,0,sizeof(v));}
	void upd(ll x){
		for(int i=60;i>=0;i--){
			if(x&mi[i]){
				if(!v[i]){
					v[i]=x;break;
				}
				else x^=v[i];
			}
		}
	}
}b[N][16];
inline void merge(basic &v1,basic v2){
	for(int i=0;i<=60;i++){
		if(v2.v[i]) v1.upd(v2.v[i]);
	}
	return;
}

ll val[N];
int pl[N][16],dep[N];
void dfs(int x,int f){
	b[x][0].upd(val[f]);pl[x][0]=f;
	for(int k=1;pl[x][k-1];k++){
		pl[x][k]=pl[pl[x][k-1]][k-1];
		b[x][k]=b[x][k-1];
		merge(b[x][k],b[pl[x][k-1]][k-1]);
	}
	dep[x]=dep[f]+1;
	for(int i=fi[x];~i;i=p[i].nxt){
		int to=p[i].to;
		if(to==f) continue;
		dfs(to,x);
	}
	return;
}
basic res;
ll query(int x,int y){
	memset(res.v,0,sizeof(res.v));
	res.upd(val[x]);res.upd(val[y]);
	if(dep[x]<dep[y]) swap(x,y);
	for(int k=15;k>=0;k--){
		if(dep[x]-mi[k]<dep[y]) continue;
		merge(res,b[x][k]);x=pl[x][k];
	}
	if(x!=y){
		for(int k=15;k>=0;k--){
			if(pl[x][k]==pl[y][k]) continue;
			merge(res,b[x][k]);merge(res,b[y][k]);
			x=pl[x][k];y=pl[y][k];
		}
		res.upd(val[pl[x][0]]);
	}
	ll ans(0);
	for(int i=60;i>=0;i--){
		//if(res.v[i]) printf("i=%d v=%lld\n",i,res.v[i]);
		if(!res.v[i]||(ans&mi[i])) continue;
		ans^=res.v[i];
		//printf("  ans=%lld\n",ans);
	}
	//putchar('\n');
	return ans;
}
int main(){
#ifndef ONLINE_JUDGE
	//freopen("a.in","r",stdin);
	//freopen("a.out","w",stdout);
#endif
	memset(fi,-1,sizeof(fi));cnt=-1;
	mi[0]=1;
	for(int i=1;i<=60;i++) mi[i]=(mi[i-1]<<1);
	n=read();m=read();
	for(int i=1;i<=n;i++) val[i]=read();
	for(int i=1;i<n;i++){
		int x=read(),y=read();
		addline(x,y);addline(y,x);
	}
	dfs(1,0);
	for(int i=1;i<=m;i++){
		int x=read(),y=read();
		printf("%lld\n",query(x,y));
	}
	return 0;
}
/*
*/