最近智商有点不在线。其实一直不在线。
先是想用$f[i][j][k][0/1/2]$表示摆了i行时有j列空着,k列有了一个炮,且当下摆了0/1/2个的状态,转移方程写的出来但是极其繁琐。于是又设法听取评讲者题解修改状态,最后的012完全可以删去。那么仍可以表示这一行那些列摆过1个,那些列摆过0个的种类。转移时分类即可。
$f[i][j][k]+=f[i-1][j][k]$ 什么都不摆
$f[i][j][k]+=(j+1)*f[i-1][j+1][k-1]$ 摆1个炮
$f[i][j][k]+=(k+1)*f[i-1][j][k+1]$ 摆1个炮
$f[i][j][k]+=(j+1)*(j+2)/2*f[i-1][j+2][k-2]$ 摆两个炮,下同
$f[i][j][k]+=(k+1)*(k+2)/2*f[i-1][j][k+2]$
$f[i][j][k]+=k*(j+1)*f[i-1][j+1][k]$
注意边界就行。
1 #include<bits/stdc++.h> 2 using namespace std; 3 typedef long long ll; 4 typedef pair<int,int> pii; 5 template<typename T>inline char MIN(T&A,T B){return A>B?A=B,1:0;} 6 template<typename T>inline char MAX(T&A,T B){return A<B?A=B,1:0;} 7 template<typename T>inline T _min(T A,T B){return A<B?A:B;} 8 template<typename T>inline T _max(T A,T B){return A>B?A:B;} 9 template<typename T>inline T read(T&x){ 10 x=0;int f=0;char c;while(!isdigit(c=getchar()))if(c==‘-‘)f=1; 11 while(isdigit(c))x=(x<<1)+(x<<3)+(c^48),c=getchar();return f?x=-x:x; 12 } 13 const int N=100+7,P=9999973; 14 int n,m,ans,f[N][N][N]; 15 16 inline void inc(int&x,int y){(x+=y)>=P?x-=P:1;} 17 18 int main(){//freopen("tmp.in","r",stdin);freopen("tmp.out","w",stdout); 19 read(n),read(m);if(m<2){printf("%d\n",1+n+(n)*(n-1)/2);return 0;} 20 f[1][m][0]=1,f[1][m-1][1]=m,f[1][m-2][2]=m*(m-1)/2; 21 for(register int i=2;i<=n;++i) 22 for(register int j=0;j<=m;++j) 23 for(register int k=0;k<=m-j;++k){ 24 inc(f[i][j][k],f[i-1][j][k]); 25 if(k)inc(f[i][j][k],(j+1)*f[i-1][j+1][k-1]%P); 26 inc(f[i][j][k],(k+1)*f[i-1][j][k+1]%P); 27 if(k>=2)inc(f[i][j][k],((j+1)*(j+2)>>1)*1ll*f[i-1][j+2][k-2]%P); 28 inc(f[i][j][k],((k+1)*(k+2)>>1)*1ll*f[i-1][j][k+2]%P); 29 inc(f[i][j][k],1ll*k*(j+1)*f[i-1][j+1][k]%P); 30 } 31 for(register int i=0;i<=m;++i)for(register int j=0;j<=m-i;++j)inc(ans,f[n][i][j]); 32 printf("%d\n",ans); 33 return 0; 34 }
原文地址:https://www.cnblogs.com/saigyouji-yuyuko/p/10358509.html
时间: 2024-10-17 11:19:35