Description
有n个整数,问从他们中取出若干个数字相与之后结果是0的有多少组。
答案比较大,输出对于 1,000,000,007 (1e9+7)取模后的结果。
n
≤
1
0
6
a
i
≤
1
0
6
n\le10^6\\a_i\le10^6
n≤106ai≤106
Solution
很容易想到设f[i]表示至少i个位置不为0的答案,我们容斥一下就能得到恰好i个位置不为0的答案,最终结果就是恰好0个位置不为0了
这里设g[x]表示二进制中含有x这个子集的元素数量,转移的时候可以像背包一样限制最高位就不会算重了
于是就可以用g[]算出f[]了
Code
#include <stdio.h>
#include <string.h>
#include <algorithm>
#define rep(i,st,ed) for (int i=st;i<=ed;++i)
#define drp(i,st,ed) for (int i=st;i>=ed;--i)
typedef long long LL;
const int MOD=1000000007;
const int N=2000005;
LL f[N],c[N],pow[N];
int read() {
int x=0,v=1; char ch=getchar();
for (;ch<'0'||ch>'9';v=(ch=='-')?(-1):(v),ch=getchar());
for (;ch<='9'&&ch>='0';x=x*10+ch-'0',ch=getchar());
return x*v;
}
void upd(LL &x,LL v) {
x+=v; (x>=MOD)?(x-=v):0;
}
int main(void) {
freopen("data.in","r",stdin);
freopen("myp.out","w",stdout);
int n=read();
rep(i,1,n) f[read()]++;
drp(i,20,0) {
rep(j,0,(1<<20)-1) {
if ((j>>i)&1) upd(f[j^(1<<i)],f[j]);
}
}
// for (int i=0;i<(1<<20);++i) printf("%lld\n", f[i]);
rep(i,0,(1<<20)-1) c[i]=c[i>>1]+(i&1);
pow[0]=1; rep(i,1,(1<<20)) pow[i]=(pow[i-1]+pow[i-1])%MOD;
LL ans=0;
rep(i,0,(1<<20)-1) {
// if (c[i]&1) upd(ans,MOD-pow[f[i]]+1);
// else upd(ans,pow[f[i]]-1);
if (c[i]&1) ans=(ans+MOD-(pow[f[i]]-1))%MOD;
else ans=(ans+pow[f[i]]-1)%MOD;
}
printf("%lld\n", ans);
return 0;
}