USACO Corn Fields

以此博客来纪念我做的第一道状压DP

题目描述
Farmer John has purchased a lush new rectangular pasture composed of M by N (1 ≤ M ≤ 12; 1 ≤ N ≤ 12) square parcels. He wants to grow some yummy corn for the cows on a number of squares. Regrettably, some of the squares are infertile and can’t be planted. Canny FJ knows that the cows dislike eating close to each other, so when choosing which squares to plant, he avoids choosing squares that are adjacent; no two chosen squares share an edge. He has not yet made the final choice as to which squares to plant.

Being a very open-minded man, Farmer John wants to consider all possible options for how to choose the squares for planting. He is so open-minded that he considers choosing no squares as a valid option! Please help Farmer John determine the number of ways he can choose the squares to plant.

农场主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的余数。

这道题的数据范围很小只有12，但是状态非常多，所以我们要开启一个新世界的大门——状态压缩。什么是状态压缩呢。。这个我们可以问问百度~针对于放不放，正反面之类的问题。我们可以把他们的状态用二进制表达出来。比如在这道题种，放我们选择用1来表达，不放选择用0来表达。这道题的状态只与当前行和上一行有关，满足无后效性。
我们可以写出状态转移方程 f【i】【j】+=f【i-1】【k】（i为行数，j，k是状态，当满足k和j合法时）；
我们知道，如果同一行内两个1相邻，一定不合法。上下两个1相邻，也不合法，且土地不肥沃也不合法。我们要对不合法情况进行预处理。土地不肥沃则一定种不了东西，土地肥沃则种不种都行，取决于土地。下面是代码：

#include<bits/stdc++.h>
using namespace std;

const int MAXN=17;
const int MOD=100000000;

int dp[MAXN][1<<MAXN],n,m;
bool ok[MAXN][1<<MAXN],a[MAXN][MAXN];//a数组维护土地肥沃，ok数组维护合法性。

bool pd(int x,int h){
//判断同一行中的合法性，如果出现了连续2个1则不合法，我们可以用右移操作不断删去最后位并记录。
    int tem=0;
    int p=x;
    while(p){
        if(!tem){
            tem=(p&1);
            p>>=1;
        }
        if(tem){
            if(p&1)return 0;
            else {
                tem=0;
                p>>=1;
            }
        }
    }
    //判断土地能不能种东西 如果种了东西则看土地肥沃度（也可以看土地肥沃度，如果不肥沃看种没种东西）
    int times=n;
    p=x;
    while(p){
        if(p&1){
            if(a[h][times]){
                p>>=1;
                times--;
            }
            else return 0;
        }
        else {
            p>>=1;
            times--;
        }
    }
    return 1;
}

bool judge(int x,int y){//判断上下是否合法，直接按位与。
    return (x&y);
}

int main(){
    memset(dp,0,sizeof(dp));
    memset(ok,0,sizeof(ok));
    scanf("%d%d",&m,&n);
    for(int i=1;i<=m;i++)
        for(int j=1;j<=n;j++){
            scanf("%d",&a[i][j]);
        }
    for(int i=1;i<=m;i++)
        for(int j=0;j<=(1<<n)-1;j++){//注意要从0开始。。0也是状态，之前就忘了。
            ok[i][j]=pd(j,i);   
        }
/*  for(int i=1;i<=m;i++)
        for(int j=0;j<=(1<<n)-1;j++){
            printf("%d ",ok[i][j]);
        }*///调试每行合法性
    for(int j=0;j<=(1<<n)-1;j++){//想一下就知道初始状态只在第一行。
        if(ok[1][j])dp[1][j]=1;
    }
    for(int i=2;i<=m;i++)
        for(int j=0;j<=(1<<n)-1;j++)
            for(int k=0;k<=(1<<n)-1;k++){
                if(ok[i][j]&&ok[i-1][k]&&!judge(j,k))dp[i][j]+=dp[i-1][k];
            }
    int ans=0;
    for(int j=0;j<=(1<<n)-1;j++){
        ans=(ans+dp[m][j])%MOD;
    }
    printf("%d",ans);
    return 0;
}
//完结