题目链接

线性基可以\(O(log^2)\)暴力合并。又是树上路径问题,考虑点分治。

对于每个点i求解 LCA(u,v)i 时的询问(u,v),只需求出这个点到其它点的线性基后,暴力合并。

LCA不能直接求啊。。树形态在变。

Get一种新点分治方法,q[x]存x这棵子树的询问,递归时Solve(son[x])。

但是处理询问是对LCAroot的(先找到这棵子树的root),即从root开始DFS,记一下每个点属于哪棵子子树;再枚举这棵子树的询问,如果LCA是root则查询,否则在同一棵子子树的话就把询问分给那棵子树。

注意这样的话改root的信息而不是x,x可能还会再访问,so还要清空x的询问。so处理x这棵子树前要存下询问后再清空,因为还会有给x的询问。。不能直接清了!

注意没有询问时剪个枝。

复杂度\(O(60n\log n+60^2q)\)。

//17604kb	5528ms
#include <cstdio>
#include <cctype>
#include <vector>
#include <cstring>
#include <algorithm>
//#define gc() getchar()
#define MAXIN 200000
#define gc() (SS==TT&&(TT=(SS=IN)+fread(IN,1,MAXIN,stdin),SS==TT)?EOF:*SS++)
#define Bit 60
typedef long long LL;
const int N=2e4+5,Qs=2e5+5; int n,Q,Enum,tmp[Qs],H[N],nxt[N<<1],to[N<<1],sz[N],bel[N],X[Qs],Y[Qs],Min,root;
bool vis[N];
LL A[N],Ans[Qs];
std::vector<int> q[N];
char IN[MAXIN],*SS=IN,*TT=IN;
struct Base
{
LL b[61];
inline void Clear() {memset(b,0,sizeof b);}
inline void Insert(LL x)
{
for(int i=Bit; ~i; --i)
if(x>>i & 1)
if(b[i]) x^=b[i];
else {b[i]=x; break;}
}
inline void Merge(const Base &x)
{
for(int i=Bit; ~i; --i)
if(x.b[i]) Insert(x.b[i]);
}
inline LL Query()
{
LL ans=0;
for(int i=Bit; ~i; --i) ans=std::max(ans,ans^b[i]);
return ans;
}
}base[N]; inline int read()
{
int now=0;register char c=gc();
for(;!isdigit(c);c=gc());
for(;isdigit(c);now=now*10+c-'0',c=gc());
return now;
}
inline LL readll()
{
LL now=0;register char c=gc();
for(;!isdigit(c);c=gc());
for(;isdigit(c);now=now*10+c-'0',c=gc());
return now;
}
inline void AddEdge(int u,int v)
{
to[++Enum]=v, nxt[Enum]=H[u], H[u]=Enum;
to[++Enum]=u, nxt[Enum]=H[v], H[v]=Enum;
}
void Get_root(int x,int f,int tot)
{
int mx=0; sz[x]=1;
for(int v,i=H[x]; i; i=nxt[i])
if(!vis[v=to[i]] && v!=f)
{
Get_root(v,x,tot), sz[x]+=sz[v];
if(sz[v]>mx) mx=sz[v];
}
mx=std::max(mx,tot-sz[x]);
if(mx<Min) Min=mx, root=x;
}
void DFS(int x,int f,int Bel)
{
bel[x]=Bel, base[x]=base[f], base[x].Insert(A[x]);
for(int i=H[x]; i; i=nxt[i])
if(to[i]!=f&&!vis[to[i]]) DFS(to[i],x,Bel);
}
void Solve(int x)
{
if(!q[x].size()) return;//!
Min=N, Get_root(x,x,sz[x]);
vis[root]=1, bel[root]=root,/*!*/ base[root].Clear(), base[root].Insert(A[root]);
for(int i=H[root]; i; i=nxt[i]) if(!vis[to[i]]) DFS(to[i],root,to[i]);
int cnt=q[x].size();
for(int i=0; i<=cnt; ++i) tmp[i]=q[x][i];//!
q[x].clear();//! 当然如果用边表清空就快了
Base b;
for(int i=0,id; i<cnt; ++i)
if(bel[X[id=tmp[i]]]==bel[Y[id]]) q[bel[X[id]]].push_back(id);
else b=base[X[id]], b.Merge(base[Y[id]]), Ans[id]=b.Query();
for(int i=H[root]; i; i=nxt[i]) if(!vis[to[i]]) Solve(to[i]);
} int main()
{
n=read(), Q=read();
for(int i=1; i<=n; ++i) A[i]=readll();
for(int i=1; i<n; ++i) AddEdge(read(),read());
for(int i=1; i<=Q; ++i)
{
X[i]=read(), Y[i]=read();
if(X[i]!=Y[i]) q[1].push_back(i);
else Ans[i]=A[X[i]];
}
sz[1]=n, Solve(1);
for(int i=1; i<=Q; ++i) printf("%lld\n",Ans[i]); return 0;
}

