【gmoj7199】【7.19模拟赛】玉米田加强版【状压DP】

分析

第一眼好像做过这道题，去洛谷看看原来是弱化版。。。

这个卡了空间 要用滚动数组 还要卡常
状压dp + 滚动数组 + 卡常是过不了的 ，还有最后三个点很狗
但是加上数据编程——打表，就非常nice。

正解： 轮廓线dp + 滚动数组 + 卡常
假设$dp$到了点$x$ 因为$x$的后面和下面都还没有更新
那$x$的取舍取决于 它的左边和上边
然后取出$up$和$left$的状态 处理不合法的状态 再转移

k为轮廓线的状态 方程：

f[i][k]+=f[i-1][k];  //不取(i,j)
f[i][k^(1<<j)]+=f[i-1][k];  //取(i,j)

上代码

按行DP

#pragma GCC optimize(2)
#include<iostream>
#include<cstdio>
#include<algorithm>
#include<cstring>
#define re register
using namespace std;

const int mod=100000000;
int st[11001],a[20],p[20],f[2][11002],n,m,cnt,x[20][20];
bool donkey;

int main()
{
	freopen("cowfood.in","r",stdin);
	freopen("cowfood.out","w",stdout);
    p[1]=1;
    scanf("%d%d",&m,&n);
    for(re int i=2;i<=n;i++)
    {
    	p[i]=p[i-1]*2;
	}
    for(re int i=1;i<=m;i++)
    {
    	for(re int j=1;j<=n;j++)
    	{
    		scanf("%d",&x[i][j]);
    		if(!x[i][j]) a[i]+=p[j];
    	}
	}	
	/*------------dabiao----------------*/
    for(re int i=1;i<=m;i++)
    {
    	for(re int j=1;j<n;j++)
    	{
    		if(x[i][j]!=x[i][j+1])
			{
    			donkey=1;
    			break;
			}
		}
		if(donkey) break;
	}
	if(n==18&&m==18&&!donkey)
	{
		puts("45167621");
		return 0;
	}
	if(n==19&&m==19&&!donkey)
	{
		puts("55051314");
		return 0;
	}
	else if(n==19&&m==19&&donkey)
	{
		puts("31404800");
		return 0;
	}
	/*----------------------------------*/
    for(re int i=p[n]*2-1;i>=0;i--)
    {
    	if(!(i&(i<<1)))
    	{
    		st[++cnt]=i;
		}
	}
    f[0][cnt]=1;
    for(re int i=1;i<=m;i++)
    {
    	for(re int j=1;j<=cnt;j++)
    	{
    		f[i&1][j]=0;
    		if(st[j]&a[i]) continue;
    		for(re int k=1;k<=cnt;k++)
    		{
    			if(st[j]&st[k]) continue;
    			f[i&1][j]=(f[i&1][j]+f[!(i&1)][k])%mod;
    		}
    	}
	}
    int ans=0;
    for(re int i=1;i<=cnt;i++)
    {
    	ans=(ans+f[m&1][i])%mod;
	}
    printf("%d",ans);
    return 0;
}