组合数学（1）

poj3252

题意：给出两个数n，m，求出其区间中round number的个数，round number即满足在二进制写法中，0的个数大于等于1的个数的正整数。

利用前缀和思想，求出（1，n），（1，m）中round number的个数，相减即可。

求法：

对于数字21，它的二进制写法为 10101，共五位，我们首先求得位数低于五位的round number有多少，例如在求四位round number时，首位一定为1，后面三位可以有三个0或两个0，该计算可通过组合数的性质进行化简。

接下来，对于五位的round number，我们从最高位依次下推，遇到1就计算出从该位以下000…到111…这一区间中round number的个数，在此之前需要记录下高位中0和1的个数。最后需要补充的是，到达最低位后对10101这个数本身的检验，因为其后无1，所以会漏下该情况。

ac代码

#include<cstdio>
#include<iostream>
#include<cstring>


using namespace std;

int bit[34];
int c[34][34];

void init() {
	c[0][0]=1;
	for(int i=1; i<=33; i++) {
		c[i][0]=1;
		for(int j=1; j<=i; j++)
			c[i][j]=c[i][j-1]*(i-j+1)/j;
	}
	return;
}

int solve(int x) {
	int len=0;
	while(x) {
		bit[++len]=x&1;
		x>>=1;
	}
	int ans=0;
	for(int i=len-1; i>1; i--) {
		if(i&1)
			ans+=((1<<(i-1))-c[i-1][(i-1)/2])/2;
		else
			ans+=(1<<(i-2));
	}
	int cnt0=0,cnt1=1;
	for(int i=len-1; i>0; i--) {
		if(bit[i]) {
			for(int j=i-1; j>=(cnt1-cnt0+i-1)/2; j--)
				ans+=c[i-1][j];
			cnt1++;
		} else {
			cnt0++;
		}
	}
	if(cnt0>=cnt1)
        ans++;
	return ans;
}

int main() {
	init();
	int n,m;
	while(~scanf("%d%d",&n,&m)&&n+m) {
		printf("%d\n",solve(m)-solve(n-1));
	}
	return 0;
}