[Codeforces1245F]Daniel and Spring Cleaning

本文探讨了在特定区间内寻找满足x+y=x⊕y条件的(x,y)对的数量,通过数位DP算法实现高效求解,适用于0≤l≤r≤10^9的大范围限制。

摘要生成于 C知道 ,由 DeepSeek-R1 满血版支持, 前往体验 >

题意

[l,r][l,r][l,r]之间有多少对(x,y)(x,y)(x,y)满足x+y=x⊕yx+y=x\oplus yx+y=xy

0≤l≤r≤1090\le l\le r\le10^90lr109


题解

注意到l,rl,rl,r的范围,又涉及到二进制位上运算,所以考虑按二进制来数位dpdpdp

f(l,r)f(l,r)f(l,r)表示满足以下条件的(x,y)(x,y)(x,y)有多少对

x+y=x⊕yx+y=x\oplus yx+y=xy

0≤x≤l0\le x\le l0xl

0≤y≤r0\le y\le r0yr

根据容斥可以得到Ans=f(r,r)−2f(l−1,r)+f(l−1,l−1)Ans=f(r,r)-2f(l-1,r)+f(l-1,l-1)Ans=f(r,r)2f(l1,r)+f(l1,l1)

根据数位dpdpdp的基本思路

dp(p,Limx,Limy)dp(p,Lim_x,Lim_y)dp(p,Limx,Limy)表示在二进制第ppp位,x,yx,yx,y是否有限制的答案

有限制表示已经选的填的数的高位是否与l,rl,rl,r相同

如果有限制这一位填的数就不能超过l,rl,rl,r在当前位的值

时间复杂度O(log⁡r)O(\log r)O(logr)

#include<bits/stdc++.h>
using namespace std;
using ll=long long;
int L,R;ll f[33][2][2];
ll dp(int p,bool Lim_x,bool Lim_y){
	if(p==-1)return 1;
	ll&g=f[p][Lim_x][Lim_y];
	if(g!=-1)return g;
	g=0;
	int Up_x=Lim_x?(L>>p)&1:1,
		Up_y=Lim_y?(R>>p)&1:1;
	for(int i=0;i<=Up_x;++i)
		for(int j=0;j<=Up_y;++j)
			if(!(i&j))
				g+=dp(p-1,Lim_x&&i==Up_x,Lim_y&&j==Up_y);
	return g;
}
inline ll Sol(int l,int r){
	if(l<0)return 0;
	memset(f,-1,sizeof f);
	L=l,R=r;
	return dp(log2(R+1)+1,1,1);
}
int main(){
	int t,l,r;
	scanf("%d",&t);
	while(t--)
		scanf("%d%d",&l,&r),
		printf("%lld\n",Sol(r,r)-2*Sol(l-1,r)+Sol(l-1,l-1));
return 0;
}
评论 3
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值