记录编号 |
463795 |
评测结果 |
AAAAAAAAAA |
题目名称 |
[HNOI 2012]集合选数 |
最终得分 |
100 |
用户昵称 |
wumingshi |
是否通过 |
通过 |
代码语言 |
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;
}