「SCOI2016」美味
题目描述
一家餐厅有 (n) 道菜,编号 (1 ldots n) ,大家对第 (i) 道菜的评价值为 (a_i :( 1 leq i leq n ))。有 (m) 位顾客,第 (i) 位顾客的期望值为 (b_i),而他的偏好值为 (x_i)。因此,第 (i) 位顾客认为第 (j) 道菜的美味度为 (b_i ext{xor} (a_j+x_i)) (( ext{xor}) 表示异或运算)。 第 (i) 位顾客希望从这些菜中挑出他认为最美味的菜,即美味值最大的菜,但由于价格等因素,他只能从第 (l_i) 道到第 (r_i) 道中选择。请你帮助他们找出最美味的菜。
解题思路 :
问题要求从区间里找一个数 (a_j) ,使得 ((a_j + x_i) ext{ xor } b_i) 最大,先考虑当 (x_i = 0) 的时候是一个经典问题该怎么做
当 (x_i = 0) 时显然直接用一个可持久化 (Trie) 树按位维护区间内所有 (a_j) 查询即可
但是有加上 (x_i) 的操作由于 (Trie) 的局限性就不那么可做了,不妨考虑上述算法的本质
(Trie) 树从高到低走的每一位本质上是对最终选的数的可行区间不断减少
设当前第 (i) 位的大小是 (2^i) ,那么如果选了它答案一定 $ geq 2^i$ ,否则答案一定 (< 2 ^i) ,因为有 (sum_{j=0}^{i-1} 2^j =2^i - 1)
所以我们可以用一个线段树维护 ([0, 2^{lim-1}]) 的权值来代替原本的 (Trie),每次走左儿子本质上是这一位选了 (0),走右儿子本质上是选了 (1)
考虑查询加上 (x_i) 后的 (a_j) 本质上是将原来的 (a_j - x_i) 代替新的 (a_j) ,只需要查询的区间更改一下就可以了
因为要回答区间的答案,所以还要把线段树可持久化一下,每次查询从高到低枚举 (b_i) 的每一位看看能否选相反的即可
/*program by mangoyang*/
#include<bits/stdc++.h>
#define inf (0x7f7f7f7f)
#define Max(a, b) ((a) > (b) ? (a) : (b))
#define Min(a, b) ((a) < (b) ? (a) : (b))
typedef long long ll;
using namespace std;
template <class T>
inline void read(T &x){
int f = 0, ch = 0; x = 0;
for(; !isdigit(ch); ch = getchar()) if(ch == '-') f = 1;
for(; isdigit(ch); ch = getchar()) x = x * 10 + ch - 48;
if(f) x = -x;
}
const int N = 1000005, len = (1 << 19) - 1;
int rt[N], n, m;
struct SegmentTree{
int sz[N*25], lc[N*25], rc[N*25], size;
inline void ins(int &u, int pr, int l, int r, int pos){
u = ++size;
lc[u] = lc[pr], rc[u] = rc[pr], sz[u] = sz[pr] + 1;
if(l == r) return;
int mid = l + r >> 1;
if(pos <= mid) ins(lc[u], lc[pr], l, mid, pos);
else ins(rc[u], rc[pr], mid + 1, r, pos);
}
inline int query(int x, int y, int l, int r, int L, int R){
if(l >= L && r <= R) return (sz[y] - sz[x] > 0);
int mid = l + r >> 1, res = 0;
if(L <= mid) res |= query(lc[x], lc[y], l, mid, L, R);
if(mid < R) res |= query(rc[x], rc[y], mid + 1, r, L, R);
return res;
}
}van;
int main(){
read(n), read(m);
for(int i = 1, x; i <= n; i++)
read(x), van.ins(rt[i], rt[i-1], 0, len, x);
while(m--){
int now, add, l = 0, r = len, res = 0, L, R;
read(now), read(add), read(L), read(R);
for(int i = 18; ~i; i--){
int s = ((now >> i) & 1) ^ 1, mid = l + r >> 1;
int nl = s ? mid + 1 : l, nr = s ? r : mid;
if(van.query(rt[L-1], rt[R], 0, len, Max(nl - add, 0), Max(nr - add, 0)))
res |= 1 << i, l = nl, r = nr;
else l = s ? l : mid + 1, r = s ? mid : r;
}
printf("%d
", res);
}
return 0;
}