题意 :
如果一个数化成二进制后零的个数比一的个数多(不算前导零),称这个数为Round Number,
问[a,b]内有多少个Round Number。(1<=a<=b<=
2×109
2
×
10
9
)
分析:
可以先求[0,a-1]内的Round Number数目,再求[0,b]内的Round Number数目.相减就是答案。
设一个数n化为二进制后存到bin数组内,位数是len;
想要求[0,n]内的Round Number数目,分两步:
1.求位数小于len的数中是Round Number的数目,这个很容易求。
假设一个数x的位数是len_x.(len_x
for(int i=1;i<bin[0]-1;i++){
for(int j=i/2+1;j<=i;j++){
ans+=c[i][j];
}
}
bin[0]就是n的位数(也就是len),第一层循环对i从1到bin[0]-2遍历,求出二进制是i+1位数的满足是Round Number数的数量。之所以是i+1,是因为当是len_x位数时,最高位一定是1,所以只有len_x-1位可以用来组合填剩余的数。这里的i就是len_x-1。
第二层循环是求当位数是i+1的时候,至少要放i/2+1个0,才能满足条件。之多放i个,因为最高位是1。
2.然后就是求位数时len的时候满足条件的数目。
方法是对n从len-1位向低位搜索,遇到0,记录当前0出现的个数。遇到1,把这一位当作零,对后面的位数组合数求解。
举个栗子:
n=11111.这里就是从第二个1开始向最后一个1搜索。
搜到第二个1的时候,把它看作0,对后面的三位组合数求解 即求出10000—–10111内满足条件的数目。
搜到第三个1的时候,把它看作0,对后面的两位组合数求解 即求出11000—–11011内满足条件的数目。
搜到第四个1的时候,把它看作0,对后面的一位组合数求解 即求出11100—–11101内满足条件的数目。
搜到第五个1的时候,把它看作0,对后面的零位组合数求解 即求出11110—–11110内满足条件的数目。
合起来正好是10000—11110内满足条件的数目。
再加上第一步求得,就是0—11110内所有满足条件的数目。
即是所有比n小的数中满足条件的数目。
由于题目要求的是[0,n]内满足条件的数目。因此求的时候对n+1求就可以了。
int zero=0;
for(int i=bin[0]-1;i>=1;i--){
if(bin[i]){
for(int j=(bin[0]+1)/2-(zero+1);j<=i-1;j++)
ans+=c[i-1][j];
}
else
zero++;
}
zero是用来统计到当前位为止0出现的次数。
j=(bin[0]+1)/2-(zero+1)是因为到当前为止已经出现过zero个0,又把当前位的这个1当做0,所以是zero+1。
本来是要求至少出现(len+1)/2个0。现在已经出现了zero+1个,
所以后面剩余的位数中,至少出现(bin[0]+1)/2-(zero+1)个0就可以了。
最后一点;组合数用帕斯卡公式打表就可以了
c[i][j]=c[i-1][j-1]+c[i-1][j]
代码:
#include<cstdio>
#include<iostream>
#include<algorithm>
#include<cstring>
#include<cmath>
using namespace std;
#define ll long long
ll a,b;
int bin[123];
ll c[33][33];
void CC()
{
memset(c,0,sizeof(c));
for(int i=0;i<=32;i++){
for(int j=0;j<=i;j++){
if(i==j||!j)
c[i][j]=1;
else
c[i][j]=c[i-1][j-1]+c[i-1][j];
}
}
return;
}
void get_bin(int x)
{
memset(bin,0,sizeof(bin));
while(x){
bin[++bin[0]]=x%2;
x/=2;
}
// for(int i=1;i<=bin[0];i++)
// cout<<bin[i]<<" ";
return;
}
ll solve(ll x)
{
ll ans=0;
get_bin(x);
for(int i=1;i<bin[0]-1;i++){
for(int j=i/2+1;j<=i;j++){
ans+=c[i][j];
}
}
int zero=0;
for(int i=bin[0]-1;i>=1;i--){
if(bin[i]){
for(int j=(bin[0]+1)/2-(zero+1);j<=i-1;j++)
ans+=c[i-1][j];
}
else
zero++;
}
return ans;
}
int main()
{
CC();
while(cin>>a>>b){
ll ans=solve(b+1)-solve(a);
cout<<ans<<endl;
}
return 0;
}