SDU_week11_E - 选做题11-1 东东与 ATM(多重背包+二进制拆分)

题目描述

一家银行计划安装一台用于提取现金的机器。
机器能够按要求的现金量发送适当的账单。
机器使用正好N种不同的面额钞票，例如D_k，k = 1,2,…,N，并且对于每种面额D_k，机器都有n_k张钞票。
例如，
N = 3，
n_1 = 10，D_1 = 100，
n_2 = 4，D_2 = 50，
n_3 = 5，D_3 = 10
表示机器有10张面额为100的钞票、4张面额为50的钞票、5张面额为10的钞票。
东东在写一个 ATM 的程序，可根据具体金额请求机器交付现金。
注意，这个程序计算程序得出的最大现金少于或等于可以根据设备的可用票据供应有效交付的现金。

Input
程序输入来自标准输入。 输入中的每个数据集代表特定交易，其格式为：Cash N n1 D1 n2 D2 … nN DN其中0 <= Cash <= 100000是所请求的现金量，0 <= N <= 10是 纸币面额的数量，0 <= nk <= 1000是Dk面额的可用纸币的数量，1 <= Dk <= 1000，k = 1，N。 输入中的数字之间可以自由出现空格。 输入数据正确。

Output
对于每组数据，程序将在下一行中将结果打印到单独一行上的标准输出中。

Sample Input
735 3 4 125 6 5 3 350
633 4 500 30 6 100 1 5 0 1
735 0
0 3 10 100 10 50 10 10

Sample Output
735
630
0
0

Hint
第一个数据集指定一笔交易，其中请求的现金金额为 735。 机器包含3种面额的纸币：4张钞票 125、6张钞票 5和3张钞票 350。 机器可以交付所需现金的确切金额。
在第二种情况下，机器的票据供应不能满足所要求的确切现金数量。 可以交付的最大现金为 630。 请注意，在机器中组合钞票以匹配交付的现金有多种可能性。
在第三种情况下，机器是空的，没有现金交付。 在第四种情况下，请求的现金金额为 0，因此机器不交付现金。

分析

观察数据范围，Cash* （N *D）很大，直接用01背包必然TLE，所以我们采用二级制拆分的方法解决问题。

优化方法是利用每个正整数都可以由二级制组成，我们可以将一个数拆分。比如7=(111)2。那么我们就可以将7拆分成100, 010, 001，也就是4，2，1。这三个数可以组成1-7之间的任意数。所以7个物品就拆分成了三个物品，第一个物品时4倍单价，第二个是2倍单价，第三个是单价。这样这三种物品是否选择就可以构成1-7个物品是否选择。

那么当我们要拆分的数是13的话怎么处理？由于13=(1101)2，我们不可能将其拆分成1000，0100，0001。因为这样有部分数据无法表示，比如2就没有办法表示。我们对这个的处理方法是将13化为7+6。其中7用上面的拆分方法可以表示1-7之间的任意数。第四个物品就是6倍的单价。如果第四个物品不选，前三个就可以构成选择1-7件物品的情况，如果第6种物品选择，就代表了选择7-13个的情况。这样13就拆成了四个物品。又如255，如果不使用二进制拆分，那么这种物品就要枚举1-255种情况分别表示第1个是否选择，第2个是否选择… … 。但是如果我们使用了二进制拆分，那么就可以拆分成8个物品。

代码

#define _ ios::sync_with_stdio(false);cin.tie(0);cout.tie(0);
#include <iostream>
#include <cstdio>
#include <cstdlib>
#include <cmath>
#include <cstring>
#include <string>
#include <climits>
#include <algorithm>
#include <queue>
#include <vector>
#include <map>
#include <stack>
#include<utility> //pair
using namespace std;

int cash,N, n[11],d[11];//请求的现金量,纸币面额的数量，数量，面额
int newbag[10009],dp[100009];

int main(){
    while(cin>>cash)
    {
        cin>>N;
        memset(dp,0,sizeof(dp));
        for (int i=1; i<=N; i++)
            cin>>n[i]>>d[i];

        int cnt=0;
        for (int i=1; i<=N; i++)
        {//二进制拆分
            int temp=n[i];
            for (int k=1; k<=temp; k<<=1)
            {
                cnt++;
                newbag[cnt]=k*d[i];
                temp-=k;
            }
            if(temp>0)
            {//13=7+6中的6
                cnt++;
                newbag[cnt]=temp*d[i];
            }
        }
        for (int i=1; i<=cnt; i++)//转化为仅考虑质量的01背包
            for (int j=cash; j>=newbag[i]; j--)//注意，j>=newbag[i]而非j>=0
                dp[j]=max(dp[j],dp[j-newbag[i]]+newbag[i]);
        cout<<dp[cash]<<endl;
    }
    return 0;
}