06普及组开心的金明程序思路

Posted 2023-04-11

请说明思路，再写出程序，非常感谢！！！
c++

用动态规划来解背包问题
在历届NOIP竞赛中，有4道初赛题和5道复赛题均涉及到背包问题，所谓的背包问题，可以描述如下：一个小偷打劫一个保险箱，发现柜子里有N类不同大小与价值的物品，但小偷只有一个容积为M的背包来装东西，背包问题就是要找出一个小偷选择所偷物品的组合，以使偷走的物品总价值最大。
如有4件物品，容积分别为： 3 4 5 8
对应的价值分别为： 4 5 7 10
小偷背包的载重量为：12
则取编号为1 2 3的物品，得到最大价值为16。
算法分析：如果采用贪心法，则先取价值最大的10，消耗了容积8，下面只能取容积为4的物品，得到价值5，这样总价值是15，这不是最优解，因此贪心法是不正确的。
采用穷举法，用一个B数组来表示取数的标记，当B[i]=0时表示第i件物品不取，当B[i]=1时表示第i件物品已取，初始化全部取0，以下算法是从后面的物品开始取起，通过B数组的取值把15种取法全部穷举出来，价值MAX初始化为0。
B[0] B[1] B[2] B[3] B[4]
0 0 0 0 0 初始化
0 0 0 0 1 取第4件物品，容积为8，不超，价值为10，将MAX替换为10
0 0 0 1 0 ｛取物品3，容积为5，不超，价值为7，不换｝
0 0 0 1 1 ｛取物品3、4，容积为13，超｝
0 0 1 0 0 ｛取物品2，容积为4，不超，价值为5，不换｝
0 0 1 0 1
0 0 1 1 0
0 0 1 1 1
……
0 1 1 1 0 ｛这是最佳方案｝
0 1 1 1 1
1 0 0 0 0 ｛当B〔0〕＝1时停止，B〔0〕称为哨兵｝
生成B数组中数据的方法如下：
fillchar(b,sizeof(b),0);
while b[0]=0 do
begin j:=n;
while b[j]=1 do dec(j);
b[j]:=1;
for i:=j+1 to n do
b[i]:=0;
end；
小结：以上每件物品只能取1件，所以取法只有0和1两种情况，我们称之为0、1背包，算法的时间复杂度为O（2N），在1秒内N只能做到20。
例1：选数（NOIP2002 初中组复赛第2题）
[问题描述]:已知 n 个整数 x1,x2,…,xn，以及一个整数 k（k＜n）。从 n 个整数中任选 k 个整数相加，可分别得到一系列的和。例如当 n=4，k＝3，4 个整数分别为 3，7，12，19 时，可得全部的组合与它们的和为：
3＋7＋12=22 3＋7＋19＝29 7＋12＋19＝38 3＋12＋19＝34。
现在，要求你计算出和为素数共有多少种。
例如上例，只有一种的和为素数：3＋7＋19＝29。
[输入]：
键盘输入，格式为：
n , k （1<=n<=20，k＜n）
x1,x2，…,xn （1<=xi<=5000000）
n[输出]：
屏幕输出，格式为：
一个整数（满足条件的种数）。
[输入输出样例]:
输入：
4 3
3 7 12 19
输出：
1
[算法分析]:本题应用背包问题中取数的方法进行穷举，在取数的过程中，当B数组中有K个1的时候将对应的K个数相加，再判断是不是素数。
主要程序段如下：
readln(n,k); sum:=0;
for i:=1 to n do read(a[i]);
fillchar(b,sizeof(b),0);
while b[0]=0 do
begin j:=n;
while b[j]=1 do dec(j);
b[j]:=1;
for i:=j+1 to n do b[i]:=0;
m:=0;
for i:=1 to n do
if b[i]=1 then m:=m+1; 统计1的个数
if m=k then
begin 计算此种取数方法得到的和S；
if S 是素数 then sum:=sum+1;
end;
end;
例2：采药(NOIP2005 初中组复赛第3题)
【问题描述】
辰辰是个天资聪颖的孩子，他的梦想是成为世界上最伟大的医师。为此，他想拜附近最有威望的医师为师。医师为了判断他的资质，给他出了一个难题。医师把他带到一个到处都是草药的山洞里对他说：“孩子，这个山洞里有一些不同的草药，采每一株都需要一些时间，每一株也有它自身的价值。我会给你一段时间，在这段时间里，你可以采到一些草药。如果你是一个聪明的孩子，你应该可以让采到的草药的总价值最大。”
如果你是辰辰，你能完成这个任务吗？

