记录编号 463795 评测结果 AAAAAAAAAA
题目名称 [HNOI 2012]集合选数 最终得分 100
用户昵称 Gravatarwumingshi 是否通过 通过
代码语言 C++ 运行时间 0.042 s
提交时间 2017-10-24 19:39:11 内存使用 2.31 MiB
显示代码纯文本
#include<cstdio>
#include<cstring>
#define N 100005
#define mod 1000000001
int a[1005][500];
int f[2][1000];
bool b[N];
int n,tot,ans=1;
inline void find(int x)
{
	tot=0,a[0][0]=1;
	for(int i=x;i<=n;i<<=1)
	{
		tot++,a[tot][0]=0;
		for(int j=i;j<=n;j*=3)
		  a[tot][0]++,b[j]=1;
		int top=(1<<a[tot][0]);a[tot][0]=0;
		for(int j=0;j<top;j++)
		  if(!(j&(j>>1))) a[tot][++a[tot][0]]=j;
	}
}
inline void dp()
{
	memset(f,0,sizeof(f));
	f[0][1]=1;
	for(int i=1;i<=tot;i++)
	{
		memset(f[i&1],0,sizeof(f[i&1]));
		for(int j=1;j<=a[i-1][0];j++)
		  for(int k=1;k<=a[i][0];k++)
			if(!(a[i-1][j]&a[i][k])) f[i&1][k]=(f[i&1][k]+f[~i&1][j])%mod;
	}
}
inline void solve()
{
	for(int i=1;i<=n;i++)
	  if(!b[i])
	  {
		find(i),dp();
		int tmp=0;
		for(int j=1;j<=a[tot][0];j++)
		  tmp=(tmp+f[tot&1][j])%mod;
		ans=1ll*ans*tmp%mod;
	  }
}
int main()
{
	freopen("bzoj_2734.in","r",stdin);
	freopen("bzoj_2734.out","w",stdout);
	scanf("%d",&n);
	solve();
	printf("%d",ans);
	return 0;
}