状压DP : [USACO06NOV]玉米田

玉米田

内存限制:128 MiB
时间限制:1000 ms
标准输入输出
题目类型:传统
评测方式:文本比较
题目描述
农场主John新买了一块长方形的新牧场,这块牧场被划分成M行N列(1 ≤ M ≤ 12; 1 ≤ N ≤ 12),每一格都是一块正方形的土地。John打算在牧场上的某几格里种上美味的草,供他的奶牛们享用。

遗憾的是,有些土地相当贫瘠,不能用来种草。并且,奶牛们喜欢独占一块草地的感觉,于是John不会选择两块相邻的土地,也就是说,没有哪两块草地有公共边。

John想知道,如果不考虑草地的总块数,那么,一共有多少种种植方案可供他选择?(当然,把新牧场完全荒废也是一种方案)

输入格式
第一行:两个整数M和N,用空格隔开。

第2到第M+1行:每行包含N个用空格隔开的整数,描述了每块土地的状态。第i+1行描述了第i行的土地,所有整数均为0或1,是1的话,表示这块土地足够肥沃,0则表示这块土地不适合种草。

输出格式
一个整数,即牧场分配总方案数除以100,000,000的余数。

样例
样例输入
2 3
1 1 1
0 1 0
样例输出
9

(此题注意方案数要除以余数)

这道题很容易想到DP,并且是状压DP。为什么呢,设状态(f_{i})为到第i行为止总的方案数,如果我们不能表示出(i-1)行的状态的话,那么第(i)行也推不出来。(第(i)行和第(i-1)行不能有土地邻接)

于是考虑状态压缩。用一个长度为(m)的二进制数(j)来表示第(i)行选择土地的情况,(j)的第(k)位为1表示在第(i)行的第(k)列选了一个土地。

于是我们便得到了状态转移方程:
(f_{i,j} = sum_{f_{i-1,k}})

那么这篇题解到这里就结束了

好吧别打我。

这道题的难点就是在于判断(j),(k),(j)(k)这两个状态是否合法。这里就体现出了位运算的优越性。

对于判断一行农田(j)是否有相邻的已选的农田,我们可以这样判断:

if( (j<<1) & j)  不合法
else 合法

(位运算这篇博客里不会详讲,这种方法的可行性请读者自己翻阅资料思考)

同理对于状态(k)我们也可以这样判断。

对于上下相接的两个状态(j)(k),我们要判断它们的每一位是否出现了都是1的情况:

if(i&j) 不合法
else 合法

(还有一个居然忘讲了)

需要注意的是,在枚举(j)的过程中,我们还要判断我们是否选择了“贫瘠”的土地,这是我使用的方法:

for(int i=1;i<=n;i++)  for(int j=1,L;j<=m;j++) {
        cin>>L;
        P[i]=(P[i]<<1)+L;
	}//将初始的农田转换成2进制
	for(int i=0;i<(1<<m);i++)
        if((!((i<<1) & i)) && ((P[1] & i) == i))
            f[1][i]=1; //初始化状态,((P[1] & i)  ==  i)就是我用来判断是否
            			//选择了贫瘠土地的方法.

(关于状压DP,我这里要说一句,熟练掌握位运算是必要条件,这一点可以通过多做题来训练)

代码

#include <iostream>
#include <cstdio>
#include <vector>
using namespace std;

#define N 15
#define LL long long
#define MOD 100000000

LL n,m,f[15][100010],ans,P[N];

int main() {
	cin>>n>>m;
	for(int i=1;i<=n;i++)  for(int j=1,L;j<=m;j++) {
        cin>>L;
        P[i]=(P[i]<<1)+L;
	}
    for(int i=0;i<(1<<m);i++)
        if((!((i<<1) & i)) && ((P[1] & i) == i))
            f[1][i]=1;
    for(int i=2;i<=n;i++)
        for(int y=0;y<(1<<m);y++)
            if((!((y<<1) & y)) && ((P[i] & y) == y))
                for(int x=0;x<(1<<m);x++)
                    if((!((x<<1) & x)) && ((P[i-1] & x) == x)) {
                        if(x & y) continue;
                        f[i][y]=(f[i][y]+f[i-1][x])%MOD;
                    }
    for(int i=0;i<(1<<m);i++)
        if((!((i<<1) & i)) && ((P[n] & i) == i))
                ans=(ans+f[n][i])%MOD;
    cout<<ans;
}

原文地址:https://www.cnblogs.com/MisakaMKT/p/11252270.html