解析
虽然使用三个log的倍增算法艹过去了
但是我们还是来聊聊正解吧
考虑点分治
对于当前的根,dfs求出联通块内每个点到当前根的线性基
一条路径的答案应该在路径出现上第一个成为根的点时统计到
具体来说,就是路径的两端点在同一个solve函数的不同子树内
通过奇怪的时间戳标记可以实现
对于询问,每个结点开个vector暴力遍历即可
由于每个结点最多被dfslog遍,所以复杂度是对的
复杂度瓶颈在dfs求线性基上,时间复杂度nlog2nnlog^2nnlog2n
然后就ok啦
但是懒得再写一遍了,所以还是贴的3log代码
代码
#include<bits/stdc++.h>
using namespace std;
#define ll long long
#define inf (n+1)
//#define debug(...) fprintf(stderr,__VA_ARGS__)
const int N=2e4+100;
const int M=2e5+10500;
const double eps=1e-5;
inline ll read(){ll x(0),f(1);char c=getchar();while(!isdigit(c)){if(c=='-')f=-1;c=getchar();}while(isdigit(c)){x=(x<<1)+(x<<3)+c-'0';c=getchar();}return x*f;
}int n,m;struct node{int to,nxt;
}p[N<<1];
int fi[N],cnt;
inline void addline(int x,int y){p[++cnt]=(node){y,fi[x]};fi[x]=cnt;return;
}ll mi[65];struct basic{ll v[62];basic(){memset(v,0,sizeof(v));}void upd(ll x){for(int i=60;i>=0;i--){if(x&mi[i]){if(!v[i]){v[i]=x;break;}else x^=v[i];}}}
}b[N][16];
inline void merge(basic &v1,basic v2){for(int i=0;i<=60;i++){if(v2.v[i]) v1.upd(v2.v[i]);}return;
}ll val[N];
int pl[N][16],dep[N];
void dfs(int x,int f){b[x][0].upd(val[f]);pl[x][0]=f;for(int k=1;pl[x][k-1];k++){pl[x][k]=pl[pl[x][k-1]][k-1];b[x][k]=b[x][k-1];merge(b[x][k],b[pl[x][k-1]][k-1]);}dep[x]=dep[f]+1;for(int i=fi[x];~i;i=p[i].nxt){int to=p[i].to;if(to==f) continue;dfs(to,x);}return;
}
basic res;
ll query(int x,int y){memset(res.v,0,sizeof(res.v));res.upd(val[x]);res.upd(val[y]);if(dep[x]<dep[y]) swap(x,y);for(int k=15;k>=0;k--){if(dep[x]-mi[k]<dep[y]) continue;merge(res,b[x][k]);x=pl[x][k];}if(x!=y){for(int k=15;k>=0;k--){if(pl[x][k]==pl[y][k]) continue;merge(res,b[x][k]);merge(res,b[y][k]);x=pl[x][k];y=pl[y][k];}res.upd(val[pl[x][0]]);}ll ans(0);for(int i=60;i>=0;i--){//if(res.v[i]) printf("i=%d v=%lld\n",i,res.v[i]);if(!res.v[i]||(ans&mi[i])) continue;ans^=res.v[i];//printf(" ans=%lld\n",ans);}//putchar('\n');return ans;
}
int main(){
#ifndef ONLINE_JUDGE//freopen("a.in","r",stdin);//freopen("a.out","w",stdout);
#endifmemset(fi,-1,sizeof(fi));cnt=-1;mi[0]=1;for(int i=1;i<=60;i++) mi[i]=(mi[i-1]<<1);n=read();m=read();for(int i=1;i<=n;i++) val[i]=read();for(int i=1;i<n;i++){int x=read(),y=read();addline(x,y);addline(y,x);}dfs(1,0);for(int i=1;i<=m;i++){int x=read(),y=read();printf("%lld\n",query(x,y));}return 0;
}
/*
*/