1. 程式人生 > >2018.10.25 uestc上天的卿學姐(計數dp)

2018.10.25 uestc上天的卿學姐(計數dp)

傳送門 看了DZYODZYO的題解之後發現自己又sbsb了啊。 直接dpdpO(2d)O(2^d)更新,O(1)O(1)查詢或者O(1)O(1)更新,O(2n)O(2^n)查詢的。 然後我就不會了233. 顯然可以利用分塊暴力的思想。 每次列舉前半段來計算當前答案,然後列舉後半段來更新dpdp陣列。 這樣效率O(n2d2)O(n*2^{\frac d 2})可以通過全部測試點。 程式碼:

#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; }