POJ 3252 Round Numbers （组合数，二进制）

题意 ：
如果一个数化成二进制后零的个数比一的个数多(不算前导零)，称这个数为Round Number，
问[a,b]内有多少个Round Number。(1<=a<=b<=$2\times10^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;
}