1801: [Ahoi2009]chess 中国象棋

#include<iostream>
#include<cstdio>
using namespace std;
inline long long c(long long x){return x*(x-1)/2;}
int n,m;
long long ans,f[101][101][101];
int main(){
scanf("%d%d",&n,&m);
f[1][0][0]=1;f[1][1][0]=m;f[1][2][0]=c(m);
for(int i=2;i<=n;i++)
for(int j=0;j<=m;j++)
for(int k=0;k<=m-j;k++){
f[i][j][k]=f[i-1][j][k];
if(j)f[i][j][k]+=f[i-1][j-1][k]*(m-j-k+1);
if(j>1)f[i][j][k]+=f[i-1][j-2][k]*c(m-j-k+2);
if(k&&j<m)f[i][j][k]+=f[i-1][j+1][k-1]*(j+1);
if(j&&k)f[i][j][k]+=f[i-1][j][k-1]*(m-j-k+1)*j;
if(k>1&&j<m-1)f[i][j][k]+=f[i-1][j+2][k-2]*c(j+2);
f[i][j][k]%=9999973;
}
for(int i=0;i<=m;i++)
for(int j=0;j<=m-i;j++)
ans+=f
[i][j];
printf("%lld",ans%9999973);
return 0;
}