洛谷3773 BZOJ4903 CTSC2017 吉夫特 数学 dp

题目链接
题意：
给你一个单调不升的数组a，求有多少个长度大于等于2的子序列满足${\prod }_{i=2}^k{C}_{b_i}^{b_{i-1}}\mathrm{m}\mathrm{o}\mathrm{d}2=1$。答案对1000000007取模。

题解：
首先题意是让你求一个相邻两数组合数膜2意义下都是1的子序列，也就是要求相邻两数组合数是奇数。我们考虑对于一个组合数$C_n^m=\frac{n!}{m!(n-m)!}$，如果它是个奇数，意味着我们能从$n!$中提取的因子2的个数和从$m!(n-m)!$中提取的因子2的个数相同。
根据组合数取模的卢卡斯定理，我们有$C_n^m\%2=C_{n\%2}^{m\%2}\ast C_{n/2}^{m/2}$。我们可以看出，这个式子相当于在对n和m进行二进制分解，从二进制的最低位到最高位依次判断n和m在该位对应的01，我们发现，如果在某一位n=0,m=1，那么$C_0^1$是不存在方案的，看作结果是0，不管其他位结果如何，答案在乘0之后肯定会变成0，所以我们可以推出，如果想要$C_n^m\%2=1$，那么需要满足的条件就是二进制下n的每一个是1的位，m对应位可以是1也可以是0，n是0的位，m对应位也一定要是0，也就是$n\&m=m$。如果像状压dp那样把一个二进制下的01串看作一个集合，那么$m\&m=m$就相当于要求m在二进制下是n的一个子集。
那么我们考虑dp，设dp[i]为数字i结尾的方案数，那么这个状态可以转移到所有i的二进制数的子集。最后对于每一个在$a_i$中出现的数，它的初值设为1，最后对答案的贡献要减1，相当于减去了长度为1只有$a_i$本身的子序列。
代码很短：

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

int n,tong[300010],a[300010];
int dp[300010],ans;
const int mod=1000000007;
int main()
{
 scanf("%d",&n);
 for(int i=1;i<=n;++i)
 {
  scanf("%d",&a[i]);
  tong[a[i]]=i;
  dp[a[i]]=1;
 }
 for(int i=1;i<=n;++i)
 {
  for(int j=(a[i]-1)&a[i];j;j=(j-1)&a[i])
  {
   if(tong[j]>i)
   dp[j]=(dp[j]+dp[a[i]])%mod;
  }
 }
 for(int i=1;i<=n;++i)
 ans=(ans+dp[a[i]]-1)%mod;
 printf("%d\n",ans);
 return 0;
}