不止代码：区间圆数（ybtoj-数位DP）

最新推荐文章于 2022-02-21 15:32:45 发布

wind__whisper

最新推荐文章于 2022-02-21 15:32:45 发布

阅读量171

点赞数 1

分类专栏：动态规划文章标签：数位DP c++

本文链接：https://blog.csdn.net/BUG_Creater_jie/article/details/116209256

版权

动态规划专栏收录该内容

61 篇文章 2 订阅

订阅专栏

本文探讨了一道编程题目，该题目要求统计特定范围内二进制数中1的个数。作者通过动态规划的方法，用dp[pos][s0]表示二进制数pos位有s0个0的数的个数，并给出详细转移方程。代码实现中，首先初始化dp数组，然后计算不同位置和0的个数对应的二进制数个数，最后解决边界条件并返回结果。文章适合对动态规划和二进制计数感兴趣的读者。

摘要由CSDN通过智能技术生成

题目描述

在这里插入图片描述

解析

一寸山河一寸血

理解万岁！

首先，这题统计[l,r]的个数，可以用[1,r]-[1,l-1]来实现

接下来就是如何统计出[1,n]的个数了
首先，用dp[pos][s0]来表示一个二进制pos位有s0个0的数的个数
我们允许前导零，那么转移就非常简单：

dp[0][0]=1;
dp[pos][num]=dp[pos-1][num]+dp[pos-1][num-1];

现在就要考虑如何求了
因为有前导零，所以要分为两部分来求
具体实现看代码吧

代码

#include<bits/stdc++.h>
#define ll long long
#define mem(a,b) memset(a,b,sizeof(a))
using namespace std;
const int N=600;
const int M=2e6+100;
int m,n;
ll dp[50][50];
void Dp(){
	dp[0][0]=1;
	for(int pos=1;pos<=31;pos++){
		dp[pos][0]=1;
		for(int num=1;num<=pos;num++){
			dp[pos][num]=dp[pos-1][num]+dp[pos-1][num-1];
		}
	}
}

int solve(int n){
	if(n==1||n==0) return 0;
	int ans=0,s[50],x=n,s0=0,s1=0,tot=0;
	while(x){
		s[++tot]=x&1;
		x /= 2;
	}
	for(int pos=tot;pos>=1;pos--){
		//最高位特判不要算 
		if(s[pos]&&tot!=pos){//有一的话：前面与n相同，pos填1 
			s0++;
			for(int i=pos-1;i>=0&&s0+i>=tot-s0-i;i--) ans+=dp[pos-1][i];
			s0--;
		}
		if(tot!=pos){//以pos作为最高位，填 1 
			for(int i=pos-1;i>=0&&i>=pos-i;i--) ans+=dp[pos-1][i];
		}
		s1+=s[pos];
		s0+=!s[pos];
	}
	if(s0>=s1) ans++;
	return ans;
}
int a,b;
int main(){
	Dp();
	scanf("%d%d",&a,&b);
	printf("%d\n",solve(b)-solve(a-1));
	return 0;
}
/*
13
100
200
1000
*/