bzoj 4568: [Scoi2016]幸运数字 倍增维护线性基

题意

给出一棵树，每个节点都有一个权值，每次询问x到y的路径上任意选若干个节点的最大异或值。
$n<=20000,m<=200000,a[i]<=2^{60}$

分析

一看到题目就知道是用线性基来维护，想到了用倍增来维护线性基，但觉得复杂度太大了，应该会超时，于是就不敢打。去网上搜了题解发现其实就是这样做就好了。
也就是每次暴力合并线性基即可。

代码

#include<iostream>
#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<algorithm>
#include<cmath>
#define N 20005
#define ll long long
using namespace std;

int n,m,last[N],cnt,sz,dep[N],fa[N][20],logn,num[N][20];
ll bas[N*20][65],a[N],bin[65];
struct edge{int to,next;}e[N*2];

int read()
{
    int x=0,f=1;char ch=getchar();
    while (ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=getchar();}
    while (ch>='0'&&ch<='9'){x=x*10+ch-'0';ch=getchar();}
    return x*f;
}

void addedge(int u,int v)
{
    e[++cnt].to=v;e[cnt].next=last[u];last[u]=cnt;
    e[++cnt].to=u;e[cnt].next=last[v];last[v]=cnt;
}

void merge(int x,int y)
{
    for (int i=60;i>=0;i--)
        if (bas[y][i])
            if (!bas[x][i]) bas[x][i]=bas[y][i];
            else
            {
                ll w=bas[y][i];
                for (int j=i;j>=0;j--)
                    if (w&bin[j])
                        if (!bas[x][j])
                        {
                            bas[x][j]=w;
                            break;
                        }
                        else w^=bas[x][j];
            }
}

void ins(ll x,int y)
{
    for (int i=60;i>=0;i--)
        if (x&bin[i])
            if (!bas[y][i])
            {
                bas[y][i]=x;
                break;
            }
            else x^=bas[y][i];
}

void dfs(int x)
{
    dep[x]=dep[fa[x][0]]+1;
    for (int i=1;i<=logn;i++)
    {
        fa[x][i]=fa[fa[x][i-1]][i-1];
        if (fa[x][i])
        {
            num[x][i]=++sz;
            merge(sz,num[x][i-1]);
            merge(sz,num[fa[x][i-1]][i-1]);
        }
    }
    for (int i=last[x];i;i=e[i].next)
    {
        if (e[i].to==fa[x][0]) continue;
        fa[e[i].to][0]=x;
        num[e[i].to][0]=++sz;
        ins(a[e[i].to],sz);
        dfs(e[i].to);
    }
}

void getlca(int x,int y)
{
    if (dep[x]<dep[y]) swap(x,y);
    for (int i=logn;i>=0;i--)
        if (dep[fa[x][i]]>=dep[y])
        {
            merge(0,num[x][i]);
            x=fa[x][i];
        }
    if (x==y)
    {
        ins(a[x],0);
        return;
    }
    for (int i=logn;i>=0;i--)
        if (fa[x][i]!=fa[y][i])
        {
            merge(0,num[x][i]);
            merge(0,num[y][i]);
            x=fa[x][i];y=fa[y][i];
        }
    merge(0,num[x][0]);
    merge(0,num[y][0]);
    ins(a[fa[x][0]],0);
}

ll getans()
{
    ll ans=0;
    for (int i=60;i>=0;i--)
        if ((ans^bas[0][i])>ans) ans^=bas[0][i];
    return ans;
}

void debug(int x)
{
    for (int i=0;i<=6;i++)
        printf("%lld ",bas[x][i]);
    cout<<endl;
}

int main()
{
    bin[0]=1;
    for (int i=1;i<=60;i++)
        bin[i]=(ll)bin[i-1]*2;
    n=read();m=read();
    for (int i=1;i<=n;i++)
        scanf("%lld",&a[i]);
    for (int i=1;i<n;i++)
    {
        int x=read(),y=read();
        addedge(x,y);
    }
    logn=log(n)/log(2);
    dfs(1);
    for (int i=1;i<=m;i++)
    {
        int x=read(),y=read();
        for (int j=0;j<=60;j++)
            bas[0][j]=0;
        getlca(x,y);
        printf("%lld\n",getans());
    }
    return 0;
}