『容斥』51nod 1407 与与与与

$\mathrm{P}\mathrm{r}\mathrm{o}\mathrm{b}\mathrm{l}\mathrm{e}\mathrm{m}$

---

$\mathrm{S}\mathrm{o}\mathrm{l}\mathrm{u}\mathrm{t}\mathrm{i}\mathrm{o}\mathrm{n}$

这道题的方向不是$O(n)$,而是 $O(v)$。

我们发现答案为：全集 $-\mathrm{a}\mathrm{n}\mathrm{d}$有至少 $1$ 位 $0$ 的 $+\mathrm{a}\mathrm{n}\mathrm{d}$ 后有至少 $2$ 位 $0$ 的…

因此我们可以根据一个数字 $\mathrm{a}\mathrm{n}\mathrm{d}$ 操作后的位数来确定容斥系数。

考虑 $\mathrm{a}\mathrm{n}\mathrm{d}$ 后位数至少 $k$ 的答案，假设有 $f_x$个数的其中一个子集为$x$，那么把这些数 $and$ 以后，这些值至少为 $x$ ，而这些方案恰好是至少有$cn{t}_x$位的贡献。

因此答案为：$$\sum _x(-1{)}^{cn{t}_x}\times (2^{f_x}-1)$$

$f_x$ 的求解方法参考这篇精美博客哦！

---

$\mathrm{C}\mathrm{o}\mathrm{d}\mathrm{e}$

#include <bits/stdc++.h>
#define int long long

using namespace std;
const int N = 2100000;
const int P = 1e9 + 7;

int n;
int power[N], f[N], cnt[N];

int read(void)
{
	int s = 0, w = 0; char c = getchar();
	while (!isdigit(c)) w |= c == '-', c = getchar();
	while (isdigit(c)) s = s*10+c-48, c = getchar();
	return w ? -s : s;
}

void work(void)
{
	int res = 0;
	memset(f, 0, sizeof f);
	for (int i=1;i<=n;++i) f[read()] ++;
	for (int j=20;j>=0;--j)
		for (int i=1e6;i>=0;--i)
			if (((i >> j) & 1) == 0)
				f[i] += f[i | 1 << j];
	for (int i=0;i<=1e6;++i)
	{
		if (cnt[i] % 2 == 0) res = (res + (power[f[i]] - 1)) % P;
		if (cnt[i] % 2 == 1) res = (res - (power[f[i]] - 1)) % P;
	}
	res = (res % P + P) % P;
	cout << res << endl;
	return;
}

signed main(void)
{
	power[0] = 1;
	for (int i=1;i<=1e6;++i) {
		power[i] = power[i-1] * 2 % P;
		cnt[i] = cnt[i >> 1] + (i & 1);
	}
	while (cin >> n) work();
	return 0;
}