【输入文件】
输入文件medic.in的第一行有两个整数T（1 <= T <= 1000）和M（1 <= M <= 100），用一个空格隔开，T代表总共能够用来采药的时间，M代表山洞里的草药的数目。接下来的M行每行包括两个在1到100之间（包括1和100）的整数，分别表示采摘某株草药的时间和这株草药的价值。
【输出文件】
输出文件medic.out包括一行，这一行只包含一个整数，表示在规定的时间内，可以采到的草药的最大总价值。
【样例输入】
70 3
71 100
69 1
1 2
【样例输出】
3
【数据规模】
对于30%的数据，M <= 10；
对于全部的数据，M <= 100。
【算法分析】本题如果采用上述方法来解，只能将M算到20，而这里M〈=100，所以只能拿30%的分数，比较好的算法是采用动态规划，为了能说清算法，现重新举一个例子，若输入：
10 3
3 4
4 5
5 6
表示背包的容量是10，有3种物品。用一个数组用来表示背包容量与其最大价值的关系，上例中设置一个数组count,用下标表示容量，初始化为0。然后按物品的顺序一一来统计此时的最大价值，每种药品对应各种背包容量时得到的最大价值为：
对于是第i件物品，背包容量为j时的最大价值Cmax(j)=MAX(Cmaxj,Pi+余下空间的最大价值Cmax(j-i物品所占的空间))，如上例中，根据物品的不断增加，各容量背包得到的最大价值不断替换：

容量 1 2 3 4 5 6 7 8 9 10
价值
序号 0 0 0 0 0 0 0 0 0 0
1 0 0 4 4 4 4 4 4 4 4
2 0 0 4 5 5 5 9 9 9 9
3 0 0 4 5 6 6 9 10 11 11

[数据结构] time,price数组分别用来存入时间和价值，count来存入背包的价值。
var
time,price:array[1..100] of longint;
t:longint; i,m,j:integer;
count:array[0..1000] of longint;
begin
assign(input,'medic.in');
assign(output,'medic.out');
reset(input);
rewrite(output);
readln(t,m);
for i:=1 to m do
readln(time[i],price[i]);
fillchar(count,sizeof(count),0);
for i:=1 to m do
for j:=t downto 1 do
begin
if (j>=time[i]) and (price[i]+count[j-time[i]]>count[j]) then
count[j]:=price[i]+count[j-time[i]];
end; j>=time[i]表示当前的容量能放入背包，price[i]+count[j-time[i]]>count[j]表示第i件物品的价值加上第i件物品对于背包容量为j时余下空间的最大价值大于当前背包容量为j时的最大价值
writeln(count[t]);
close(input);
close(output);
end.
例3：开心的金明（NOIP2006 初中组复赛第2题）
题目较长，省略，本题与例2相比，在比较时要先将价值乘以一个数，其余一样，但要注意的是：本题N的范围是<=26，如果用0、1背包穷举法在1秒内只能过10个点中的8个点。
总结：采用动态规划的时间复杂度为O(n*m)，范围比穷举法大多了，但也有弱点，当数据不是整数时，就不可使用；如果还要求出具体的取法，那也相当麻烦。 参考技术A DP,0-1背包问题.
把总钱数看作背包可容重量.
把希望购买物品的个数看作要放进背包的物品的数量.
把物品的价格看作物品的重量.
把物品的价格与物品的重要度的乘积看作物品的价值.
i:前i个物品 j:背包剩余可容重量 w:第i个物品的重量 r:第i个物品的价值
i=0 f[i,j]=0
i>0 f[i,j]=maxf[i-1,j],f[i-1,j-w]+r

program happy(input,output);
var
i,j,m,n:longint;
w,r:array[1..25] of longint;
a:array[0..1,-10000..100000] of longint;
begin
assign(input,'happy.in');reset(input);
assign(output,'happy.out');rewrite(output);
read(n,m);
for i:=1 to m do
begin
read(w,r);
r:=r*w;
end;
for i:=1 to m do
begin
a[1]:=a[0];
for j:=w to n do
if a[0,j-w]+r>a[1,j] then a[1,j]:=a[0,j-w]+r;
a[0]:=a[1];
end;
writeln(a[1,n]);
close(input);close(output);
end.

