bzoj2734: [HNOI2012]集合选数 压状dp

写出如下矩阵

1 3  9  27...

2 6 18 54...

4 12 36 108...

发现最多有11列。。。

我们在其中选取一些数，相邻的不能选择

然后就可以状压求方案数了，但是5没有出现，同样5的倍数也没有出现，7也如此。。

应该记录哪些数字出现过，没出现过就作为矩阵的第一个元素，最后把若干个矩阵的方案相乘

#include <algorithm>
#include <cstring>
#include <cstdio>
#include <cmath>
#include <queue>
#include <iostream>
#include <vector>
using namespace std;
#define mod 1000000001
typedef long long sint;
sint tot=1;
int bin[20],n,a[20][20],b[20],f[20][2048];
bool mark[100005];
int work(int x)
{
    memset(b,0,sizeof(b));
	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]*2;
		else 
            a[i][1]=n+1;
    }
    for(int i=1;i<=18;i++)
		for(int j=2;j<=11;j++)
        {
            if(a[i][j-1]*3<=n)
				a[i][j]=a[i][j-1]*3;
			else 
                a[i][j]=n+1;
        }
			
	for(int i=1;i<=18;i++)
		for(int j=1;j<=11;j++)
			if(a[i][j]<=n)
			{
				b[i]+=bin[j-1];
				mark[a[i][j]]=1;
			}
    memset(f,0,sizeof(f));
	f[0][0]=1;
	for(int i=0;i<18;i++)
		for(int x=0;x<=b[i];x++)
		{
			if(f[i][x])
				for(int y=0;y<=b[i+1];y++)
					if(((x&y)==0)&&((y&(y>>1))==0))
						f[i+1][y]=(f[i][x]+f[i+1][y])%mod;
		}
	return f[18][0];
}
int main()
{
    scanf("%d",&n);
    bin[0]=1;
    for(int i=1;i<20;i++)
        bin[i]=bin[i-1]<<1;
	for(int i=1;i<=n;i++)
		if(!mark[i])
            tot=(tot*work(i))%mod;
	printf("%d",tot);
    return 0;
}