poj 2411 Mondriaan's Dream 【dp】

题目:poj 2411 Mondriaan's Dream


题意:给出一个n*m的矩阵,让你用1*2的矩阵铺满,然后问你最多由多少种不同的方案。


分析:这是一个比较经典的题目,网上各种牛B写法一大堆。题解也是

我们可以定义状态:dp【i】【st】:在第 i 行状态为 st 的时候的最大方案数、

然后转移方程:dp【i】【st】 = sum (dp【i-1】【ss】)

即所有的当前行都是由上一行合法的状态转移而来。

而状态的合法性由两种铺法得到,第一种横放,注意要求前一行全满,然后竖放,上一行为空,可以留空。


AC代码:

#include <cstdio>
#include <cstring>
#include <string>
#include <iostream>
#include <algorithm>
#include <vector>
using namespace std;
const long long N = 15;
long long dp[N][1<<N];
long long ans[N][N];
long long n,m;
bool solve(long long st)
{
    long long tmp=0;
    for(long long i=0; i<m; i++)
    {
        if(st&(1<<i))
            tmp++;
        else
        {
            if(tmp%2)
                return false;
            tmp=0;
        }
    }
    if(tmp%2)
        return false;
    return true;
}
bool cmp(long long st,long long up)
{
    for(long long i=0; i<m; i++)
    {
        if(st&(1<<i))
        {
            if(up&(1<<i))
            {
                if(i==m-1 || !(st&(1<<(i+1))) || !(up&(1<<(i+1))))
                    return false;
                else
                    i++;
            }//否则的话竖放
        }
        else
        {
            if(!(up&(1<<i)))
                return false;
        }
    }
    return true;
}
int main()
{
    memset(ans,0,sizeof(ans));
    while(~scanf("%lld%lld",&n,&m))
    {
        if(n==0 && m==0)
            break;
        if((m*n)%2)
        {
            printf("0\n");
            continue;
        }
        if(m>n)
            swap(m,n);
        if(ans[n][m])
        {
            printf("%lld\n",ans[n][m]);
            continue;
        }
        memset(dp,0,sizeof(dp));
        long long len = 1<<m;
        for(long long i=0; i<len; i++)
        {
            if(solve(i))
                dp[1][i]=1;
        }
        for(long long i=2; i<=n; i++)
        {
            for(long long st=0; st<len; st++) //当前行
            {
                for(long long k = 0; k<len; k++) //上一行
                {
                    if(cmp(st,k))
                        dp[i][st]+=dp[i-1][k];
                }
            }
        }
        printf("%lld\n",dp[n][len-1]);
        ans[n][m] = dp[n][len-1];
    }
    return 0;
}


  • 1
    点赞
  • 1
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
题目描述 给出一个$n\times m$的矩阵,每个位置上有一个非负整数,代表这个位置的海拔高度。一开始时,有一个人站在其中一个位置上。这个人可以向上、下、左、右四个方向移动,但是只能移动到海拔高度比当前位置低或者相等的位置上。一次移动只能移动一个单位长度。定义一个位置为“山顶”,当且仅当从这个位置开始移动,可以一直走到海拔高度比它低的位置上。请问,这个矩阵中最多有多少个“山顶”? 输入格式 第一行两个整数,分别表示$n$和$m$。 接下来$n$行,每行$m$个整数,表示整个矩阵。 输出格式 输出一个整数,表示最多有多少个“山顶”。 样例输入 4 4 3 2 1 4 2 3 4 3 5 6 7 8 4 5 6 7 样例输出 5 算法1 (递归dp) $O(nm)$ 对于这道题,我们可以使用递归DP来解决,用$f(i,j)$表示以$(i,j)$为起点的路径最大长度,那么最后的答案就是所有$f(i,j)$中的最大值。 状态转移方程如下: $$ f(i,j)=\max f(x,y)+1(x,y)是(i,j)的下一个满足条件的位置 $$ 注意:这里的状态转移方程中的$x,y$是在枚举四个方向时得到的下一个位置,即: - 向上:$(i-1,j)$ - 向下:$(i+1,j)$ - 向左:$(i,j-1)$ - 向右:$(i,j+1)$ 实现过程中需要注意以下几点: - 每个点都需要搜一遍,因此需要用双重for循环来枚举每个起点; - 对于已经搜索过的点,需要用一个数组$vis$来记录,防止重复搜索; - 在进行状态转移时,需要判断移动后的点是否满足条件。 时间复杂度 状态数为$O(nm)$,每个状态转移的时间复杂度为$O(1)$,因此总时间复杂度为$O(nm)$。 参考文献 C++ 代码 算法2 (动态规划) $O(nm)$ 动态规划的思路与递归DP类似,只不过转移方程和实现方式有所不同。 状态转移方程如下: $$ f(i,j)=\max f(x,y)+1(x,y)是(i,j)的下一个满足条件的位置 $$ 注意:这里的状态转移方程中的$x,y$是在枚举四个方向时得到的下一个位置,即: - 向上:$(i-1,j)$ - 向下:$(i+1,j)$ - 向左:$(i,j-1)$ - 向右:$(i,j+1)$ 实现过程中需要注意以下几点: - 每个点都需要搜一遍,因此需要用双重for循环来枚举每个起点; - 对于已经搜索过的点,需要用一个数组$vis$来记录,防止重复搜索; - 在进行状态转移时,需要判断移动后的点是否满足条件。 时间复杂度 状态数为$O(nm)$,每个状态转移的时间复杂度为$O(1)$,因此总时间复杂度为$O(nm)$。 参考文献 C++ 代码

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值