POJ_1742 Coins(DP)

题目链接

题解:

这道题和杭电上的2844是一样的,但是poj上面的数据可能要强一些。在杭电上可能模板多重背包就可以过,但是这里就要变通一下,今天再做到,配合着杭电2844,把这道题的做法好好整理了一下,一共三种方法dp。


方法一:多重背包记录状态

因为题目只要求得到有多少个可达,所以可以用bool表示dp数组,变化递推式,最后压线AC。

递推式变化为 dp[i] |= dp[i-cost];

代码实现:

#include <iostream>
#include <cstdio>
#include <cstdlib>
#include <cstring>

using namespace std;

const int MAX_N = 110;
const int MAX_M = 100010;

int N,M;
int res;
int A[MAX_N];
int C[MAX_N];
bool dp[MAX_M];
void ZeroOnePack(int cost,int weight);
void CompletePack(int cost,int weight);
void MultiplePack(int cost,int weight,int amount);
int main(){
    while( scanf("%d%d",&N,&M) != EOF ){
        if( N == 0 && M == 0 ){
            break;
        }
        for( int i = 0; i < N; i++ ){
            scanf("%d",&A[i]);
        }
        for( int i = 0; i < N; i++ ){
            scanf("%d",&C[i]);
        }
        res = 0;
        memset(dp,0,sizeof(dp));
        dp[0] = 1;
        for( int i = 0; i < N; i++ ){
            MultiplePack(A[i],A[i],C[i]);
        }
        for( int i = 1; i <= M; i++ ){
            //表示容量为i的背包装满
            if( dp[i] ){
                res++;
            }
        }
        printf("%d\n",res);
    }
    return 0;
}

void MultiplePack(int cost,int weight,int amount){
    if( cost*amount >= M ){
        CompletePack(cost,weight);
        return ;
    }
    int k = 1;
    while( k < amount ){
        ZeroOnePack(k*cost,k*weight);
        amount -= k;
        k *= 2;
    }
    ZeroOnePack(amount*cost,amount*weight);
    return ;
}

void CompletePack(int cost,int weight){
    for( int i = cost; i <= M; i++ ){
        dp[i] |= dp[i-cost];
        //dp[i] = max(dp[i],dp[i-cost]+cost);
    }
    return ;
}

void ZeroOnePack(int cost,int weight){
    for( int i = M; i >= cost; i-- ){
        dp[i] |= dp[i-cost];
        //dp[i] = max(dp[i],dp[i-cost]+cost);
    }
    return ;
}



方法二:记录可达状态,dp钱币剩余数目

钱币i从0~N-1遍历,对于每一种钱币,若在这之前dp[j]已经可达,赋值dp[j] = C[i],然后记录剩余数目dp,排除不可达的情况。

代码实现:

#include <iostream>
#include <cstdio>
#include <cstdlib>
#include <cstring>

using namespace std;

const int MAX_M = 100010;
const int MAX_N = 110;

int N,M;
int res;
int A[MAX_N];
int C[MAX_N];
int dp[MAX_M];
int main(){
    while( scanf("%d%d",&N,&M) != EOF ){
        if( N == 0 && M == 0 ){
            break;
        }
        for( int i = 0; i < N; i++ ){
            scanf("%d",&A[i]);
        }
        for( int i = 0; i < N; i++ ){
            scanf("%d",&C[i]);
        }
        memset(dp,-1,sizeof(dp));
        dp[0] = 0;
        res = 0;
        for( int i = 0; i < N; i++ ){
            for( int j = 0; j <= M; j++ ){
                if( dp[j] >= 0 ){
                    dp[j] = C[i];
                }
                else if( j < A[i] || dp[j-A[i]] < 0 ){
                    dp[j] = -1;
                }
                else{
                    dp[j] = dp[j-A[i]]-1;
                }
            }
        }
        for( int i = 1; i <= M; i++ ){
            if( dp[i] >= 0 ){
                res++;
            }
        }
        printf("%d\n",res);
    }
    return 0;
}



方法三:记录可达状态,对于每一种钱币dp使用数目

设置vis数组,每次更新dp数组,在之前可达的基础上dp,dp数组记录使用数目。

代码实现:

#include <iostream>
#include <cstdio>
#include <cstdlib>
#include <cstring>

using namespace std;

const int MAX_M = 100010;
const int MAX_N = 110;

int N,M;
int res;
int A[MAX_N];
int C[MAX_N];
int dp[MAX_M];
int vis[MAX_M];
int main(){
    while( scanf("%d%d",&N,&M) != EOF ){
        if( N == 0 && M == 0 ){
            break;
        }
        for( int i = 0; i < N; i++ ){
            scanf("%d",&A[i]);
        }
        for( int i = 0; i < N; i++ ){
            scanf("%d",&C[i]);
        }
        memset(vis,0,sizeof(vis));
        vis[0] = 1;res = 0;
        for( int i = 0; i < N; i++ ){
            memset(dp,0,sizeof(dp));
            for( int j = A[i]; j <= M; j++ ){
                if( vis[j-A[i]] && !vis[j] ){
                    if( dp[j-A[i]] < C[i] ){
                        res++;
                        dp[j] = dp[j-A[i]]+1;
                        vis[j] = 1;
                    }
                }
            }
        }
        printf("%d\n",res);
    }
    return 0;
}



  • 0
    点赞
  • 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、付费专栏及课程。

余额充值