POJ 3252 Round Numbersm

题意：一个整数的二进制表示中，如果0的个数比1的相等或更多，则称这个整数为round number. 现在给出范围i，j 要求输出[i,j]之间的round number 的个数。

i,j在2,000,000,000范围内。

规模太大，离线打表不太现实，不过思路并不难想其实。。我们计算m及其之前的round number的个数，记为cal(m)，那么输入m和n时，答案就是cal(m)-cal(n-1)

cal(m)的计算方式是，先得到m到底有几个二进制位，并同时记下有几个0有几个1（这个在后面会用到），假设有x位，那么一切少于x位的二进制数转成十进制一定是小于m的，那么对于y位二进制数（y<x），符合要求的round number 的个数即C(y,y/2)+C(y,y/2+1)+……+C(y,y)，即从二进制位中抽y/2个，y/2+1个，...y个分别作为数字0的位置，得到的数自然是小于m的round number.这样求组合数之和即可。

另外还需要考虑同样有x个二进制位，小于m，且也是个round number的数，这些稍微不好考虑一些，从高位开始考虑，碰到一个0，就将0的数目加一，碰到一个1，就暂时将其抹成0（这确保得到的新数缺失小于原数），同时0的个数加一，这时再考虑这一位及其之后0的个数多于1 的情况数目。。同样是组合数，做完这部分处理之后再把之前暂时抹成0的那一位抹回成1。

组合数是用预先打表的，这个用c[i][j]=c[i-1][j-1]+c[i-1][j] 递推即可，非常简单。

code:

#include<iostream>
#include<cstdio>
#include<cmath>
#include<cstring>
#include<string>
#include<stack>
#include<queue>
#include<algorithm>
#include<cstdlib>

#define maxn 201314
#define inf 0x3f3f3f3f
#define LL long long

using namespace std;
int s,e;
int ans=0;
int bit[33];
int c[33][33];

void init()
{
    for(int i=0; i<33; i++)c[i][0]=1;
    for(int i=1; i<33; i++)
    {
        for(int j=1; j<=i; j++)
        {
            c[i][j]=c[i-1][j-1]+c[i-1][j];
        }
    }
}
int cal(int n)//calculate the ans under or equals n
{
    if(n<=1)return 0;
    int len=0;
    int num1,num0;
    num1=num0=0;
    while(n!=0)//from low to high,get bits
    {
        if(n&1)
        {
            bit[len++]=1;
            num1++;
        }
        else
        {
            bit[len++]=0;
            num0++;
        }
        n>>=1;
    }
    //for(int i=0;i<len;i++)cout<<bit[i];
    //cout<<endl;
    int ans=0;
    for(int i=1; i<len; i++) //the ans of the numbers which length<len
    {
        if(i%2==1)
        {
            int tmp=(i+1)/2;
            for(int k=tmp; k<=i-1; k++)
            {
                ans+=c[i-1][k];
            }
        }
        else
        {
            int tmp=(i/2);
            for(int k=tmp; k<=i-1; k++)
            {
                ans+=c[i-1][k];
            }
        }

    }
    if(num0>=num1)ans++;
    num0=0;
    num1=1;
    for(int i=len-2;i>=0;i--)
    {
        if(bit[i]==1)
        {
            for(int j=i;j>=0;j--)
            {
                if(num0+j+1>=num1+i-j)
                    ans+=c[i][j];
            }
            num1++;
        }
        else num0++;
    }
    return ans;

}
int main()
{
    init();
    while(scanf("%d%d",&s,&e)!=EOF)
    {
        printf("%d\n",cal(e)-cal(s-1));
    }
}