BZOJ.4516.[SCOI2016]幸运数字(线性基 点分治)的更多相关文章

  1. BZOJ 4568: [Scoi2016]幸运数字 [线性基 倍增]

    4568: [Scoi2016]幸运数字 题意:一颗带点权的树,求树上两点间异或值最大子集的异或值 显然要用线性基 可以用倍增的思想,维护每个点向上\(2^j\)个祖先这些点的线性基,求lca的时候合 ...

  2. BZOJ 4568 [Scoi2016]幸运数字 ——线性基 倍增

    [题目分析] 考虑异或的最大值,维护线性基就可以了. 但是有多次的询问,树剖或者倍增都可以. 想了想树剖动辄数百行的代码. 算了,我还是写倍增吧. 注:被位运算和大于号的优先级坑了一次,QaQ [代码 ...

  3. 洛谷P3292 [SCOI2016]幸运数字 线性基+倍增

    P3292 [SCOI2016]幸运数字 传送门 题目描述 A 国共有 n 座城市,这些城市由 n-1 条道路相连,使得任意两座城市可以互达,且路径唯一.每座城市都有一个幸运数字,以纪念碑的形式矗立在 ...

  4. 洛谷P3292 [SCOI2016] 幸运数字 [线性基,倍增]

    题目传送门 幸运数字 题目描述 A 国共有 n 座城市,这些城市由 n-1 条道路相连,使得任意两座城市可以互达,且路径唯一.每座城市都有一个幸运数字,以纪念碑的形式矗立在这座城市的正中心,作为城市的 ...

  5. bzoj4568 [Scoi2016]幸运数字 线性基+树链剖分

    A 国共有 n 座城市,这些城市由 n-1 条道路相连,使得任意两座城市可以互达,且路径唯一.每座城市都有一个 幸运数字,以纪念碑的形式矗立在这座城市的正中心,作为城市的象征.一些旅行者希望游览 A ...

  6. P3292 [SCOI2016]幸运数字 线性基

    正解:线性基+倍增 解题报告: 先放下传送门QAQ 然后这题,其实没什么太大的技术含量,,,?就几个知识点套在一起,除了代码长以外没任何意义,主要因为想复习下线性基的题目所以还是写下,,, 随便写下思 ...

  7. BZOJ4568: [Scoi2016]幸运数字(线性基 倍增)

    题意 题目链接 Sol 线性基是可以合并的 倍增维护一下 然后就做完了?? 喵喵喵? // luogu-judger-enable-o2 #include<bits/stdc++.h> # ...

  8. [SCOI2016]幸运数字 线性基

    题面 题面 题解 题面意思非常明确:求树上一条链的最大异或和. 我们用倍增的思想. 将这条链分成2部分:x ---> lca , lca ---> y 分别求出这2个部分的线性基,然后合并 ...

  9. P3292 [SCOI2016]幸运数字 [线性基+倍增]

    线性基+倍增 // by Isaunoya #include <bits/stdc++.h> using namespace std; #define rep(i, x, y) for ( ...

随机推荐

  1. 关于在函数中使用Array.prototype.slice.call而不是直接用slice

    arguments是每个函数在运行的时候自动获得的一个近似数组的对象(除了length外没有其他属性),这个arguments对象其实并不是Array,所以没有slice方法. Array.proto ...

  2. python collection 中的队列

    认识中的队列 在以前的认知里,队列是先进先出,就是一头进,一头出,Queue.而无意间看到了deque 双向队列. 即从该队列的头或者尾部都能插入和移除元素.而起时间复杂度竟然是一样的!O(1),是不 ...

  3. spring-boog-测试打桩-Mockito

    Mockito用于测试时进行打桩处理:通过它可以指定某个类的某个方法在什么情况下返回什么样的值. 例如:测试 controller时,依赖 service,这个时候就可以假设当调用 service 某 ...

  4. Linux触摸屏驱动测试程序范例【转】

    转自:http://blog.sina.com.cn/s/blog_4b4b54da0102viyl.html 转载2015-05-09 16:28:27 标签:androiditlinux 触摸屏驱 ...

  5. URIEncoding与useBodyEncodingForURI 在tomcat中文乱码处理上的区别

    大家知道tomcat5.0开始,对网页的中文字符的post或者get,经常会出现乱码现象. 具体是因为Tomcat默认是按ISO-8859-1进行URL解码,ISO-8859-1并未包括中文字符,这样 ...

  6. 十一、springboot之web开发之Filter

    我们常常在项目中会使用filters用于录调用日志.排除有XSS威胁的字符.执行权限验证等等.Spring Boot自动添加了OrderedCharacterEncodingFilter和Hidden ...

  7. cvpr densnet论文

  8. 洛谷P2017晕牛

    传送门啦 这个题没有想象中复杂. 我们先有向边建立,因为我们无法改变有向边的方向. 建立完之后跑一边拓扑排序,我们按拓扑序小的指向大的就好了. 解释: 我们知道如果一个点在另一个点顺序的后面的话,如果 ...

  9. Laravel 程序架构设计思路:使用动作类

    当我们谈论到应用程序的架构的时候,经常会问到一个经典的问题,那就是"这段代码应该放在哪里比较好". 因为 Laravel 是一个相当灵活的框架,所以要回答这个问题其实没那么容易.我 ...

  10. Django数据库数据表操作

    建立表单 django通过设置类来快速建表,打开models.py 例: from __future__ import unicode_literals from django.db import m ...