解法(背包DP问题)
(下为转)
其实感觉 像此题这种类型的并不属于dp范畴 虽然程序看起来使用的是递推这一过程,但总不能说开个二重循环就是dp吧 如果只从求解上来讲(不考虑数据值的范围), 只有枚举这唯一途径, 而此题的复杂度为O(20^20), 大约等于 10^26, 微机是很难在短时间内求解的。 而此题还有另一个限制, 就是数据值的范围, 虽然状态数那么多, 但是它们的范围有限这个时候, 状态产生大量重复, 于是, 程序可以优化为压缩这些重复状态, 从而减少开销 在实现的过程中, 确实划分了阶段, 但这个过程更像是枚举, 优化过后的枚举 在朴素枚举当中, 阶段数与复杂度呈指数关系, 而在此题这种特殊情况下, 每个阶段数的最大不同状态数是有限的, 压缩以后就只是常数关系了 所以那 'dp'方程更像是一个'压缩状态'过程, 并没有求最优这一意思, 最后输出是状态数, 并不是最优解 再看一些常见的dp问题, 它们的意义很明确, 一定是求某种最优的结果, 其中普遍存在max,min 最为关键的一点, 凡是能够被dp的问题, 都是p类问题, 它们的复杂度都为多项式表达式 而该题类似的问题, 它们似乎往往都是npc问题, 它们很难有特别有效的算法 像本题这种就是通过狭小的数据范围来进行优化, 但仍然逃不脱枚举这一本质, 如果稍微加大数据范围, 您看看还有更好的方法么?
这道题目是按照平衡度来考虑的,即挂完砝码后整个天平的值,从而可以划分为逐次挂砝码的阶段,则此次挂砝码和上一次的过程就存在了重叠子问题。运用DP就可以
有效的减少状态数。
我们使用DP[i][j]表示挂第i个砝码时平衡度为j的状态数。显然i的范围是1-20,对j来说,砝码重量是1-25,且不同砝码重量两两不同,而挂钩的范围是-15-+15,最多有20个砝码,故取值范围为(5+6+...+25)*15,即-4500~+4500,由此可以确定需要开出的table大小。
代码如下:
1 #include<iostream> 2 using namespace std; 3 int hook[20]; 4 int weight[20]; 5 int balance[22][9002]; 6 7 int main() 8 { 9 int hC,wG,i,j,k,value; 10 cin>>hC>>wG; 11 memset(balance,0,sizeof(balance)); 12 for(i=0;i<hC;i++) 13 { 14 cin>>hook[i]; 15 } 16 for(i=0;i<wG;i++) 17 { 18 cin>>weight[i]; 19 } 20 for(i=0;i<hC;i++) 21 { 22 balance[0][weight[0]*hook[i]+4500]++; 23 } 24 for(i=1;i<wG;i++) 25 { 26 for(j=0;j<9002;j++) 27 { 28 if(balance[i-1][j]!=0) 29 { 30 for(k=0;k<hC;k++) 31 { 32 value=weight[i]*hook[k]; 33 balance[i][j+value]+=balance[i-1][j]; 34 } 35 } 36 } 37 } 38 cout<<balance[wG-1][0+4500]<<endl; 39 return 0; 40 }
这种方法的复杂度为 O(C*G*9000)