[DP] BZOJ1801: [Ahoi2009]chess 中国象棋

最新推荐文章于 2018-10-31 14:06:22 发布

Lynstery

最新推荐文章于 2018-10-31 14:06:22 发布

阅读量298

点赞数

分类专栏： --------------DP--------------

本文链接：https://blog.csdn.net/CHHNZ/article/details/78413766

版权

--------------DP-------------- 专栏收录该内容

54 篇文章 0 订阅

订阅专栏

大水题。显然限制就是同行同列个数不能超过2。 $f_{i,j,k}$ 表示前 $i$ 行，放了 $0$ 个的列数是 $j$ ，放了 $1$ 个的列数是 k <script id="MathJax-Element-95" type="math/tex">k</script> 。就好了…

#include<cstdio>
#include<algorithm>
using namespace std;
const int maxn=105,MOD=9999973;
typedef long long LL;
int n,m;
LL f[maxn][maxn][maxn],ans;
int main(){
    freopen("bzoj1801.in","r",stdin);
    freopen("bzoj1801.out","w",stdout);
    scanf("%d%d",&n,&m);
    f[0][m][0]=1;
    for(int i=0;i<=n-1;i++)
     for(int j=0;j<=m;j++)
      for(int k=0;j+k<=m;k++) if(f[i][j][k]){
        (f[i+1][j][k]+=f[i][j][k])%=MOD;
        if(j) (f[i+1][j-1][k+1]+=f[i][j][k]*j%MOD)%=MOD;
        if(k) (f[i+1][j][k-1]+=f[i][j][k]*k%MOD)%=MOD;
        if(j>=2) (f[i+1][j-2][k+2]+=f[i][j][k]*j*(j-1)/2%MOD)%=MOD;
        if(k>=2) (f[i+1][j][k-2]+=f[i][j][k]*k*(k-1)/2%MOD)%=MOD;
        if(j&&k) (f[i+1][j-1][k]+=f[i][j][k]*j*k%MOD)%=MOD;
      }
    for(int j=0;j<=m;j++)
     for(int k=0;j+k<=m;k++) (ans+=f[n][j][k])%=MOD;
    printf("%lld\n",ans);
    return 0; 
}