bzoj2734：[HNOI2012]集合选数（状压DP）

　　菜菜的喵喵题~

　　化序列为矩阵！化腐朽为神奇！左上角为1，往右每次*3，往下每次*2，这样子就把问题转化成了在矩阵里选不相邻的数有几种可能。

　　举个矩阵的例子

　　1 3 9 27
　　2 6 18 54
　　4 12 36 108

　　这样最多11列，最多17行，那么方案数就可以用状压了。　

　　但是我们会发现，矩阵里没有5，所以我们要把5作为左上角再算一次答案，最后把所有矩阵的答案用乘法原理乘起来就好

#include<iostream>
#include<cstring>
#include<cstdlib>
#include<cstdio>
#include<cmath>
#include<algorithm>
#define MOD(x) ((x)>=mod?(x)-mod:(x))
using namespace std;
const int maxn=,mod=1e9+;
int n,ans,cnth;
int f[][<<],cntl[];
bool v[maxn];
inline void read(int &k)
{
    int f=;k=;char c=getchar();
    while(c<''||c>'')c=='-'&&(f=-),c=getchar();
    while(c<=''&&c>='')k=k*+c-'',c=getchar();
    k*=f;
}
inline int find(int x)
{
    cnth=;memset(cntl,,sizeof(cntl));
    for(int i=,fir=x;fir<=n;i++,fir*=,cnth++)
    for(int j=,sec=fir;sec<=n;j++,sec*=,cntl[i]++)v[sec]=;
    f[][]=;
    for(int i=;i<=cnth;i++)for(int j=;j<=(<<cntl[])-;j++)f[i][j]=;
    for(int i=;i<=cnth;i++)
    for(int j=;j<(<<cntl[i]);j++)
    if(!(j&(j>>)))
    for(int k=;k<(<<cntl[i-]);k++)
    if(!(k&(k>>)))if(!(j&k))f[i][j]=MOD(f[i][j]+f[i-][k]);
    int sum=;
    for(int i=;i<=(<<cntl[cnth])-;i++)sum=MOD(sum+f[cnth][i]);
    return sum;
}
int main()
{
    read(n);ans=;
    for(int i=;i<=n;i++)
    if(!v[i])ans=1ll*ans*find(i)%mod;
    printf("%d\n",ans);
}