实验名称
回溯法求解0-1背包。
实验目的
回溯法求解0-1背包。
实验原理
给定任意几组数据,利用回溯算法的限界与剪枝的思想,求解0-1背包问题并输出答案。
实验步骤
① 将用例数据从文件中读取,初始化数据;
② 将物品按照单位重量的价值从大到小排序;
③ 从第0层开始进行递归函数Backtrack,
当i>n时,算法搜索至叶子结点,得到一个新的物品装包方案,此时算法适时更新当前的最优价值和最优解;
当i<n时,当前扩展结点位于排列树的第(i-1)层,此时算法选择下一个要安排的物品,以深度优先方式递归的对相应的子树进行搜索,对不满足上界约束的结点,则剪去相应的子树;
④当整个解空间树遍历完成则输出最优值与最优解。
时间复杂度分析
该问题为回溯法的子集树问题,对每一个物品都有放入与不放两种决策,因此遍历解空间树的时间复杂度未O(2^n)。
实验心得
通过这次实验,我回顾了回溯算法的基本原理,为更难更深层次的问题的回溯法求解打下基础。
#include <iostream>
#include <stdio.h>
#include<algorithm>
#include <time.h>
#include <windows.h>
#include<fstream>
using namespace std;
ifstream ifile("input.txt");
ofstream ofile("output.txt");
int n;//物品数量
double c;//背包容量
double cw = 0.0;//当前背包重量 current weight
double cp = 0.0;//当前背包中物品总价值 current value
double bestp = 0.0;//当前最优价值best price
int x[100];//设置是否装入,为1的时候表示选择该组数据装入,为0的表示不选择该组数据
int bestx[100];//记录最优解
struct WP{
int id; //编号
double pre; //单价
double v; //价值
double w; //重量
};
WP wp[1000];
double bound(int i);
//回溯函数
void backtrack(int i)
{ //i用来指示到达的层数(第几步,从0开始),同时也指示当前选择完几个物品
if(i>=n) //递归结束的判定条件
{
if (bestp<cp){
for(int j=0;j<n;j++){
bestx[j]=x[j];
}
bestp = cp;
}
return;
}
//如若左子节点可行,则直接搜索左子树;
//对于右子树,先计算上界函数,以判断是否将其减去
if(cw+wp[i].w<=c)//将物品i放入背包,搜索左子树
{
cw+=wp[i].w;//同步更新当前背包的重量
cp+=wp[i].v;//同步更新当前背包的总价值
x[i]=1;
backtrack(i+1);//深度搜索进入下一层
cw-=wp[i].w;//回溯复原
cp-=wp[i].v;//回溯复原
x[i]=0;
}
if(bound(i+1)>bestp)//如若符合条件则搜索右子树
backtrack(i+1);
}
//计算上界函数,功能为剪枝
double bound(int i)
{ //判断当前背包的总价值cp+剩余容量可容纳的最大价值<=当前最优价值
double leftw= c-cw;//剩余背包容量
double b = cp;//记录当前背包的总价值cp,最后求上界
//以物品单位重量价值递减次序装入物品
while(i<n && wp[i].w<=leftw)
{
leftw-=wp[i].w;
b+=wp[i].v;
i++;
}
//装满背包
if(i<n)
b+=wp[i].pre*leftw;
return b;//返回计算出的上界
}
bool cmp(WP a,WP b)
{
return a.pre>b.pre;
}
void init(){
ifile>>n>>c;
for(int i=0;i<n;i++){
ifile>>wp[i].w>>wp[i].v;
wp[i].pre=wp[i].v/wp[i].w; //计算单位价值(单位重量的物品价值)
wp[i].id=i;
}
sort(wp,wp+n,cmp);
cw=0;cp=0;bestp=0;
}
int main()
{
LARGE_INTEGER frequency;
double v,beginoftime,endoftime,dt,t;//时钟的频率,起始,实际时间差,实际时间
//样例
init();
backtrack(0);
ofile<<bestp<<endl;
for(int i=0;i<n;i++)
{
if(bestx[i]==1)
ofile<<wp[i].id<<' ';
}
//10
ifile.close();
ifile.open("10.txt");
cout<<"10 SISE:\n";
init();
QueryPerformanceFrequency(&frequency);//获得时钟频率
v=(double)frequency.QuadPart;//取得频率*/
QueryPerformanceCounter(&frequency);
beginoftime=frequency.QuadPart;//获得初始值
backtrack(0);
QueryPerformanceCounter(&frequency);
endoftime=frequency.QuadPart;//获得终止值
dt=(double)(endoftime-beginoftime);//差值
cout<<bestp;
t=dt/v;//差值除以频率得到时间
cout<<endl<<"Time is "<<t*1000<<"ms\n";
//100
ifile.close();
ifile.open("100.txt");
cout<<"100 SISE:\n";
init();
QueryPerformanceFrequency(&frequency);//获得时钟频率
v=(double)frequency.QuadPart;//取得频率*/
QueryPerformanceCounter(&frequency);
beginoftime=frequency.QuadPart;//获得初始值
backtrack(0);
QueryPerformanceCounter(&frequency);
endoftime=frequency.QuadPart;//获得终止值
dt=(double)(endoftime-beginoftime);//差值
cout<<bestp;
t=dt/v;//差值除以频率得到时间
cout<<endl<<"Time is "<<t*1000<<"ms\n";
//1000
ifile.close();
ifile.open("1000.txt");
cout<<"1000 SISE:\n";
init();
QueryPerformanceFrequency(&frequency);//获得时钟频率
v=(double)frequency.QuadPart;//取得频率*/
QueryPerformanceCounter(&frequency);
beginoftime=frequency.QuadPart;//获得初始值
backtrack(0);
QueryPerformanceCounter(&frequency);
endoftime=frequency.QuadPart;//获得终止值
dt=(double)(endoftime-beginoftime);//差值
cout<<bestp;
t=dt/v;//差值除以频率得到时间
cout<<endl<<"Time is "<<t*1000<<"ms\n";
return 0;
}
#include <iostream>
#include <time.h>
#include <fstream>
//***本程序用于产生随机数据
using namespace std;//RAND_MAX=32767
int w[10000]={0};
int v[10000]={0};
int c,n;
int Size[4]={10,100,10000};//文件大小
void print(ofstream &outfile,int n,int c)//输出到文件
{
outfile<<n<<" "<<c<<endl;
for(int i=0;i<n;i++)
{
outfile<<w[i]<<' '<<v[i]<<endl;
}
}
int main()
{
int n=0;
//ifstream
ofstream out_10("10.txt"),out_100("100.txt"),out_10000("10000.txt");//输入代查找数据;
srand(time(NULL));//时间种子
//生成测试文件
for(int i=0;i<Size[0];i++)
{
w[i]=rand()%10+1;//得到[0,10)内的数
v[i]=rand()%200+1;
}
print(out_10,Size[0],20);
for(int i=0;i<Size[1];i++)
{
w[i]=rand()%10+1;
v[i]=rand()%200+1;
}
print(out_100,Size[1],20);
for(int i=0;i<Size[2];i++)
{
w[i]=rand()%10+1;
v[i]=rand()%200+1;
}
print(out_10000,Size[2],20);
}