POJ3252 Round Numbers

POJ3252 http://poj.org/problem?id=3252

题目大意：

找出小于等于N的且二进制表示中0的数目大于等于1的数目的数字(Round Numbers ,RN)的个数

 

解题思路：

参考

首先考虑，以23(10111 B)为例。

首先计算出它的二进制表示长度，即5.

首先考虑比它长度小的，（长度更小当然数字也更小）。

四位：第四位首先首位是1（否则无法单独构成），此时cnt1 = 1，接着考虑后面3位，要保证0的个数大于或等于1的个数，即3位中必须要有2位(C[3][2])或者3位均为0（C[3][3]); 即C[3][2] + C[3][3];

三位：C[2][2];

两位：C[1][1];

一位：0

接下来考虑与23同等长度且<=情况：

23 = 10111B 所以前两位都限定为1,0（不然无法满足）

第三位如果为0 即后面两位只能有一个1或者没有1 即0的个数为C[2][2] + C[2][1];

然后将第三位恢复为1，考虑第四位，第四位如果为0，则最后一位必定为0 C[1][1];

然后将第四位恢复为1，考虑第五位，最后不存在；

所有情况即为综上。

 

下面开始推公式。

长度为L 的数字（即考虑后L - 1）所包括的RN数：

(参考：kuangbin https://www.cnblogs.com/kuangbin/archive/2012/08/22/2651730.html)

L - 1为偶数
sum = C[l - 1][(l - 1) / 2 + 1] + C[l - 1][(l - 1) / 2 + 2] +...+C[l-1][l-1]
    =  (2^(l-1)-C[l - 1][(l - 1) / 2] / 2;
L - 1是奇数
sum = (2^(l-1)) / 2;

 

 

 

 

AC代码：

#include <iostream>
#include <queue>
#include <algorithm>
#include <cstdio>
using namespace std;
int C[33][33];
void init() {
    C[0][0] = 1;//坑
    C[1][0] = 1;
    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 - 1] + C[i - 1][j];//YANGHUI TRIANGLE
        }
        C[i][i] = 1;
    }
}
int bits[33];
int solve(int n) {
    if(n <= 1) return 0;
    int len = 0;
    while(n) {
        if(n&1) bits[len++] = 1;
        else bits[len++] = 0;
        n >>= 1;
    }
    int ans = 0;
    for(int i = len - 1; i > 0; i--) {
        if(i % 2==0) ans +=((1 << (i - 1))) / 2;
        else ans += ((1 << (i - 1)) - C[i - 1][(i - 1)/2])/ 2;
    }
    int cnt0 = 0,cnt1 = 0;
    for(int i = 0;i < len; i++) {
        if(bits[i] == 0) cnt0++;
        else cnt1++;
    }
    if(cnt0 >= cnt1) ans++; // n
    cnt0 = 0;
    cnt1 = 1;
    for(int i = len - 2; i >= 0; i--) {
        if(bits[i] == 1) {
            for(int j = i; j >= 0 && j + cnt0 + 1 >= i - j + cnt1;j--) //{
                ans+=C[i][j];
                cnt1++;
            //}
        }
        else cnt0 ++;
    }
    return ans;
}
int main() {
    init();
    int a,b;
    while(~scanf("%d%d",&a,&b)) {
        printf("%d\n",solve(b) - solve(a - 1));
    }

    return 0;

}