轮廓线DP基础题。
不得不说轮廓线DP的思想真的很精妙,但是轮廓线状态的表示方法我也是看了好久才理解。
轮廓线的状态不是按照一般状压DP不太一样,个人觉得这篇博客讲得比较简短清楚,这篇博客虽然具体但是不太清晰,再结合一下代码应该就能理解了。
话说我之前学的时候明明看到一篇很好的博客啊,现在怎么死活找不着了。。
#include<cstdio>
#include<iostream>
#include<map>
#include<algorithm>
#include<cstring>
#define LL long long
using namespace std;
const int N=12;
int n,m;
LL f[2][1<<N];
void read(int &x){
char ch=getchar();x=0;int w=1;
for(;ch<'0'||ch>'9';ch=getchar()) if (ch=='-') w=-1;
for(;ch>='0'&&ch<='9';ch=getchar()) x=(x<<3)+(x<<1)+ch-'0';
x*=w;
}
int main(){
for(;;){
read(n);read(m);
if (n==0&&m==0) break;
int cur=1;
//memset(f,0,sizeof f);
f[0][(1<<m)-1]=1;
for(int i=1;i<=n;i++)
for(int j=1;j<=m;j++){
cur^=1;
memset(f[cur^1],0,sizeof f[cur^1]);
for(int k=0;k<(1<<m);k++){
if (k&(1<<m-1)) f[cur^1][(k<<1)^(1<<m)]+=f[cur][k];
if (i>1 && !(k&(1<<m-1))) f[cur^1][k<<1|1]+=f[cur][k];
if (j>1 && (k&(1<<m-1)) && !(k&1)) f[cur^1][(k<<1)^(1<<m)|3]+=f[cur][k];
}
}
printf("%lld\n",f[cur^1][(1<<m)-1]);
}
return 0;
}