POJ 3252 Round Numbers （组合计数）

  题目：POJ 3252  Round Numbers

  题目大意：二进制0的个数>=二进制1的个数为圆数，问一个较大范围里面圆数的个数。

  解题思路：之前这块书没好好复习，纯自己想的，基本上来讲花了我一下午...对于一个二进制有tot位的数，分两步操作，第一步：计算二进制位小于tot的圆数的个数，假设n位，那么，第n位肯定要是1，这样剩下的n-1个位置上要选1+ceil((n-2)/2)个0或者更多。然后就是组合计数加预处理解决问题。之后从最低位枚举，假设为m位，当这个位置的数为0的时候continue,因为从这个位置向下找不到更小的数（下一位在之前已经枚举过了），若为1，那么类分为二，若更高位值的1比0多，假设多x个，那么就有要从m-1个里面挑选出来x+ceil((m-1-x)/2)个0或者更多；若0比1多，又类分为2，若多很多，比剩下的数都多，那么就是选几个都可以，或者假设多y个，需要在m-1个数里面选(m-1-y)/2个数。之后就是代码的实现。

  ac代码：

#include<iostream>
#include<cstring>
#include<cstdio>
#include<cmath>
#include<algorithm>
using namespace std;
long long c[33][33];
long long sum[33];
long long gcd(long long a,long long b){
    return b==0?a:gcd(b,a%b);
}
void getc(){
    long long i,j;
    for (i=0;i<=32;i++){
        for (j=0;j<=i;j++){
            if (!j||j==i){
                c[i][j]=1;
            }
            else c[i][j]=c[i-1][j]+c[i-1][j-1];
        }
    }
    sum[1]=1;
    for (i=2;i<=32;i++){
        long long x=i-1;
        long long y=1+(ceil)((double)(i-2)/2);
        for (j=y,sum[i]=sum[i-1];j<=x;j++){
            sum[i]+=c[x][j];
        }
    }
}
long long solve(long long fun){
    long long a[33],i,j;
    long long tot=0;
    long long sum0=0;
    long long sum1=0;
    long long xx=fun;
    while (xx){
        if (xx%2==0)sum0++;
        else sum1++;
        a[++tot]=xx&1;
        xx>>=1;
    }
    long long ans;
    if (sum0>=tot/2)ans=1+sum[tot-1];
    else ans=sum[tot-1];
    for (i=1;i<=tot-1;i++){
        if (a[i]==1){
            sum1--;
            sum0++;
            long long x=i-1;
            long long yy=sum1-sum0;
            long long y;
            if (yy>=0){
                y=yy+(ceil)((double)(x-yy)/2);
            }
            else {
                if (-yy>=x){
                    y=0;
                }
                else {
                    y=(ceil)((double)(x+yy)/2);
                }
            }
            for (j=y;j<=x;j++){
                ans+=c[x][j];
            }
            sum0--;
        }
        else sum0--;
    }
    return ans;
}
int main(){
    getc();
    long long i,j,n,m,l,r;
    while (scanf("%lld%lld",&l,&r)!=EOF){
        cout<<solve(r)-solve(l-1)<<endl;
    }
}