bzoj4571&&SCOI2016美味

Posted

tags:

篇首语:本文由小常识网(cha138.com)小编为大家整理,主要介绍了bzoj4571&&SCOI2016美味相关的知识,希望对你有一定的参考价值。

4571: [Scoi2016]美味

Time Limit: 30 Sec  Memory Limit: 256 MB
Submit: 656  Solved: 350
[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

 

三步走

●维护xor最大值的话用trie树可以搞一搞

●维护区间xor最大值可以用可持久化trie树

●维护区间xor并支持加减的就是可持久化权值线段树了,这就是本题的算法,用贪心+主席树实现

 

 

贪心?肯定先确定最高位使得xor后此位为1,由此得到

策略1:从二进制高位向低位贪心,尽量满足 b的i位与(a[j]+x)的i位不同

 

 

怎么判断第i位是否可以不同?可以假设(a[j]+x)的i位与b的i位不同,来确定所选的数的第i位,再判定假设是否合法。

合法的条件是什么?记已确定的数为ans,只用查找区间内是否存在数这样的满足条件

1.前面位和ans相同

2.此位满足已确定的数

 

不太容易懂,举一个例子:

假设b的当前位为1,我们应尽量使a[j]+x当前位为0,前面位确定为ans,ans第i位设定为0,查找区间内是否存在 在[ans,ans+(1<<i)-1]之内的数,如果有,则此位可以为0

 

 

技术分享

判断区间的数的个数,权值线段树

但此时我们注意到查找是双关键字的,一个是区间一个是权值,所以用主席树来完成操作

总结得出

策略2:用主席树判断区间数的存在性问题来确定当前位究竟是什么

 

到了这里,此题也就差不多了,但有几个细节需要注意

●主席树空间开够(我就是这里错了好久--为什么dev不报RE而是乱搞数组啊!!)

●由于我们查询的数是a[j]+x,而树中插入的是a[],所以查询区间统一减去x并考虑区间是否超界

●最后注意,我们的到的ans是选出来的数,所以输出时要xor b

 

没有了-------------------------------------------------------------------

 

 

 

#include<cstdio>
#include<algorithm>
#include<iostream>
#include<cstring>
#define ll long long
#define N 200005
#define mx (1<<18)-1
using namespace std;
int n,m,a[N],rt[N],ls[N<<5],rs[N<<5],sum[N<<5],sz;
void insert(int pre,int &u,int l,int r,int pos){
	u=++sz;ls[u]=ls[pre];rs[u]=rs[pre];sum[u]=sum[pre]+1;
	if(l==r)return;int mid=(l+r)>>1;
	if(pos<=mid)insert(ls[pre],ls[u],l,mid,pos);
	else insert(rs[pre],rs[u],mid+1,r,pos);
}

int query(int pre,int u,int l,int r,int L,int R){
	if(L<=l&&r<=R)return sum[u]-sum[pre];
	int mid=(l+r)>>1;int t=0;
	if(L<=mid)t+=query(ls[pre],ls[u],l,mid,L,R);
	if(R>mid)t+=query(rs[pre],rs[u],mid+1,r,L,R);
	return t;
}

int main(){
	scanf("%d%d",&n,&m);
	for(int i=1;i<=n;i++)scanf("%d",&a[i]);
	for(int i=1;i<=n;i++)
	insert(rt[i-1],rt[i],0,mx,a[i]);
	for(int i=1;i<=m;i++){
		int b,x,l,r;
		scanf("%d%d%d%d",&b,&x,&l,&r);
		int ans=0;
		for(int j=17;j>=0;j--){
			if(b&(1<<j)){
				int L=max(ans-x,0),R=ans+(1<<j)-x-1;
				if(R<0||!query(rt[l-1],rt[r],0,mx,L,R))ans^=(1<<j);
			}
			else{
				ans^=(1<<j);
				int L=max(ans-x,0),R=ans+(1<<j)-x-1;
				if(R<0||!query(rt[l-1],rt[r],0,mx,L,R))ans^=(1<<j);
			}
		}
		printf("%d\\n",ans^b);
	}
	return 0;
}

 

  

 

以上是关于bzoj4571&&SCOI2016美味的主要内容,如果未能解决你的问题,请参考以下文章

[BZOJ4571][Scoi2016]美味 贪心+主席树

BZOJ4571[Scoi2016]美味 主席树

bzoj 4571 [Scoi2016]美味——主席树

BZOJ4571: [Scoi2016]美味

bzoj4571 [Scoi2016]美味

[BZOJ 4571][Scoi2016]美味(主席树)