不是很懂这类题有什么意思。
可以正着做也可以反着做,反着做状态要少一点,常数小一点。
尝试状压DP,发现复杂度是 O ( n k 2 k ) O(nk2^k) O(nk2k),感觉很方,测了一下极限数据,发现非0状态只有 7000 7000 7000多个,于是只存有用状态即可。
n n n甚至可以开到 500 500 500。
如果 k k k再小一点,有效状态更少,可以上矩阵快速幂, n n n甚至可以开 1 e 9 1e9 1e9直接用转移矩阵的特征多项式艹过去。
好的我会记住把这个改成毒瘤题的
代码:
#include<bits/stdc++.h> #define ll long long #define re register #define cs const using std::cerr; using std::cout; cs int mod=998244353; inline int add(int a,int b){a+=b-mod;return a+(a>>31&mod);} inline int dec(int a,int b){a-=b;return a+(a>>31&mod);} inline int mul(int a,int b){ll r=(ll)a*b;return r>=mod?r%mod:r;} inline int po(int a,int b){ int r=1;for(;b;b>>=1,a=mul(a,a)) if(b&1)r=mul(r,a);return r; } inline void Inc(int &a,int b){a+=b-mod;a+=a>>31&mod;} inline void Dec(int &a,int b){a-=b;a+=a>>31&mod;} inline void Mul(int &a,int b){a=mul(a,b);} cs int N=1<<21|7; int n,k,lim,l,tmp,S; struct array{ int vis[N],val[N]; int st[N],tp; int &operator[](int o){ if(!vis[o]){ vis[o]=true; val[o]=0; st[++tp]=o; }return val[o]; } inline void clear(){ while(tp){ vis[st[tp]]=0; val[st[tp]]=0; --tp; } } }f; int st[N],val[N],tp; #undef zxyoi signed main(){ #ifdef zxyoi freopen("sequence.in","r",stdin); #else #ifndef ONLINE_JUDGE freopen("sequence.in","r",stdin);freopen("sequence.out","w",stdout); #endif #endif scanf("%d%d%d",&n,&k,&l);tmp=l%mod; lim=std::min(k-1,l);S=1<<k; l=l>k?(l-k)%mod:0;f[1]=1; for(int re i=1;i<=n;++i){ tp=f.tp; for(int re j=1;j<=tp;++j){ st[j]=f.st[j]; val[j]=f.val[st[j]]; }f.clear(); for(int re j=1;j<=tp;++j){ int s=st[j],vl=val[j]; if(l)Inc(f[s],mul(vl,l)); for(int re k=0;k<=lim;++k)if(!(S&(s<<k))) Inc(f[(s|(s<<k))&(S-1)],vl); } }int ans=0; for(int re i=1;i<=f.tp;++i)Inc(ans,f.val[f.st[i]]); cout<<dec(po(tmp+1,n),ans)<<"\n"; return 0; }