内存限制:512 MiB 时间限制:1000 ms

题目描述

给定长度为n的数列X={x1,x2,...,xn}和长度为m的数列Y={y1,y2,...,ym},令矩阵A中第i行第j列的值Aij=xi xor yj,每次询问给定矩形区域i∈[u,d],j∈[l,r],找出第k大的Aij。

输入格式

第一行包含两个正整数n,m,分别表示两个数列的长度

第二行包含n个非负整数xi

第三行包含m个非负整数yj

第四行包含一个正整数p,表示询问次数

随后p行,每行均包含5个正整数,用来描述一次询问,每行包含五个正整数u,d,l,r,k,含义如题意所述。

输出格式

共p行,每行包含一个非负整数,表示此次询问的答案。

样例

样例输入

3 3
1 2 4
7 6 5
3
1 2 1 2 2
1 2 1 3 4
2 3 2 3 4

样例输出

6
5
1

数据范围与提示

对于100%的数据,0<=Xi,Yj<2^31,

1<=u<=d<=n<=1000,

1<=l<=r<=m<=300000,

1<=k<=(d-u+1)*(r-l+1),

1<=p<=500


区间有关异或的问题,都可以用trie树解决

本题的$m$,显然用可持久化trie维护

对于每个询问,考虑枚举$[u,d]$,在区间$[l,r]$上的trie中查询。

每次从答案的最高位开始,判断区间内该位为$1$的A_{i,j}的个数,如果$>=k$则答案的该位则为$1$。

对每个枚举的$i(u<=i<=d)$,维护$L,R$两个指针表示当前询问的区间,每次判断完更新

#include<iostream>
#include<cstdio>
#include<cstring>
using namespace std;
int read(){
char c=getchar(); int x=;
while(c<''||c>'') c=getchar();
while(''<=c&&c<='') x=x*+c-,c=getchar();
return x;
}
#define N 1005
#define M 300005
int n,m,a[N],L[N],R[N];
int u,rt[M],nx[M*][],sum[M*];
void ins(int o,int p,int d,int v){
int k=(v>>d)&;
nx[o][!k]=nx[p][!k];
sum[nx[o][k]=++u]=sum[nx[p][k]]+;
if(d) ins(nx[o][k],nx[p][k],d-,v);
}
int ask(int l,int r,int k,int d){
if(d<) return ;
int tt=,b=<<d;
for(int i=l;i<=r;++i){
if(a[i]&b) tt+=sum[nx[R[i]][]]-sum[nx[L[i]][]];
else tt+=sum[nx[R[i]][]]-sum[nx[L[i]][]];
}
if(tt>=k){//该位为1的个数>=k,说明答案的该位为1
for(int i=l;i<=r;++i){
if(a[i]&b) L[i]=nx[L[i]][],R[i]=nx[R[i]][];
else L[i]=nx[L[i]][],R[i]=nx[R[i]][];
}
return ask(l,r,k,d-)|b;
}else{
for(int i=l;i<=r;++i){
if(a[i]&b) L[i]=nx[L[i]][],R[i]=nx[R[i]][];
else L[i]=nx[L[i]][],R[i]=nx[R[i]][];
}
return ask(l,r,k-tt,d-);
}
}
int main(){
n=read(); m=read();
for(int i=;i<=n;++i) a[i]=read();
for(int i=;i<=m;++i) rt[i]=++u,ins(rt[i],rt[i-],,read());
for(int Q=read(),l1,l2,r1,r2,k;Q;--Q){
l1=read(),l2=read(),r1=read(),r2=read(),k=read();
for(int i=l1;i<=l2;++i) L[i]=rt[r1-],R[i]=rt[r2];
printf("%d\n",ask(l1,l2,k,));
}return ;
}

