第十一周作业——东东与ATM

东东与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 的程序，可根据具体金额请求机器交付现金。注意，这个程序计算程序得出的最大现金少于或等于可以根据设备的可用票据供应有效交付的现金。
　　
二、输入
　　程序输入来自标准输入。 输入中的每个数据集代表特定交易，其格式为：Cash N n1 D1 n2 D2 … nN DN其中0 <= Cash <= 100000是所请求的现金量，0 <= N <= 10是 纸币面额的数量，0 <= nk <= 1000是Dk面额的可用纸币的数量，1 <= Dk <= 1000，k = 1，N。 输入中的数字之间可以自由出现空格。 输入数据正确。
　　
三、输出
　　对于每组数据，程序将在下一行中将结果打印到单独一行上的标准输出中。

四、样例输入输出
Ｉｎｐｕｔ

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

Ｏｕｔｐｕｔ

735
630
0
0

五、解题思路
　　多重背包问题：二进制拆分+0-1背包问题+滚动数组
　　首先对所有面额纸币数量进行二进制拆分，将其拆分为多个面额的结果。二进制拆分的方法：举例说明，例如对7（111）进行拆分，将其拆分为（001）（010）（101），通过筛选这三组能够完全覆盖所有决策。对13（1101）进行拆分，将其拆分为7（001）（010）（100）和6（110），即可完成覆盖。
　　使用0-1背包方法对拆分的结果进行处理，0-1背包的转移方程为：fi,j=max(f[i-1][j] ,f[i-1][j-wi]+vi)。代表的状态分别是前i-1件物品放入容量为j的背包，和前i-1件物品放入容量为j-w[i]的背包。其时间复杂度和空间复杂度均为O(NV）。
　　为了进一步优化其空间复杂度，发现f[i][j]的更新只用到了f[i-1][j]的信息，所以采用滚动数组的方法将其优化成一维，第二层循环更改为逆序，因为顺序可能会使得同一个拆分结果被多次选择。
　　
六、样例代码

正确样例（二进制拆分做法）：

#include<stdio.h>
#include<iostream>
using namespace std;
int N,V;
int v[100010],vv[100010];
int w[20],ww[10010],C[20];
int cnt;
int f[100010];

void split()
{
	cnt=0;
	for(int i=1;i<=N;i++)
	{
		int t=C[i];
		for(int k=1;k<=t;k<<=1)
		{
			cnt++;
			ww[cnt]=k*w[i];
			t=t-k;
		}
		if(t>0)
		{
			cnt++;
			ww[cnt]=t*w[i];
		}
	}
}

int main()
{
	while(scanf("%d",&V)!=EOF)
	{
		scanf("%d",&N);
		for(int i=1;i<=N;i++)
     	{
	    	scanf("%d %d",&C[i],&w[i]);
    	}
    	
    	split();
    	for(int i=0;i<=V;i++) f[i]=0;

    	for(int i=1;i<=cnt;i++)
    	{
    		for(int j=V;j>0;j--)
    		{
    			if(j-ww[i]>=0)
    			f[j]=max(f[j],f[j-ww[i]]+ww[i]);
			}
		}
    	printf("%d\n",f[V]);
	}
	return 0;
}

错误样例（超时）：

//非二进制拆分 
#include<iostream>
using namespace std;

int f[15][100010];
int c[15];
int w[15];
int cash,n;

int main()
{
	while(scanf("%d",&cash)!=EOF)
	{
    	cin>>n;
    	for(int i=1;i<=n;i++)
    	{
	    	cin>>c[i]>>w[i];
	    }
    	for(int i=0;i<=cash;i++) f[0][i]==0;
    	for(int i=1;i<=n;i++) 
        {
        	for(int j=1;j<=cash;j++)
        	{
    	    	f[i][j]==f[i-1][j];
    	    	for(int k=0;k<=c[i];k++)
    	    	{
    			if(j-k*w[i]>=0)
    		    	f[i][j]=max(f[i][j],f[i-1][j-k*w[i]]+k*w[i]);
		    	}
	    	}
    	}
     	cout<<f[n][cash]<<endl;		
	}
}