[CTSC2017]吉夫特

CTSC2017

一句话题意：

• 这道题是我写的第三遍了，真是一点长进也没有……甚至还看了点题解。
• 你要选择一个不上升子序列，并且使得${\Pi }_{i=2}{C}_{a_{bi}}^{a_{bi-1}}mod2>0$。因为在$mod$ 2的意义下，只有0和1这两个值。所以对于任意的$C_{a_{bi}}^{a_{bi-1}}=1$。我们如何保证这个必须等于1呢？这样看好像看不出来什么东西。尝试把阶乘写出来消一消，还是没什么思路。这个$mod$ 2到底有什么特殊的作用呢？
• Lucas定理啊！咋没想到，哎。一个组合数C，在$mod$ 2的意义下，套用Lucas定理，相当于两个数在二进制下的每一位求一个组合数的连乘。那么每一位的组合数都得为1，即$a_{bi-1}$的每一位都不大于$a_{bi}$。相当于$a_{bi-1}$是$a_{bi}$的子集。
• 所以循环的时候，枚举当前值的子集，递推方案数就好了。
• 实际操作中，题目的组合数是反着的，所以你需要倒序遍历。还有就是注意这个序列的长度必须不小于2。

Coding

#include<bits/stdc++.h>
using namespace std;
const int N=1e6+10;
const int mod=1e9+7;
int n,a[N],f[N],flag[N];
int main(){
	scanf("%d",&n);
	for(int i=1;i<=n;++i) scanf("%d",&a[i]);
	for(int i=n;i>0;--i){
		flag[a[i]]=1;
		for(int j=(a[i]-1)&a[i];j;j=(j-1)&a[i]){
			if(flag[j]&&!f[j]) f[a[i]]=(f[a[i]]+1)%mod;
			else if(flag[j]) f[a[i]]=(f[a[i]]+f[j]+1)%mod;
		}
	}
	int ans=0;
	for(int i=1;i<=n;++i) ans=(ans+f[a[i]])%mod;
	cout<<ans<<endl;
}