回溯法求解0-1背包

实验名称
回溯法求解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);
}