[BZOJ2728][HNOI2012]与非（并查集+数位DP）

$$\mbox{NOT }A=A\mbox{ NAND }A$$
$$A\mbox{ AND }B=\mbox{NOT}(A\mbox{ NAND }B)$$
$$A\mbox{ OR }B=(\mbox{NOT }A)\mbox{NAND}(\mbox{NOT }B)$$
$$A\mbox{ XOR }B=((\mbox{NOT }A)\mbox{AND }B)\mbox{OR}(A\mbox{ AND}(\mbox{NOT B}))$$
从上面 $4$ 个式子得出一个神奇的性质—— $\mbox{NAND}$ 可以表示出所有的位运算！
但这样是否就可以表示出所有 $[0,2^K)$ 内的数呢？答案是否定的。
考虑到如果存在 $i,j\in[1,K]$ ，使得 每个 $A_x$ （ $1\leq x\leq n$ ） 都满足 $A_x$ 的第 $i$ 位和第 $j$ 位相等，那么不管怎样位运算，位运算结果的第 $i$ 位和第 $j$ 位总是相等的。所以用一个并查集维护位之间的相等关系，如果找到上面这样条件的 $i$ 和 $j$ ，就把 $i$ 和 $j$ 所在集合并起来。
这样，在区间 $[0,2^K)$ 内能计算出的数的个数就是 $2^{集合个数}$ 。
但我们要计算的是 $[L,R]$ 而不仅仅是 $[0,2^K)$ ，所以还要进行 数位 DP。
先预处理出 $pre_i$ ：是否存在一个 $j>i$ 使得 $j$ 和 $i$ 在同一个集合内，有则为 $j$ ，否则为 $0$ 。
先把询问拆成 $[0,R]$ 和 $[0,L-1]$ ，设现在要求 $[0,m]$ ， $d(x)$ 为 $m$ 第 $x$ 位的值。
定义状态 $f[i][0/1]$ 表示从第 $K$ 位到第 $i$ 位，小于 / 等于 $m$ 的第 $K$ 位到第 $i$ 位的个数。边界 $f[K][...]$ 直接处理，然后从 $K$ 到 $1$ 转移：
$$f[i][0]= \begin{cases} f[i+1][0] & pre_i\neq 0\\ f[i+1][0]\times 2 & pre_i=0 \end{cases}$$
$$f[i][1]= \begin{cases} 0 & pre_i\neq 0\mbox{ AND }d(i)\neq d(pre_i)\\ f[i+1][1] & \mbox{ELSE} \end{cases}$$
此外，如果 $d(i)=1$ ，还有一个转移：
$$f[i][0]+= \begin{cases} 0 & pre_i\neq 0\mbox{ AND }d(i)=d(pre_i)\\ f[i+1][1]&\mbox{ELSE} \end{cases}$$
最后答案 $f[1][0]+f[1][1]$ 。
代码：

#include <cmath>
#include <cstdio>
#include <cstring>
#include <iostream>
#include <algorithm>
#define For(i, a, b) for (i = a; i <= b; i++)
#define Rof(i, a, b) for (i = a; i >= b; i--)
using namespace std;
typedef long long ll;
const int N = 1005, M = 65;
int n, K, fa[M]; ll L, R, a[N], f[M][3], orz[M];
int cx(int x) {
    if (fa[x] != x) fa[x] = cx(fa[x]); return fa[x];
}
void zm(int x, int y) {
    int ix = cx(x), iy = cx(y); if (ix != iy) fa[iy] = ix;
}
int cmp(int x, int y) {
    if (x < y) return 0; if (x > y) return 2; return 1;
}
ll solve(ll m) {
    if (m == -1) return 0; int i; memset(f, 0, sizeof(f));
    f[K][cmp(0, m >> K - 1)]++; f[K][cmp(1, m >> K - 1)]++;
    Rof(i, K - 1, 1) {
        f[i][0] = f[i + 1][0] * (orz[i] ? 1 : 2);
        int x = (m >> i - 1) & 1, y = orz[i] ? (m >> orz[i] - 1) & 1 : -1;
        f[i][1] = f[i + 1][1] * (orz[i] ? x == y : 1);
        if (x) f[i][0] += f[i + 1][1] * (orz[i] ? x != y : 1);
    }
    return f[1][0] + f[1][1];
}
int main() {
    int i, j, k; cin >> n >> K >> L >> R;
    L = max(L, 0ll); R = min(R, (1ll << K) - 1);
    For (i, 1, K) fa[i] = i; For (i, 1, n) scanf("%lld", &a[i]);
    For (i, 1, K) For (j, i + 1, K) {
        bool flag = 1; For (k, 1, n)
            if (((a[k] >> i - 1) & 1) != ((a[k] >> j - 1) & 1))
                {flag = 0; break;}
        if (flag) zm(i, j);
    }
    For (i, 1, K) For (j, i + 1, K) if (cx(i) == cx(j)) orz[i] = j;
    cout << solve(R) - solve(L - 1) << endl;
    return 0;
}