【容斥&&状压】 CF449D Jzzhu and Numbers

题目描述：

---

首先我们考虑这么一个问题：
求$n$个数中，任意$j$，使得$a[j]\&i=i$，有多少个数满足。
设$f_i$表示与值为$i$的数的个数，首先可以得出以下转移：
$$f_i=\sum _{j\notin i}{f}_{i|(1<<j)}$$

但是这样转移很显然会有重复。
一般的转移式这样的：

	for (int i = Max; i >= 0; i--)
	  for (int j = 0; j <= 19; j++)
	    if (((i>>j)&1) == 0)
	      f[i] = (f[i] + f[i|(1<<j)]*x)%P;

但是我们神奇的发现，如果把$i$和$j$的循环翻一下，就可以避免重复，而是计算子集的答案了。

	for (int j = 0; j <= 19; j++)
	  for (int i = Max-1; i >= 0; i--)
	    if (((i>>j)&1) == 0)
	      f[i] = (f[i] + f[i|(1<<j)]*x)%P;

因此就是一个结论：

这是从我孙女的博客上摘下来的。
所以，对于这道题，做法也十分显然了

我们先用上面的方法做一遍，求出以$i$为子集的数的个数。
而后，从这些数中任意挑几个数，相与之后也一定以$i$为子集。
所以方案数是：
$$g[i]=2^{f[i]}-1$$

在类比之前的问题，只需要把加好变为减号即可。

---

#include<bits/stdc++.h>
using namespace std;

#define int long long
const int N = 2e6+10,P=1e9+7;
int n,a[N];
int f[1<<20];
int Max;

void dp(int x){
	for (int j = 0; j <= 19; j++)
	  for (int i = Max-1; i >= 0; i--)
	    if (((i>>j)&1) == 0)
	      f[i] = (f[i] + f[i|(1<<j)]*x)%P;
}

int power(int x,int y){
	int sum = 1;
	while (y){
		if (y&1) sum = (sum * x) % P; 
		y>>=1 , x = (x*x)%P;
	}
	return sum;
}

signed main(){
	scanf("%lld",&n);
	for (int i = 1; i <= n; i++) scanf("%lld",&a[i]) , f[a[i]]++;
	Max = (1<<20) - 1;
	dp(1);
	for (int i = 0; i <= Max; i++) f[i] = power(2,f[i])-1;
	dp(-1);
	printf("%lld",(f[0]+P)%P);
}