BZOJ 1662: [Usaco2006 Nov]Round Numbers 圆环数(数位DP+恶心细节)

BZOJ 1662: [Usaco2006 Nov]Round Numbers 圆环数

Time Limit: 5 Sec  Memory Limit: 64 MB

Description

正如你所知,奶牛们没有手指以至于不能玩“石头剪刀布”来任意地决定例如谁先挤奶的顺序。她们甚至也不能通过仍硬币的方式。 所以她们通过"round number"竞赛的方式。第一头牛选取一个整数,小于20亿。第二头牛也这样选取一个整数。如果这两个数都是 "round numbers",那么第一头牛获胜,否则第二头牛获胜。 如果一个正整数N的二进制表示中,0的个数大于或等于1的个数,那么N就被称为 "round number" 。例如,整数9,二进制表示是1001,1001 有两个'0'和两个'1'; 因此,9是一个round number。26 的二进制表示是 11010 ; 由于它有2个'0'和 3个'1',所以它不是round number。 很明显,奶牛们会花费很大精力去转换进制,从而确定谁是胜者。 Bessie 想要作弊,而且认为只要她能够知道在一个指定区间范围内的"round numbers"个数。 帮助她写一个程序,能够告诉她在一个闭区间中有多少round numbers。区间是 [start, finish],包含这两个数。 (1 <= Start < Finish <= 2,000,000,000)

Input

* Line 1: 两个用空格分开的整数,分别表示Start 和 Finish。

Output

* Line 1: Start..Finish范围内round numbers的个数

Sample Input

2 12

Sample Output

6

题解
其实我真的懒得动写题解。。。
不过这题网上做法千奇百怪。3维~5维DP都有,还有用记忆化搜索的。我实在看不下去了。这题其实两维DP+几个循环就足够了!
就是个恶心版数位DP!!!
我们设dp[i][j]为转成二进制后 前导可以为0 长度为i的串,其中零的个数为j的二进制串总个数。
很显然得到状态转移方程:dp[i][j]=dp[i-1][j]+dp[i-1][j-1]
dp数组可以预处理得到。。。
sum(i,j)函数可以累计DP数组求出长度为i且零比一至少多j个的二进制串总个数
那么如何计算ST-ED范围内的答案?
用1...ED的答案减去1...ST-1的答案即可。(solve函数传入x即求1...二进制_x的答案)
比如计算1..110100的答案:
其中几部分是1-1,10-11,100-111,...,10000-11111
最后一部分是求100000-110100的答案
可以拆成100000-101111和110000-110011分别求解,这样就很好解决了,怎么求解还是很明显的。
由此可见求100000-110100的答案之类最后一部分的基本思想是:
从左往右扫到第i位时 当前1比0多d个 如果是此位为0跳过 如果是1则累计ans+=sum(i-1,d-1)。
大体上就是这样,不过细节恶心!
首先,若读入的数为0直接跳过,不然会无限RE。。。
其次,若读入的数转成二进制,0的个数>=1的个数则最后ans++,说明这个数本身合法。
另外还有一些起始点,边界等细节问题。
下面给出参考代码:(其实还是很短的2333)
#include<bits/stdc++.h>
using namespace std;

int n,m; 
int dp[35][35];

int sum(int l,int x)
{
	int s=0;
	for(int i=0;i<=l;i++)if(2*i>=l+x)s+=dp[l][i];
	return s;
}

int solve(int x)
{
	int len=0,ans=0,d=1,t=x,v[35];
	while(t)
	{
		v[++len]=t&1;
		t>>=1;
	}
	for(int i=1;i<len;i++)ans+=sum(i-1,1);
	len--;
	while(len)
	{
		if(v[len])ans+=sum(len-1,d-1),d++; else d--;
		len--;
	}
	if(d<=0)ans++;
	return ans; 
}

int main()
{
	for(int i=1;i<=32;i++){dp[i-1][0]=1;for(int j=1;j<=i;j++)dp[i][j]+=dp[i-1][j]+dp[i-1][j-1];}
	scanf("%d%d",&n,&m);
	if(n>2)printf("%d\n",solve(m)-solve(n-1)); else printf("%d\n",solve(m));
	return 0;
}

  

 
posted @ 2017-08-28 00:35  winmt  阅读(259)  评论(0编辑  收藏  举报
当前时间: