1. 程式人生 > >4806 炮 [dp]

4806 炮 [dp]

題意:MxN的棋盤中最多放多少個炮,使得互相不能被打。
dp解決,M,N的範圍為100, 不考慮狀壓。
狀態:dp[i][j][k]表示當前行為i,前i-1行有j列沒有炮,有k列有1個炮的方案數。
決策:
if(j-1>=0) dp[i][j][k] = (dp[i][j][k] + dp[i-1][j-1][k]*(m-(j-1)-k))%p; ///在沒炮的地方放一個
if(j+1 <= m && k >= 1) dp[i][j][k] = (dp[i][j][k] + dp[i-1][j+1][k-1]*(j+1))%p; ///在一個炮的地方放一個
if(j >= 2) dp[i][j][k] = (dp[i][j][k] + dp[i-1][j-2][k]*C(m-k-(j-2)))%p; ///在沒炮的地方放兩個
if(k >= 1 && j>= 1) dp[i][j][k] = (dp[i][j][k] + dp[i-1][j][k-1](m-j-(k-1))

(j))%p; ///在沒炮的地方和一個炮的地方分別放一個
if(j+2 <= m && k >= 2) dp[i][j][k] = (dp[i][j][k] + dp[i-1][j+2][k-2]*C(j+2))%p; ///在兩個一個炮的地方分別放一個炮

邊界條件:dp[0][0][0] = 1, else = 0

#include<iostream>
#include<string>
#include<cstdio>
#include<cstring>
#include<bitset>
#include<algorithm> #include<map> #include<set> #include<queue> #include<vector> #include<cstdlib> #include<list> #include<stack> #include<cmath> using namespace std; typedef long long LL; const int p = 999983; int n, m; LL dp[105][105][105], ans; LL C(LL x) { return
x*(x-1)/2; } int main() { ios::sync_with_stdio(false); #ifndef ONLINE_JUDGE freopen("input.txt", "r", stdin); #endif // ONLINE_JUDGE while(cin >> n >> m) { memset(dp, 0, sizeof(dp)); dp[0][0][0] = 1; for(int i = 1; i <= n; i++) for(int j = 0; j <= m; j++) for(int k = 0; j+k <= m; k++) { dp[i][j][k] += dp[i-1][j][k] %= p;///不放 if(j-1>=0) dp[i][j][k] = (dp[i][j][k] + dp[i-1][j-1][k]*(m-(j-1)-k))%p; ///在沒炮的地方放一個 if(j+1 <= m && k >= 1) dp[i][j][k] = (dp[i][j][k] + dp[i-1][j+1][k-1]*(j+1))%p; ///在一個炮的地方放一個 if(j >= 2) dp[i][j][k] = (dp[i][j][k] + dp[i-1][j-2][k]*C(m-k-(j-2)))%p; ///在沒炮的地方放兩個 if(k >= 1 && j>= 1) dp[i][j][k] = (dp[i][j][k] + dp[i-1][j][k-1]*(m-j-(k-1))*(j))%p; ///在沒炮的地方和一個炮的地方分別放一個 if(j+2 <= m && k >= 2) dp[i][j][k] = (dp[i][j][k] + dp[i-1][j+2][k-2]*C(j+2))%p; ///在兩個一個炮的地方分別放一個炮 } LL ans = 0; for(int i = 0; i <= m; i++) for(int j = 0; j+i <= m; j++) ans = (ans+dp[n][i][j])%p; cout << ans << endl; } return 0; }