[BZOJ4762]最小集合/[JZOJ5151]幻梦终醒

题目大意

给定$n$个数$a_i$，你要从中选出一个非空子集，使得这个子集$\mathrm{AND}$和是$0$，并且这个子集的任意一个非空子集的$\mathrm{AND}$和都不是$0$。

$1\le n\le10^3,0\le a_i<2^{10}$

---

题目分析

又是一道好题。
在范老师@AwD的博客看懂的，大家可以过去膜拜一下。
思路是这样的：令$S_a$为$S$删去$a$之后形成的集合，令$f(S)$表示$S$集合中所有元素的$\mathrm{AND}$和。
那么答案就是

$$\sum_S[(f(S)=0)\land \operatorname*{\LARGE\land}_{a\in S}(f(S_a)\not=0)]$$
这个看起来不太好算，我们对后面的部分容斥一下，变成
$$\sum_S\sum_{S'\subseteq S}[(f(S)=0)\land \operatorname*{\LARGE\land}_{a\in S'}(f(S_a)=0)]\left(-1\right)^{|S'|}$$
其实就是
$$\sum_S\sum_{S'\subseteq S}[(f(S)=0)\land \left(\operatorname*{\LARGE\lor}_{a\in S'}f(S_a)=0\right)]\left(-1\right)^{|S'|}$$
这个怎么计算呢？令 $f_{i,k,j}$表示做到了第 $i$个位置，$f(S)=k$， $\operatorname*{\LARGE\lor}_{a\in S'}f(S_a)=j$的方案数（注意这个是乘上容斥系数的）。
显然必须要满足 $k\subseteq j$，也就是说状态数只有 $O(n3^{10})$。
转移的话，假设第 $i+1$个数是 $d$，那么有三种转移：
$\bullet\ $不选 $d$：$f_{i+1,k,j}\leftarrow f_{i,k,j}$
$\bullet\ S$中选择了 $d$，但是$S'$中没有选择： $f_{i+1,k\land d,j\land d}\leftarrow f_{i,k,j}$
$\bullet\ S$以及 $S'$中都选择了 $d$：$f_{i+1,k\land d,(j\land d)\lor k}\leftarrow-f_{i,k,j}$
初始状态 $f_{0,2^{10}-1,2^{10}-1}=1$，目标状态 $f_{n,0,0}$。
时间复杂度 $\mathrm O(n3^{10})$。

---

代码实现

#include <iostream>
#include <cctype>
#include <cstdio>

using namespace std;

int read()
{
    int x=0,f=1;
    char ch=getchar();
    while (!isdigit(ch)) f=ch=='-'?-1:f,ch=getchar();
    while (isdigit(ch)) x=x*10+ch-'0',ch=getchar();
    return x*f;
}

const int P=1000000007;
const int N=1005;
const int S=1024;

int f[2][S][S];
int a[N];
int n;

void update(int &x,int y){x=(x+y)%P;}

int dp()
{
    f[0][S-1][S-1]=1;
    for (int i=0;i<n;++i)
    {
        int cur=i&1,nxt=cur^1,x=a[i+1];
        for (int j=0;j<S;++j)
        {
            for (int k=j;k;k=(k-1)&j) f[nxt][k][j]=0;
            f[nxt][0][j]=0;
        }
        for (int j=0;j<S;++j)
        {
            for (int k=j;k;k=(k-1)&j) update(f[nxt][k][j],f[cur][k][j]),update(f[nxt][k&x][j&x],f[cur][k][j]),update(f[nxt][k&x][k|j&x],P-f[cur][k][j]);
            update(f[nxt][0][j],f[cur][0][j]);
        }
    }
    return f[n&1][0][0];
}

int main()
{
    freopen("never.in","r",stdin),freopen("never.out","w",stdout);
    n=read();
    for (int i=1;i<=n;++i) a[i]=read();
    printf("%d\n",dp());
    fclose(stdin),fclose(stdout);
    return 0;
}