某 SCOI 模拟赛 T1 松鼠串门（squirrel）【线性基 贪心】

题意

给定一棵 $n$ 个节点的带权有根树，$q$ 次操作：

• 给某点添加一个儿子；
• 对于某个点到其祖先的一段路径上的点构成的集合，求一个子集使得子集中的点的点权异或和最大。

$n,q\le 5\times 10^5$，权值 $<2^{64}$。

题解

修改操作显然可以提到最前面。

• 对于树上的每个点，保存其到根的路径上所有点的点权构成的线性基，并对每个向量记录提供这个向量对应点的深度；
• 在对线性基从高位到低位进行插入操作时，若某位上的向量非空，应当优先用深度大的替换深度小的，替换后拿被换掉的那个继续尝试往下插入；因为如果深度小的在那一段路径上，深度大的肯定也在那一段路径上；
• 查询时，仍然从高位到低位依次贪心，但是只考虑深度合法的向量。

时间复杂度：$O(nw)$（$w=64$）

代码：

#include<bits/stdc++.h>
using namespace std;
#define ll long long
#define ull unsigned long long
ll getint(){
	ll ans=0,f=1;
	char c=getchar();
	while(c<'0'||c>'9'){
		if(c=='-')f=-1;
		c=getchar();
	}
	while(c>='0'&&c<='9'){
		ans=ans*10-'0'+c;
		c=getchar();
	}
	return ans*f;
}
ull getuint(){
	ull ans=0;
	char c=getchar();
	while(c<'0'||c>'9')c=getchar();
	while(c>='0'&&c<='9'){
		ans=ans*10-'0'+c;
		c=getchar();
	}
	return ans;
}
void putuint(ull x){
	if(x>=10)putuint(x/10);
	putchar(x%10+'0');
}
const int N=1e6+10,A=64;

int n,dep[N],fa[N];
ull a[N];

ull v[N][A];
int d[N][A];
void insert(ull x,int dx,ull *a,int *d){
	for(int i=A-1;i>=0&&x;--i){
		if((x>>i)&1){
			if(a[i]){
				if(dx<d[i])x^=a[i];
				else swap(x,a[i]),swap(dx,d[i]),x^=a[i];
			}
			else{
				a[i]=x;d[i]=dx;
				return;
			}
		}
	}
}
ull query(int dx,ull *a,int *d){
	ull ans=0;
	for(int i=A-1;i>=0;--i)
		if(d[i]>=dx&&(ans^a[i])>ans)ans^=a[i];
	return ans;
}
void ss(){
	for(int i=1;i<=n;i++){
		dep[i]=dep[fa[i]]+1;
		memcpy(v[i],v[fa[i]],sizeof(v[i]));
		memcpy(d[i],d[fa[i]],sizeof(d[i]));
		insert(a[i],dep[i],v[i],d[i]);
	}
}

int qu[N],qd[N],cnt=0;
int main(){
	n=getint();
	int m=getint();
	for(int i=1;i<=n;i++){
		fa[i]=getint();
		a[i]=getuint();
	}
	for(int i=0;i<m;i++){
		int tp=getint();
		if(tp==0){
			qu[cnt]=getint();
			qd[cnt]=getint();
			++cnt;
		}else{
			++n;
			fa[n]=getint();
			a[n]=getuint();
		}
	}
	ss();
	for(int i=0;i<cnt;i++){
		putuint( query(dep[qu[i]]-qd[i],v[qu[i]],d[qu[i]]) );
		putchar('\n');
	}
	return 0;
}