4571: [Scoi2016]美味

4571: [Scoi2016]美味

Time Limit: 30 Sec   Memory Limit: 256 MB
Submit: 374   Solved: 192
[ Submit][ Status][ Discuss]

Description

一家餐厅有 n 道菜,编号 1...n ,大家对第 i 道菜的评价值为 ai(1≤i≤n)。有 m 位顾客,第 i 位顾客的期
望值为 bi,而他的偏好值为 xi 。因此,第 i 位顾客认为第 j 道菜的美味度为 bi XOR (aj+xi),XOR 表示异或
运算。第 i 位顾客希望从这些菜中挑出他认为最美味的菜,即美味值最大的菜,但由于价格等因素,他只能从第 
li 道到第 ri 道中选择。请你帮助他们找出最美味的菜。

Input

第1行,两个整数,n,m,表示菜品数和顾客数。
第2行,n个整数,a1,a2,...,an,表示每道菜的评价值。
第3至m+2行,每行4个整数,b,x,l,r,表示该位顾客的期望值,偏好值,和可以选择菜品区间。
1≤n≤2×10^5,0≤ai,bi,xi<10^5,1≤li≤ri≤n(1≤i≤m);1≤m≤10^5

Output

 输出 m 行,每行 1 个整数,ymax ,表示该位顾客选择的最美味的菜的美味值。

Sample Input

4 4
1 2 3 4
1 4 1 4
2 3 2 3
3 2 3 3
4 1 2 4

Sample Output

9
7
6
7

HINT

Source

[ Submit][ Status][ Discuss]



多次查询区间内的数加上一个数以后xor一个数的最大值。。。。。。。。。。

考虑按高位到低位贪心,什么时候该位的答案能为1?

首先要满足前面高位的限制,然后调整该位取的数,使得答案的这一位为1

这样后面的区间有个取值范围(000000...~11111...),满足这个范围的数字都是可行的

询问要求先加上一个数字,那就把这个区间左右端点都减去这个数字

那问题就可以转换成询问一个区间内对于某个值域,是否存在属于这个值域的数字

主席树就能解决了。

#include<iostream>
#include<cstring>
#include<vector>
#include<queue>
#include<algorithm>
#include<cmath>
#include<cstdio>
#include<bitset>
using namespace std;

const int maxn = 2E5 + 20;
const int T = 20;

int n,m,cur = 1,cnt,a[maxn],A[maxn],rt[maxn],lc[maxn*T],rc[maxn*T],c[maxn*T];

bool query(int o1,int o2,int l,int r,int ql,int qr)
{
	if (ql <= l && r <= qr) return c[o2] - c[o1] > 0;
	int mid = (l + r) >> 1; bool ret = 0;
	if (ql <= mid) ret = query(lc[o1],lc[o2],l,mid,ql,qr);
	if (ret) return ret;
	if (qr > mid) ret = query(rc[o1],rc[o2],mid+1,r,ql,qr);
	return ret;
}

bool Query(int l,int r,int A,int B)
{
	if (B < a[1] || a[n] < A) return 0;
	int pa = lower_bound(a + 1,a + cur + 1,A) - a;
	int pb = lower_bound(a + 1,a + cur + 1,B) - a;
	if (a[pb] > B) --pb; if (pa > pb) return 0;
	return query(rt[l-1],rt[r],1,cur,pa,pb);
}

int Insert(int o,int l,int r,int pos)
{
	int ret = ++cnt;
	c[ret] = c[o] + 1;
	if (l == r) return ret;
	int mid = (l + r) >> 1;
	if (pos <= mid) rc[ret] = rc[o],lc[ret] = Insert(lc[o],l,mid,pos);
	else lc[ret] = lc[o],rc[ret] = Insert(rc[o],mid+1,r,pos);
	return ret;
}

int getint()
{
	char ch = getchar(); int ret = 0;
	while (ch < '0' || '9' < ch) ch = getchar();
	while ('0' <= ch && ch <= '9')
		ret = ret*10 + ch - '0',ch = getchar();
	return ret;
}

int main()
{
	#ifdef DMC
		freopen("DMC.txt","r",stdin);
	#endif
	
	n = getint(); m = getint();
	for (int i = 1; i <= n; i++)
		A[i] = getint(),a[i] = A[i];
	sort(a + 1,a + n + 1);
	for (int i = 2; i <= n; i++)
		if (a[i] != a[i-1]) a[++cur] = a[i];
	a[++cur] = ~0U>>1;
	for (int i = 1; i <= n; i++)
	{
		A[i] = lower_bound(a + 1,a + cur + 1,A[i]) - a;
		rt[i] = Insert(rt[i-1],1,cur,A[i]);
	}
	
	while (m--)
	{
		int b,x,l,r,Ans = 0;
		b = getint(); x = getint();
		l = getint(); r = getint();
		for (int i = 20; i >= 0; i--)
		{
			int k = Ans^b>>i<<i; k ^= 1<<i;
			if (Query(l,r,k-x,(k|(1<<i)-1)-x)) Ans |= 1<<i;
		}
		printf("%d\n",Ans);
	}
	return 0;
}

  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值