题目
某大佬的题解
以下内容改编自大佬的题解

题解：

dp[i][j][k]表示在前i行中有j列已经有2个，有k列已经有1个，那么0的个数就是m-j-k个
不放：直接加上dp[i-1][j][k]
放一个：1.把一列0变成1；2.把一列1变成2
放两个：1.把两列0变成1；2.把两列1变成2；3.把一列0变成1 再把另外一列1变成2
关于组合数：
放一个：直接乘要改变的数的个数（0或1）
放两个：前两个直接C（n，2） （n为可选个数）
最后一个根据乘法原理，先乘上原来0的个数，再乘上原来1的个数

标程：

#include<bits/stdc++.h>
using namespace std;
const int M=9999973;
long long t,f[103][103][103];
int i,j,k,n,m;
int main(){
    cin>>n>>m;
    f[0][0][0]=1;
    for (i=1;i<=n;i++)
        for (j=0;j<=m;j++)
            for (k=0;k<=m-j;k++){
                t=f[i-1][j][k];//不放
                if (j) t+=f[i-1][j-1][k+1]*(k+1);//放1个，1变2
                if (k) t+=f[i-1][j][k-1]*(m-j-k+1);//放1个，0变1
                if (j && k) t+=f[i-1][j-1][k]*k*(m-j-k+1);//放两个，0变1,1变2
                if (j>1) t+=f[i-1][j-2][k+2]*(k+2)*(k+1)/2;//放两个，两列1变2
                if (k>1) t+=f[i-1][j][k-2]*(m-j-k+2)*(m-j-k+1)/2;//放两个，两列0变1
                f[i][j][k]=t%M;
            }
    t=0;
    for (i=0;i<=m;i++)
        for (j=0;j<=m-i;j++) t=(t+f[n][i][j])%M;
    cout<<t;
}