BZOJ 4568 幸运数字（在线倍增法+线性基）

Description

$A$国共有$n$座城市，这些城市由$n-1$条道路相连，使得任意两座城市可以互达，且路径唯一。每座城市都有一个幸运数字，以纪念碑的形式矗立在这座城市的正中心，作为城市的象征。一些旅行者希望游览$A$国。旅行者计划乘飞机降落在$x$号城市，沿着$x$号城市到$y$号城市之间那条唯一的路径游览，最终从$y$ 城市起飞离开$A$国。在经过每一座城市时，游览者就会有机会与这座城市的幸运数字拍照，从而将这份幸运保存到自己身上。然而，幸运是不能简单叠加的，这一点游览者也十分清楚。他们迷信着幸运数字是以异或的方式保留在自己身上的。例如，游览者拍了 $3$ 张照片，幸运值分别是 $5，7，11$，那么最终保留在自己身上的幸运值就是$9（5 xor 7 xor 11）$。有些聪明的游览者发现，只要选择性地进行拍照，便能获得更大的幸运值。例如在上述三个幸运值中，只选择$5$和$11$，可以保留的幸运值为$14$ 。现在，一些游览者找到了聪明的你，希望你帮他们计算出在他们的行程安排中可以保留的最大幸运值是多少。

Input

第一行包含$2$个正整数$n ，q$，分别表示城市的数量和旅行者数量。第二行包含$n$个非负整数，其中第$i$个整

数 $G_i$ 表示$i$号城市的幸运值。随后$n-1$行，每行包含两个正整数$x ，y$，表示$x$号城市和$y$号城市之间有一

条道路相连。随后$q$行，每行包含两个正整数$x ，y$，表示这名旅行者的旅行计划是从$x$号城市到$y$号城市。

$N\le 20000,Q\le 200000,G_i\le 2^{60}$

Output

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

Sample Input

4 2

11 5 7 9

1 2

1 3

1 4

2 3

1 4

Sample Output

14

11

Solution

一条路径上的最大幸运值在求出这条路径上边权集里的线性基后从高位到低位贪心取即可，问题在于如果求路径上的线性基，考虑倍增，求出从$i$点到其第$2^j$层祖先路径上边权集的线性基，线性基合并复杂度$O(60^2)$，预处理时间复杂度$O(3600nlogn)$，单次查询时间复杂度$O(3600logn)$

Code

#include<cstdio>
#include<iostream>
#include<cstring>
#include<algorithm>
#include<cmath>
#include<vector>
#include<queue>
#include<map>
#include<set>
#include<ctime>
using namespace std;
typedef long long ll;
const int maxn=20005;
void Add(ll *b,ll a)
{
    for(int i=60;i>=0;i--)
        if((a>>i)&1)
        {
            if(!b[i])
            {
                b[i]=a;
                return ;
            }
            else a^=b[i];
        }
}
void Unite(ll *x,ll *y)
{
    for(int i=0;i<=60;i++)
        if(y[i])Add(x,y[i]);
}
int n,q,fa[maxn][16],dep[maxn];
ll base[maxn][16][61],ans[61];
vector<int>g[maxn];
void dfs(int u,int f)
{
    dep[u]=dep[f]+1;
    fa[u][0]=f;
    for(int i=0;i<g[u].size();i++)
    {
        int v=g[u][i];
        if(v!=f)dfs(v,u);
    }
}
void LCA(int u,int v)
{
    if(dep[u]<dep[v])swap(u,v);
    for(int i=15;i>=0;i--)
        if(dep[fa[u][i]]>=dep[v])
        {
            Unite(ans,base[u][i]);
            u=fa[u][i];
        }
    if(u==v)
    {
        Unite(ans,base[u][0]);
        return ;
    }
    for(int i=15;i>=0;i--)
        if(fa[u][i]!=fa[v][i])
        {
            Unite(ans,base[u][i]);u=fa[u][i];
            Unite(ans,base[v][i]);v=fa[v][i];
        }
    Unite(ans,base[u][0]);u=fa[u][0];
    Unite(ans,base[v][0]);v=fa[v][0];
    Unite(ans,base[u][0]);
}
int main()
{
    scanf("%d%d",&n,&q);
    for(int i=1;i<=n;i++)
    {
        ll a;
        scanf("%lld",&a);
        Add(base[i][0],a);
    }
    for(int i=1;i<n;i++)
    {
        int u,v;
        scanf("%d%d",&u,&v);
        g[u].push_back(v),g[v].push_back(u);
    }
    dfs(1,0);
    for(int j=1;j<=15;j++)
        for(int i=1;i<=n;i++)
        {
            fa[i][j]=fa[fa[i][j-1]][j-1];
            for(int k=0;k<=60;k++)base[i][j][k]=base[i][j-1][k];
            Unite(base[i][j],base[fa[i][j-1]][j-1]);
        }
    while(q--)
    {
        int u,v;
        scanf("%d%d",&u,&v);
        memset(ans,0,sizeof(ans));
        LCA(u,v);
        ll res=0;
        for(int i=60;i>=0;i--)res=max(res,res^ans[i]);
        printf("%lld\n",res);
    }
    return 0;
}