【CTSC2017】【BZOJ4903】吉夫特 卢卡斯定理 DP

时间:2023-03-09 17:27:49
【CTSC2017】【BZOJ4903】吉夫特 卢卡斯定理 DP

题目描述

  给你一个长度为\(n\)的数列\(a\),求有多少个长度\(\geq 2\)的不上升子序列\(a_{b_1},a_{b_2},\ldots,a_{b_k}\)满足

\[\prod_{i=2}^k\binom{a_{b_{i-1}}}{a_{b_i}}\mod 2>0
\]

  答案对\({10}^9+7\)取模。

  \(n\leq211985,a_i\leq 233333\)

  \(\forall i\neq j,a_i\neq a_j\)

题解

  水题。

  先忽略长度\(\geq 2\)这个条件。

  根据卢卡斯定理,有\(a_{b_i}|a_{b_{i-1}}\)。

  从前往后DP。

  设\(f_i\)为前面那部分,最后一个数是\(i\)的方案数。

  转移直接枚举\(a_i|j\),让\(f_{a_i}+=f_j\)。

  时间复杂度:枚举子集的复杂度,\(O(3^{\log \max_{i=1}^na_i})\)

  p.s. gift在德语中的意思是毒。

代码

#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;
const int p=1000000007;
int f[1000010];
int main()
{
int n;
scanf("%d",&n);
int i,x,j;
int ans=0;
for(i=1;i<=n;i++)
{
scanf("%d",&x);
f[x]=1;
for(j=(x+1)|x;j<=233333;j=(j+1)|x)
f[x]=(f[x]+f[j])%p;
ans=(ans+f[x])%p;
}
ans=(ans-n+p)%p;
printf("%d\n",ans);
return 0;
}