0-1背包问题的分支限界法源代码

合集下载
  1. 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
  2. 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
  3. 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
{
int i;
n=nn;
c=cc;
p=new int[n];
w=new int[n];
M=new int[n];
for(i=0;i<n;i++)
{
p[i]=pp[i];
w[i]=ww[i];
M[i]=i; //用M数组记录大小顺序关系
}
front=new node[1];
front->next=NULL;
cin>>w[i];
cout<<endl;
cout<<"请输入这"<<n<<"个物品的价值P:"<<endl;
for(i=0;i<n;i++)
cin>>p[i];
Knap knbag(p,w,c,n);
knbag.solvebag();
getch();
return 0;
}
//#include "stdafx.h"
#include<iostream>
#include<cstdio>
#include<conio.h>
#include<iomanip>
using namespace std;
int *x;
struct node //结点表结点数据结构
{
node *parent;//父结点指针
{
minl=1.0*p[i]/w[i];
k=0;
for(j=1;j<=n-i;j++)
{
if(minl<1.0*p[j]/w[j])
{
minl=1.0*p[j]/w[j];
swap(p[k],p[j]);
swap(w[k],w[j]);
swap(M[k],M[j]);
k=j;
}
}
}
}
Knap::Knap(int *pp,int *ww,int cc,int nn)
{
ubb=Bound(i,aa->cw+w[i],aa->cp+p[i]);
addnode(nnoder(aa,1,ubb));//将结点加入到优先队列中
}
ubb=ubb=Bound(i,aa->cw,aa->cp);
if(ubb>lbestp)
addnode(nnoder(aa,0,ubb));
}
void Knap::display()
{
int i;
cout<<endl;
cout<<"当前最优价值为:"<<lbestp<<endl;
for(i=n;i>=1;i--)
{
x[M[i-1]]=bestp->bag;
bestp=bestp->parent;
}
cout<<"变量值x = ";
for(i=1;i<=n;i++)
node *next; //后继结点指针
int level;//结点的层
int bag;//节点的解
int cw;//当前背包装载量
int cp;//当前背包价值
float ub; //结点的上界值
};
//类Knap中的数据记录解空间树中的结点信息,以减少参数传递及递归调用所需的栈空间
class Knap
i=aa->level;
//当叶子节点处的解>最优解时,更新最优解
if(i==n-1)
{
if(aa->cw+w[i]<=c&&(long)(aa->cp+p[i])>lbestp)
{
lbestp=aa->cp+p[i];
bestp=nnoder(aa,1,(float)lbestp);//将一个新的结点插入到子集树和优先队列中
float b=(float)cp; //价值上界
//以物品单位重量价值减序装填剩余容量
while (i<n&&w[i]<=cleft)
{
cleft-=w[i];
b+=p[i];
i++;
}
//装填剩余容量装满背包
if (i<n) b+=1.0*p[i]/w[i]*cleft;
return b;
}
//计算最优值和变量值
{
if(p->ub<ub){no->next=p;next1->next=no;break;}
next1=p;
p=p->next;
}
if(p==NULL){next1->next=no;}
}
//计算结点所相应价值的上界
float Knap::Bound(int i,int cw,int cp)
{
int cleft=c-cw; //剩余容量
~Knap();
float Bound(int i,int cw,int cp);//计算上界限
node *nnoder(node *pa,int ba,float uub);//生成一个结点ba=1生成左节点ba=0生成右节点
void addnode(node *nod);//向队列中添加活结点
void deletenode(node *nod);//将结点从队列中删除
cout<<endl;
cout<<"请输入背包容量C = ";
cin>>c;
cout<<endl;
x=new int[n]; //变量值
p=new int[n]; //物品价值
w=new int[n]; //物品重量
cout<<"请分别输入这"<<n<<"个物品的重量W:"<<endl;
for(i=0;i<n;i++)
}
display();
}
int main()
{
int c,n;
int i=0;
int *p;
int *w;
cout<<endl;
cout<<"|***********分支限界法解0-1背包问题***********| "<<endl;
cout<<endl;
cout<<"请输入物品数量n = ";
cin>>n;
lbestp=0;
bestp=new node[1];
bestp=NULL;
Sort();
}
Knap::~Knap()
{
delete []first;
delete []front;
delete []bestp;
delete []p;
delete []w;
}
//取上限最大结点
node *Knap::nextnode()
{
private:
struct node *front, //队列队首
*bestp,*first; //解结点、根结点
int *p,*w,n,c,*M;//背包价值、重量、物品数、背包容量、记录大小顺序关系
long lbestp;//背包容量最优解
public:
void Sort();
Knap(int *pp,int *ww,int cc,int nn);
{
node *p=front->next;
front->next=p->next;
return(p);
}
//将一个新的结点插入到子集树和优先队列中
node * Knap::nnoder(struct node *pa,int ba,float uub)
{//生成一个新结点
node * nodell=new(node);
nodell->parent=pa;
nodell->next=NULL;
nodell->level=(pa->level)+1;
nodell->bag=ba;
nodell->ub=uub;
if(ba==1)
{
nodell->cw=pa->cw+w[pa->level];
nodell->cp=pa->cp+p[pa->level] ;
struct node *nextnode(); //取下一个节点
void display(); //输出结果
void solvebag(); //背包问题求解
};
//按物品单位重量的价值排序
void Knap::Sort()
{
int i,j,k,kkl;
float minl;
for(i=1;i<n;i++)
}
else
{
nodell->cw=pa->cw;
nodell->cp=pa->cp;
}
return(nodell);
}
//将结点加入优先队列
void Knap::addnode(node *no)
{
node *p=front->next,*next1=front;
float ub=no->ub;
while(p!=NULL)
}
if((long)(aa->cp)>lbestp)
{
lbestp=aa->cp;
bestp=nnoder(aa,0,(float)lbestp);
}
}
//非叶子结点,递归调用Bound函数计算上界
if(i<n-1)
{
if(aa->cw+w[i]<=c&&Bound(i+1,aa->cw+w[i],aa->cp+p[i])>(float)lbestp)
first->level=0; //用level记录结点的层
first->cw=0;
first->cp=0;
first->bag=0;
ubb=Bound(0,0,0);
first->ub=ubb;
front->next=first;
while(front->next!=NULL)t;x[i-1];
cout<<endl;
}
//背包问题求解
void Knap::solvebag()
{
int i;
float ubb;
node *aa; //记录上限最大结点
first=new node[1]; //根结点
first->parent=NULL;
first->next=NULL;
相关文档
最新文档