bzoj 4568 [Scoi2016]幸运数字 倍增+线性基

首先我们可以先预处理倍增链上的线性基，需要合并两个线性基时暴力合并就行了。bzoj时限60s 可以卡过去。

#include<cstdio>
#include<iostream>
#include<cstring>
#include<vector>
#include<algorithm>
#define ll long long
using namespace std;
const int maxn=22010;
vector<int> f[maxn];
int n,m,x,y,cnt,fa[maxn][19],deep[maxn];
ll val[maxn],F[maxn][19][65],ans[maxn];
void merge1(ll *a,ll b)
{
    for(int i=60;~i;i--) if(b>>i){
        if(a[i]) b^=a[i]; 
        else{
            a[i]=b;return ; 
        }  
    } 
}
void merge2(ll *a,ll *b)
{
    for(int i=60;~i;i--) if(b[i]) merge1(a,b[i]);
}
void dfs1(int u,int fath,int dep)
{
    fa[u][0]=fath;deep[u]=dep+1;
    for(int i=0;i<f[u].size();i++)
    {
        int v=f[u][i];
        if(v==fath) continue; 
        dfs1(v,u,dep+1);
    }
}
void pre()
{
    for(int j=1;j<=16;j++)
        for(int i=1;i<=n;i++){
            fa[i][j]=fa[fa[i][j-1]][j-1];
            memcpy(F[i][j],F[i][j-1],sizeof(F[i][j-1]));
            merge2(F[i][j],F[fa[i][j-1]][j-1]);
        }
}

void lca(int x,int y)
{
    if(deep[x]<deep[y]) swap(x,y);
    for(int i=16;~i;i--)
        if(deep[fa[x][i]]>=deep[y]){
            merge2(ans,F[x][i]);
            x=fa[x][i];
        }
    if(x==y){
        merge2(ans,F[x][0]);
        return ;
    }
    for(int i=16;~i;i--)
        if(fa[x][i]!=fa[y][i]){
            merge2(ans,F[x][i]);merge2(ans,F[y][i]);
            x=fa[x][i];y=fa[y][i];
        }
    merge2(ans,F[x][0]),merge2(ans,F[y][0]);
    merge2(ans,F[fa[x][0]][0]);
}

int main()
{
    scanf("%d%d",&n,&m);
    for(int i=1;i<=n;i++) 
    {
        scanf("%lld",&val[i]);merge1(F[i][0],val[i]);
    }
    for(int i=1;i<n;i++)
    {
        scanf("%d%d",&x,&y);
        f[x].push_back(y);f[y].push_back(x);
    }
    dfs1(1,0,0);
    pre();
    for(int i=1;i<=m;i++)
    {
        memset(ans,0,sizeof(ans));
        scanf("%d%d",&x,&y);
        lca(x,y);
        ll Sum=0;
        for(int j=60;~j;j--) Sum=max(Sum,Sum^ans[j]);
        printf("%lld\n",Sum);
    }
    return 0;
}