（bzoj 2734 [HNOI2012]集合选数）<状压DP>

题目

Time Limit: 10 Sec
Memory Limit: 128 MB

Description

《集合论与图论》这门课程有一道作业题，要求同学们求出{1, 2, 3, 4, 5}的所有满足以 下条件的子集：若 x 在该子集中，则 2x 和 3x 不能在该子集中。同学们不喜欢这种具有枚举性 质的题目，于是把它变成了以下问题：对于任意一个正整数 n≤100000，如何求出{1, 2,…, n} 的满足上述约束条件的子集的个数（只需输出对 1,000,000,001 取模的结果），现在这个问题就 交给你了。

Input

只有一行，其中有一个正整数 n，30%的数据满足 n≤20。

Output

仅包含一个正整数，表示{1, 2,…, n}有多少个满足上述约束条件 的子集。

Sample Input

4

Sample Output

8

样例解释

有8 个集合满足要求，分别是空集，{1}，{1，4}，{2}，{2，3}，{3}，{3，4}，{4}。

题解

比如选1，那么不能选的数有：2,3,4,6,8,9…
打表就是这个样子：
1 3 9 27…
2 6 18 54…
4 12 36 108…
8 24…
其中每个数向右是它的三倍，向下是它的两倍
也就是选了它就不能选和它同行和同列的数
需要求出方案数
是不是有点状压的意思了？
此外还需要一个数组记录选1时未考虑过的数（表格中未出现），并以这些点再次进行状压DP，再记下未考虑的点，然后继续更新
每次答案相乘取模
不用担心爆空间或时间，因为n≤100000，所以表最多也没几层

代码

// by spli
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<iostream>
#define LL long long
using namespace std;

int mod=1000000001;
int n;
LL ans=1;
int a[21][21],f[21][75000];
bool vis[100010];
int line[20];

LL cnt(int x){
    int len;
//  memset(f,0,sizeof(f));  //亲测这样会TLE  QAQ
    a[1][1]=x;
    for(int i=2;i<=18;++i)
        if((a[i-1][1]*2)<=n){
            a[i][1]=a[i-1][1]<<1;
            vis[a[i][1]]=1;
        }
        else{
            len=i-1;
            break;
        }
    for(int i=1;i<=len;++i)
        for(int j=2;j<=12;++j)
            if(a[i][j-1]*3<=n){
                a[i][j]=a[i][j-1]*3;
                vis[a[i][j]]=1;
            }
            else{
                line[i]=j-1;
                break;
            }
    line[0]=1;
    for(int i=1;i<=len;++i)
        for(int j=0;j<(1<<line[i]);++j)
            f[i][j]=0;
    f[0][0]=1;
    line[len+1]=0;
    f[len+1][0]=0;
    for(int i=0;i<=len;++i)
        for(int j=0;j<(1<<line[i]);++j)
            if(f[i][j]){
                if(j&(j>>1)) continue;
                for(int k=0;k<(1<<line[i+1]);++k){
                    if(k&j) continue;
                    if(k&(k>>1)) continue;
                    f[i+1][k]=(f[i+1][k]+f[i][j])%mod;
            }
        }
    return (LL)f[len+1][0]%mod;
}


int main(){
    scanf("%d",&n);
    for(int i=1;i<=n;++i)
        if(!vis[i])
            ans=(ans*cnt(i))%mod;
    cout<<ans%mod<<endl;
    return 0;
}