当前位置: 代码迷 >> 综合 >> [Usaco2006 Nov]Corn Fields牧场的安排 壮压DP
  详细解决方案

[Usaco2006 Nov]Corn Fields牧场的安排 壮压DP

热度:2   发布时间:2024-01-13 17:23:03.0

看到第一眼就发觉是壮压DP

然后就三进制枚举子集吧。

这题真是壮压入门好题。。。


对于dp[i][j] 表示第i行,j状态下前i行的分配方案数。

那么dp[i][j]肯定是从i-1行转过来的

那么由于不能挨着放,那么我们肯定是枚举i - 1行状态时不能包含j的任何一位。

那么只要令k = ((1 << n) - 1) ^ j,k中肯定就不包含j的位了

是这样枚举k的子集

int sub = k;

do

{

    sub = k& (sub - 1);

}while(sub != k);

然后对每个子集,判断合法性,然后相加即可。

#include <iostream>
#include <cstdio>
#include <cstring>
#include <vector>
#include <algorithm>
#define MAXN 1005
#define INF 1000000000
using namespace std;
int dp[13][1 << 13];
int n, m;
int st[13];
int mod = 100000000;
bool ok(int s, int pos)
{if((s | st[pos]) > st[pos]) return false;for(int i = 0; i < n; i++)if(s & (1 << i)){if(s & (1 << (i + 1))) return false;}return true;
}
int main()
{int x;scanf("%d%d", &m, &n);for(int i = 1; i <= m; i++){for(int j = 0; j < n; j++){scanf("%d", &x);if(x) st[i] |= (1 << j);}}dp[0][0] = 1;for(int i = 1; i <= m; i++){for(int k = 0; k < (1 << n); k++){int s = ((1 << n) - 1) ^ k;if(ok(k, i)){//printf("%d %d\n", k, s);dp[i][k] += dp[i - 1][0];for(int j = s; j; j = s & (j - 1)){if(ok(j, i - 1))dp[i][k] = (dp[i][k] + dp[i - 1][j]) % mod;}}}}int ans = 0;for(int i = 0; i < (1 << n); i++)ans = (ans + dp[m][i]) % mod;printf("%d\n", ans);return 0;
}