【bzoj4568】[Scoi2016]幸运数字 线性基+树链剖分

什么鬼东西呀？

就是把两个线性基合并的过程放在线段树上，倍增可做，反正我写的树剖，听说学校里某人在某次考试的考场上的O(nlogn)算法被O(nlog^2n)的树剖秒成渣了。

异或没打括号，调了好久。

#include<cstdio>
#include<cstring>
#include<cstdlib>
#include<cmath>
#include<algorithm>
#include<iostream>
#define maxn 20010

using namespace std;

struct yts1
{
	long long c[61];
}ans;

void insert(yts1 &a,long long x)
{
	for (int i=60;i>=0;i--)
	  if (x&(1ll<<i))
	  {
	  	if (!a.c[i]) {a.c[i]=x;break;}
	  	else x^=a.c[i];
	  }
}

void merge(yts1 &a,yts1 &b,yts1 &c)
{
	yts1 d;
	for (int i=60;i>=0;i--) d.c[i]=b.c[i];
	for (int i=60;i>=0;i--) if (c.c[i]) insert(d,c.c[i]);
	for (int i=60;i>=0;i--) a.c[i]=d.c[i];
}

long long calc()
{
	long long Ans=0;
	for (int i=60;i>=0;i--) if ((Ans^ans.c[i])>Ans) Ans^=ans.c[i];
	return Ans;
}

struct yts
{
	int l,r;
	yts1 w;
}t[4*maxn];

int head[maxn],to[2*maxn],next[2*maxn],fa[20][maxn],g[maxn],dep[maxn],size[maxn],dd[maxn];
int rank[maxn],seq[maxn],e[maxn],st[maxn];
long long w[maxn];
bool vis[maxn];
int n,m,num,T,tot;

void addedge(int x,int y)
{
	num++;to[num]=y;next[num]=st[x];st[x]=num;
}

void dfs(int x)
{
	e[++tot]=x;size[x]=1;
	for (int p=st[x];p;p=next[p])
	  if (to[p]!=g[x])
	  {
	  	fa[0][to[p]]=g[to[p]]=x;dep[to[p]]=dep[x]+1;
	  	dfs(to[p]);
	  	if (size[to[p]]>size[dd[x]]) dd[x]=to[p];
	  	size[x]+=size[to[p]];
	  }
}

void build(int i,int l,int r)
{
	t[i].l=l;t[i].r=r;
	if (l==r)
	{
		memset(t[i].w.c,0,sizeof(t[i].w.c));
		insert(t[i].w,w[seq[l]]);
		return;
	}
	int mid=(l+r)/2;
	build(i<<1,l,mid);build(i<<1|1,mid+1,r);
	merge(t[i].w,t[i<<1].w,t[i<<1|1].w);
}

void query(int i,int l,int r)
{
	if (l<=t[i].l && t[i].r<=r) {merge(ans,ans,t[i].w);return;}
	int mid=(t[i].l+t[i].r)/2;
	if (l<=mid) query(i<<1,l,r);
	if (mid<r) query(i<<1|1,l,r);
}

int go_up(int x,int d)
{
	if (d==0) return x;
	for (int i=0;i<=16;i++)
	  if (d&(1<<i)) x=fa[i][x];
	return x;
}

int LCA(int x,int y)
{
	if (dep[x]>dep[y]) x=go_up(x,dep[x]-dep[y]);
	else if (dep[y]>dep[x]) y=go_up(y,dep[y]-dep[x]);
	if (x==y) return x;
	for (int i=16;i>=0;i--)
	  if (fa[i][x]!=fa[i][y]) x=fa[i][x],y=fa[i][y];
	return g[x];
}

int main()
{
	scanf("%d%d",&n,&T);
	for (int i=1;i<=n;i++) scanf("%lld",&w[i]);
	for (int i=1;i<n;i++)
	{
		int x,y;
		scanf("%d%d",&x,&y);
		addedge(x,y);addedge(y,x);
	}
	dfs(1);
	int qwer=0;
	for (int i=1;i<=n;i++)
	  if (!vis[e[i]])
	  {
	  	int k=e[i];
	  	while (k)
	  	{
	  		seq[++qwer]=k;vis[k]=1;head[k]=e[i];k=dd[k];
	  	}
	  }
	for (int i=1;i<=n;i++) rank[seq[i]]=i;
	for (int j=1;j<=17;j++)
	  for (int i=1;i<=n;i++)
	    fa[j][i]=fa[j-1][fa[j-1][i]];
	build(1,1,n);
	while (T--)
	{
		int x,y;
		scanf("%d%d",&x,&y);
		if (dep[x]<dep[y]) swap(x,y);
		memset(ans.c,0,sizeof(ans.c));
		int z=LCA(x,y);
		while (dep[head[x]]>dep[z])
		{
			query(1,rank[head[x]],rank[x]);x=g[head[x]];
		}
		query(1,rank[z],rank[x]);
		if (y!=z)
		{
			int w=go_up(y,dep[y]-dep[z]-1);
			while (dep[head[y]]>dep[w])
			{
				query(1,rank[head[y]],rank[y]);y=g[head[y]];
			}
			query(1,rank[w],rank[y]);
		}
		printf("%lld\n",calc());
	}
	return 0;
}