[BZOJ4903][UOJ300][Ctsc2017]吉夫特（DP+优化）

显然，要满足

$$\prod_{i = 2}^k \binom{a_{b_{i-1}}}{a_{b_i}} \bmod 2= \binom{a_{b_1}}{a_{b_2}} \times \binom{a_{b_2}}{a_{b_3}} \times \cdots \times \binom{a_{b_{k-1}}}{a_{b_k}} \bmod 2 > 0$$
当且仅当上面乘式的每一项都是奇数。
根据Lucas定理可以推出， $\binom{n}{m}$为奇数，当且仅当对于任何一个 $i$，都有：
如果$n$的第 $i$个二进制位为$0$，那么 $m$的第$i$个二进制位不能为 $1$。
如果把一个二进制数看作一个集合（第$i$位为 $1$表示集合里有$i$，否则表示集合里没有 $i$），
那么可以得出：
$\binom{n}{m}\equiv1(\mod 2)\Leftrightarrow m$是 $n$的子集
定义状态$f[i]$表示以 $i$为终点的合法子序列个数。
加一个统计数组$sum$， $sum[s]$表示当前满足 $s$是$a_i$的子集的 $f_i$之和。
然后每算出一个 $f[i]$之后枚举 $a_i$的子集加入 $sum$。
由于所有的 $a_i$互不相同，所以复杂度为 $O(3^{\log a_i})$。
代码：

#include <cmath>
#include <cstdio>
#include <cstring>
#include <iostream>
#include <algorithm>
using namespace std;
inline int read() {
    int res = 0; bool bo = 0; char c;
    while (((c = getchar()) < '0' || c > '9') && c != '-');
    if (c == '-') bo = 1; else res = c - 48;
    while ((c = getchar()) >= '0' && c <= '9')
        res = (res << 3) + (res << 1) + (c - 48);
    return bo ? ~res + 1 : res;
}
const int N = 3e5 + 5, MX = 1e9 + 7;
int n, a[N], cnt[N], f[N], ans;
int main() {
    int i, j; n = read(); for (i = 1; i <= n; i++) a[i] = read();
    for (i = 1; i <= n; i++) {
        ans = (ans + (f[i] = cnt[a[i]])) % MX; f[i] = (f[i] + 1) % MX;
        for (j = a[i]; j; j = (j - 1) & a[i]) cnt[j] = (cnt[j] + f[i]) % MX;
    }
    cout << ans << endl;
    return 0;
}