[bzoj4568] [loj#2013] [Scoi2016] 幸运数字
Description
\(A\) 国共有 \(n\) 座城市,这些城市由 \(n-1\) 条道路相连,使得任意两座城市可以互达,且路径唯一。每座城市都有一个幸运数字,以纪念碑的形式矗立在这座城市的正中心,作为城市的象征。一些旅行者希望游览 \(A\) 国。旅行者计划乘飞机降落在 \(x\) 号城市,沿着 \(x\) 号城市到 \(y\) 号城市之间那条唯一的路径游览,最终从 \(y\) 城市起飞离开 \(A\) 国。在经过每一座城市时,游览者就会有机会与这座城市的幸运数字拍照,从而将这份幸运保存到自己身上。然而,幸运是不能简单叠加的,这一点游览者也十分清楚。他们迷信着幸运数字是以异或的方式保留在自己身上的。例如,游览者拍了 3 张照片,幸运值分别是 5,7,11,那么最终保留在自己身上的幸运值就是 9(5 \(xor\) 7 \(xor\) 11)。有些聪明的游览者发现,只要选择性地进行拍照,便能获得更大的幸运值。例如在上述三个幸运值中,只选择 5 和 11 ,可以保留的幸运值为 14 。现在,一些游览者找到了聪明的你,希望你帮他们计算出在他们的行程安排中可以保留的最大幸运值是多少。
Input
第一行包含 \(2\) 个正整数 \(n\) ,\(q\),分别表示城市的数量和旅行者数量。第二行包含 \(n\) 个非负整数,其中第 \(i\) 个整数 \(Gi\) 表示 \(i\) 号城市的幸运值。随后 \(n-1\) 行,每行包含两个正整数 \(x\) ,\(y\),表示 \(x\) 号城市和 \(y\) 号城市之间有一条道路相连。随后 \(q\) 行,每行包含两个正整数 \(x\) ,\(y\),表示这名旅行者的旅行计划是从 \(x\) 号城市到 \(y\) 号城市。
\(N \leq20000,Q \leq 200000,Gi \leq 2^60\)
Output
输出需要包含 \(q\) 行,每行包含 \(1\) 个非负整数,表示这名旅行者可以保留的最大幸运值。
Sample Input
4 2
11 5 7 9
1 2
1 3
1 4
2 3
1 4
Sample Output
14
11
想法
要求一堆数中的最大异或和要用线性基。
那么这道题我们就要把一条路径上的所有店的值扔到线性基里去求最大。
怎么扔进去呢?倍增 ,维护每个点网上跳 \(2^i\) 个点这条链上的线性基
查询的时候 \(lca\) 的时候合并一下线性基就可以了。
没有卡常,极其感动。
代码
#include<cstdio>
#include<iostream>
#include<algorithm>
#include<cstring>
using namespace std;
const int N = 20005;
typedef long long ll;
ll lread(){
ll x=0;
char ch=getchar();
while(!isdigit(ch)) ch=getchar();
while(isdigit(ch)) x=x*10+ch-'0',ch=getchar();
return x;
}
int read(){
int x=0;
char ch=getchar();
while(!isdigit(ch)) ch=getchar();
while(isdigit(ch)) x=x*10+ch-'0',ch=getchar();
return x;
}
int n;
ll G[N];
struct node{
int v;
node *nxt;
}pool[N*2],*h[N];
int cnt;
void addedge(int u,int v){
node *p=&pool[++cnt],*q=&pool[++cnt];
p->v=v;p->nxt=h[u];h[u]=p;
q->v=u;q->nxt=h[v];h[v]=q;
}
int f[N][16],dep[N];
ll g[N][16][65];
void ins(ll *A,ll x) {
for(int i=60;i>=0;i--)
if(x&(1ll<<i)){
if(!A[i]) { A[i]=x; return; }
x^=A[i];
}
}
void merge(ll *A,ll *B) { for(int i=60;i>=0;i--) if(B[i]) ins(A,B[i]); }
ll MAX(ll *A) {
ll ret=0;
for(int i=60;i>=0;i--) ret=max(ret,ret^A[i]);
return ret;
}
void dfs(int u){
int v;
for(node *p=h[u];p;p=p->nxt)
if(!dep[v=p->v]){
dep[v]=dep[u]+1;
f[v][0]=u; ins(g[v][0],G[u]);
for(int j=1;j<16;j++){
f[v][j]=f[f[v][j-1]][j-1];
memcpy(g[v][j],g[v][j-1],sizeof(g[v][j]));
merge(g[v][j],g[f[v][j-1]][j-1]);
}
dfs(v);
}
}
ll lca(int x,int y){
ll c[65]; memset(c,0,sizeof(c));
ins(c,G[x]); ins(c,G[y]);
if(dep[x]<dep[y]) swap(x,y);
for(int i=15;i>=0;i--)
if(dep[f[x][i]]>=dep[y]) merge(c,g[x][i]),x=f[x][i];
if(x==y) return MAX(c);
for(int i=15;i>=0;i--)
if(f[x][i]!=f[y][i]){
merge(c,g[x][i]); x=f[x][i];
merge(c,g[y][i]); y=f[y][i];
}
merge(c,g[x][0]);
return MAX(c);
}
int main()
{
int Q;
n=read(); Q=read();
for(int i=1;i<=n;i++) G[i]=lread();
for(int i=1;i<n;i++) addedge(read(),read());
dep[1]=1; dfs(1);
while(Q--)
printf("%lld\n",lca(read(),read()));
return 0;
}
[bzoj4568] [loj#2013] [Scoi2016] 幸运数字的更多相关文章
- [BZOJ4568][SCOI2016]幸运数字(倍增LCA,点分治+线性基)
4568: [Scoi2016]幸运数字 Time Limit: 60 Sec Memory Limit: 256 MBSubmit: 2131 Solved: 865[Submit][Statu ...
- 【BZOJ4568】[Scoi2016]幸运数字 倍增+线性基
[BZOJ4568][Scoi2016]幸运数字 Description A 国共有 n 座城市,这些城市由 n-1 条道路相连,使得任意两座城市可以互达,且路径唯一.每座城市都有一个幸运数字,以纪念 ...
- [BZOJ4568][Scoi2016]幸运数字 倍增+线性基
4568: [Scoi2016]幸运数字 Time Limit: 60 Sec Memory Limit: 256 MBSubmit: 1791 Solved: 685[Submit][Statu ...
- bzoj4568: [Scoi2016]幸运数字(LCA+线性基)
4568: [Scoi2016]幸运数字 题目:传送门 题解: 好题!!! 之前就看过,当时说是要用线性基...就没学 填坑填坑: %%%线性基 && 神犇 主要还是对于线性基的运用和 ...
- BZOJ 4568: [Scoi2016]幸运数字 [线性基 倍增]
4568: [Scoi2016]幸运数字 题意:一颗带点权的树,求树上两点间异或值最大子集的异或值 显然要用线性基 可以用倍增的思想,维护每个点向上\(2^j\)个祖先这些点的线性基,求lca的时候合 ...
- [SCOI2016]幸运数字 树链剖分,线性基
[SCOI2016]幸运数字 LG传送门 为了快乐,我们用树剖写这题. 强行树剖,线段树上每个结点维护一个线性基,每次查询暴力合并. 瞎分析一波复杂度:树剖两点之间\(\log n\)条重链,每条重链 ...
- bzoj 4568: [Scoi2016]幸运数字
4568: [Scoi2016]幸运数字 Time Limit: 60 Sec Memory Limit: 256 MBSubmit: 848 Solved: 336[Submit][Status ...
- [洛谷P3292] [SCOI2016]幸运数字
洛谷题目链接:[SCOI2016]幸运数字 题目描述 A 国共有 n 座城市,这些城市由 n-1 条道路相连,使得任意两座城市可以互达,且路径唯一.每座城市都有一个幸运数字,以纪念碑的形式矗立在这座城 ...
- 【BZOJ 4568】 4568: [Scoi2016]幸运数字 (线性基+树链剖分+线段树)
4568: [Scoi2016]幸运数字 Description A 国共有 n 座城市,这些城市由 n-1 条道路相连,使得任意两座城市可以互达,且路径唯一.每座城市都有一个 幸运数字,以纪念碑的形 ...
随机推荐
- 2018-11-17-win10-uwp-在-xaml-让-TextBlock-换行
title author date CreateTime categories win10 uwp 在 xaml 让 TextBlock 换行 lindexi 2018-11-17 16:2:29 + ...
- How to fix nuget Unrecognized license type MIT when pack
When I packaging license within the nupkg, I will using License to replace licentUrl. I using this c ...
- Dart常见类型转换 Int String Double
int -> string age.toString() string -> int int.parse('100'); String -> double var onePointO ...
- 洛谷$P4045\ [JSOI2009]$密码 $dp$+$AC$自动机
正解:$dp$+$AC$自动机+搜索 解题报告: 传送门$QwQ$ 首先显然先建个$AC$自动机,然后考虑设$f_{i,j,k}$表示长度为$i$,现在在$AC$自动机的第$j$个位置,已经表示出来的 ...
- JAVA字节码文件之常量池
一.常量池的内容 一个java类中定义的很多信息都是由常量池来维护和描述的,可以将常量池看作是class文件的资源仓库,比如java类中定义的方法与变量信息.常量池中主要存储两类常量:字面量(文本字符 ...
- 大数据学习之路-hdfs
1.什么是hadoop hadoop中有3个核心组件: 分布式文件系统:HDFS —— 实现将文件分布式存储在很多的服务器上 分布式运算编程框架:MAPREDUCE —— 实现在很多机器上分布式并行运 ...
- 信息熵为什么要定义成-Σp*log(p)?
信息熵为什么要定义成-Σp*log(p)? 再解释信息熵之前,需要先来说说什么是信息量. 信息量是对信息的度量,单位一般用bit. 信息论之父克劳德·艾尔伍德·香农(Claude Elwood Sha ...
- 简单聊一聊JS中的循环引用及问题
本文主要从 JS 中为什么会出现循环引用,垃圾回收策略中引用计数为什么有很大的问题,以及循环引用时的对象在使用 JSON.stringify 时为什么会报错,怎样解决这个问题简单谈谈自己的一些理解. ...
- (01)大话设计模式-简单工厂-java实现
1.运算接口 public interface Operation { public double getResult(double NumberA , double NumberB); } 2.加减 ...
- 手算CRC及其实现
前言: 这篇文章主要讲的是如何手算CRC以及运用CRC,更侧重方法的步骤,对原理方面不做探讨. 方法也是按照我个人理解的来,所以并不专业~ 一些搬过来的代码我也修改了一下下 如果想了解原理的可参考资料 ...