1. 程式人生 > >【BZOJ 1801】【AHOI 2009】中國象棋(遞推DP)

【BZOJ 1801】【AHOI 2009】中國象棋(遞推DP)

(誰告訴我是狀壓的????)

一行/列最多能放三個炮

30分是讓你爆搜的

50分是讓你狀壓的

假設有8列

dp[i][j]表示到第i行且前面列的狀態為j (沒有/有一個/有兩個炮)

然後刷表

100分其實就是在50分的基礎上修改了一點

因為發現 並不需要準確的知道之前擺放情況到底是什麼樣子的

於是dp[i][j][k]表示放了前i行,有j列是有1個棋子,有k列有兩個棋子

同樣的刷表

#include<bits/stdc++.h>
#define N 105
#define ll long long
using namespace std;
const int mod=9999973;
ll n,m,dp[N][N][N];
int main()
{
    cin>>n>>m;
    dp[0][0][0]=1;
    for(int i=0;i<=n-1;i++)
    {
        for(int j=0;j<=m;j++)
        {
            for(int k=0;k+j<=m;k++)
            {
                dp[i+1][j][k]=(dp[i+1][j][k]+dp[i][j][k])%mod;
                //放一個在空的那列
                if(j+1<=m) dp[i+1][j+1][k]=(dp[i+1][j+1][k]+dp[i][j][k]*(m-j-k))%mod; 
                //放一個在放了一個的那列
                if(k+1<=m) dp[i+1][j-1][k+1]=(dp[i+1][j-1][k+1]+dp[i][j][k]*j)%mod;
                //放兩個在兩個放了一個的那兩列
                if(k+2<=m) dp[i+1][j-2][k+2]=(dp[i+1][j-2][k+2]+dp[i][j][k]*(j*(j-1)/2))%mod; 
                //放兩個在兩列空的
                if(j+2<=m) dp[i+1][j+2][k]=(dp[i+1][j+2][k]+dp[i][j][k]*((m-j-k)*(m-j-k-1)/2))%mod;
                //放兩個在一個空的 一個放了一個的
                if(k+1<=m) dp[i+1][j][k+1]=(dp[i+1][j][k+1]+dp[i][j][k]*j*(m-j-k))%mod; 
            }
        }
    }
    ll ans=0;
    for(int i=0;i<=m;i++)
        for(int j=0;j+i<=m;j++)
            ans=(ans+dp[n][i][j])%mod;
    cout<<ans;
    return 0;
}