POJ3254 状态压缩+dp


#include <iostream>
#include <stdio.h>
#include <algorithm>
#include <string.h>
///我不懂,你为什么老是把数组啊下标啊什么的写错?for(i<cur),你用cur去判断?无语,脑子呢
///因为地图的行数和列数都小于30(int最大是2^31-1,所以每一行或者每一列能压缩到一个int数即可),且每个点只有两个选择(要么为0要么为1),所以每行或者每列可以压缩成一个状态(这里选择行),一行之中可以选择的方法是有限的,行与行之间可以枚举所有方法,对每一行可以使用的方法做一下记忆,然后继续遍历
const int mod =100000000;
int n,m;
int dp[15][1<<15];//dp[i][j] 是 当第i行为第j种状态时,第一行到第i行到目前为止共有多少种放置方法;最后结果及,遍历第m(最后一行)行的每一种状态,连续相加得到总方法数。
int type[1<<15];//
int line[15];//用一个十进制数,存储当前地图中,各行的状态
int Judge1(int x){
    return x&(x<<1);
}
int Judge2(int x,int row){
    return type[x]&(~line[row]);///注意,用type[x],而非x
}
int main(){
    while(~scanf("%d%d",&m,&n)&&m+n){
        int number;
        for(int i=1;i<=m;i++){
            line[i]=0;
            for(int j=1;j<=n;j++){
                scanf("%d",&number);
                if(number){
                    line[i]+=(1<<(n-j));
                }
            }
        }
        int total=1<<n;
        int cur=0;
        for(int i=0;i<total;i++){
            if(!Judge1(i)){
                type[cur++]=i;///第一次剪枝,去掉一行之中相邻是1的方法
                //cout<<type[cur-1]<<endl;
            }
        }
        memset(dp,0,sizeof(dp));
        for(int i=0;i<cur;i++){
            if(!Judge2(i,1)){
                dp[1][i]=1;///去掉与地图不符的方法
            }
        }
        for(int i=2;i<=m;i++){
            for(int j=0;j<cur;j++){//当前行的放置状态
                if(Judge2(j,i)) continue;///(当前行的放置状态是否和地图冲突)
                for(int k=0;k<cur;k++){//上一行的放置状态
                    if(Judge2(k,i-1)) continue;///(判断枚举的上一行状态与地图是否冲突!!之前这个忘掉了)(另外,不用担心需不需要判断上上行与地图冲突,因为如果上上行冲突,那么上行dp就是0,那么这里加上也没事。)
                    if(type[j]&type[k]) continue;///(上一行的放置状态是否和当前行的放置状态是否冲突)
                    dp[i][j]+=dp[i-1][k];
                    //cout<<dp[i][j]<<endl;
                    //dp[i][j]=dp[i][j]%mod;
                }
            }
        }
        int ans=0;
        for(int i=0;i<cur;i++){
            //cout<<dp[m][cur]<<endl;
            ans+=dp[m][i];
            ans%=mod;
        }
        cout<<ans<<endl;
    }
}


  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 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、付费专栏及课程。

余额充值