51nod1407 与与与与 容斥+DP

Description

---

有n个整数，问从他们中取出若干个数字相与之后结果是0的有多少组。

答案比较大，输出对于 1,000,000,007 (1e9+7)取模后的结果。
$$\begin{gathered} n\le 10^6 \\ {a}_i\le 10^6 \end{gathered}$$

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;
}