题解 luogu P1879 【[USACO06NOV]玉米田Corn Fields】

轮廓线dp

记录当前的轮廓线，根据当前轮廓线进行状态转移，理论时间复杂度比普通状压dp快。

---

		k3
k1	k2	*

* 如图，假设当前转移到位置，轮廓线状态为二进制数k1k2k3k4，1表示第i列已选，0表示第i列未选，转移时用位运算解压出k3与k2,作出决策，并维护每一列的轮廓线即可

---

代码如下

#include
#include
#include
#define max(x,y) ((x)>(y) ? (x) : (y))
#define min(x,y) ((x)<(y) ? (x) : (y))
#define Mod 100000000
#define N 18
using namespace std;
//运用滚动数组优化
int d[2][65546],n,m,a[N][N],cur;

int main()
{
    scanf("%d%d",&n,&m);
    for(int i=0;ifor(int j=0;jscanf("%d",&a[i][j]);
    }
    memset(d,0,sizeof d);
    cur=0;
    d[0][0]=1;
    for(int i=0;ifor(int j=0;j1;
            memset(d[cur],0,sizeof d[cur]);
            for(int k=0;k<(1<int up=(1<int left=0;//解压出左边和上面
                if(j>0) left=(1<<(j-1))&k;
                if(!i && up) continue;
                if(!j && left) continue;
                if(up && left) {//上面选了，左边选了，这个只能不选
                    d[cur][k^(1<1<1][k])%Mod;
                    continue;
                }
                if(up) {//上面选了，这个只能不选
                    d[cur][k^(1<1<1][k])%Mod;
                    continue;
                }
                if(left) {//左边选了，这个只能不选
                    d[cur][k]=(d[cur][k]+d[cur^1][k])%Mod;
                    continue;
                }
                //障碍不能选
                if(!a[i][j]) {
                    d[cur][k]=(d[cur][k]+d[cur^1][k])%Mod;
                    continue;
                }
                //除去以上状态，考虑选和不选两种情况
                d[cur][k]=(d[cur][k]+d[cur^1][k])%Mod;
                d[cur][k^(1<1<1][k])%Mod;
            }
        }
    }
    int ans=0;
    for(int i=0;i<(1<printf("%d\n",ans);
    return 0;
}