有趣的异或

描述

这是一个有趣的题，给你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

4
1 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");
    }
}