BZOJ4568 [Scoi2016]幸运数字

BZOJ4568 [Scoi2016]幸运数字

---

题目描述

传送门

题目分析

这个题，求树上链的(XOR)最大值，可以不选某些点。

考虑到线性基可以用很快的速度求出(XOR)的最大值，我们首先可以确定使用线性基来搞定这个题的询问。

然后可以考虑用树剖解决，发现一个问题，使用树剖时，由于一边跳一遍计算线性基，还要合并，复杂度比较高，虽然的确可以通过本题，但是明显可以通过随便增大数据就卡死了。倍增貌似拥有更好的常数，但是，空间利用太低，容易(MLE)。

发现每次询问只询问链，那么我们可以使用LCT点分治来解决这个题。

具体方法就是，我们点分治计算每个子树的时候，把询问中lca是这个子树的根的直接求出来然后合并，不是的就把询问传下去，这样复杂度就是(O(Nlog^2N))的，非常优秀可以通过。

是代码呢

#include <bits/stdc++.h>
using namespace std;
const int MAXN=2e5+7;
#define ll long long
const ll inf=1ll<<61;
ll n,m,a[MAXN],p[MAXN],ans[MAXN],size[MAXN],wson[MAXN],vis[MAXN],U[MAXN],V[MAXN],bel[MAXN],tq[MAXN],sum,root;
bool use[MAXN];
inline ll read()
{
    ll x=0,c=1;
    char ch=' ';
    while((ch<'0'||ch>'9')&&ch!='-')ch=getchar();
    while(ch=='-')c*=-1,ch=getchar();
    while(ch>='0'&&ch<='9')x=x*10+ch-'0',ch=getchar();
    return x*c;
}
vector<int> edge[MAXN],q[MAXN];
struct T{
	ll b[65];
	inline void clear() {memset(b,0,sizeof(b));}
	inline void insert(ll x){
		for(int i=60;~i;i--)
			if(x>>i&1){
				if(b[i]) x^=b[i];
				else {b[i]=x;break;}
			}
	}
	inline void merge(const T &x){
		for(int i=60;~i;i--){
			if(x.b[i]) insert(x.b[i]);
		}
	}
	inline ll query(){
		ll ans=0;
		for(int i=60;~i;i--){
			if((ans^b[i])>ans) ans^=b[i];
		}
		return ans;
	}
}base[MAXN];
inline void get_root(int u,int fa)
{
	size[u]=1;wson[u]=0;
	for(int i=0;i<edge[u].size();i++){
		int v=edge[u][i];
		if(v==fa||vis[v]) continue;
		get_root(v,u);
		size[u]+=size[v];
		if(wson[u]<size[v]) wson[u]=size[v];
	}
	if(sum-size[u]>wson[u]) wson[u]=sum-size[u];
	if(wson[u]<wson[root]) root=u;
}

inline void dfs(int u,int fa,int tt){
	base[u]=base[fa]; bel[u]=tt;base[u].insert(a[u]);
	for(int i=0;i<edge[u].size();i++){
		int v=edge[u][i];
		if(v!=fa&&!vis[v]) dfs(v,u,tt);
	}
}

inline void solve(int x){
	if(!q[x].size()) return;
	wson[0]=inf;sum=size[x];root=0;
	get_root(x,0);
	int u=root;
	vis[root]=1;bel[root]=root;
	base[root].clear();base[root].insert(a[root]);
	for(int i=0;i<edge[root].size();i++){
		int v=edge[root][i];
		if(!vis[v]) dfs(v,root,v);
	}
	int tot=q[x].size();
	for(int i=0;i<tot;i++) tq[i]=q[x][i];
	q[x].clear();
	T tmp;ll id;
	for(int i=0;i<tot;i++){
		id=tq[i];
		if(bel[U[id]]==bel[V[id]])
			q[bel[U[id]]].push_back(id);
		else tmp=base[U[id]],tmp.merge(base[V[id]]),ans[id]=tmp.query(),use[id]=1;
	}
	for(int i=0;i<edge[u].size();i++){
		int v=edge[u][i];
		if(vis[v]) continue;
		solve(v);
	}
}
int main()
{
	n=read();m=read();
	for(int i=1;i<=n;i++) a[i]=read();
	for(int i=1;i<n;i++){
		int u=read(),v=read();
		edge[u].push_back(v);
		edge[v].push_back(u);
	}
	for(int i=1;i<=m;i++){
		U[i]=read();V[i]=read();
		if(U[i]==V[i]) ans[i]=a[U[i]];
		else q[1].push_back(i);
	}
	size[1]=n;solve(1);
	for(int i=1;i<=m;i++) printf("%lld
", ans[i]);
}