解析

雖然使用三個log的倍增算法艸過去了
但是我們還是來聊聊正解吧

考慮點分治
對于當前的根，dfs求出聯通塊內每個點到當前根的線性基
一條路徑的答案應該在路徑出現上第一個成為根的點時統計到
具體來說，就是路徑的兩端點在同一個solve函數的不同子樹內
通過奇怪的時間戳標記可以實現
對于詢問，每個結點開個vector暴力遍歷即可
由于每個結點最多被dfslog遍，所以復雜度是對的
復雜度瓶頸在dfs求線性基上，時間復雜度$nlo{g}^{2}n$
然后就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; } /* */

總結

以上是生活随笔為你收集整理的YBTOJ洛谷P3292：幸运数字（线性基、点分治/倍增）的全部內容，希望文章能夠幫你解決所遇到的問題。

如果覺得生活随笔網站內容還不錯，歡迎將生活随笔推薦給好友。