2018.10.25 uestc上天的卿学姐（计数dp）

传送门

看了DZYODZYODZYO的题解之后发现自己又sbsbsb了啊。

直接dpdpdp是O(2d)O(2^d)O(2d)更新，O(1)O(1)O(1)查询或者O(1)O(1)O(1)更新，O(2n)O(2^n)O(2n)查询的。

然后我就不会了233.

显然可以利用分块暴力的思想。

每次枚举前半段来计算当前答案，然后枚举后半段来更新dpdpdp数组。

这样效率O(n∗2d2)O(n*2^{\frac d 2})O(n∗22d​)可以通过全部测试点。

代码：

#include<bits/stdc++.h>
using namespace std;
inline int read(){
	int ans=0;
	char ch=getchar();
	while(!isdigit(ch))ch=getchar();
	while(isdigit(ch))ans=(ans<<3)+(ans<<1)+(ch^48),ch=getchar();
	return ans;
}
inline int read_2(){
	int ans=0;
	char ch=getchar();
	while(!isdigit(ch))ch=getchar();
	while(isdigit(ch))ans=(ans<<1)+(ch^48),ch=getchar();
	return ans;
}
const int mod=1e9+7,lim=1<<8;
int f[lim+5][lim+5],n,d,ans=0;
int main(){
	n=read(),d=read();
	for(int i=1;i<=n;++i){
		int x=read_2(),fi=x>>8,se=x^(fi<<8),stat,sum=1;
		for(stat=fi;;stat=fi&(stat-1)){
			(sum+=f[stat][se])%=mod;
			if(!stat)break;
		}
		for(stat=se;stat<lim;stat=se|(stat+1))(f[fi][stat]+=sum)%=mod;
	}
	for(int i=0;i<lim;++i)(ans+=f[i][lim-1])%=mod;
	cout<<ans;
	return 0;
}