Wall Painting ( 位运算 + 数学 )

HDU4810

#include<bits/stdc++.h>
#define LL long long
using namespace std;
const int maxn=1e3+5;
const int mod=1e6+3;

LL fac[maxn],inv[maxn],a[maxn],ant[maxn];
LL qpow(LL x,int p){
    LL ans=1;
    while(p){
        if(p&1) ans=ans*x%mod;
        x=x*x%mod;
        p>>=1;
    }return ans;
}
void init(){
    fac[0]=1;
    for(int i=1;i<=1e3;i++) fac[i]=fac[i-1]*i%mod;
    inv[1000]=qpow(fac[1000],mod-2);
    for(int i=999;i>=0;i--) inv[i]=inv[i+1]*(i+1)%mod;
}
LL C(int n,int m){
    if(m>n) return 0;
    if(m==0) return 1;
    return fac[n]*inv[m]%mod*inv[n-m]%mod;
}
int main(){
    init();
    int n;
    while(~scanf("%d",&n)){
        for(int i=1;i<=n;i++) ant[i]=0;
        for(int i=1;i<=n;i++) scanf("%lld",&a[i]);
        for(int i=0;i<=62;i++){
            int cnt=0;
            for(int j=1;j<=n;j++){
                cnt+=((a[j]>>i)&1);
            }
            for(int j=1;j<=cnt;j+=2){
                for(int k=0;k<=n-cnt;k++){
                    ant[j+k]=(ant[j+k]+C(cnt,j)*C(n-cnt,k)%mod*(1LL<<i)%mod)%mod;
                }
            }
        }
        for(int i=1;i<=n;i++) printf("%lld%c",ant[i],(i==n)?'\n':' ');
    }
    return 0;
}