bzoj4806炮

题目链接

分析

和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;
}