zoukankan      html  css  js  c++  java
  • dp完全背包问题,hdu1114

    http://acm.hdu.edu.cn/showproblem.php?pid=1114

    这个题目关键在于两点,一个是完全背包问题,一个是背包是否装满

    完全背包  就是每个物品都可以无限放入,所以每个dp[j]放的都是目前状态比较后的结果,化为一维数组,为dp[v]=min(dp[v],dp[v-wi]+ci);

    循环记得是顺序的,这样能保证这一个放入的一定是之前所已知的,因为,是目前的状态

    题目
    有N种物品和一个容量为V的背包,每种物品都有无限件可用。第i种物品的费用是c[i],价值是w[i]。求解将哪些物品装入背包可使这些物品的费用总和不超过背包容量,且价值总和最大。

    基本思路
    这个问题非常类似于01背包问题,所不同的是每种物品有无限件。也就是从每种物品的角度考虑,与它相关的策略已并非取或不取两种,而是有取0件、取1件、取2件……等很多种。如果仍然按照解01背包时的思路,令f[i][v]表示前i种物品恰放入一个容量为v的背包的最大权值。仍然可以按照每种物品不同的策略写出状态转移方程,像这样:
        f[i][v]=max{f[i-1][v-k*c[i]]+k*w[i]|0<=k*c[i]<=v}
    这跟01背包问题一样有O(N*V)个状态需要求解,但求解每个状态的时间已经不是常数了,求解状态f[i][v]的时间是O(v/c[i]),总的复杂度是超过O(VN)的。
    将01背包问题的基本思路加以改进,得到了这样一个清晰的方法。这说明01背包问题的方程的确是很重要,可以推及其它类型的背包问题。但我们还是试图改进这个复杂度。

    一个简单有效的优化
    完全背包问题有一个很简单有效的优化,是这样的:若两件物品i、j满足c[i]<=c[j]且w[i]& gt;=w[j],则将物品j去掉,不用考虑。这个优化的正确性显然:任何情况下都可将价值小费用高得j换成物美价廉的i,得到至少不会更差的方案。对于随机生成的数据,这个方法往往会大大减少物品的件数,从而加快速度。然而这个并不能改善最坏情况的复杂度,因为有可能特别设计的数据可以一件物品也去不掉。
    这个优化可以简单的O(N^2)地实现,一般都可以承受。另外,针对背包问题而言,比较不错的一种方法是:首先将费用大于V的物品去掉,然后使用类似计数排序的做法,计算出费用相同的物品中价值最高的是哪个,可以O(V+N)地完成这个优化。这个不太重要的过程就不给出伪代码了,希望你能独立思考写出伪代码或程序。

    转化为01背包问题求解
    既然01背包问题是最基本的背包问题,那么我们可以考虑把完全背包问题转化为01背包问题来解。最简单的想法是,考虑到第i种物品最多选V/c[i]件,于是可以把第i种物品转化为V/c[i]件费用及价值均不变的物品,然后求解这个01背包问题。这样完全没有改进基本思路的时间复杂度,但这毕竟给了我们将完全背包问题转化为01背包问题的思路:将一种物品拆成多件物品。
    更高效的转化方法是:把第i种物品拆成费用为c[i]*2^k、价值为w[i]*2^k的若干件物品,其中k满足c[i]*2^k<=V。这是二进制的思想,因为不管最优策略选几件第i种物品,总可以表示成若干个2^k件物品的和。这样把每种物品拆成O(log(V/c[i]))件物品,是一个很大的改进。

    O(VN)的算法
    这个算法使用一维数组,先看伪代码:
    for i=1..N
        for v=0..V
            f[v]=max{f[v],f[v-cost]+weight}
    你会发现,这个伪代码与P01的伪代码只有 v的循环次序不同而已。为什么这样一改就可行呢?首先想想为什么P01中要按照v=V..0的逆序来循环。这是因为要保证第i次循环中的状态f[i] [v]是由状态f[i-1][v-c[i]]递推而来。换句话说,这正是为了保证每件物品只选一次,保证在考虑“选入第i件物品”这件策略时,依据的是一个绝无已经选入第i件物品的子结果f[i-1][v-c[i]]。而现在完全背包的特点恰是每种物品可选无限件,所以在考虑“加选一件第i种物品”这种策略时,却正需要一个可能已选入第i种物品的子结果f[i][v-c[i]],所以就可以并且必须采用v=0..V的顺序循环。这就是这个简单的程序为何成立的道理。
    这个算法也可以以另外的思路得出。例如,基本思路中的状态转移方程可以等价地变形成这种形式:
        f[i][v]=max{f[i-1][v],f[i][v-c[i]]+w[i]}
    将这个方程用一维数组实现,便得到了上面的伪代码。
    最后抽象出处理一件完全背包类物品的过程伪代码,以后会用到:
    procedure CompletePack(cost,weight)
        for v=cost..V
            f[v]=max{f[v],f[v-c[i]]+w[i]}

    总结
    完全背包问题也是一个相当基础的背包问题,它有两个状态转移方程,分别在“基本思路”以及“O(VN)的算法“的小节中给出。希望你能够对这两个状态转移方程都仔细地体会,不仅记住,也要弄明白它们是怎么得出来的,最好能够自己想一种得到这些方程的方法。事实上,对每一道动态规划题目都思考其方程的意义以及如何得来,是加深对动态规划的理解、提高动态规划功力的好方法。

    背包是否装满,除了dp[0]=0;其他的都按情况赋为无穷大,或者无穷小,

    如果是第一种问法,要求恰好装满背包,那么在初始化时除了f[0]为0其它f[1..V]均设为-∞,这样就可以保证最终得到的f[N]是一种恰好装满背包的最优解。

    如果并没有要求必须把背包装满,而是只希望价格尽量大,初始化时应该将f[0..V]全部设为0。

    为什么呢?可以这样理解:初始化的f数组事实上就是在没有任何物品可以放入背包时的合法状态。如果要求背包恰好装满,那么此时只有容量为0的背包可能被价值为0的nothing“恰好装满”,其它容量的背包均没有合法的解,属于未定义的状态,它们的值就都应该是-∞了。如果背包并非必须被装满,那么任何容量的背包都有一个合法解“什么都不装”,这个解的价值为0,所以初始时状态的值也就全部为0了。

     1 #include<stdio.h>
     2 int min(int a,int b)
     3 {
     4     return a>b?b:a;
     5 }
     6 int main()
     7 {
     8     int t,e,f,n,v[555],w[555],i,j,dp[55500];
     9     scanf("%d",&t);
    10     while(t--)
    11     {
    12        scanf("%d%d%d",&e,&f,&n);
    13        f-=e;
    14        for(i=1;i<=n;i++)
    15        scanf("%d%d",&v[i],&w[i]);
    16        dp[0]=0;//把0件物品装满0体积的最小价值为0 
    17        for(i=1;i<=f;i++)dp[i]=9999999;
    18        for(i=1;i<=n;i++)
    19          for(j=w[i];j<=f;j++)
    20            dp[j]=min(dp[j],dp[j-w[i]]+v[i]);
    21        if(dp[f]==9999999)//不是原来的无穷大,说明装满了,否则就没装满 
    22        printf("This is impossible.
    ");
    23        else 
    24        printf("The minimum amount of money in the piggy-bank is %d.
    ",dp[f]);
    25     }
    26     return 0;
    27 }
    View Code
  • 相关阅读:
    系统的buffer 最小和最大值
    SpringCloud Alibaba微服务实战十二
    socket.recv(fd, buf, size, flag) &&setsockopt(fd, SOL_SOCKET, SO_RCVBUF, value, valuesize)
    源码阅读:Masonry(三)—— MASViewAttribute
    源码阅读:Masonry(三)—— MASViewAttribute
    源码阅读:Masonry(三)—— MASViewAttribute
    长江存储推全新3D NAND架构 挑战三星存储
    长江存储推全新3D NAND架构 挑战三星存储
    长江存储推全新3D NAND架构 挑战三星存储
    长江存储推全新3D NAND架构 挑战三星存储
  • 原文地址:https://www.cnblogs.com/huzhenbo113/p/3246820.html
Copyright © 2011-2022 走看看