题目大意,求一个数列的方案数,这个数列满足,只由1到k之间的数字组成,而且当ai=aj时,i和j的差值至少为ai
之前那个50分的是因为压得方式不同
设f[i,i1,i2...ik]表示长度为i,数字1上一次出现位置在第i位的前i1位,2上次在i的前i2位……有多少种方案
当i1>=1时,i1=0。当i2>=2时,i2=0……
通过实践发现,k=7时,合法的i1,i2……ik状态最多900多种
现在来考虑转移,枚举第i位是1..k,如果第i位可以填i,就是一种转移。
通过实践发现,k=7时,合法的转移最多3000多种!
这样,我们可以预处理出f[i]到f[i+1]的每种转移,然后dp
k=7时时间复杂度o(3000*n),可以承受
显然我们需要两次dfs一次暴力转移一次暴力初始值
状态显然用7进制压缩即可
合法的状态满足没有相同数字【意外的没想到结果状态数一直是5040= =写了一天的我满脸大汗
/* ***********************************************
Author :BPM136
Created Time :2016-3-18 17:12:25
File Name :E:\BPM136\BPM\DP\SSL2418GG.cpp
************************************************ */
#include<stdio.h>
#include<string.h>
#include<iostream>
#include<algorithm>
#include<vector>
#include<queue>
#include<set>
#include<map>
#include<bitset>
#include<string>
#include<math.h>
#include<stdlib.h>
#include<time.h>
#define LL long long
#define fo(i,a,b) for(int i=a;i<=b;i++)
#define dfo(i,a,b) for(int i=a;i>=b;i--)
#define efo(i,x) for(int i=last[x];i!=0;i=e[i].next)
using namespace std;
LL read()
{
LL d=0,f=1;char s=getchar();
while(s<48||s>57){if(s==45)f=-1;s=getchar();}
while(s>=48&&s<=57){d=d*10+s-48;s=getchar();}
return d*f;
}
#define inf 910109
int f[2][900000];
struct S
{
int x,y;
}s[20000];
int n,m,num=0;
int tt[10],c[10],aa[10],ma=0,mi=0;
int status[10000];
bool tot[900000];
int po7[9];
void add(int x,int y)
{
//cout<<x<<' '<<y<<' '<<num+1<<endl;
if(tot[x]==0)
{
status[++status[0]]=x;
tot[x]=1;
}
s[++num].x=x;s[num].y=y;
mi=min(mi,min(x,y));
ma=max(ma,max(x,y));
}
void check()
{
cout<<"CCCC:";
fo(i,1,m-1)cout<<c[i]<<' ';
cout<<endl;
cout<<"AAAA:";
fo(i,1,m)cout<<aa[i]<<' ';
cout<<endl;
}
void dfs(int k)
{
if(k==m+1)
{
int x=0,y=0;
memset(tt,0,sizeof(tt));
fo(i,1,m)tt[c[i]]++;
fo(i,1,m)
{
if(tt[i]>1)return ;
x+=po7[i-1]*c[i];
if(c[i]!=0)y+=((c[i]+1)%i)*po7[i-1];
}
fo(i,1,m)
if(c[i]==0)
{
if(i==1)add(x,y);else
add(x,y+po7[i-1]);
}
return ;
}
fo(i,0,k-1)
{
c[k]=i;
dfs(k+1);
c[k]=0;
}
}
void dfs2(int k)
{
if(k==m)
{
memset(aa,0,sizeof(aa));
for(int i=m-1,j=2;i>=1;i--,j++)
{
if(c[i]>=j&&aa[c[i]]==0)
{
aa[c[i]]=j-1;
}
}
// check();
int x=0;
fo(i,1,m)x+=po7[i-1]*aa[i];
f[1-m%2][x]++;
// cout<<x<<endl;
return ;
}
bitset<10>a;a.reset();c[k]=1;
dfs2(k+1);c[k]=0;
for(int i=k-1,j=2;j<=m;i--,j++)
{
if(i>0)a[c[i]]=1;
if(a[j]==0)
{
c[k]=j;
dfs2(k+1);
c[k]=0;
}
}
}
void prework()
{
po7[0]=1;
fo(i,1,8)po7[i]=po7[i-1]*7;
memset(tot,0,sizeof(tot));
memset(c,0,sizeof(c));
memset(tot,0,sizeof(tot));
dfs(1);
memset(f,0,sizeof(f));
memset(c,0,sizeof(c));
dfs2(1);
}
int main()
{
//freopen(".in","r",stdin);
//freopen(".out","w",stdout);
n=read(),m=read();
prework();
fo(i,m,n)
{
int t=i%2;
fo(j,1,status[0])f[t][status[j]]=0;
fo(j,1,num)
{
f[t][s[j].y]+=f[1-t][s[j].x];
if(f[t][s[j].y]>=inf)
{
f[t][s[j].y]%=inf;
}
}
}
int ans=0,t=n%2;
fo(i,1,status[0])
{
ans+=f[t][status[i]];
if(ans>=inf)ans%=inf;
}
cout<<ans<<endl;
return 0;
}