Poj 2411 Mondriaan's Dream 状态压缩(方案数)

来源:互联网 发布:两个excel筛选重复数据 编辑:程序博客网 时间:2024/06/08 15:52

题目链接:http://poj.org/problem?id=2411

题目大意:给你一块N*M的地M,让你用1*2的地砖去铺满,求总的方案数。

解题思路:可以用深搜得到每行的合法状态,但本人弱菜,以后学习补上,直接暴力合法状态,具体解释见代码。

代码如下:

#include<cstdio>#include<cstring>#include<algorithm>using namespace std;int n,m;long long  dp[12][1<<12],sta[1<<12],vis[12];void init(){    memset(dp,0,sizeof(dp));    memset(sta,0,sizeof(sta));    memset(vis,0,sizeof(vis));}void solve(){    int num=(1<<m),i,j,k,p,ct=0;    for(i=0;i<num;i++)   //筛选一行的合法状态数,其中1代表横放,则两个1相邻才合法    {        memset(vis,0,sizeof(vis));        for(j=1; j<m; j++)           //j-1,故j从 1开始            if( (i&(1<<j)) &&  (i&(1<<(j-1))) && !vis[j-1])                vis[j-1]=vis[j]=1;        for(j=0;j<m;j++)            if( (i&(1<<j)) && !vis[j])            break;        if(j==m) sta[ct++]=i;    }    for(i=0;i<ct;i++)            //初始化,求方案数初始化为1        dp[0][ sta[i] ]=1;    for(i=1;i<n;i++)         //行数    {        for(j=0;j<num;j++)    //枚举i-1行状态        {            for(k=0;k<ct;k++)   //枚举第i行合法状态            {                int t=0;                for(p=m-1;p>=0;p--)  //i行与i-1行状态对比,得到第i行状态                {                    //上下同为1,表示上下均横放,同为0 ,表示上下均竖放,均合法,且该位置已填满为1;                    if( (j&(1<<p))==(sta[k]&(1<<p)) ) t=(t<<1)+1;                    //上面为1,下面为0,表示上面横放,下面竖放,说明该位置未填为0;                    else if( ((j&(1<<p))) &&  !(sta[k]&(1<<p)) )  t=(t<<1);                    //上为0,下为1,说明下面横放,上面未放满,不合法,跳出                    else break;                }                if(p==-1)                      dp[i][t]+=dp[i-1][j];            }        }    }    printf("%I64d\n",dp[n-1][num-1]);//最终状态为第n-1行铺满,即dp[n-1][num-1]}int main(){    while(scanf("%d%d",&n,&m)&&(n+m))    {        if(n<m) swap(n,m);        init();        solve();    }    return 0;}


0 0
原创粉丝点击