BZOJ4568 [Scoi2016]幸运数字 【点分治 + 线性基】

题目链接

BZOJ4568

题解

选任意个数异或和最大，使用线性基

线性基插入(O(logn))，合并(O(log^2n))

我们要求树上两点间异或和最大值，由于合并是(O(log^2n))的，我们尽量只合并一次
那就采用点分治
每次求出到分治重心的线性基，将过分治重心的询问的两个线性基合并即可
复杂度(O(60^2q + 60nlogn))

#include<algorithm>
#include<iostream>
#include<cstring>
#include<cstdio>
#include<cmath>
#include<map>
#define Redge(u) for (int k = h[u],to; k; k = ed[k].nxt)
#define REP(i,n) for (int i = 1; i <= (n); i++)
#define mp(a,b) make_pair<int,int>(a,b)
#define cls(s) memset(s,0,sizeof(s))
#define cp pair<int,int>
#define LL long long int
using namespace std;
const int maxn = 20005,maxm = 200005,INF = 1000000000;
inline LL read(){
	LL out = 0,flag = 1; char c = getchar();
	while (c < 48 || c > 57){if (c == '-') flag = -1; c = getchar();}
	while (c >= 48 && c <= 57){out = (out << 3) + (out << 1) + c - 48; c = getchar();}
	return out * flag;
}
struct Bit{
	LL A[60];
	void init(){for (int i = 59; ~i; i--) A[i] = 0;}
	void copy(LL* B){for (int i = 59; ~i; i--) A[i] = B[i];}
	void ins(LL x){
		for (int i = 59; ~i; i--)
			if ((x >> i) & 1){
				if (A[i]) x ^= A[i];
				else {A[i] = x; break;}
			}
	}
	LL ask(){
		LL re = 0;
		for (int i = 59; ~i; i--) if ((re ^ A[i]) > re) re ^= A[i];
		return re;
	}
}B[maxn],T;
int h[maxn],ne = 1;
struct EDGE{int to,nxt;}ed[maxn << 1];
inline void build(int u,int v){
	ed[++ne] = (EDGE){v,h[u]}; h[u] = ne;
	ed[++ne] = (EDGE){u,h[v]}; h[v] = ne;
}
int g[maxn],nxt[30 * maxm],tq[30 * maxn],cnt;
int n,q,x[maxm],y[maxm],vis[maxn],pos[maxn],Vis[maxn],now;
LL G[maxn],ans[maxm];
void Add(int u,int v){
	nxt[++cnt] = g[u]; tq[cnt] = v; g[u] = cnt;
}
int F[maxn],siz[maxn],fa[maxn],N,rt;
void getrt(int u){
	F[u] = 0; siz[u] = 1;
	Redge(u) if (!vis[to = ed[k].to] && to != fa[u]){
		fa[to] = u; getrt(to);
		siz[u] += siz[to];
		F[u] = max(F[u],siz[to]);
	}
	F[u] = max(F[u],N - siz[u]);
	if (F[u] < F[rt]) rt = u;
}
void dfs(int u,int R){
	pos[u] = R; siz[u] = 1; Vis[u] = now;
	B[u].copy(B[fa[u]].A); B[u].ins(G[u]);
	Redge(u) if (!vis[to = ed[k].to] && to != fa[u]){
		fa[to] = u; dfs(to,R);
		siz[u] += siz[to];
	}
}
void solve(int u){
	F[rt = 0] = INF; N = siz[u]; getrt(u);
	//printf("u%d  rt%d
",u,rt);
	pos[rt] = rt; vis[rt] = true; siz[u] = 1; Vis[rt] = ++now;
	B[rt].init(); B[rt].ins(G[rt]);
	Redge(rt) if (!vis[to = ed[k].to]){
		fa[to] = rt; dfs(to,to);
		siz[u] += siz[to];
	}
	for (int k = g[u],i,a,b; k; k = nxt[k]){
		i = tq[k]; a = x[i]; b = y[i];
		if (Vis[a] != now || Vis[b] != now) continue;
		if (pos[a] == pos[b]) Add(pos[a],i);
		else {
			T.copy(B[a].A);
			for (int j = 59; ~j; j--)
				if (B[b].A[j]) T.ins(B[b].A[j]);
			ans[i] = T.ask();
		}
	}
	Redge(rt) if (!vis[to = ed[k].to]){
		solve(to);
	}
}
int main(){
	n = read(); q = read();
	for (int i = 1; i <= n; i++) G[i] = read();
	for (int i = 1; i < n; i++) build(read(),read());
	for (int i = 1; i <= q; i++){
		x[i] = read(); y[i] = read();
		if (x[i] == y[i]) ans[i] = G[x[i]];
		else Add(1,i);
	}
	siz[1] = n; solve(1);
	for (int i = 1; i <= q; i++)
		printf("%lld
",ans[i]);
	return 0;
}