#2013. 「SCOI2016」幸运数字

题目描述

A 国共有 n nn 座城市,这些城市由 n−1 n - 1n−1 条道路相连,使得任意两座城市可以互达,且路径唯一。每座城市都有一个幸运数字,以纪念碑的形式矗立在这座城市的正中心,作为城市的象征。一些旅行者希望游览 A 国。旅行者计划乘飞机降落在 x xx 号城市,沿着 x xx 号城市到 y yy 号城市之间那条唯一的路径游览,最终从 y yy 城市起飞离开 A 国。

在经过每一座城市时,游览者就会有机会与这座城市的幸运数字拍照,从而将这份幸运保存到自己身上。然而,幸运是不能简单叠加的,这一点游览者也十分清楚。他们迷信着幸运数字是以异或的方式保留在自己身上的。例如,游览者拍了 3 33 张照片,幸运值分别是 5 55、7 77、11 1111,那么最终保留在自己身上的幸运值就是 9 99(5xor7xor11)。

有些聪明的游览者发现,只要选择性地进行拍照,便能获得更大的幸运值。例如在上述三个幸运值中,只选择 5 55 和 11 1111 ,可以保留的幸运值为 14 1414 。现在,一些游览者找到了聪明的你,希望你帮他们计算出在他们的行程安排中可以保留的最大幸运值是多少。

输入格式

第一行包含两个正整数 n nn、q qq,分别表示城市的数量和旅行者数量。
第二行包含 n nn 个非负整数,其中第 i ii 个整数 Gi G_iG​i​​ 表示 i ii 号城市的幸运值。随后 n−1 n - 1n−1 行,每行包含两个正整数 x xx、y yy,表示 x xx 号城市和 y yy 号城市之间有一条道路相连。
随后 q qq 行,每行包含两个正整数 x xx、y yy,表示这名旅行者的旅行计划是从 x xx 号城市到 y yy 号城市。

输出格式

输出需要包含 q qq 行,每行包含 1 11 个非负整数,表示这名旅行者可以保留的最大幸运值。

样例

样例输入

4 2
11 5 7 9
1 2
1 3
1 4
2 3
1 4

样例输出

14
11

数据范围与提示

N≤20000,Q≤200000,Gi≤260 N \leq 20000, Q \leq 200000, G_i \leq 2 ^ {60}N≤20000,Q≤200000,G​i​​≤2​60​​

#include<iostream>
#include<cstdio>
#include<cstring>
#include<cmath>
#define maxn 20010
using namespace std;
int n,q,num,head[maxn],fa[maxn][],dep[maxn];
long long st[maxn][][],t[],lb[],ans;
struct node{int to,pre;}e[maxn*];
long long qread(){
long long i=,j=;
char ch=getchar();
while(ch<''||ch>''){if(ch=='-')j=-;ch=getchar();}
while(ch<=''&&ch>=''){i=i*+ch-'';ch=getchar();}
return i*j;
}
void Insert(int from,int to){
e[++num].to=to;
e[num].pre=head[from];
head[from]=num;
}
void dfs(int x){
for(int i=;dep[x]-(<<i)>;i++)
fa[x][i]=fa[fa[x][i-]][i-];
for(int i=head[x];i;i=e[i].pre){
int to=e[i].to;
if(to==fa[x][])continue;
dep[to]=dep[x]+;fa[to][]=x;
dfs(to);
}
}
void merge(long long a[],long long c[]){
for(int i=;i<=;i++)t[i]=a[i];
for(int i=;i<=;i++)
for(int j=;j>=;j--){
if(!t[i])break;
if(t[i]>>j&){
if(!c[j]){c[j]=t[i];break;}
else t[i]^=c[j];
}
}
}
int lca(int x,int y){
if(dep[x]<dep[y])swap(x,y);
for(int i=;i>=;i--)
if(fa[x][i]&&dep[fa[x][i]]>dep[y])
x=fa[x][i];
if(dep[x]>dep[y])x=fa[x][];
if(x==y)return x;
for(int i=;i>=;i--)
if(fa[x][i]&&fa[y][i]&&fa[x][i]!=fa[y][i])
x=fa[x][i],y=fa[y][i];
return fa[x][];
}
int cal(int x,int goal){
if(dep[x]==goal)return x;
for(int i=;i>=;i--)
if(fa[x][i]&&dep[fa[x][i]]>goal)
x=fa[x][i];
return fa[x][];
}
int main(){
scanf("%d%d",&n,&q);
for(int i=;i<=n;i++){
long long x=qread();
for(int j=;j>=;j--)
if(x>>j&){
st[i][][j]=x;
break;
}
}
for(int i=;i<n;i++){
int x=qread(),y=qread();
Insert(x,y);Insert(y,x);
}
dep[]=;
dfs();
for(int i=;<<i<=n;i++)
for(int j=;j<=n;j++)
if(dep[j]-(<<i)>=)
merge(st[j][i-],st[j][i]),
merge(st[cal(j,dep[j]-(<<i-))][i-],st[j][i]);
int x,y,z,p1,p2;
while(q--){
memset(lb,,sizeof(lb));
scanf("%d%d",&x,&y);
ans=;z=lca(x,y);
p1=log2(dep[x]-dep[z]+);
merge(st[x][p1],lb);
merge(st[cal(x,dep[z]+(<<p1)-)][p1],lb);
p2=log2(dep[y]-dep[z]+);
merge(st[y][p2],lb);
merge(st[cal(y,dep[z]+(<<p2)-)][p2],lb);
for(int i=;i>=;i--)
if(!(ans>>i&))ans^=lb[i];
printf("%lld\n",ans);
}
return ;
}

