【51Nod1407】与与与与

有n个整数，问从他们中取出若干个数字相与之后结果是0的有多少组。
答案比较大，输出对于 1,000,000,007 (1e9+7)取模后的结果。

Input
第一行输入一个整数n。(1<=n<=1,000,000).
第二行有n个整数a[0],a[1],a[2],…a[n-1]，以空格分开.(0<=a[i]<=1,000,000)
Output
对于每一组数据，输出一个整数。
Input示例
3
2 3 3
4
0 1 2 3
Output示例
0
10

题解
http://blog.csdn.net/bahuia/article/details/55806426

代码

#include<bits/stdc++.h>
#define mod 1000000007
#define inf 1000000005
#define pa pair<int,int>
typedef long long ll;
using namespace std;
inline int read()
{
    int x=0,f=1;char ch=getchar();
    while (ch<'0'||ch>'9'){if (ch=='-')f=-1;ch=getchar();}
    while (ch>='0'&&ch<='9'){x=x*10+ch-'0';ch=getchar();}
    return x*f;
}
int n,dp[1000005],ans,fac[1000005];
bool g[1000005];
int main()
{
    n=read();fac[0]=1;
    for (int i=1;i<=1000000;i++) fac[i]=(fac[i-1]<<1)%mod;
    for (int i=1;i<=n;i++)
    {
        int x=read();
        dp[x]++;
    }
    for (int j=20;j>=0;j--)
        for (int i=1000000;i>=0;i--)
            if (i&(1<<j)) (dp[i^(1<<j)]+=dp[i])%=mod;
    ans=fac[dp[0]]-1;
    for (int i=1;i<=1000000;i++)
    {
        g[i]=(g[i>>1]+i&1)&1;
        if (g[i]) (ans-=fac[dp[i]]-1)%=mod;
        else (ans+=fac[dp[i]]-1)%=mod;
    }
    printf("%d",(ans+mod)%mod);
    return 0;
}