HDU 4084 Campus Design 状压dp

原创 2013年12月02日 13:45:23

题意:在n*m的格子里面放1*1和1*2的砖块,要求不能放在格子为0的地方,并且使用1*1的数目在[c,d]之间,求方案总数。

想了一天也没有想到有效的转移,其实是陷入了一种错误的想法中,转移应该是一块一块的,不应该是一层一层的,终于在看过了别人的代码以后明白了。

又跟学长交流了一下,经过学长的指导,发现是自己想轴了,一层一层的转移也是可行的,比较一下两个方式的复杂度,一块一块的转移是O(n*m*d*2^m)1437 ms,一层一层的转移是O(n*d*2^m*dfs一层的复杂度)6406 ms,试了一下,也是可以AC的。

一块一块的转移

#include <iostream>  
#include <stdio.h>  
#include <string.h>  
  
using namespace std;  
const int mod=1000000007;  
char s[110][12];  
int n,m,c,d;  
int dp[2][1100][25];  
int now,next;  
void update(int nextp,int p,int nextk,int k)  
{  
    dp[next][nextp][nextk]+=dp[now][p][k];  
    if (dp[next][nextp][nextk]>=mod)   
        dp[next][nextp][nextk]-=mod;  
}  
int main()  
{  
   // freopen("in.txt","r",stdin);  
    while (~scanf("%d%d%d%d",&n,&m,&c,&d))  
    {  
        for (int i=1;i<=n;i++)  
            scanf("%s",s[i]);  
        memset(dp,0,sizeof(dp));  
        dp[0][(1<<m)-1][0]=1;  
  
        now=0;  
  
        for (int i=1;i<=n;i++)  
            for (int j=0;j<m;j++)  
        {  
            next=now^1;  
            //状态表示的是当前格子的前m块格子是否放砖的情况;  
            if (s[i][j]=='0')  
            {  
                for (int k=0;k<=d;k++)  
                    for (int p=(1<<(m-1));p<(1<<m);p++)//格子是'0',什么也不放,并且必须保证这个格子的上面那个格子放过了;  
                    if (dp[now][p][k])  
                    update(((p<<1)-(1<<m))^1,p,k,k);  
            }  
            else//格子是'1';  
            {  
                for (int k=0;k<=d;k++)  
                for (int p=(1<<(m-1));p<(1<<m);p++)//什么也不放,并且必须保证这个格子的上面那个格子放过了;  
                if (dp[now][p][k])  
                update(((p<<1)-(1<<m)),p,k,k);  
  
                for (int k=0;k<=d;k++)  
                for (int p=(1<<(m-1));p<(1<<m);p++)//放1*1,并且必须保证这个格子的上面那个格子放过了;  
                if (dp[now][p][k])  
                update(((p<<1)-(1<<m))^1,p,k+1,k);  
  
                if (j!=0)  
                {  
                    for (int k=0;k<=d;k++)  
                    for (int p=(1<<(m-1));p<(1<<m);p++)//横放1*2,并且必须保证这个格子的上面那个格子放过了,前面的格子是空的;  
                    if (!(p&1)&&dp[now][p][k]) update(((p<<1)-(1<<m))^3,p,k,k);  
                }  
  
                for (int k=0;k<=d;k++)  
                for (int p=0;p<(1<<(m-1));p++)//竖放1*2,并且必须保证这个格子的上面的格子是空的;  
                if (dp[now][p][k])  
                update((p<<1)^1,p,k,k);  
            }  
            memset(dp[now],0,sizeof(dp[now]));  
            now^=1;  
        }  
        int ans=0;  
  
        for (int k=c;k<=d;k++)  
        {  
            ans+=dp[now][(1<<m)-1][k];  
            if (ans>=mod) ans-=mod;  
        }  
  
    }  
  
    return 0;  
} 



一层一层的转移:


#include <iostream>  
#include <stdio.h>  
#include <string.h>  
  
