题目链接:http://acm.hdu.edu.cn/showproblem.php?pid=4810
思路:先把每个数字按位分离出来,存放1的个数,那么每位0的个数为n - 1的个数,然后利用组合数学和异或的原理,枚举奇数个1的情况,然后利用乘法和加法计数原理累加出来的就是该位的答案,最后乘上改为对应的数值最后加起来就是答案
代码:
#include <stdio.h>
#include <string.h>
const __int64 MOD = 1000003;
const int N = 1005;
int n, num, sum[32];
__int64 C[N][N], mi[32];
void tra(int num) {
int sn = 0;
while (num) {
sum[sn++] += num % 2;
num /= 2;
}
}
__int64 cal(int day) {
__int64 ans = 0;
for (int i = 0; i < 32; i++) {
for (int j = 1; j <= day && j <= sum[i]; j += 2) {
if (sum[i] < j || n - sum[i] < day - j) continue;
ans = (ans + C[sum[i]][j] * C[n - sum[i]][day - j] % MOD * mi[i] % MOD) % MOD;
}
}
return ans;
}
int main() {
mi[0] = 1;
for (int i = 1; i < 32; i++)
mi[i] = mi[i - 1] * 2;
for (int i = 0; i < N; i++) {
C[i][0] = C[i][i] = 1;
for (int j = 1; j < i; j++)
C[i][j] = (C[i - 1][j] + C[i - 1][j - 1]) % MOD;
}
while (~scanf("%d", &n)) {
memset(sum, 0, sizeof(sum));
for (int i = 0; i < n; i++) {
scanf("%d", &num);
tra(num);
}
for (int i = 1; i < n; i++) {
printf("%I64d ", cal(i));
}
printf("%I64d\n", cal(n));
}
return 0;
}HDU 4810 Wall Painting(组合数学),布布扣,bubuko.com
原文:http://blog.csdn.net/accelerator_/article/details/26005923