[CF1245F]Daniel and Spring Cleaning

最新推荐文章于 2023-08-31 19:11:51 发布

StaroForgin

最新推荐文章于 2023-08-31 19:11:51 发布

阅读量120

点赞数 1

分类专栏： ------动态规划------ # 数位dp

本文链接：https://blog.csdn.net/Tan_tan_tann/article/details/107935504

版权

------动态规划------ 同时被 2 个专栏收录

174 篇文章 0 订阅

订阅专栏

数位dp

11 篇文章 0 订阅

订阅专栏

Daniel and Spring Cleaning

题解

很容易发现，原式等价于 $\sum_{a=l}^{r}\sum_{b=l}^{r}[a \& b= 0]$ 。

再差分一下， $\sum^{r}\sum^{r}[a\&b =0]- 2\sum^{l-1}\sum^{r}[a\& b=0]+\sum^{l-1}\sum^{l-1}[a\& b=0]$ 。

令 $calc(x,y)= \sum^{x}\sum^{y}[a\& b= 0]$ ，接下来，来考虑一下 calc 函数该怎么求。

由于是二进制，很容易想到数位dp，令 $dp_{i,j,k}$ 为在转化为二进制后，的第位为，的第位为时的满足条件的方案数。

dp的转换也十分好想，判断是否达到上界即可。

源码

#include<bits/stdc++.h>
using namespace std;
typedef long long LL;
typedef pair<int,int> pii;
int t,l,r;LL dp[40][2][2];
LL calc(LL a,LL b){
	if(a<0||b<0)return 0;
	memset(dp,0,sizeof(dp));dp[32][1][1]=1;
	for(int i=31;i>=0;i--)
		for(int j=0;j<2;j++)
			for(int k=0;k<2;k++){
				bool nxtj=(j&&!((a>>i)&1LL)),nxtk=(k&&!((b>>i)&1LL));
				dp[i][nxtj][nxtk]+=dp[i+1][j][k];
				if(((a>>i)&1)||!j)dp[i][j][nxtk]+=dp[i+1][j][k];
				if(((b>>i)&1)||!k)dp[i][nxtj][k]+=dp[i+1][j][k];
			}
	return dp[0][0][0]+dp[0][1][0]+dp[0][0][1]+dp[0][1][1];
}
signed main(){
	scanf("%d",&t);
	while(t--){
		scanf("%d %d",&l,&r);l--;
		printf("%lld\n",calc(r,r)-calc(l,r)*2ll+calc(l,l));
	}
	return 0;
}