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

来源：互联网发布：js缩略图点击放大编辑：程序博客网时间：2024/04/29 07:47

给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;}

0 0