[HDOJ 4810][Vjudge 48316] Wall Painting [动态规划]

给n个数，问从中选k个，取异或，问所有不同的选择方案得到的异或值的和为多少

数据范围：n不超过1000

动态规划，定义状态dp[i][j][k]为前i个数中选j个，第k位异或值为1的方案数

#include <cstdio>
#include <cstring>

const int mod=1000003;

int dp[2][1001][32];
int c[1001][1001]={0};
int a[1001];
int n;

int main() {
	int i,j,k;
	c[0][0]=1;
	for (i=0;i<1001;i++) {
		c[i][0]=c[i][i]=1;
		for (j=1;j<i;j++) 
			c[i][j]=(c[i-1][j-1]+c[i-1][j])%mod;
	}
	while (scanf("%d",&n)!=EOF) {
		memset(dp,0,sizeof(dp));
		for (i=1;i<=n;i++) scanf("%d",&a[i]);
		for (i=1;i<=n;i++) {
			for (j=1;j<=i;j++) {
				for (k=0;k<32;k++) {
					if (a[i]&1<<k) dp[i&1][j][k]=c[i-1][j-1]-dp[i&1^1][j-1][k]+mod;
					else dp[i&1][j][k]=dp[i&1^1][j-1][k];
					dp[i&1][j][k]=(dp[i&1][j][k]+dp[i&1^1][j][k])%mod;
				}
			}
		}
		for (i=1;i<=n;i++) {
			int ans=0;
			for (k=0;k<32;k++) ans=(ans+(long long)dp[n&1][i][k]*(1<<k))%mod;
			if (i==1) printf("%d",ans);
			else printf(" %d",ans);
		}
		printf("\n");
	}
	return 0;
}