DP（大容量背包的一种做法）—— POJ 1786

最新推荐文章于 2020-05-13 13:46:30 发布

FeBr2

最新推荐文章于 2020-05-13 13:46:30 发布

阅读量3.1k

点赞数 1

分类专栏： ACM算法(题解)：动态规划 ——背包文章标签： dp

本文链接：https://blog.csdn.net/FeBr2/article/details/55095045

版权

ACM算法(题解)：同时被 3 个专栏收录

237 篇文章 3 订阅

订阅专栏

动态规划

51 篇文章 0 订阅

订阅专栏

——背包

7 篇文章 0 订阅

订阅专栏

题目链接： http://acm.hdu.edu.cn/showproblem.php?pid=1786
题意：给出一个背包容量最大为1E9的背包，给出N（N<10）种物品，每种物品给出质量（不超过100），数量无限，每拿一次物品需要3s时间，求填满背包所需要的最小时间
分析：这题看起来是一道普通的背包，但是背包容量太大，所以不可能直接将背包容量作为状态来做，我们需要另外找一个合适的量来作为状态。
设这N种物品所构成的总重为 $W_s$ ，那么有：
$W_s = W_1 * n_1 + W_2 * n_2 + ... + W_N * n_N$
( 将 $W_1..W_N$ 从小到大排序 )
因为题目要求拿物品的次数最小，所以我们要缩小 $n_i$ 的大小，如果 $n_i>=W_N (i<N)$ 我们就可以将它用 $W_N$ 替换掉一部分，以保证 $n_i<W_N (i<N)$ ，这样，拿物品的总次数就可以最小。然后，上式可以转换为：
$W_s = W_1 * n_1 + W_2 * n_2 + ... + W_N * n_N <= (N-1)*W_N*W_N + W_N*n_N$
我们可以发现， $W_s - W_N*n_N <=(N-1)*W_N*W_N (最大为80000)$ ，所以我们就可以用 N中物品组成的总重量减去质量最大的物品组成的重量，即前N-1种物品组成的总重量来表示一个状态。而只要前N-1种物品的状态确定了，第N种物品的数量也可以求出。
状态：
设DP[j]表示前 N-1 种物品构成质量 j 最小需要拿多少次
转移方程：

if(j-A[i] == 0) DP[j] = 1;
eles
{
    if(j-A[i]>0 && DP[j-A[i]!=INF)
    {
        DP[i] = MIN(DP[i], DP[j-A[i]]+1);
    }       
}

上述DP方程帮助我们求出了前N-1中物品可能构成的重量状态最少需要拿的次数，然后我们还需要结合第N中物品拿的次数从中选取一个最小值：

int k = M%A[0] //求出若全用A[0]来填背包，背包最终会余下k质量不能被填满
while(k<len)   //len为前N-1种物品在最小次数情况下可能构成的最大质量
{
    if(DP[k]!=INF)
    {
        ans = MIN(ans, DP[k]+(M-k)/A[0])
    }
    k+=A[0]
}

AC 代码：

/*************************************************************************
    > File Name: test.cpp
    > Author: Akira 
    > Mail: qaq.febr2.qaq@gmail.com 
 ************************************************************************/

#include <iostream>
#include <cstdio>
#include <cstring>
#include <string>
#include <cstdlib>
#include <algorithm>
#include <queue>
#include <stack>
#include <map>
#include <cmath>
#include <vector>
#include <set>
#include <list>
#include <ctime>
#include <climits>
typedef long long LL;
typedef unsigned long long ULL;
typedef long double LD;
#define MST(a,b) memset(a,b,sizeof(a))
#define CLR(a) MST(a,0)
#define Sqr(a) ((a)*(a))
using namespace std;

#define MaxN 100000
#define MaxM MaxN*10
#define INF 5E10
#define bug cout<<88888888<<endl;
#define MIN(x,y) (x<y?x:y)
#define MAX(x,y) (x>y?x:y)

template<typename _> inline void scan(_& t)
{
    int c;
    while((c = getchar()) < '0' || c > '9');
    t = c - '0';
    while((c = getchar()) >= '0' && c <= '9') t = t * 10 + c - '0';
}
template<typename _> inline void print(_ x)
{
    int len = 0, p[20];
    if(x < 0) putchar('-'), x = -x;
    while(x) p[++len] = x % 10, x /= 10;
    if(!len) p[++len] = 0;
    while(len) putchar(p[len--] + '0');
}
int N,M;
int A[10];
int DP[80008];
void solve()
{
    int len = 8*A[0]*A[0];
    for(int i=1;i<=len;i++) DP[i] = INT_MAX;
    DP[0] = 0;
    for(int i=1;i<=len;i++)
    {
        for(int j=1;j<N;j++)
        {
            if(i-A[j] == 0) DP[i] = 1;
            else if(i-A[j]>0 && DP[i-A[j]] != INT_MAX) DP[i] = MIN(DP[i], DP[i-A[j]] + 1);
        }
    }

    int ans = INT_MAX;
    int k = M%A[0];
    while(k<len)
    {
        if(DP[k]!=INT_MAX) 
            ans = MIN(ans, DP[k]+(M-k)/A[0]);
        k+=A[0];
    }
    ans = ans*3;
    printf("%d\n", ans);
}  

int main()
{
    while(cin >> N >> M && N && M)
    {
        for(int i=0;i<N;i++)
        {
            cin >> A[i] ;
        }
        sort(A,A+N,greater<int>());
        solve();
    }
    system("pause");
}