bzoj4806 炮

这道题其实是同bzoj1801的，
只是模数不同
用dp[i][j][k]表示第i行有j列有没有子，k列一个子
再枚举一下是从那列转移来的就行了

#include <cstdio>
#include <cstring>
#include <iostream>
#include <algorithm>
using namespace std;
const int mod=999983;
inline int read()
{
    int x=0,f=1,ch=getchar();
    while(ch<'0'||ch>'9'){if(ch=='-'){f=-1;}ch=getchar();}
    while(ch>='0'&&ch<='9'){x=x*10+ch-'0';ch=getchar();}
    return x*f;
}
int n,m;
int dp[105][105][105];
int C[105][3];
int dfs(int x,int cnt1,int cnt2)
{
    if(x>n) return 1;
    if(dp[x][cnt1][cnt2]) return dp[x][cnt1][cnt2];
    dp[x][cnt1][cnt2]=dfs(x+1,cnt1,cnt2);
    if(cnt1>1) dp[x][cnt1][cnt2]=(dp[x][cnt1][cnt2]+1ll*dfs(x+1,cnt1-2,cnt2+2)*C[cnt1][2]%mod)%mod;
    if(cnt1>0) dp[x][cnt1][cnt2]=(dp[x][cnt1][cnt2]+1ll*dfs(x+1,cnt1-1,cnt2+1)*C[cnt1][1]%mod)%mod;
    if(cnt2>1) dp[x][cnt1][cnt2]=(dp[x][cnt1][cnt2]+1ll*dfs(x+1,cnt1,cnt2-2)*C[cnt2][2]%mod)%mod;
    if(cnt2>0) dp[x][cnt1][cnt2]=(dp[x][cnt1][cnt2]+1ll*dfs(x+1,cnt1,cnt2-1)*C[cnt2][1]%mod)%mod;
    if(cnt1&&cnt2) dp[x][cnt1][cnt2]=(dp[x][cnt1][cnt2]+1ll*dfs(x+1,cnt1-1,cnt2)*C[cnt1][1]*C[cnt2][1])%mod;
    return dp[x][cnt1][cnt2];
}
int main()
{
    n=read(),m=read();
    for(int i=0;i<=m;i++) C[i][0]=1;
    for(int i=1;i<=m;i++) C[i][1]=i;
    for(int i=2;i<=m;i++) C[i][2]=(C[i-1][1]+C[i-1][2])%mod;
    printf("%d\n",dfs(1,m,0));
    /*for(int i=1;i<=n;i++)
        for(int j=0;j<=m;j++)
            for(int k=0;j+k<=m;k++)
                cout<<i<<"\t"<<j<<"\t"<<k<<"\t"<<dp[i][j][k]<<endl;
                */
    return 0;
}