题目大意:给定n,求集合S={1,2,3,...,n}有多少子集满足对于任意集合中任意两个数x和y,x≠2y并且x≠3y
原题解见 http://www.cnblogs.com/Randolph87/p/3677786.html
对于n以内任意与6互质的整数x,我们列出一个矩阵
x 3x 9x 27x ...
2x 6x 18x 54x ...
4x 12x 36x 108x ...
...
向下是*2,向右是*3,这样这个矩阵的任意两个相邻的数都不能出现在同一集合中
方案数状压DP即可 最后把所有x的矩阵方案数用乘法原理乘在一起即可
很巧妙的一道题
#include<cstdio>
#include<cstring>
#include<iostream>
#include<algorithm>
#define Mo 1000000001
using namespace std;
int n,f[2][1<<12];
bool usable[1<<12];
long long ans=1;
int State_Compression_DP(int now)
{
int i,j,s1,s2,last=0,re=0;
memset(f,0,sizeof f);
f[1][0]=1;
for(i=0;now*(1<<i)<=n;i++)
{
int temp;
memset(f+(i&1),0,sizeof(f)>>1 );
for(j=0,temp=1;now*(1<<i)*temp<=n;j++,temp*=3);
for(s1=0;s1<1<<last;s1++)
if(usable[s1])
for(s2=0;s2<1<<j;s2++)
if(usable[s2])
if( (s1&s2)==0 )
f[i&1][s2]+=f[~i&1][s1],f[i&1][s2]%=Mo;
last=j;
}
for(j=0;j<1<<last;j++)
re+=f[~i&1][j],re%=Mo;
return re;
}
int main()
{
int i;
cin>>n;
for(i=0;i<1<<12;i++)
if( !(i>>1&i) && !(i<<1&i) )
usable[i]=1;
for(i=1;i<=n;i++)
if(i%2&&i%3)
ans*=State_Compression_DP(i),ans%=Mo;
printf("%d\n",(int)ans);
return 0;
}