[SCOI2016]幸运数字

题目描述

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 个整数 Gi 表示 i 号城市的幸运值。随后 n-1 行，每行包含两个正整数 x ，y，表示 x 号城市和 y 号城市之间有一条道路相连。随后 q 行，每行包含两个正整数 x ，y，表示这名旅行者的旅行计划是从 x 号城市到 y 号城市。N<=20000,Q<=200000,Gi<=2^60

输出格式：

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

输入输出样例

输入样例#1：

4 2
11 5 7 9
1 2
1 3
1 4
2 3
1 4

输出样例#1：

14 
11

我们每次找到重心，处理与重心相关的路径。

处理时我们将重心到每个节点这一段的线性基存起来，然后找到所有经过重心的路径

在遍历以重心G为根的一个子树过程中，找到与这棵子树中节点u有关的询问(u,v)，判断是否在之前遍历过的以重心为根的其他子树中出现过，如果出现过，我们可以将G->u和G->v的线性基合并

找到合并后的线性基中的最大值就好了。对于合并，直接暴力拆开一个线性基，一个一个插入到另一个线性基中

处理完这棵子树之后，再遍历一遍，打上访问过的标记。这样保证所有询问都只计算过一次。

注意，处理完这个重心之后，清空标记时，不要memset，直接再遍历一遍就好了，快得多。

值得注意的是，要单独考虑与G有关的询问如(G,v)

#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<cmath>
using namespace std;
typedef long long lol;
struct Node
{
  int next,to,dis;
}edge[40001],edgeq[400001];
int num,numq,head[20001],headq[200001];
int size[20001],maxsize[20001],minsize,root,n,q;
bool vis[20001],judge[20001];
lol ans[200001],pw[64],val[20001];
struct BASE
{
  lol a[64];
  void clear()
  {
    for (int i=0;i<=62;i++)
      a[i]=0;
  }
  void insert(lol x)
  {
    for (int i=62;i>=0;i--)
      {
    if (x&pw[i])
      {
        if (a[i]==0)
          {a[i]=x;break;}
        else x^=a[i];
      }
      }
  }
  lol getmax()
  {
    lol x=0;
    for (int i=62;i>=0;i--)
      {
    if ((x^a[i])>x)
      x^=a[i];
      }
    return x;
  } 
  void merge(BASE b)
  {
    for (int i=0;i<=62;i++)
      insert(b.a[i]);
  }
  void copy(BASE b)
  {
    for (int i=0;i<=62;i++)
      a[i]=b.a[i];
  }
}base[20002];
void add(int u,int v)
{
  num++;
  edge[num].next=head[u];
  head[u]=num;
  edge[num].to=v;
}
void addq(int u,int v,int d)
{
  numq++;
  edgeq[numq].next=headq[u];
  headq[u]=numq;
  edgeq[numq].to=v;
  edgeq[numq].dis=d;
}
void get_size(int x,int fa)
{
  int i;
  size[x]=1;
  maxsize[x]=0;
  for (i=head[x];i;i=edge[i].next)
    {
      int v=edge[i].to;
      if (vis[v]==0&&v!=fa)
    {
      get_size(v,x);
      size[x]+=size[v];
      maxsize[x]=max(maxsize[x],size[v]);
    }
    }
}
void get_root(int r,int x,int fa)
{
  int i;
  maxsize[x]=max(maxsize[x],size[r]-size[x]);
  if (maxsize[x]<minsize)
    {
      root=x;
      minsize=maxsize[x];
    }
  for (i=head[x];i;i=edge[i].next)
    {
      int v=edge[i].to;
      if (vis[v]==0&&v!=fa)
    {
      get_root(r,v,x);
    }
    }
}
void get_ans(int r,int x,int fa)
{int i;
  base[x].copy(base[fa]);
  base[x].insert(val[x]);
  for (i=headq[x];i;i=edgeq[i].next)
    {
      int v=edgeq[i].to;
      if (judge[v])
    {
      BASE tmp;
      tmp.copy(base[x]);
      tmp.merge(base[v]);
      ans[edgeq[i].dis]=tmp.getmax();
    }
      else if (v==r)
    {
      ans[edgeq[i].dis]=base[x].getmax();
    }
    }
  for (i=head[x];i;i=edge[i].next)
    {
      int v=edge[i].to;
      if (vis[v]==0&&v!=fa)
    {
      get_ans(r,v,x);
    }
    }
}
void get_update(int x,int fa)
{
  judge[x]=!judge[x];
  for (int i=head[x];i;i=edge[i].next)
    {
      int v=edge[i].to;
      if (vis[v]==0&&v!=fa)
    {
      get_update(v,x);
    }
    }
}
void doit(int x)
{
  minsize=2e9;
  get_size(x,0);
  get_root(x,x,0);
  vis[root]=1;
  base[root].clear();
  base[root].insert(val[root]);
  for (int i=head[root];i;i=edge[i].next)
    {
      int v=edge[i].to;
      if (vis[v]==0)
    {
      get_ans(root,v,root);
      get_update(v,root);
    }
    }
    for (int i=head[root];i;i=edge[i].next)
    {
      int v=edge[i].to;
      if (vis[v]==0)
    {
      get_update(v,root);
    }
    }
 for (int i=head[root];i;i=edge[i].next)
    {
      int v=edge[i].to;
      if (vis[v]==0)
    {
      doit(v);
    }
    }
}
int main()
{int i,u,v;
  cin>>n>>q;
  pw[0]=1;
  for (i=1;i<=62;i++)
    pw[i]=pw[i-1]*2;
  for (i=1;i<=n;i++)
    {
      scanf("%lld",&val[i]);
    }
  for (i=1;i<=n-1;i++)
    {
      scanf("%d%d",&u,&v);
      add(u,v);add(v,u);
    }
  for (i=1;i<=q;i++)
    {
      scanf("%d%d",&u,&v);
      if (u!=v)
    addq(u,v,i),addq(v,u,i);
      else ans[i]=val[u];
    }
  doit(1);
  for (i=1;i<=q;i++)
    {
      printf("%lld
",ans[i]);
    }
}