有趣的异或

描述

这是一个有趣的题，给你n个数，聪明的大家都知道从中挑出k个数的方案数是多少？（ps:当然不会考你这么简单的题）。

现在我们按照如下规则产生A(k)：

1.从n个数中挑出k个数

2.将这k个数作异或运算（S=a1^a2^..ak）

3.将每一种可能的方案产生的S求和

输入

多组测试数据，以EOF结束程序。
每组测试数据第一行包含一个数n（1 <= n <= 1000）.
第二行包含n个数，以空格隔开。保证每个数都在int范围内。

输出

对于每组数据输出一行n个数，以空格隔开，依次表示A(1),A(2)...A(n),行末没有空格。
需要注意的是，输出的结果可能很大，需要对1000003（10^6+3）取余。

样例输入1

41 2 10 1

样例输出1

14 36 30 8

思路：

对于每个数字拆解为二进制，如果二进制上某位置有贡献度，那必然是取了奇数个1，这样排列组合一下

#include<bits/stdc++.h>using namespace std;typedef long long ll;ll C[1005][1005];const ll mod=1000003;ll digit[1005];ll a[1005];int main(){    for(int i=0;i<=1000;i++)    {        C[i][0]=1;        for(int j=1;j<=i;j++)        {            C[i][j]=C[i-1][j-1]+C[i-1][j];            C[i][j]%=mod;        }    }    int n;    while(scanf("%d",&n)!=EOF)    {        memset(digit,0,sizeof(digit));        int big=0;        for(int i=1;i<=n;i++)        {            int tot=0;            scanf("%lld",&a[i]);            while(a[i])            {                if(a[i]&1)                {                    digit[tot]++;                }                a[i]>>=1;                tot++;            }            big=max(big,tot);        }        for(int i=1;i<=n;i++)        {            ll ans=0;            for(int j=0;j<=big;j++)            {                for(int k=1;k<=i&&k<=digit[j];k+=2)                {                    if(i-k<=n-digit[j])                    {                        ans+=((C[digit[j]][k]*C[n-digit[j]][i-k])%mod)*pow(2,j);                        ans%=mod;                    }                }            }            if(i!=1)                printf(" ");            printf("%lld",ans);        }        printf("\n");    }}