[SCOI2016]美味-题解

题意见题面。


首先,考虑异或最大可以按照二进制位重高位到低位贪心,所以我们对于一个询问的 b i b_i bi,假如当前答案要确定的这一位上 b i b_i bi 0 0 0,那么我们只需要权值在 ( ⋯ 10000 ) 2 ∼ ( ⋯ 11111 ) 2 (\cdots10000)_2\sim(\cdots11111)_2 (10000)2(11111)2的菜品(也就是这一位是 1 1 1的范围,省略号表示前面已经确定了的),不断的修改左右范围,一位一位确定即可,复杂度为 l o g v logv logv的。

由于询问还要加上 x i x_i xi,所以我们反过来思考,那么真正的查询区间就是原来的 L ∼ R L\sim R LR变成了 L − x ∼ R − x L-x\sim R-x LxRx(注意和0取max)。

所以每一位我们只需查询在当前的可选取的菜品内有没有 a i a_i ai在这个范围内的,这个可以用主席树在 l o g n logn logn的时间内实现。

那么总复杂度就是 O ( n l o g n + m l o g 2 n ) O(nlogn+mlog^2n) O(nlogn+mlog2n),那么这个题就解决了。

#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;
const int M=2e5+10,Log=20;
const int UP=2e5+1;
int sum[M*Log],ls[M*Log],rs[M*Log],tot,root[M];
int n,m;
void update(int pre,int &o,int l,int r,int p){
	o=++tot;ls[o]=ls[pre];rs[o]=rs[pre];sum[o]=sum[pre];
	if(l==r){++sum[o];return;}
	int mid=l+r>>1;
	if(p<=mid) update(ls[pre],ls[o],l,mid,p);
	else update(rs[pre],rs[o],mid+1,r,p);
	sum[o]=sum[ls[o]]+sum[rs[o]];
}
int query(int o,int l,int r,int L,int R){
	if(L<=l&&r<=R) return sum[o];
	int mid=l+r>>1;
	if(L>mid) return query(rs[o],mid+1,r,L,R);
	if(R<=mid) return query(ls[o],l,mid,L,R);
	return query(ls[o],l,mid,L,R)+query(rs[o],mid+1,r,L,R);
}
int query_area(int l,int r,int vl,int vr){
	vr=min(vr,UP);vl=min(vl,UP);
	return query(root[r],0,UP,vl,vr)-query(root[l-1],0,UP,vl,vr);
}
int b,x,l,r;
int main(){
	scanf("%d%d",&n,&m);
	for(int i=1;i<=n;i++){
		scanf("%d",&x);
		update(root[i-1],root[i],0,UP,x);
	}
	while(m--){
		scanf("%d%d%d%d",&b,&x,&l,&r);
		int ans=0;
		for(int i=19,low=0,up=(1<<20)-1;i>=0;i--){
			if((b>>i)&1){
				up^=(1<<i);
				if(!query_area(l,r,max(0,low-x),max(0,up-x)))up^=(1<<i),low^=(1<<i);
				else ans|=(1<<i);
			}else{
				low^=(1<<i);
				if(!query_area(l,r,max(0,low-x),max(0,up-x)))low^=(1<<i),up^=(1<<i);
				else ans|=(1<<i);
			}
		}
		printf("%d\n",ans);
	}
	return 0;
}
### 关于 SCOI2009 WINDY 数的解法 #### 定义与问题描述 WINDY数是指对于任意两个相邻位置上的数字,它们之间的差至少为\(2\)。给定正整数区间\([L, R]\),计算该范围内有多少个WINDY数。 #### 动态规划方法解析 为了高效解决这个问题,可以采用动态规划的方法来处理。定义状态`dp[i][j]`表示长度为`i`且最高位是`j`的WINDY数的数量[^3]。 - **初始化** 对于单个数字的情况(即只有一位),显然每一位都可以单独构成一个合法的WINDY数,因此有: ```cpp dp[1][d] = 1; // d ∈ {0, 1,...,9} ``` - **状态转移方程** 当考虑多位数时,如果当前位选择了某个特定数值,则下一位的选择会受到限制——它必须满足与前一位相差不小于2的要求。具体来说就是当上一高位为`pre`时,当前位置可选范围取决于`pre`的具体取值: - 如果`pre >= 2`, 则可以选择`{0... pre-2}` - 否则只能从剩余的有效集合中选取 这样就可以通过遍历所有可能的状态来进行状态间的转换并累加结果。 - **边界条件处理** 特殊情况下需要注意的是,在实际应用过程中还需要考虑到给出区间的上下限约束。可以通过逐位比较的方式判断是否越界,并据此调整有效状态空间大小。 ```cpp // 计算不超过num的最大windy数数量 int calc(int num){ int f[15], g[15]; memset(f, 0, sizeof(f)); string s = to_string(num); n = s.size(); for (char c : s) { a[++len] = c - '0'; } // 初始化f数组 for (int i=0;i<=9;++i)f[1][i]=1; // DP过程省略... return sum; } long long solve(long long L,long long R){ return calc(R)-calc(L-1); } ``` 此代码片段展示了如何利用预处理好的`dp`表快速查询指定范围内的WINDY数总量。其中`solve()`函数用于返回闭区间\[L,R\]内符合条件的总数;而辅助函数`calc()`负责根据传入参数构建相应的状态序列并最终得出答案。
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包

打赏作者

VictoryCzt

你的鼓励将是我创作的最大动力

¥1 ¥2 ¥4 ¥6 ¥10 ¥20
扫码支付:¥1
获取中
扫码支付

您的余额不足,请更换扫码支付或充值

打赏作者

实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

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

余额充值