BZOJ 2734: [HNOI2012]集合选数

惊了 妙啊 我好菜啊

膜了题解之后 忍不住惊叹的蒟蒻决定也要写篇blog

写出这样一个矩阵
x 3x 9x 27x….
2x 6x 18x 54x…
4x 12x 36x 108x…
.
.

这个矩阵其实很小 当x=1时 才达到17*11的级别（后面更小。。
然后 不难发现这个矩阵相邻的数不能选 就可以状压DP算这个矩阵的方案数了！

那么x我们只需要每次取没出现过的最小值即可 1,5,7….
这些矩阵之间显然是互不影响的，所以把每个矩阵的答案乘起来就是最终答案了

代码不难打。一打就A了QAQ

#include<bits/stdc++.h>using namespace std;const int mod=1000000001,N=1e5+5,M=18;char B[1<<14],*S=B,*T=B;#define g (S==T&&(T=(S=B)+fread(B,1,1<<14,stdin),S==T)?-1:*S++)inline int read(){    int x=0,f=1; char ch=g;    while(ch<'0' || ch>'9'){if(ch=='-')f=-1; ch=g;}    while(ch>='0' && ch<='9'){x=(x<<1)+(x<<3)+ch-'0'; ch=g;}    return x*f;}int n,ans,a[M][15],d[M],b[M],f[M][(1<<15)+5];bool v[N];inline void add(int &x,int y){x=(x+y)%mod;}void solve(int x){    int m,i,j,k; a[1][1]=x;    for(m=1;a[m][1]<=n;++m){        a[m+1][1]=a[m][1]<<1 , b[m]=0;        for(i=1;a[m][i]<=n;++i)            b[m]+=d[i], a[m][i+1]=a[m][i]*3, v[a[m][i]]=1;        for(i=0;i<=b[m];++i)f[m][i]=0;    }    --m;    f[0][0]=1; b[0]=1;    for(i=0;i<m;++i)        for(j=0;j<=b[i];++j) if( !( j&(j<<1) ) && !( j&(j>>1) ) )            for(k=0;k<=b[i+1];++k) if( !( k&(k<<1) ) && !( k&(k>>1) ) )                if( !(j&k) ) add(f[i+1][k],f[i][j]);    int s=0;    for(i=0;i<=b[m];++i) add(s,f[m][i]);    ans=(1ll*ans*s)%mod;}int main(){    n=read(); ans=d[1]=1;    for(int i=2;i<M;++i) d[i]=d[i-1]<<1;    for(int i=1;i<=n;++i) if(!v[i]) solve(i);    printf("%d\n",ans);    return 0;}