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

看到第一眼就发觉是壮压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 1000000000using 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;}