poj 3252 组合数

    主要考察组合数知识，初始化的时候参考公式

首先先推个公式，就是长度为len的Round Numbers的个数。

     长度为len,第一位肯定是1了。

     那么后面剩下 len-1位。

     如果len-1是偶数。

     那么  C（len-1,(len-1)/2+1)+C（len-1,(len-1)/2+2)+````C(len-1,len-1)

=   ( 2^(len-1)-C(len-1,(len-1)/2) )/2;

    如果len是奇数

   那么就是 (  2^(len-1) )/2

#include<cstdio>
#include<iostream>
#include<algorithm>
#include<cstring>
#include<cmath>
#include<queue>
using namespace std;
const int maxn=35;
int m,t;
int c[maxn][maxn];
int bit[maxn];
void init()     //初始化组合数
{
    c[0][0]=c[1][0]=c[1][1]=1;
    for(int i=2;i<=33;i++)
    {
        c[i][0]=1;
        for(int j=1;j<i;j++)
            c[i][j]=c[i-1][j]+c[i-1][j-1];
        c[i][i]=1;
    }
}
int calc(int n) //求小于等于n的round数的和
{
    if(n<=1)    return 0;
    int len=0;
    while(n>0)  //求出该数的二进制表示
    {
        if(n&1) bit[len++]=1;
        else bit[len++]=0;
        n>>=1;
    }
    int ans=0;
    for(int i=len-1;i>0;i--)    //求出长度小于len的round数的和
    {
        if(i%2==0)ans+=((1<<(i-1)))/2;
        else  ans+=((1<<(i-1))-c[i-1][(i-1)/2])/2;
    }
    int c1=0,c0=0;
    for(int i=0;i<len;i++)
    {
        if(bit[i]==0)    c0++;
        else c1++;
    }
    if(c0>=c1)  ans++;      //该数本身是round数
    c1=1;c0=0;
    for(int i=len-2;i>=0;i--)    //注意计数是从零开始的
    {
        if(bit[i]==1)   //该位可以变成零
        {
            for(int j=i;j>=0&&j+c0+1>=i-j+c1;j--)    ans+=c[i][j];   //j是准备添加的0的数目
            c1++;
        }
        else c0++;
    }
    return ans;
}
int main()
{
    int i,j,k;
    init();
    int a,b;
    //freopen("1.in","r",stdin);
    while(scanf("%d%d",&a,&b)!=EOF)
    {
        printf("%d
",calc(b)-calc(a-1));
    }
    return 0;
}