其他分享
首页 > 其他分享> > luogu P3292 [SCOI2016]幸运数字

luogu P3292 [SCOI2016]幸运数字

作者:互联网

题面传送门
树上问题没有修改,考虑倍增。
然而发现如果暴力合并线性基要\(O(nlognlog^2w+qlognlog^2w)\)铁定过不去。
考虑优化。
可以发现线性基有一个很优秀的性质就是两个有重叠部分的线性基合并是没有影响的。
那么可以用st表的思虑拿来合并了,复杂度降到\(O(nlognlog^2w+qlog^2w)\)几乎不用怎么卡就过了。
代码实现:

#include<cstdio>
#include<cstring>
#define ll long long
using namespace std;
int lg[100039],n,m,k,x,y,z,fa[100039][39],d[100039],lcas,now;
ll w[100039];
struct base{
	ll a[69];int i;
	base(){memset(a,0,sizeof(a));}
	inline void get(ll  x){
		for(i=61;i>=0;i--){
			//if(!x) break;
			if(x&(1ll<<i)){
				if(!a[i]){a[i]=x;break;}
				else x^=a[i];
			}
		}
	}
	inline void merge(base b){for(int i=61;i>=0;i--)if(b.a[i]) get(b.a[i]);}
	inline ll find(){
		ll ans=0;
		for(i=61;i>=0;i--) if((ans^a[i])>ans) ans^=a[i];
		return ans;
	}
}f[20039][20],tmp;
struct yyy{int to,z;};
struct ljb{
	int head,h[100039];yyy f[100039];
	inline void add(int x,int y){f[++head]=(yyy){y,h[x]};h[x]=head;}
}s;
inline void dfs(int x,int last){
	yyy tmp;d[x]=d[last]+1;fa[x][0]=last;f[x][0].get(w[x]);
	for(int i=1;i<=lg[d[x]]+1;i++) fa[x][i]=fa[fa[x][i-1]][i-1],f[x][i]=f[x][i-1],f[x][i].merge(f[fa[x][i-1]][i-1]);
	for(int cur=s.h[x];cur;cur=tmp.z){
		tmp=s.f[cur];
		if(tmp.to^last) dfs(tmp.to,x);
	}
}
inline int kth(int x,int k){
	for(int i=20;i>=0;i--)if(k&(1<<i)) x=fa[x][i];
	return x;
}
inline void swap(int &x,int &y){x^=y^=x^=y;}
inline int lca(int x,int y){
	if(d[x]<d[y]) swap(x,y);
	while(d[x]!=d[y]) x=fa[x][lg[d[x]-d[y]]];
	if(x==y) return x;
	for(int i=20;i>=0;i--)if(fa[x][i]!=fa[y][i]) x=fa[x][i],y=fa[y][i];
	return fa[x][0];
}
int main(){
	freopen("1.in","r",stdin);
	register int i;
	scanf("%d%d",&n,&m);
	for(i=2;i<=n;i++) lg[i]=lg[i/2]+1;
	for(i=1;i<=n;i++)scanf("%lld",&w[i]);
	for(i=1;i<n;i++) scanf("%d%d",&x,&y),s.add(x,y),s.add(y,x);
	dfs(1,0);
	for(i=1;i<=m;i++){
		scanf("%d%d",&x,&y);
		lcas=lca(x,y);
		now=lg[d[y]-d[lcas]+1];tmp=f[y][now];
		tmp.merge(f[kth(y,d[y]-d[lcas]+1-(1<<now))][now]);
		now=lg[d[x]-d[lcas]+1];tmp.merge(f[x][now]);
		tmp.merge(f[kth(x,d[x]-d[lcas]+1-(1<<now))][now]);
		printf("%lld\n",tmp.find());
	}
}

标签:--,int,luogu,ll,ans,P3292,fa,100039,SCOI2016
来源: https://www.cnblogs.com/275307894a/p/14319200.html