using namespace std;  
const int mod=1000000007;  
char s[110][12];  
int n,m,c,d;  
int dp[110][1100][25];  
void dfs(int i,int j,int k,int num,int p,int kk)  
{  
    if (num>=m)//达到一种状态;  
    {  
        if(k+kk<=d)  
        {  
            dp[i+1][p][k+kk]+=dp[i][j][k];  
            if (dp[i+1][p][k+kk]>=mod) dp[i+1][p][k+kk]-=mod;  
        }  
        return;  
    }  
    if (j&(1<<num))//当前位置已放;  
    {  
        if (s[i+1][num]=='0')  
        dfs(i,j,k,num+1,p+(1<<num),kk);  
        else dfs(i,j,k,num+1,p,kk);  
        return;  
    }  
    if (s[i+1][num]=='0') //放1*1  
        dfs(i,j,k,num+1,p+(1<<num),kk+1);  
        else dfs(i,j,k,num+1,p,kk+1);  
  
    if (num<m-1&&!(j&(1<<(num+1))))//横放1*2;  
    {  
        int pp=p;  
        if (s[i+1][num]=='0') pp+=(1<<num);  
        if (s[i+1][num+1]=='0') pp+=(1<<(num+1));  
        dfs(i,j,k,num+2,pp,kk);  
    }  
    if (s[i+1][num]=='1')//竖放1*2;  
       dfs(i,j,k,num+1,p+(1<<num),kk);  
  
    return;  
}  
int main()  
{  
   // freopen("in.txt","r",stdin);  
    while (~scanf("%d%d%d%d",&n,&m,&c,&d))  
    {  
        for (int i=1;i<=n;i++)  
            scanf("%s",s[i]);  
  
            for (int i=0;i<m;i++)  
                s[n+1][i]='0';  
  
        memset(dp,0,sizeof(dp));  
        int sum=(1<<m)-1;  
  
        int p=0;  
        for (int i=0;i<m;i++)//把最后一行的下面一行填满;  
            if (s[1][i]=='0') p+=(1<<i);  
  
        dp[1][p][0]=1;  
  
        for (int i=1;i<=n;i++)  
        {  
            for (int j=0;j<=sum;j++)  
            for (int k=0;k<=d;k++)  
            if (dp[i][j][k])  
            {  
                dfs(i,j,k,0,0,0);  
            }  
        }  
  
        int ans=0;  
        for(int i=c;i<=d;i++)  
        {  
            ans+=dp[n+1][sum][i];  
            if (ans>=mod) ans-=mod;  
  
        }  
        printf("%d\n",ans);  
    }  
    return 0;  
  
  
}




相关文章推荐

HDU 4804 Campus Design(状压DP)

Campus Design Time Limit: 15000/8000 MS (Java/Others)    Memory Limit: 32768/32768 K (Java/Others) ...
  • Yasola
  • Yasola
  • 2017年05月10日 22:40
  • 164

HDU 4804 Campus Design 轮廓线DP

原题 Description Nanjing University of Scienceand Technology is celebrating its 60th anniversary. In...
  • SongBn
  • SongBn
  • 2016年08月05日 18:46
  • 185

hdu 4804 Campus Design 轮廓线dp

hdu 4804 Campus Design 轮廓线dp 题意: 给出一个n*m的01矩阵,其中要求把矩阵里面的1用1*1或1*2的砖块铺满,矩阵里面的0为障碍物,问使用1*1的砖块数>=c ...
  • whai362
  • whai362
  • 2015年07月09日 15:53
  • 668

hdu 4804 Campus Design (2013 南京 轮廓线dp)

hdu 4804 Campus Design (2013 南京 ) 水轮廓线dp,有一个类似的题目 Uva 11270 Tiling Dominoes 本题是在其基础上加了一些东西而已, ...
  • guognib
  • guognib
  • 2013年11月30日 20:23
  • 2381

HDU - 4804:Campus Design

HDU - 4804:Campus Design Nanjing University of Science and Technology is celebrating its 60th anniv...

HDU4804 Campus Design

题意:给一张地图,可以放

pku3229 The Best Travel Design 状压DP

题意太过于坑爹,数据也迷之坑,反正坑的不行。重边啊题目说是整数给小数啊什么鬼都出来了= = 嘛直接状压,伪TSP问题,floyed做完prework一切好办 Problem: 3229 User...
  • BPM136
  • BPM136
  • 2016年03月24日 16:23
  • 191

HDU DP 题集

  • 2013年04月24日 18:31
  • 724KB
  • 下载

Medianet Campus QoS Design 4.0

  • 2014年01月09日 00:14
  • 2.8MB
  • 下载

HDU 4539 郑厂长系列故事――排兵布阵 状压dp

题目链接:https://vjudge.net/contest/159644#problem/H题意:注意 是曼哈顿距离恰好是2,恰好!!!题解:与poj1185差不多的… 滚动数组,否则MLE,是...
  • yxg_123
  • yxg_123
  • 2017年04月21日 22:08
  • 151
内容举报
返回顶部
收藏助手
不良信息举报
您举报文章:HDU 4084 Campus Design 状压dp
举报原因:
原因补充:

(最多只允许输入30个字)