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}。
题解:
这题不看题解真的是想不到啊,因为这个集合不能让2x和3x共存,所以(dalao们)就想到构造出一个如下的矩阵:
x 3x 9x …
2x 6x 18x …
4x 12x 36x …
…
这样对于每个x只要在矩阵中出现一次即可,就有了若干个矩阵,每个矩阵之间互不干预,所以只要将所有矩阵的合法选取方案数乘起来即可,我们看这个矩阵的行最大不超过17,所以我们直接硬上状压DP就行了,F[i][j]表示第i行状态为j的方案数,在转移的时候扔掉两个1相邻的情况,再判断一下上下有相同的1的情况即可。
#include<cstdio>
#include<cstdlib>
#include<iostream>
#include<iomanip>
#include<ctime>
#include<cmath>
#include<algorithm>
#include<cstring>
#include<string>
using namespace std;
const int mod=1e9+1;
bool vis[300000];
int a[20][20];
long long f[20][200000];
int top[20];
int n;
long long calc(int x)
{
a[1][1]=x;
int i;
for(i=1;;i++)
{
int j;
for(j=1;;j++)
{
vis[a[i][j]]=true;
a[i][j+1]=a[i][j]*3;
if(a[i][j+1]>n) break;
}
top[i]=j;
a[i+1][1]=a[i][1]*2;
if(a[i+1][1]>n) break;
}
top[i+1]=0;
int temp=i;
for(int i=1;i<=temp;i++)
for(int j=0;j<(1<<top[i]);j++)
f[i][j]=0;
f[temp+1][0]=0;
f[0][0]=1;
top[0]=0;
for(i=0;i<=temp;i++)
{
for(int j=0;j<(1<<(top[i]));j++)
{
if(f[i][j]==0) continue;
if((j&(j<<1))) continue;
for(int k=0;k<(1<<(top[i+1]));k++)
{
if(j&k) continue;
(f[i+1][k]+=f[i][j])%=mod;
}
}
}
return f[temp+1][0];
}
int main()
{
scanf("%d",&n);
long long ans=1;
for(int i=1;i<=n;i++)
if(!vis[i])
(ans*=calc(i))%=mod;
printf("%lld\n",ans);
return 0;
}