(参见经典0-1背包-采药) 参考技术B 把0-1背包问题一改就行了 只是把加的 重量 改成 价值与重要度的乘积 ！
program happy;
var
f:array[0..25,0..30000] of longint;
i,j:integer;
n,m:longint;
v,p:array[1..25] of longint;
begin
readln(m,n);
for i:=1 to n do
read(v[i],p[i]);
fillchar(f,sizeof(f),0);
for i:=1 to n do
for j:=1 to m do
begin
f[i,j]:=f[i-1,j];
if (j>=v[i])and(f[i,j]<f[i-1,j-v[i]]+p[i]*v[i]) then
f[i,j]:=f[i-1,j-v[i]]+p[i]*v[i];
end;
writeln(f[n,m]);
end. 参考技术C 典型的01背包问题，时间复杂度为O(Nm)。
一种做法是使用滚动数组，每次可以从前往后更新，也可以从后往前更新；另一种是只使用一个一维数组，每次只能从后往前更新。
在表示状态的时候，可以用f[i]表示正好花了i元，能够得到的最大总和；也可以用f[i]表示花了不超过i元，能够得到的最大总和。
程序1：滚动数组+表示正好
var
n,m,i,j,ans:longint;
a:array[0..25,0..2] of longint;
dp,dp2:array[0..30000] of longint;
begin
assign(input,'happy.in'); reset(input);
assign(output,'happy.out'); rewrite(output);
fillchar(a,sizeof(a),0);
read(n,m);
for i:=1 to m do begin
read(a[i,0],a[i,1]); 输入价格和重要程度
a[i,2]:=a[i,0]*a[i,1]; 计算乘积
end;
fillchar(dp2,sizeof(dp2),$FF); 初始时dp[i]=-1，i>0，表示正好花i元得不到任何总和
dp2[0]:=0; 只有dp[0]=0，因为正好花了0元，得到的总和为0
for i:=1 to m do begin
dp:=dp2; 以后用dp来推dp2
for j:=0 to n do 枚举花了j元钱
if (dp[j]>-1) and (j+a[i,0]<=n) and (dp2[j+a[i,0]]<dp[j]+a[i,2]) then j元钱的方案存在，加上这次的价格没有超过可用的n元，并且这个方案得到的总和更大，就更新
dp2[j+a[i,0]]:=dp[j]+a[i,2];
end;
ans:=-1; 求dp2中的最大值
for i:=0 to n do
if dp2[i]>ans then ans:=dp2[i];
writeln(ans);
close(input); close(output);
end.

程序2：单个一维数组+表示不超过
var
n,m,i,j:longint;
a:array[0..25,0..2] of longint;
dp:array[0..30000] of longint;
begin
assign(input,'happy.in'); reset(input);
assign(output,'happy.out'); rewrite(output);
fillchar(a,sizeof(a),0);
read(n,m);
for i:=1 to m do begin
read(a[i,0],a[i,1]); 输入价格和重要程度
a[i,2]:=a[i,0]*a[i,1]; 计算乘积
end;
fillchar(dp,sizeof(dp),0); 初始时dp[i]=0，表示花不超过i元可以得到0总和
for i:=1 to m do
for j:=n downto 0 do 从后往前枚举花了j元钱
if (j+a[i,0]<=n) and (dp[j+a[i,0]]<dp[j]+a[i,2]) then 加上这次的价格没有超过可用的n元，并且这个方案得到的总和更大，就更新
dp[j+a[i,0]]:=dp[j]+a[i,2];
writeln(dp[n]); 输出最后花不超过n元能得到的最大总和
close(input); close(output);
end.

2006开心的金明

题目描述 Description

金明今天很开心，家里购置的新房就要领钥匙了，新房里有一间他自己专用的很宽敞的房间。更让他高兴的是，妈妈昨天对他说：“你的房间需要购买哪些物品，怎么布置，你说了算，只要不超过N元钱就行”。今天一早金明就开始做预算，但是他想买的东西太多了，肯定会超过妈妈限定的N元。于是，他把每件物品规定了一个重要度，分为5等：用整数1~5表示，第5等最重要。他还从因特网上查到了每件物品的价格（都是整数元）。他希望在不超过N元（可以等于N元）的前提下，使每件物品的价格与重要度的乘积的总和最大。

设第j件物品的价格为v[j]，重要度为w[j]，共选中了k件物品，编号依次为j1，j2，……，jk，则所求的总和为：

v[j1]*w[j1]+v[j2]*w[j2]+…+v[jk]*w[jk]。（其中*为乘号）

请你帮助金明设计一个满足要求的购物单。

输入描述 Input Description

输入的第1行，为两个正整数，用一个空格隔开：N m

（其中N（<30000）表示总钱数，m（<25）为希望购买物品的个数。）

从第2行到第m+1行，第j行给出了编号为j-1的物品的基本数据，每行有2个非负整数

v p

（其中v表示该物品的价格(v<=10000)，p表示该物品的重要度(1~5)）

输出描述 Output Description

输出只有一个正整数，为不超过总钱数的物品的价格与重要度乘积的总和的最大值（<100000000）。

样例输入 Sample Input

1000 5

800 2

400 5

300 5

400 3

200 2

样例输出 Sample Output

3900

数据范围及提示 Data Size & Hint

 

 

 

题解：

动归背包。

这道题目可以用搜索做（数据水，东西只有25个），不过容易超时，所以建议用的背包做。

用f【j】表示j元钱能卖多少，f【j】=max（f【j】，f【j-a【i】】+b【i】*a【i】），答案f【m】。

var n,m,i,j:longint;

    a,b,f:array[0..30001]of longint;

function max(a,b:longint):longint;

 begin

  if a>b then exit(a);

  exit(b);

 end;

begin

 readln(m,n);

 for i:=1 to n do read(a[i],b[i]);

 for i:=1 to n do b[i]:=b[i]*a[i];

 for i:=1 to n do

  for j:=m downto a[i] do

   f[j]:=max(f[j],f[j-a[i]]+b[i]);

 write(f[m]);

end.