洛谷2051 [AHOI2009]中国象棋

题目链接

题意概述：n行m列棋盘放若干个棋子每行每列最多两个求方案总数，答案对9999973取模。

可以比较容易看出这是个dp，设f[i][j][k]表示前i行j列放1个棋子k列放2个棋子的方案总数。转移时分类讨论+计数就好了。分类讨论比较麻烦需要仔细考虑一下。

吐槽：我把循环里的m写成n还有50，查了半天（真的半天）取模QAQ

#include<cstring>
#include<iostream>
#include<cctype>
#include<cstdio>
#include<algorithm>
#define now f[i][j][k]
#define ll unsigned long long
using namespace std;
inline int read()
{
    register int X=0;
    register char ch=0;
    bool flag=0;
    for(; !isdigit(ch); ch=getchar()) if(ch=='-') flag=1;
    for(; isdigit(ch); ch=getchar()) X=(X<<3)+(X<<1)+ch-'0';
    return (flag ? -X : X);
}
int N=101;
const ll mod=9999973;
int n,m;
ll f[201][201][201],ans;
int C(int x)
{
    return x*(x-1)/2;
}
int main()
{
    n=read(),m=read(),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(now)
                {
                    (f[i+1][j][k]+=now)%=mod;
                    if(m-j-k > 0) (f[i+1][j+1][k]+=now*(m-j-k))%=mod;
                    if(j > 0) (f[i+1][j-1][k+1]+=now*j)%=mod;
                    if(m-j-k > 1) (f[i+1][j+2][k]+=now*(C(m-j-k)%mod))%=mod;
                    if(j > 0 && m-j-k > 0) (f[i+1][j][k+1]+=now*j%mod*(m-j-k))%=mod;
                    if(j > 1) (f[i+1][j-2][k+2]+=now*(C(j)%mod))%=mod;
                }
            }
    for(int i=0; i<=m; i++)
        for(int j=0; j+i<=m; j++)
            (ans+=f[n][i][j])%=mod;
    printf("%llu\n",ans);
}