luogu P2051 [AHOI2009]中国象棋

Posted Excim

tags:

篇首语:本文由小常识网(cha138.com)小编为大家整理,主要介绍了luogu P2051 [AHOI2009]中国象棋相关的知识,希望对你有一定的参考价值。

原题链接:https://www.luogu.org/problem/show?pid=2051

本题是一道DP,思维难度主要是在多种状态的考虑上。

对于每一行,最多放两个炮,所以,对于每一行,只会有放一个炮,放两个炮或者不放炮的三种情况

以f[i][j][k]表示前i行中,有j列有1个炮,k列有两个炮,则有m-j-k列没有炮,a[i]记录C(i,2)。

1.不放炮:直接将本行的方案数加到下一行即可。

2.放一个炮,放在一个空列上。要注意空行有m-j-k个,所以要乘上m-j-k再加到下一行

3.放一个炮,放在一个有一个炮的列上。

尤其要注意,此时可能没有任何一列只有一个炮的,所以要特判能否进行此次转移。

4.放两个炮,都放在空列上。空列有(m-j-k)个,而从中选两个放炮,就有C(m-i-j,2)种选法。

5.放两个炮,分别放在一个空列上和一个只有一个炮的列上,其实相当于在一个空列上放了两个炮,空列有(m-j-k)个,只有一个炮的列有j个。

6.放两个炮,都放在只有一个炮的列上,首先特判j>=2时才进行转移,有C(j,2)种选法。

由此,我们就能得到六个转移方程(未特判):

 

f[i+1][j][k]+=f[i][j][k];
f[i+1][j+1][k]+=((f[i][j][k])*(m-k-j))%c;
f[i+1][j-1][k+1]+=(f[i][j][k]*j)%c;
f[i+1][j+2][k]+=(f[i][j][k]*a[m-j-k])%c;
f[i+1][j][k+1]+=(f[i][j][k]*(m-j-k)*j)%c;
f[i+1][j-2][k+2]+=(f[i][j][k]*a[j])%c;

 

转移完成后统计答案,f[i][j][k]中j+k<=m的方案符合要求,加入答案中(j+k>m的情况显然是不合实际的)

 

#include<cstdio>
const int c=9999973; 
long long ans,f[105][105][105];//前i行j 1,k 2 m-i-j 0
int n,m,a[105];
int main()
{
    scanf("%d %d",&n,&m);
    for(int i=2;i<=100;i++) a[i]=(i*(i-1))/2;
    f[0][0][0]=1;
    for(int i=0;i<n;i++)
    {
        for(int j=0;j<=m;j++)
        {
            for(int k=0;k<=m;k++)
            {
                f[i][j][k]%=c;
                f[i+1][j][k]+=f[i][j][k];
                f[i+1][j+1][k]+=((f[i][j][k])*(m-k-j))%c;
                if(j>0) f[i+1][j-1][k+1]+=(f[i][j][k]*j)%c;
                f[i+1][j+2][k]+=(f[i][j][k]*a[m-j-k])%c;
                f[i+1][j][k+1]+=(f[i][j][k]*(m-j-k)*j)%c;
                if(j>1) f[i+1][j-2][k+2]+=(f[i][j][k]*a[j])%c;
            }
        }
    }
    for(int i=0;i<=m;i++)
    {
        for(int j=0;i+j<=m;j++) ans=(ans+f[n][i][j])%c;
    }
    printf("%lld",ans);
    return 0;
}

 

以上是关于luogu P2051 [AHOI2009]中国象棋的主要内容,如果未能解决你的问题,请参考以下文章

luogu P2051 [AHOI2009]中国象棋

luogu P2051 [AHOI2009]中国象棋

luogu P2051 [AHOI2009]中国象棋 dp 状态压缩+容斥

Luogu P2051 [AHOI2009]中国象棋

Luogu P2051 [AHOI2009]中国象棋 //DP

Luogu P2051 [AHOI2009]中国象棋 | dp