luogu P2051 [AHOI2009]中国象棋

luogu P2051 [AHOI2009]中国象棋
简单dp?
f[i][j][k]表示dp了前i行，有j列放了1个棋子，有k列放了2个棋子的方案数
转移自己手写
（太长了）
注意：转移时记得排列组合
f[0][0][0]初始化为1！
std:

#include<bits/stdc++.h>
using namespace std;
#define mod 9999973 
const int N=105;
typedef long long ll;
int n,m;
ll f[N][N][N];//前i行放了1个的有j列，放了2个的有k列情况总数 
ll ans=0;
ll C(int x,int y){
    if(y==1){
        return (ll)x;
    }
    else{
        return (ll)x*(x-1)/2;
    }
}
int main(){
    memset(f,0,sizeof(f));
    f[0][0][0]=1;
    scanf("%d%d",&n,&m);
    for(int i=0;i<n;i++){
        for(int j=0;j<=m;j++){
            for(int k=0;k+j<=m;k++){
                //printf("%d %d %d %d\n",i,j,k,f[i][j][k]);
                if(k+j>m) break;
                int l=m-j-k;//0个的列数
                //放0个的方案数
                f[i+1][j][k]+=f[i][j][k];
                f[i+1][j][k]%=mod;
                //放一个
                //再有1个的地方
                if(j>0) {
                    f[i+1][j-1][k+1]+=f[i][j][k]*C(j,1);
                    f[i+1][j-1][k+1]%=mod;
                }
                //0个
                if(l>0) {
                    f[i+1][j+1][k]+=f[i][j][k]*C(l,1);
                    f[i+1][j+1][k]%=mod;
                } 
                //在不同列（全0） 
                if(l>1) {
                    f[i+1][j+2][k]+=f[i][j][k]*C(l,2);
                    f[i+1][j+2][k]%=mod;
                }
                //一1一0
                if(j>0&&l>0) {
                    f[i+1][j][k+1]+=f[i][j][k]*C(l,1)*C(j,1);
                    f[i+1][j][k+1]%=mod;
                }
                //全1
                if(j>1) {
                    f[i+1][j-2][k+2]+=f[i][j][k]*C(j,2); 
                    f[i+1][j-2][k+2]%=mod;
                }
            } 
        } 
    }
    for(int i=0;i<=m;i++){
        for(int j=0;j+i<=m;j++){
            ans+=f[n][i][j];
            ans%=mod;
        }
    }
    printf("%lld\n",ans);
    return 0;
}