分析
第一眼好像做过这道题,去洛谷看看原来是弱化版。。。
这个卡了空间 要用滚动数组 还要卡常
状压dp + 滚动数组 + 卡常是过不了的 ,还有最后三个点很狗
但是加上数据编程——打表,就非常nice。
正解: 轮廓线dp + 滚动数组 + 卡常
假设
d
p
dp
dp到了点
x
x
x 因为
x
x
x的后面和下面都还没有更新
那
x
x
x的取舍取决于 它的左边和上边
然后取出
u
p
up
up和
l
e
f
t
left
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;
}