bzoj 4806: 炮dp

Posted lokiii

tags:

篇首语:本文由小常识网(cha138.com)小编为大家整理,主要介绍了bzoj 4806: 炮dp相关的知识,希望对你有一定的参考价值。

同1801
注意到一行只能放012个炮,我们只需要知道列的状态,不用状压行
所以设f[i][j][k]表示前i行有j列有1个炮,有k列有2个炮的方案数
然后分情况讨论转移就行了

#include<cstdio>
#include<iostream>
using namespace std;
const int N=105,mod=999983;
long long n,m,f[N][N][N],ans;
int main()
{
    scanf("%lld%lld",&n,&m);
    f[0][0][0]=1;
    for(int i=1;i<=n;i++)
        for(int j=0;j<=m;j++)
            for(int k=0;k<=m-j;k++)
            {
                f[i][j][k]=f[i-1][j][k];
                if(j>=1)
                    f[i][j][k]=(f[i][j][k]+(f[i-1][j-1][k]*(m-j+1ll-k))%mod)%mod;
                if(k>=1&&j+1<=m)
                    f[i][j][k]=(f[i][j][k]+(f[i-1][j+1][k-1]*(j+1ll))%mod)%mod;
                if(j>=2)
                    f[i][j][k]=(f[i][j][k]+(f[i-1][j-2][k]*((m-j+2ll-k)*(m-j+1ll-k)/2)%mod)%mod)%mod;
                if(k>=2&&j+2<=m)
                    f[i][j][k]=(f[i][j][k]+(f[i-1][j+2][k-2]*((long long)(j+2ll)*(long long)(j+1ll)/2)%mod)%mod)%mod;
                if(j>=1&&k>=1)
                    f[i][j][k]=(f[i][j][k]+(f[i-1][j][k-1]*(long long)j%mod*(m-j-k+1)%mod)%mod)%mod;
                if(i==n)
                    ans=(ans+f[i][j][k])%mod;
            }
    printf("%lld
",ans);
    return 0;
}



以上是关于bzoj 4806: 炮dp的主要内容,如果未能解决你的问题,请参考以下文章

bzoj4806 炮

BZOJ 4806 - 4809 象棋四题

bzoj4806~bzoj4809 象棋四连发 DP-高精度-匈牙利算法-dfs

bzoj1801: [Ahoi2009]chess 中国象棋 [dp]

BZOJ1801 [Ahoi2009]chess 中国象棋 dp

炮(棋盘DP)