POJ3252 Round Numbers

组合数学的一道题，题意是求区间内二进制下数‘0’不少于‘1’的个数。

可以看出转化为二进制时，只要考虑后面几位中取几个‘0’符合要求就行。。。这样的话就是组合数的求解，可以发现一个规律：不同长度的二进制数位要取的个数为杨辉三角的一半，而且当杨辉三角该项项数为奇数时，中间那位数不取。。将所有情况保存在二维数组里取用就行。

#include<iostream>
#include<cstdio>
#include<string>
#include<cstring>
#include<cmath>
using namespace std;
long long sum,num,ans;
int a[33],b[33];
int Cn[32][32];
void YH()//构建杨辉三角
{
    int i,j;
    for(i=0;i<32;++i)
    {
        Cn[0][i]=Cn[1][i]=0;
    }
    Cn[0][0]=Cn[1][1]=1;
    for(i=2;i<32;++i)
    {
        Cn[i][0]=0;
        for(j=1;j<32;++j)
        {
            Cn[i][j]=Cn[i-1][j-1]+Cn[i-1][j];
        }
    }
}
int main()
{
    int n,m,i,j,p=0,q=0,k;
    YH();
    scanf("%d%d",&n,&m);
    m++;
    while(m>0)//前后区间分别转化为二进制
    {
        b[q]=m%2;
        m=m/2;
        q++;
    }
    while(n>0)
    {
        a[p]=n%2;
        n=n/2;
        p++;
    }
    sum=0;
    for(i=1;i<p;++i)
    {
        for(j=(i/2);j>0;--j)
        {
            sum+=Cn[i][j];
        }
    }
    k=0;
    for(i=p-2; i>=0; --i)
    {
        if(a[i]==1)//对后面空的位数进行排列组合
        {
            for(k++,j=(p/2)-k;j>=0;--j)
            {
                sum+=Cn[i+1][j+1];
            }
        }
    }
    ans=sum;
    sum=0;
    for(i=1;i<q;++i)
    {
        for(j=(i/2);j>0;--j)
        {
            sum+=Cn[i][j];
        }
    }
    k=0;
    for(i=q-2; i>=0; --i)
    {
        if(b[i]==1)
        {
            for(k++,j=(q/2)-k;j>=0;--j)
            {
                sum+=Cn[i+1][j+1];
            }
        }
    }
    num=sum;
    cout<<num-ans<<endl;
    return 0;
}