POJ 2411 Mondriaan's Dream // 状压dp

POJ 2411 Mondriaan's Dream // 状压dp 

http://poj.org/problem?id=2411

突然觉得这题有丶意思，拿来写写看。（是个典型例题，用2*1砖头铺满矩形的方案数）

将一列压缩成一个二进制，就是一个状态，总共2^w*h个状态，直接dp+dfs。时间复杂度O（2^w*h）。

转移的时候，已知前一行满，将当前列填满，这样转移就行了。

h或w过大时，可以用滚动数组转移。

（好像还可以轮廓线解  有时间看看轮廓线dp）

#include<iostream>
#include<cstdio>
#include<algorithm>
#include<cstdlib>
#include<cstring>
using namespace std;
#define LL long long
int n,m;
LL dp[3000][15];/*前一行满，当前行的二进制,填补当前行,转移到后一行*/
void dfs(int now,int l,int h,int st,LL sum){
	if(h==n+1) {   //跳出
		dp[st][l+1]+=sum;
		return;
	}
	if(((now>>(h-1))&1)==1) { //当前位置已经被前面的占领了，那么直接跳过这个位置，到下一行
		dfs(now,l,h+1,st,sum);
		return ;
	}
	if(h==n) {//最后一个只能横着放
		dfs(now,l,h+1,st+(1<<(h-1)),sum);
		return;
	}
	dfs(now,l,h+1,st+(1<<(h-1)),sum); //横着放
	if(h<=n-1&&((now>>h)&1)==0) //如果能竖着放，竖着放然后转移
		dfs(now,l,h+2,st,sum);
}
int main(){
	while(cin>>n>>m){
		memset(dp,0,sizeof(dp));
		dp[0][1]=1;
		if(n==0) break;
		for(int i=1;i<=m;i++){ 
			for(int j=0;j<=(1<<n)-1;j++){ //对每个状态更新后面的状态
				dfs(j,i,1,0,dp[j][i]);
			}
		}
		cout<<dp[0][m+1]<<endl;
	}
	return 0;
}