第十二周算法心得体会

回溯算法

回溯算法的基本思想

有“万能解题法”之称，用它可以系统地搜索一个问题的所有解或者任意解。

适用情况

①需要找出问题的解集

②要求确定满足某些约束条件的最佳解或所有解

回溯算法解决问题的3个步骤

①针对所给问题，定义问题的解空间

②确定易于搜索的解空间结构，一般为树的形式

③以深度优先方式搜索解空间，并在搜索过程中利用剪枝函数避免无效地搜索。

子集树：当给定的问题是从n个元素的集合中找出满足某种性质的子集或组合时，相应的解空间树为子集树。

排列树：当给定的问题是确定n个元素满足某种性质的排列或顺序时，相应的解空间树为排列树。

约束函数在扩展节点处剪去不满足约束条件的子树。

限界函数剪去不能得到最优解的子树。

基本的回溯法正确求解0-1背包问题

代码实现

#include <stdio.h>
#include <stdlib.h>
const int N=10;//物品数量
const int bag=N*2;//背包承重量
int max_value=0;//记录能获得的最大价值
int backpack(int a[],int v[],int w[],int n,int t)
{
	int i,max_v=0,max_w=0;
	if(t>n-1)
	{
		printf("找到一个方案");
		for(i=0;i<n;i++)
		{
			printf("%d",a[i]);
			//此处应根据a[]中记录的0和1进行背包价值和重量计算
			max_v=max_v+a[i]*v[i];
			max_w=max_w+a[i]*v[i];
		}
		if(max_w>bag) printf("\n该方案总价%d,总重%d，超重！不可行！！！\n",max_v,max_w);
		else 
		{			
			printf("\n该方案总价%d,总重%d，可行\n",max_v,max_w);
			//此处应考虑是否更新max_value的值
			if(max_value<max_v){
				max_value=max_v;
			}
		}
    }
    else
    {
        for(i=1;i>=0;i--)
        {
            a[t]=i;//记录当前物品选或不选
			backpack(a,v,w,n,t+1);//递归到下一层
        }
    }
	return max_value;//返回最大价值
}
int main()
{
	int a[N],v[N],w[N],i;
	printf("背包最大承重%d公斤\n",bag);
	for(i=0;i<N;i++)//随机生成价值和重量
	{
		v[i]=rand()%5+1;
		w[i]=rand()%5+1;
		printf("物品%d,价值%d,重量%d\n",i,v[i],w[i]);
	}
	printf("能获得的最大价值为：%d\n",backpack(a,v,w,N,0));
}

带约束剪枝的回溯法正确求解0-1背包问题

代码实现

#include <stdio.h>
#include <stdlib.h>
const int N=10;//物品数量
const int bag=N;//背包承重量
int max_value=0;//记录能获得的最大价值
bool constraint(int a[],int w[],int t)//约束剪枝判断
{
	int i,max_w=0;
	for(i=0;i<t;i++)
	{
		max_w=max_w+a[i]*w[i];//累加计算第t层节点的重量
		
	}	
	if(max_w<=bag)//加上判断条件)//判断是否超重
		return true;
	else
		return false;
}
int backpack(int a[],int v[],int w[],int n,int t)
{
	int i,max_v=0,max_w=0;
	if(t>n-1)
	{
		printf("找到一个方案");
		for(i=0;i<n;i++)
		{
			printf("%d",a[i]);
			//根据a[]中记录的0和1进行背包价值和重量计算
			max_v=max_v+a[i]*v[i];
			max_w=max_w+a[i]*w[i];		
		}
		if(max_w>bag) printf("\n该方案总价%d,总重%d，超重！不可行！！！\n",max_v,max_w);
		else 
		{
			printf("\n该方案总价%d,总重%d，可行\n",max_v,max_w);
			//此处应考虑是否更新max_value的值
			if(max_value<max_v){
				max_value=max_v;
			}
		}
    }
    else
    {
        for(i=1;i>=0;i--)
        {
            a[t]=i;//记录当前物品选或不选
			if(constraint(a,w,t+1))
				backpack(a,v,w,n,t+1);//递归到下一层
			else
				printf("在第%d层发生了一次剪枝\n",t);
        }
    }
	return max_value;
}
int main()
{
	int a[N],v[N],w[N],i;
	printf("背包最大承重%d公斤\n",bag);
	for(i=0;i<N;i++)//随机生成价值和重量
	{
		v[i]=rand()%5+1;
		w[i]=rand()%5+1;
		printf("物品%d,价值%d,重量%d\n",i,v[i],w[i]);
	}
	printf("能获得的最大价值为：%d\n",backpack(a,v,w,N,0));
}