NKOI 集合选数

问题描述

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

输入格式

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

输出格式

仅包含一个正整数，表示{1, 2,…, n}有多少个满足上述约束条件 的子集。
这道题是老师评讲的，但是代码实现出现了问题
千万别用memset！千万别用memset！千万别用memset！每次用都超时
另外这道题虽然是对1000000001取模，看似在int范围内，但是其实在取模前那一步运算可能爆int。所以数据范围一定是要比取模的范围大的，至于大多少依照题目而定

#include<stdio.h>
#include<bits/stdc++.h>
using namespace std;
long long book[100005],l[20],ans=1,f[20][500000];
int main()
{
int n,i,j,k,hei;
scanf("%d",&n);
for(i=1;i<=n;i++)
{
if(!book[i])
{
for(j=0;j<=19;j++)
l[j]=0;
hei=0;
for(j=i;j<=n;j*=2)
{
hei++;
for(k=j;k<=n;k*=3)
{
book[k]=1;
l[hei]++;
}
}
for(k=1;k<=hei;k++)
for(j=0;j<=(1<<l[k])-1;j++)
{
f[k][j]=0;
if(!(j&(j>>1)))
{
int s=(((1<<l[k-1])-1)^j);
int s0;
for(s0=s;s0;s0=(s0-1)&s)
{
f[k][j]=(f[k][j]+f[k-1][s0])%1000000001;
}
f[k][j]=(f[k][j]+f[k-1][0])%1000000001;
if(k==1)
f[k][j]=1;
}
}
long long tot=0;
for(k=0;k<=(1<<l[hei])-1;k++)
tot=(tot+f[hei][k])%1000000001;
ans=(ans*tot)%1000000001;
}
}
printf("%lld",ans);
return 0;
}