题目链接
分析
和Ahoi2009中国象棋一样,注意取mod的不同,Ahoi2009中国象棋题解
上代码
#include<bits/stdc++.h>
using namespace std;
long long dp[103][103][103],mod=999983,ans=0,n,m;
int C(int x){
return x*(x-1)/2;
}
int main(){
int n,m;
scanf("%d%d",&n,&m);
memset(dp,0,sizeof(dp));
dp[0][0][0]=1;
for(int i=1;i<=n;i++)
for(int j=0;j<=m;j++)
for(int k=0;j+k<=m;k++){
dp[i][j][k]=(dp[i][j][k]+dp[i-1][j][k])%mod;
if(j-1>=0) dp[i][j][k]=(dp[i][j][k]+dp[i-1][j-1][k]*(m-j+1-k))%mod;
if(j+1<=m&&k-1>=0) dp[i][j][k]=(dp[i][j][k] + dp[i-1][j+1][k-1]*(j+1))%mod;
if(j-2>=0) dp[i][j][k]=(dp[i][j][k]+dp[i-1][j-2][k]*C(m-k-j+2))%mod;
if(k-1>=0) dp[i][j][k]=(dp[i][j][k]+dp[i-1][j][k-1]*(m-j-k+1)*j)%mod;
if(j+2<=m&&k-2>=0) dp[i][j][k]=(dp[i][j][k]+dp[i-1][j+2][k-2]*C(j+2))%mod;
}
for(int i =0;i<=m;i++) for(int j=0;j+i<=m;j++) ans=(ans+dp[n][i][j])%mod;
printf("%lld\n",ans);
return 0;
}