BZOJ 1801 中国象棋（DP）

int n,m;
i64 f[N][N][N];


void up(i64 &x,i64 y)
{
    x+=y;
    x%=mod;
}


i64 C(int x)
{
    return x*(x-1)/2;
}


int main()
{
    RD(n,m);
    f[1][0][0]=1;
    f[1][1][0]=m;
    f[1][2][0]=m*(m-1)/2;
    int i,j,k;


    for(i=2;i<=n;i++)
    {
        for(j=0;j<=m;j++) for(k=0;k<=m;k++)
        {
            if(j+k>m) continue;


            f[i][j][k]=f[i-1][j][k];
            
            if(j>0) up(f[i][j][k],f[i-1][j-1][k]*(m-j-k+1));
            if(k>0) up(f[i][j][k],f[i-1][j+1][k-1]*(j+1));
            
            if(j>=2) up(f[i][j][k],f[i-1][j-2][k]*C(m-j-k+2));
            if(k>=1) up(f[i][j][k],f[i-1][j][k-1]*(m-j-k+1)*j);
            if(k>=2) up(f[i][j][k],f[i-1][j+2][k-2]*C(j+2));
        }
    }
    i64 ans=0;
    for(j=0;j<=m;j++) for(k=0;k<=m;k++)
    {
        up(ans,f[n][j][k]);
    }
    PR(ans);
}