1. 程式人生 > >poj2411 Mondriaan's Dream(狀壓dp)

poj2411 Mondriaan's Dream(狀壓dp)

用1*2的小木塊覆蓋n*m的格子,問方案數。n,m<=10.狀壓dp,用二進位制數表示一行的狀態,逐行轉移。橫放用11表示,豎放用01表示。a[s]存的是上一行狀態為s,能轉移到的所有狀態。初始值f[0][tot]=1,表示第0行全為1的方案有1個。最後答案就是f[n][tot]。複雜度不會算啊。。。O(T*2^n*n*2^n)???

#include <cstdio>
#include <cstring>
#include <vector>
using namespace std;
#define ll long long
#define N 12
#define inf 0x3f3f3f3f
inline int read(){ int x=0,f=1;char ch=getchar(); while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=getchar();} while(ch>='0'&&ch<='9') x=x*10+ch-'0',ch=getchar(); return x*f; } int tot; int n,m,bin[N]; ll f[N][1<<11]; vector<int>a[1<<11]; void dfs(int s,int pos,int
ss){ if(pos>m){a[s].push_back(ss);return;} //橫放1*2 if(s&bin[pos-1]&&s&bin[pos]) dfs(s,pos+2,ss|bin[pos-1]|bin[pos]); //豎放2*1 if(s&bin[pos-1]) dfs(s,pos+1,ss); else dfs(s,pos+1,ss|bin[pos-1]); } int main(){ // freopen("a.in","r",stdin); bin[0
]=1; for(int i=1;i<=11;++i) bin[i]=bin[i-1]<<1; while(1){ n=read();m=read();if(!n) break;tot=bin[m]-1; for(int i=0;i<=tot;++i) a[i].clear(),dfs(i,1,0); memset(f,0,sizeof(f));f[0][bin[m]-1]=1;ll ans=0; for(int i=0;i<=n-1;++i) for(int s=0;s<=tot;++s){ if(!f[i][s]) continue; for(int j=0;j<a[s].size();++j) f[i+1][a[s][j]]+=f[i][s]; } printf("%lld\n",f[n][tot]); }return 0; }