bzoj4568: [Scoi2016]幸运数字 线性基 倍增

洛谷T了，mmp

还是bzoj时限良心，虽然三天两头爆炸

算是简单的题吧，构造出每个点logn个的线性基表

感觉线性基最nb的就是只有log个，所以可以做的很暴力

合并就是拆开再并

#include <bits/stdc++.h>
using namespace std;
struct Bas
{
    long long a[61];
    Bas()
    {
        for(int i=0;i<=60;i++)
            a[i]=0;
    }
    void add(long long x)
    {
        for(int i=60;i>=0;i--)
        if((x>>i-1)&1)
        if(a[i]) x^=a[i];
        else
        {
            a[i]=x;
            return;
        }
    }
    long long que()
    {
        long long ret=0;
        for(int i=60;i>=0;i--)
        if((ret^a[i])>ret)
            ret^=a[i];
        return ret;
    }
} ba[20001][16];
int n,m,E,p,q;
long long w[20001];
int fa[20001][16],dep[20001];
int fir[20001],nex[40001],to[40001];
Bas merge(Bas a,Bas b)
{
    for(int i=0;i<=60;i++)
    if(a.a[i]) b.add(a.a[i]); 
    return b;
}
void add(int x,int y)
{
    to[++E]=y;nex[E]=fir[x];fir[x]=E;
}
void build(int now,int ft)
{
    fa[now][0]=ft;dep[now]=dep[ft]+1;
    for(int i=0;fa[now][i];i++)
        fa[now][i+1]=fa[fa[now][i]][i],
        ba[now][i+1]=merge(ba[now][i],ba[fa[now][i]][i]);
    for(int i=fir[now];i;i=nex[i])
    if(to[i]!=ft)
        build(to[i],now);
}
int main()
{
    scanf("%d%d",&n,&m);
    for(int i=1;i<=n;i++)
    {
        scanf("%lld",&w[i]);
        ba[i][0].add(w[i]);
    }
    for(int i=1;i<n;i++)
        scanf("%d%d",&p,&q),add(p,q),add(q,p);
    build(1,0);
    for(int i=1;i<=m;i++)
    {
        scanf("%d%d",&p,&q);
        Bas ret=Bas();
        if(dep[p]<dep[q]) swap(p,q);
        for(int i=15;dep[p]>dep[q];i--)
        if(dep[p]-dep[q]>>i)
        {
            ret=merge(ret,ba[p][i]);
            p=fa[p][i];
        }
        for(int i=15;i>=0;i--)
        if(fa[p][i]!=fa[q][i])
            ret=merge(merge(ba[p][i],ba[q][i]),ret),
            p=fa[p][i],q=fa[q][i];
        if(p!=q)
        {
            ret=merge(merge(ba[p][0],ba[q][0]),ret);
            p=fa[p][0];q=fa[q][0];
        }
        ret=merge(ba[p][0],ret);
        printf("%lld
",ret.que());
    }
    return 0;
}