可持续化Tried的应用1:18湖北省赛(I. Five Day Couple)
程序员文章站
2024-03-03 17:44:46
...
题意:数组a[1……n], 给m个区间询问:问a[l……r]中一个数异或b的最大值。数值范围1~1e9。
遇到’在一组数中找一个数与数b的异或值最大’这类问题,一般都是将数组建01字典树,因为高位与b对应的数位互异,则异或值大于所有低位之和,这是最优策略:保证高位为1.这道题因为限定了数组的区间,我们需要将思路转化到‘前缀和’上,当前位上与b对应的数位互异的数字的个数,sum[R]-sum[L-1],判断是否为0!
构建可持续化01Tried
策略:我们将组建n+1棵01字典树,插入a[i]在第T[i-1]棵树的基础上建立第T[i]棵树,为什么要在T[i-1]棵树上?因为我们的第T[i]棵01树维护的插入前i个数字的结果,所以我们在查询时仅需要同时遍历T[l-1]、T[r]上的路径!
/***********************************
*** 静态主席树 (不带修改操作)
*** 类比划分树
*** 建树 过程
*** 询问[l,r] 中第k小数
***********************************/
#include <iostream>
#include <map>
#include <cstdio>
#include <string.h>
#include <algorithm>
using namespace std;
typedef long long llt;
const int N = 1e5+77;
const int dep = 29;//最多层数
int a[N];
int L[N*64],R[N*64],sum[N*64];
int T[N];
int toUsed;
void init(){
toUsed = 0;
T[0] = ++toUsed;
L[T[0]] = R[T[0]] = sum[T[0]] = 0;//将根0的左右儿子全标为0!
L[0] = R[0] = sum[0] =0;
}
void Insert(int x,int pos,int &rt,int pre){
rt = ++toUsed;
/*pre==0,即为第一个*/
L[rt] = L[pre];
R[rt] = R[pre];
sum[rt] = sum[pre]+1;
if(pos==-1) return ;//即插入完毕!
/*插入右1节点*/
if(x&(1<<pos)) Insert(x,pos-1,R[rt],R[pre]);
else Insert(x,pos-1,L[rt],L[pre]);
//cout<<"rt: "<<rt<<" "<<"L[rt]: "<<L[rt]<<" "<<"R[rt]: "<<R[rt]<<" "<<"sum[rt]: "<<sum[rt]<<endl;
}
int query(int l,int r,int b,int pos){//主席树中编号T[L-1],T[R]
if(pos==-1) return 0;
int ans = 0;
//cout<<"l: "<<l<<" "<<"r: "<<r<<endl;
if(b&(1<<pos)){//x的pos位为1时,
int num = sum[L[r]]-sum[L[l]];//区间内有num个pos位为0的数在
if(num<=0) {
ans += 0;
ans += query(R[l],R[r],b,pos-1);
}
else {
ans += (1<<pos);
ans += query(L[l],L[r],b,pos-1);
}
}else{//x的pos位为0时,
int num = sum[R[r]]-sum[R[l]];//区间内有num个pos位为1的数
if(num<=0){
ans += 0;
ans += query(L[l],L[r],b,pos-1);
}else{
ans += (1<<pos);
ans += query(R[l],R[r],b,pos-1);
}
}
return ans;
}
int main(){
int n;
scanf("%d",&n);
init();
for(int i=1;i<=n;++i) {
scanf("%d",&a[i]);
Insert(a[i],dep,T[i],T[i-1]);
}
int m ;
scanf("%d",&m);
while(m--){
int b,l,r;
scanf("%d%d%d",&b,&l,&r);
printf("%d\n",query(T[l-1],T[r],b,dep));
}
return 0;
}