码迷,mamicode.com
首页 > 其他好文 > 详细

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

时间:2019-07-24 09:18:16      阅读:97      评论:0      收藏:0      [点我收藏+]

标签:条件   忽略   out   存在   状态   位置   col   中国   span   

hhy给我推的这道题,我不会,我去康了题解.....

题解第一个人说自己20min K一道省选DP,还一次AC

我想着哪个人这么强,往上一翻发现是__stdcall哥哥......也许这就是大佬的世界吧


50pts:

首先可以想到的是状压DP,因为一行一列要最多不超过2个棋子,

所以用3进制表示整个棋盘每列的放了几个,dp[棋盘状态],边界为dp[0]=1

刷表法比填表法好写DP方程,转移时,在0或者1的位置上加1构成新的棋盘状态,

最后统计所有棋盘状态就行了。


100pts:

在50分的DP中我们反思,什么地方浪费了空间?我们全部状压,存了每条列的顺序和状态

但是在这道题中顺序对答案不造成影响,于是我们可以只存在这个状态的列有多少条,而忽略顺序

于是正确的DP状态就来了,dp[摆到x行][有i列摆了1个][有j列摆了2个],

刷表转移,加法原理转移到下一个状态,乘法原理把当前状态和放法相乘,

最后统计答案

刚开始写填表写挂了,然后写刷表,刷表写完了改了填表

刷表法:

#include<bits/stdc++.h>
#define mod 9999973

using namespace std ;

const int MAXN = 110,MAXM = 110;
long long dp[MAXN][MAXM][MAXM],n,m;

int main(){
    cin>>n>>m;
    dp[0][0][0]=1;
    for(int x=0;x<n;x++){
        for(int i=0;i<=m;i++){
            for(int j=0;i+j<=m;j++){
                if(dp[x][i][j]){
                    (dp[x+1][i][j] += dp[x][i][j]) %= mod;
                    if(m-i-j>0)(dp[x+1][i+1][j] += dp[x][i][j] * (m-i-j)) %= mod;
                    if(i>=1)(dp[x+1][i-1][j+1] += dp[x][i][j] * i) %= mod;
                    if(m-i-j-1>0)(dp[x+1][i+2][j] += dp[x][i][j] * (m-i-j) * (m-i-j-1) / 2) %= mod;
                    if(m-i-j>=0 && i>=1)(dp[x+1][i][j+1] += dp[x][i][j] * (m-i-j) * i) %= mod;
                    if(i>=2)(dp[x+1][i-2][j+2] += dp[x][i][j] * (i-1) * i / 2) %= mod;
                }
            }
        }
    }
    long long ans = 0;
    for(int i=0;i<=m;i++){
        for(int j=0;j<=m;j++){
            if(i+j<=m){
                ans += dp[n][i][j];
                ans %= mod;
            }
        }
    }
    cout<<ans<<endl;
    return 0;
}

填表法:

#include<bits/stdc++.h>
#define mod 9999973

using namespace std ;

const int MAXN = 110,MAXM = 110;
long long dp[MAXN][MAXM][MAXM],n,m;

int main(){
    cin>>n>>m;
    dp[0][0][0]=1;
    for(int x=1;x<=n;x++){
        for(int i=0;i<=m;i++){
            for(int j=0;i+j<=m;j++){
                dp[x][i][j] += dp[x-1][i][j],dp[x][i][j] %= mod;//什么也不做 
                if(i>=1) dp[x][i][j] += dp[x-1][i-1][j] * (m-i-j+1),dp[x][i][j] %= mod;//放一个在一列0上;
                if(j>=1) dp[x][i][j] += dp[x-1][i+1][j-1] * (i+1),dp[x][i][j] %= mod;//放一个在一列1上;
                if(i>=2) dp[x][i][j] += dp[x-1][i-2][j] * (m-i-j+2) * (m-i-j+1) / 2,dp[x][i][j] %= mod;//放两个在2列0上
                if(j>=2) dp[x][i][j] += dp[x-1][i+2][j-2] * (i+2) * (i+1) / 2,dp[x][i][j] %= mod;//放两个在2列1上 
                if(j>=1 && i>=1) dp[x][i][j] += dp[x-1][i][j-1] * (m-i-j+1) * i,dp[x][i][j] %= mod;//放2个,一个在1,一个在2上 
            }
        }
    }
    long long ans = 0;
    for(int i=0;i<=m;i++){
        for(int j=0;i+j<=m;j++){
            ans += dp[n][i][j];
            ans %= mod;
        }
    }
    cout<<ans<<endl;
    return 0;
}

 状态设计是一个巧活,状态是获得答案的最简条件时,这个DP才是最好的。

TAG:SIN_XIII ⑨

 

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

标签:条件   忽略   out   存在   状态   位置   col   中国   span   

原文地址:https://www.cnblogs.com/SINXIII/p/11235592.html

(0)
(0)
   
举报
评论 一句话评论(0
登录后才能评论!
© 2014 mamicode.com 版权所有  联系我们:gaon5@hotmail.com
迷上了代码!