[Luogu P3292] [BZOJ 4568] [SCOI2016]幸运数字

洛谷传送门

BZOJ传送门

题目描述

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

一些旅行者希望游览 A 国。旅行者计划乘飞机降落在 $x$ 号城市，沿着 $x$ 号城市到 $y$ 号城市之间那条唯一的路径游览，最终从 $y$ 城市起飞离开 A 国。在经过每一座城市时，游览者就会有机会与这座城市的幸运数字拍照，从而将这份幸运保存到自己身上。然而，幸运是不能简单叠加的，这一点游览者也十分清楚。他们迷信着幸运数字是以异或的方式保留在自己身上的。

例如，游览者拍了 $3$ 张照片，幸运值分别是 $5$ ， $7$ ， $11$ ，那么最终保留在自己身上的幸运值就是 $9$ （ $5\ xor\ 7\ xor\ 11$ ）。

有些聪明的游览者发现，只要选择性地进行拍照，便能获得更大的幸运值。例如在上述三个幸运值中，只选择 $5$ 和 $11$ ，可以保留的幸运值为 $14$ 。现在，一些游览者找到了聪明的你，希望你帮他们计算出在他们的行程安排中可以保留的最大幸运值是多少。

输入输出格式

输入格式：

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

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

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

随后 $q$ 行，每行包含两个正整数 $x$ ， $y$ ，表示这名旅行者的旅行计划是从 $x$ 号城市到 $y$ 号城市。 $N\leq 20000,Q\leq 200000,G_i\leq 2^{60}$

输出格式：

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

输入输出样例

输入样例#1：

输出样例#1：

14 
11

解题分析

显然这玩意就是个树上的线性基，直接倍增维护貌似是不满的 $O(Nlog^2(G)log(N))$ ，然而玄学 $AC$ 了？？？

貌似还有点分的 $O(Nlog(G)log(N))$ 的做法，以后填坑吧。

代码如下：

#include <cstdio>
#include <cstring>
#include <cmath>
#include <cctype>
#include <algorithm>
#include <cstdlib>
#define R register
#define IN inline
#define W while
#define gc getchar()
#define MX 20050
#define ll long long
template <class T>
IN void in(T &x)
{
	x = 0; R char c = gc;
	for (; !isdigit(c); c = gc);
	for (;  isdigit(c); c = gc)
	x = (x << 1) + (x << 3) + c - 48;
}
struct Edge {int to, nex;} edge[MX << 1];
ll base[MX][16][63], val[MX], ans[63], res;
int fa[MX][16], head[MX], dep[MX], lg[MX];
int dot, q, cnt;
IN void add(R int from, R int to)
{edge[++cnt] = {to, head[from]}, head[from] = cnt;}
IN void insert(ll *dat, ll val)
{
	for (R int i = 60; ~i; --i)
	{
		if(val & (1ll << i))
		{
			if(!dat[i]) return dat[i] = val, void();
			val ^= dat[i];
		}
	}
}
IN void merge(ll *down, ll *up)
{
	for (R int i = 60; ~i; --i)
	if(up[i]) insert(down, up[i]);
}
void DFS(R int now)
{
	R int tar, par;
	for (R int i = head[now]; i; i = edge[i].nex)
	{
		if(edge[i].to == fa[now][0]) continue;
		tar = edge[i].to;
		fa[tar][0] = now;
		dep[tar] = dep[now] + 1;
		insert(base[tar][0], val[now]);
		for (R int j = 1; j <= lg[dep[tar]]; ++j)
		{
			std::memcpy(base[tar][j], base[tar][j - 1], sizeof(base[tar][j]));
			par = fa[tar][j - 1]; if(!par) break;
			merge(base[tar][j], base[par][j - 1]);
			fa[tar][j] = fa[par][j - 1];
		}
		DFS(tar);
	}
}
IN void getans(R int x, R int y)
{
	if(dep[x] < dep[y]) std::swap(x, y);
	R int del = dep[x] - dep[y];
	insert(ans, val[x]); insert(ans, val[y]);
	for (R int i = 0; i <= lg[del]; ++i)
	if(del & (1 << i)) merge(ans, base[x][i]), x = fa[x][i];
	if(x == y) return;
	for (R int i = lg[dep[x]]; ~i; --i)
	{
		if(fa[x][i] != fa[y][i])
		merge(ans, base[x][i]), merge(ans, base[y][i]), x = fa[x][i], y = fa[y][i];
	}
	merge(ans, base[x][0]);
}
int main(void)
{
	int a, b;
	in(dot), in(q);
	for (R int i = 2; i <= dot; ++i) lg[i] = lg[i >> 1] + 1;
	for (R int i = 1; i <= dot; ++i) in(val[i]);
	for (R int i = 1; i < dot; ++i) in(a), in(b), add(a, b);
	DFS(1);
	W (q--)
	{
		in(a), in(b);
		std::memset(ans, 0, sizeof(ans));
		res = 0; getans(a, b);
		for (R int i = 61; ~i; --i) if((res ^ ans[i]) > res) res ^= ans[i];
		printf("%lld\n", res);
	}
}