bzoj2734: [HNOI2012]集合选数

2734: [HNOI2012]集合选数

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}。

HINT

 

Source

day2

 

1 3  9  27…

2 6 18 54…

4 12 36 108…

以未出现过的数作为矩阵的第一个数。那么矩阵互相之间的取数不会产生影响，ans=每个矩阵的方案数相乘

写出这样的矩阵，发现横向不超过18行，纵向不超过11列，每个数不能与上下左右的数同时出现，直接将列压缩成状态进行dp就行了；

细节看代码；

 

#include<cstdio>
#include<cstring>
#include<iostream>
#include<algorithm>
#define MOD 1000000001
using namespace std;
 
int n,m,mark[100008];
long long ans,dp[20][2049],a[28][20],b[28],ex[28];
 
long long cal(int x){
    memset(b,0,sizeof(b));
    memset(a,0,sizeof(a));
    a[1][1]=x;
    for(int i=2;i<=18;i++){
        a[i][1]=a[i-1][1]*2;
    }
    for(int i=1;i<=18;i++)
        for(int j=2;j<=11;j++)
            a[i][j]=a[i][j-1]*3;
    for(int i=1;i<=18;i++)
        for(int j=1;j<=11;j++)
            if(a[i][j]<=n){
                b[i]+=ex[j-1];
                mark[a[i][j]]=1;
            }
    memset(dp,0,sizeof(dp));
    dp[0][0]=1;
    for(int i=0;i<18;i++)
        for(int j=0;j<=b[i];j++)
            if(dp[i][j]){
                for(int k=0;k<=b[i+1];k++)
                    if((j&k)==0&&(j&(j>>1))==0)
                        dp[i+1][k]=(dp[i+1][k]+dp[i][j])%MOD;
            }
    return dp[18][0];
}
 
int main(){
    scanf("%d",&n);
    ans=1;
    ex[0]=1;
    for(int i=1;i<=20;i++) ex[i]=ex[i-1]*2;
    for(int i=1;i<=n;i++){
        if(mark[i]==0){
            ans=(ans*cal(i))%MOD;
        }
    }
    printf("%lld",ans);
}