实验四回溯算法和分支限界法.docxVIP

实验四回溯算法和分支限界法 0-1背包问题 、实验冃的 1、 掌握0?1背包问题的回溯算法; 2、 进一步掌握回溯算法。 二、实验内容 给定n和物品和一人背包，物品i的重量是wi,其价值为vi,问如何选择装 入背包的物品，使得装入背包的物品的总价值最大？ 三、实验步骤 1、代码 // HS_ALGcpp : Defines the entry point for the console application. // #includeiostream #includealgorithm using namespace std; // 物体结构体 typedef struct { float w; 〃物品重量 float p; 〃物品价值 float v; 〃背包体积 int id; 〃物品个数 } OBJECT; 〃比较两物品体积bool cmp(OBJECT a, OBJECT b){ 〃比较两物品体积 return a.vb.v; int n, bool} float knapsack_back(OBJECT ob[], float M, int n, bool x[]){ 〃回溯法 int i,k; float w_cur, p_total, p_cur, w_est, p_est; bool *y = new bool[n+l]; // 计算物体的价值重量比 for(i=0; i=n; i++){ ob[i].v = ob[i].p/ob[i].w; y[i] = false; } // 按照物体的价值重量比降序排列 sort(ob, ob+n, cmp); // 初始化当前背包中的价值、重量 w_cur = p_cur 二 p_total = 0; k = 0; while(k=0){ w_est = w_cur; p_est = p_cur; // 沿当前分支可能取得的最大价值 for( i=k; in; i++) { w_est += ob[i].w; if(w_estM){ p_est += ob[i].p; }else{ p_est += ((M-w_est+ob[i].w)/ob[i].w)*ob[i]?p; break; } // 估计值大于上界 if(p_estp_total){ for(i=k; in; i++){ if(w_cur+ob| i]? wv 二 M) { // 可装入第i个物体 w_cur = w_cur + ob[i].w; p_cur = p_cur + ob[i].p; y[i] = true; }else{ // 不能装入第i个物体 y[i] = false; break; } } if(i=n){ // n个物体已经全部装入 if(p_curp_total){ // 更新当前上限 p_total = p_cur; k = n; // 保存可能的解 for(i=0; in; i++){ x[i] = y[i]； } } }else{ // 继续装入物体 k = i+l; } }else{ // 估计值小于上界时 while((i=0)(! y [i]))i—; // 沿着右分支结点方 向回溯直到左分支结点 //到达根结点算法结束 //到达根结点算法结束 //修改当前值 else{ w_cur -二 ob[i].w; p_cur -= ob[i].p; y[i] = false; k = i+1; 搜索右分支子树 } } } //delete y; return p_total; } int main(){ int n; float m; cout?H请输入背包载重 cin?m; cout?n请输入物品个数： cin?n; OBJECT* ob 二 new OBJECT[n]; { cout?HW输入物品的重量、价格:H?endl; for(int i=0; in; i++){ cin?ob[i]. w?ob[i] .p; ob[i].id = i+1; } bool* x = new boolfn]; float v = knapsack_back(ob, m, n, x); cout?n最优方案：n?endl; for(int i=0; in; i++){ if(x[i])cout?n[n?ob[i].idn] n; } cout?endl; cout?n物品总价值为：n?v?endl; } return 0; } 2、结果 执行成功. p MC:\ Documents and SettingsAdministrator廩面\算法分析与设计\程序文件\HW 请输入物品个数^ 8 请输入物品的重量、价格: 2 12 5 56 78 25 23 45 78 12 10 23 1