HDU4810:Wall Painting

题目大意

传送门
求n个数里面，取i个数异或的所有组合的和，i取1~n

思路

我们把每个数转换为2进制(因为异或)，我们可以求出有多少个数第i位是1，我们设为$bi{t}_i$,我们假设我们在这$bi{t}_i$中选k个数出来异或，那么，就有$C(bi{t}_i,k)\times C(n-bi{t}_i,t-k)$种方案,其中t是我们这位取$0$和$1$的数量总和，这个就表示从$bi{t}_i$个1中选k个数出来，从剩下的$n-bi{t}_i$中选出$t-k$个0出来。然后，如果$k$是奇数的话，它是有贡献的，每个方案有$2^i$,否则$k$是偶数的话，说明这一位就是0，就是没有贡献，然后我们就可以$O(n\times logn\times logn)$解决这道题了，需要注意的是，建议用杨辉三角预处理组合数，如果用逆元还像我一样脑残写了快速幂，你就会$T$到飞起，最后注意开long long

代码

#include <cstdio>
#include <vector>
#include <cstring>
#include <iostream>
using namespace std;

#define Int register int
#define int long long
#define mod 1000003
#define MAXN 1005

int n;
int bit[45];
int a[MAXN],fac[MAXN],ans[MAXN];

int C[MAXN][MAXN];

void init()
{
	C[0][0] = 1;
	for (Int i = 1;i <= 1000;++ i)
	{
		C[i][0] = C[i][i] = 1;
		for (Int j = 1;j < i;++ j)
			C[i][j] = (C[i - 1][j - 1] + C[i - 1][j]) % mod;
	}
}

void read (int &x)
{
	x = 0;char c = getchar();int f = 1;
	while (c < '0' || c > '9'){if (c == '-') f = -f;c = getchar();}
	while (c >= '0' && c <= '9'){x = (x << 3) + (x << 1) + c - '0';c = getchar();}
	x *= f;return ;
}

void write (int x)
{
	if (x < 0){x = -x;putchar ('-');}
	if (x > 9) write (x / 10);
	putchar (x % 10 + '0');
}

signed main()
{
	init();
	while (scanf ("%lld",&n) != EOF)
	{
		memset (a,0,sizeof (a));
		memset (bit,0,sizeof (bit));
		memset (ans,0,sizeof (ans));
		for (Int i = 1;i <= n;++ i) read (a[i]);
		for (Int i = 1;i <= n;++ i)
		{
			for (Int j = 30;j >= 0;-- j)
				if ((a[i] >> j) & 1)
					bit[j] ++;
		}
		for (Int i = 1;i <= n;++ i)
		{
			for (Int j = 0;j <= 30;++ j)
			{
				int val = (1ll << j);
				for (Int k = 1;k <= bit[j] && k <= i;k += 2)
					ans[i] = (ans[i] + 1ll* C [bit[j]][k] % mod * C [n - bit[j]][i - k] % mod * val % mod) % mod;
			}
		}
		for (Int i = 1;i < n;++ i)
			write (ans[i] % mod),putchar (' ');
		write (ans[n]),putchar ('\n');
	}
	return 0;
}