poj 3252 Round Numbers （组合数）

题目：http://poj.org/problem?id=3252

题意：求区间Start..Finish 之间 二进制0比1多或相等的数字的个数

借鉴一下别人的思路吧。。。

思路：记f(start, finish)为[start, finish]里Round Number的个数,那么要求f(start,finish),只需用f(0,finish)-f(0,start-1)即可,则问题转化为给定x,求出f(0,x).
假定x=(10101101),其长度为8位.而[0,x]中的数可分为二进制长度小于8位的和二进制长度等于8位的.

首先看二进制长度小于8位的,即求出长度在[0,7]区间内的Round Number个数.对于长度为len的二进制（最高位一定是1）,记其Round Number个数为R(len),可按照len的奇偶性分两种情况.
1. len=2k+1时,去掉最高位的1,剩下2k位里至少要有k+1个0,用C(n,m)表示n个位置选m个位置的方法.
有R(len)=C(2k,k+1)+C(2k,k+2)+…+C(2k,2k)=1/2*(2^(2k)-C(2k,k))
2. 同理可得,len=2k时,R(len)=1/2*(2^(2k-1))

接下来看二进制长度为8位的.
首先判断x本身,然后对于x的二进制,若将其中除了最高位的1以外的任意一个或多个1变为0,得到的数一定小于x,那么就能通过这种方法得到二进制位数和x相等且比x小的Round Number个数了.
最后将两部分结果求和即可.

精度方面,2000000000<2*1024*1024*1024=2^31,故用31位表示数组,又第一位总为1,所以组合数只用求到30,C(30,k) (0<=k<=30) <= 2^30,故用int即可.

#include <iostream>
#include<cstdio>
#include<cmath>
using namespace std;
int c[35][35];
int bir[35];
int init()
{
    int i,j;
    for(i=0;i<35;i++)
    {
        c[i][0]=c[i][i]=1;
    }
    for(i=2;i<35;i++)
    {
        for(j=1;j<i;j++)
        c[i][j]=c[i-1][j-1]+c[i-1][j];
    }
    return 0;
}
int chuli(__int64 x)
{
    if(x<=1)
    return 0;
    __int64  u;
    int i,j,k,n;
    int res=0;
    u=x;
    n=0;
    int num[2]={0};
    int t=0;
    while(u)
    {
        t=u%2;
        bir[n]=t;
        num[t]++;
        n++;
        u/=2;
    }
    for(i=2;i<n;i++)
    {
        if(i%2==1)
        {
            k=(i-1)/2;
            res+=(pow(2,(2*k))-c[2*k][k])/2;
        }
        else
        {
            k=i/2;
            res+=pow(2,(2*k)-2);
        }
    }
    if(num[0]>=num[1])
    {
        res++;
    }
    int n1=1;
    int n0=0;
    for(i=n-2;i>=0;i--)
    {
        if(bir[i])
        {
            for(j=i;j>=0&&j+n0+1>=i-j+n1;j--)
            {
                res+=c[i][j];
            }
            n1++;
        }
        else n0++;
    }
    return res;
}
int main()
{
    __int64 a,b;
    init();
    scanf("%I64d %I64d",&a,&b);
    printf("%d\n",chuli(b)-chuli(a-1));
    return 0;
}