bzoj 4903: [Ctsc2017]吉夫特【lucas+狀壓dp】
阿新 • • 發佈:2018-10-25
二進制 clu 根據 main namespace ios ctsc ace 狀壓dp
首先根據lucas,
\[
C_n^m\%2=C_{n\%2}^{m\%2}*C_{n/2}^{m/2}
\]
讓這個式子的結果為計數的情況只有n&m==m,因為m的每一個為1的二進制位都需要n中這一位為1,否則結果就是0
所以枚舉子集,設f[i]為以i開頭的合法子序列個數,dp的時候枚舉子集從後往前dp即可
#include<iostream> #include<cstdio> using namespace std; const int N=300005,mod=1e9+7; int n,a[N],p[N],f[N],ans; int read() { int r=0,f=1; char p=getchar(); while(p>‘9‘||p<‘0‘) { if(p==‘-‘) f=-1; p=getchar(); } while(p>=‘0‘&&p<=‘9‘) { r=r*10+p-48; p=getchar(); } return r*f; } void jia(int &x,int y) { x+=y; x>=mod?x-=mod:0; } int main() { n=read(); for(int i=1;i<=n;i++) { a[i]=read(); p[a[i]]=i; } for(int i=1;i<=233333;i++) if(p[i]) { for(int j=i;j;j=(j-1)&i) if(p[j]>p[i]) jia(f[i],f[j]); jia(ans,f[i]++); } printf("%d\n",ans); return 0; }
bzoj 4903: [Ctsc2017]吉夫特【lucas+狀壓dp】