bzoj4103 [Thu Summer Camp 2015]异或运算(可持久化trie)的更多相关文章

  1. [BZOJ4103][Thu Summer Camp 2015]异或运算 可持久化Trie树

    4103: [Thu Summer Camp 2015]异或运算 Time Limit: 20 Sec  Memory Limit: 512 MB Description 给定长度为n的数列X={x1 ...

  2. 【bzoj4103】[Thu Summer Camp 2015]异或运算 可持久化trie树

    Description 给定长度为n的数列X={x1,x2,...,xn}和长度为m的数列Y={y1,y2,...,ym},令矩阵A中第i行第j列的值Aij=xi xor yj,每次询问给定矩形区域i ...

  3. BZOJ 4103: [Thu Summer Camp 2015]异或运算 可持久化trie

    开始想了一个二分+可持久化trie验证,比正解多一个 log 仔细思考,你发现你可以直接按位枚举,然后在可持久化 trie 上二分就好了. code: #include <bits/stdc++ ...

  4. BZOJ4103 [Thu Summer Camp 2015]异或运算 【可持久化trie树】

    题目链接 BZOJ4103 题解 一眼看过去是二维结构,实则未然需要树套树之类的数据结构 区域异或和,就一定是可持久化\(trie\)树 观察数据,\(m\)非常大,而\(n\)和\(p\)比较小,甚 ...

  5. bzoj4103: [Thu Summer Camp 2015]异或运算

    对于每个询问暴力枚举x~y,然后在Trie去找第k大,开始我写了个二分答案然后算比当前答案大的个数,打了个第10个点的表就跑出19s+比bzoj垫底还慢4s+ 然而不用二分,直接1000个点一起在树上 ...

  6. 【BZOJ 4103】 [Thu Summer Camp 2015]异或运算 可持久化01Trie

    我们观察数据:树套树 PASS    主席树 PASS  一层一个Trie PASS 再看,异或!我们就把目光暂时定在01Tire然后我们发现,我们可以带着一堆点在01Trie上行走,因为O(n*q* ...

  7. 【BZOJ 4103】 4103: [Thu Summer Camp 2015]异或运算 (可持久化Trie)

    4103: [Thu Summer Camp 2015]异或运算 Time Limit: 20 Sec  Memory Limit: 512 MBSubmit: 474  Solved: 258 De ...

  8. [BZOJ 4103] [Thu Summer Camp 2015] 异或运算 【可持久化Trie】

    题目链接:BZOJ - 4103 题目分析 THUSC滚粗之后一直没有写这道题,从来没写过可持久化Trie,发现其实和可持久化线段树都是一样的.嗯,有些东西就是明白得太晚. 首先Orz ZYF-ZYF ...

  9. 【BZOJ 4104】【Thu Summer Camp 2015】解密运算

    http://www.lydsy.com/JudgeOnline/problem.php?id=4104 网上题解满天飞,我也懒得写了 #include<cstdio> #include& ...

随机推荐

  1. Memcached在Asp.net下的应用

    Memcached 是一个高性能的分布式内存对象缓存系统,用于动态Web应用以减轻数据库负载.它通过在内存中缓存数据和对象来减少读取数据库的次数,从而提高动态.数据库驱动网站的速度.Memcached ...

  2. C# 向上取整数

    PageCount = personInfodb.get_count(selected_id); //数据总数 PageCount = (int)Math.Ceiling(PageCount / (P ...

  3. 查看 XML 文件

    在所有现代浏览器中,均能够查看原始的 XML 文件. 不要指望 XML 文件会直接显示为 HTML 页面. 查看 XML 文件 查看这个 XML 文件:note.xml 打开 XML 文件 - XML ...

  4. Sublime Text 3 程序运行后中文显示乱码的解决方案

    问题描述: 问题分析: 新建python文件,输入下面代码,测试发现发送给python解释器时,默认为cp936的编码方式,而py文件为UTF-8的编码方式,cp936和UTF-8的编码方式对于英文字 ...

  5. HDU1575--Tr A(矩阵快速幂)

    Time Limit: 1000/1000 MS (Java/Others) Memory Limit: 32768/32768 K (Java/Others) Total Submission(s) ...

  6. ORA-01034、ORA-27101、ORA-00119错误

    ORA-01034ORA-27101Linux-x86_64 Error: 2:No such file or directory. 或者是 ORA-00119: invalid specificat ...

  7. Hbase数据备份&&容灾方案

    Hbase数据备份&&容灾方案 标签(空格分隔): Hbase 一.Distcp 在使用distcp命令copy hdfs文件的方式实现备份时,需要禁用备份表确保copy时该表没有数据 ...

  8. sqli-labs(36)

    0X01发现又是‘’被过滤了 ?id=%df%%20and%=%23 错误 ?id=1%df%27%20and%201=1%23 正确   存在注入 0X01爆数据库 ?id=-%df%%20unio ...

  9. [BZOJ1195]:[HNOI2006]最短母串(AC自动机+BFS)

    题目传送门 题目描述 给定n个字符串(S1,S2,…,Sn),要求找到一个最短的字符串T,使得这n个字符串(S1,S2,…,Sn)都是T的子串. 输入格式 第一行是一个正整数n,表示给定的字符串的个数 ...

  10. centos7安装nvidia驱动

    1. disable UEFI security boot! 2.添加 ELRepo 源: Import the public key: rpm --import https://www.elrepo ...