loj #2013. 「SCOI2016」幸运数字的更多相关文章

  1. loj#2013. 「SCOI2016」幸运数字 点分治/线性基

    题目链接 loj#2013. 「SCOI2016」幸运数字 题解 和树上路径有管...点分治吧 把询问挂到点上 求出重心后,求出重心到每个点路径上的数的线性基 对于重心为lca的合并寻味,否则标记下传 ...

  2. LOJ #2013「SCOI2016」幸运数字

    时限为什么这么大啊 明摆着放多$ log$的做法过啊$QAQ$ LOJ #2013 题意 有$ Q$次询问,每次询问树上一条链,点有点权,你需要选择一些链上的点使得异或和尽量大 点数$ \leq 2* ...

  3. 【LOJ】 #2013. 「SCOI2016」幸运数字

    题解 最大异或和,明显是个线性基 然而还有那么多路径--那就树分治,反正点数看起来很少,就是为了让人乘上一个60的常数嘛 把一个树的点分树记录下来,然后看看询问的两个点彼此相同的最后一个父亲是谁,把这 ...

  4. AC日记——「SCOI2016」幸运数字 LiBreOJ 2013

    「SCOI2016」幸运数字 思路: 线性基: 代码: #include <bits/stdc++.h> using namespace std; #define maxn 20005 # ...

  5. 「洛谷3292」「BZOJ4568」「SCOI2016」幸运数字【倍增LCA+线性基+合并】

    [bzoj数据下载地址]不要谢我 先讲一下窝是怎么错的... \(MLE\)是因为数组开小了.. 看到异或和最大,那么就会想到用线性基. 如果不会线性基的可以参考一下我的学习笔记:「线性基」学习笔记a ...

  6. loj2013 「SCOI2016」幸运数字

    点分治+线性基 (为了这六个字窝调了一下午一晚上QAQ #include <iostream> #include <cstring> #include <cstdio&g ...

  7. loj#2015. 「SCOI2016」妖怪 凸函数/三分

    题目链接 loj#2015. 「SCOI2016」妖怪 题解 对于每一项展开 的到\(atk+\frac{dnf}{b}a + dnf + \frac{atk}{a} b\) 令$T = \frac{ ...

  8. loj#2016. 「SCOI2016」美味

    题目链接 loj#2016. 「SCOI2016」美味 题解 对于不带x的怎么做....可持久化trie树 对于带x,和trie树一样贪心 对于答案的二进制位,从高往低位贪心, 二进制可以表示所有的数 ...

  9. loj#2012. 「SCOI2016」背单词

    题目链接 loj#2012. 「SCOI2016」背单词 题解 题面描述有点不清楚. 考虑贪心 type1的花费一定不会是优的,不考虑, 所以先把后缀填进去,对于反串建trie树, 先填父亲再填儿子, ...

随机推荐

  1. kafka集群配置和java编写生产者消费者操作例子

    kafka 安装 修改配置文件 java操作kafka kafka kafka的操作相对来说简单很多 安装 下载kafka http://kafka.apache.org/downloads tar ...

  2. vc6.0出现“cannot add new member”解决办法

    在编辑一个工程的时候,添加一个按钮,对这个按钮产生消息函数 双击按钮之后出现 cannot add new member 网上说删除.clw文件 我发现并没有这个后缀的文件 以下是我的操作办法:   ...

  3. DRF - 序列化组件(GET/PUT/DELETE接口设计)、视图优化组件

    一.序列化组件 基于上篇随笔的表结构 , 通过序列化组件的ModelSerializer设计如下三个接口 : GET 127.0.0.1:8000/books/{id} # 获取一条数据,返回值:{} ...

  4. [ShaderStaff] Sprite Outline外轮廓效果

    操作系统:Windows8.1 显卡:Nivida GTX965M 开发工具:Unity5.3.8f1 Unity提供了2D Object Sprite对象,但是没有提供外轮廓Outline效果的支持 ...

  5. error:crosses initialization of ...的解决办法

    switch(c) { case 0x01: int temp = a + b; .... break; case 0x02: break; default:break; } 此时会报如题所示错误 原 ...

  6. laravel数据迁移(创建错误列不能创建)

    创建数据表的命令 php artisan make:migration create_users_table 执行这个迁移的命令, php artisan migrate 其实感觉就像简单的方法创建数 ...

  7. JS事件冒泡和事件捕获的详解

    在学校,听老师讲解事件冒泡和事件捕获机制的时候跟听天书一样,只依稀记得IE使用的是事件冒泡,其他浏览器则是事件捕获.当时的我,把它当成IE浏览器兼容问题,所以没有深究(IE8以下版本的浏览器已基本退出 ...

  8. selenium2 定位 窗体切换等等 (二)

    定位用的html素材有两个 demo.html <html> <head> <title>UI Automation Testing</title> & ...

  9. opennebula 一些问t题讨论

    ou_ian - June 8th, 2011 11:13 am非常想请教一下如何设置VM的ip为DHCP? 我们的VM都需要通过eth0 (通过公司的DHCP获得ip)来连接到公司的网络(Inter ...

  10. 【转】Provisional headers are shown

    在chrome开发者工具的 Network 面板中,某些请求头后面会跟着下面这行文字: Provisional headers are shown 这种请求实际上根本没有产生,对应的请求头当然也不应该 ...