简单思路:dp
用f[i][j][k]表示前i行有j列放了一个子,有k列放了两个子
转移,只要枚举每一行的放置情况即可,有:不放、在空列放一个,在有一个子的列放一个,在两个空列各放一个,在两个一子列各放一个,在空列和一子列各放一个。
枚举ijk即可,最后求和
注意:逆元只有在进行除法前需要取模时才需要,否则不需要。
代码如下:
#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;
const long long mod=9999973;
int n,m;
long long f[105][105][105];
int main(){
scanf("%d%d",&n,&m);
memset(f,0,sizeof(f));
f[0][0][0]=1;
for(int i=0;i<n;i++){
for(int j=0;j<=m;j++){
for(int k=0;(k+j)<=m;k++){
if(f[i][j][k]){
f[i+1][j][k]=(f[i+1][j][k]+f[i][j][k])%mod;
if(m-j-k>=1){
f[i+1][j+1][k]=(f[i+1][j+1][k]+f[i][j][k]*(m-j-k))%mod;
}
if(j>=1){
f[i+1][j-1][k+1]=(f[i+1][j-1][k+1]+f[i][j][k]*j)%mod;
}
if(j>=1 && m-j-k>=1){
f[i+1][j][k+1]=(f[i+1][j][k+1]+f[i][j][k]*j*(m-j-k))%mod;
}
if(m-j-k>=2){
f[i+1][j+2][k]=(f[i+1][j+2][k]+f[i][j][k]*(m-j-k)*(m-j-k-1)/2)%mod;
}
if(j>=2){
f[i+1][j-2][k+2]=(f[i+1][j-2][k+2]+f[i][j][k]*(j-1)*j/2)%mod;
}
}
}
}
}
long long ans=0;
for(long long j=0;j<=m;j++){
for(long long k=0;k+j<=m;k++){
ans=(ans+f[n][j][k])%mod;
}
}
printf("%lld",ans);
return 0;
}