解析

虽然使用三个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;
}
/*
*/

YBTOJ洛谷P3292：幸运数字（线性基、点分治/倍增）相关推荐

BZOJ4568: [Scoi2016]幸运数字(线性基倍增)
题意题目链接 Sol 线性基是可以合并的倍增维护一下然后就做完了?? 喵喵喵? // luogu-judger-enable-o2 #include<bits/stdc++.h> # ...
YBTOJ洛谷P4068：数字配对（网络流）
文章目录题目描述解析代码题目描述有n种数字,第 i 种数字是 ai,有 bi个,权值是 ci. 若两个数字 ai,aj 满足, ai 是 aj 的倍数,且 ai/aj 是一个质数,那么这两个 ...
YBTOJ洛谷P4331：数字序列（左偏树）
文章目录题目描述数据范围解析代码题目描述数据范围 n<=1e6n<=1e6n<=1e6 解析先考虑简单情况如果原数列是单调递增的,显然应该使bi=aib_i=a_ib ...
[洛谷P3292] [SCOI2016]幸运数字
洛谷题目链接:[SCOI2016]幸运数字题目描述 A 国共有 n 座城市,这些城市由 n-1 条道路相连,使得任意两座城市可以互达,且路径唯一.每座城市都有一个幸运数字,以纪念碑的形式矗立在这座城 ...
洛谷试炼场 P1553 数字反转（升级版）题解
洛谷试炼场 P1553 数字反转(升级版)题解 [c] 题目描述给定一个数,请将该数各个位上数字反转得到一个新数. 这次与NOIp2011普及组第一题不同的是:这个数可以是小数,分数,百分数,整数. ...
洛谷P4390 [BOI2007]Mokia 摩基亚题解
洛谷P4390 [BOI2007]Mokia 摩基亚题解题目链接:P4390 [BOI2007]Mokia 摩基亚题意:摩尔瓦多的移动电话公司摩基亚(Mokia)设计出了一种新的用户定位系统.和 ...
洛谷P3292 [SCOI2016]幸运数字（倍增+线性基）
传送门不知道线性基是什么东西的可以看看蒟蒻的总结第一眼:这不会是个倍增LCA暴力合并线性基吧-- 打了一发--A了? 所以这真的是个暴力倍增LCA合并线性基么-- ps:据某大佬说其实可以离线之后 ...
幸运数字（洛谷-P3292）
题目描述 A 国共有 n 座城市,这些城市由 n-1 条道路相连,使得任意两座城市可以互达,且路径唯一.每座城市都有一个幸运数字,以纪念碑的形式矗立在这座城市的正中心,作为城市的象征. 一些旅行者希望 ...
YBTOJ洛谷P4869：出现位置（线性基）
解析关键结论: 若 nnn 个数组成的线性基大小为 SSS,则其子集异或组成的结果有 2S2^S2S 种,且每种结果都有 2n−S2^{n-S}2n−S 种方案. 证明:考虑 n−Sn-Sn−S 个 ...

YBTOJ洛谷P3292：幸运数字（线性基、点分治/倍增）

解析

代码

YBTOJ洛谷P3292：幸运数字（线性基、点分治/倍增）相关推荐

最新文章

热门文章