天天看点

算法设计与分析_贪心算法①

贪心算法

①在求最优解问题的过程中,依据某种贪心标准,从问题的初始状态出发,直接去求每一步的最优解,通过若干次的贪心选择,最终得出整个问题的最优解,这种求解方法就是贪心算法。

②从贪心算法的定义可以看出,贪心算法不是从整体上考虑问题,它所做出的选择只是在某种意义上的局部最优解,而由问题自身的特性决定了该题运用贪心算法可以得到最优解。

③如果一个问题可以同时用几种方法解决,贪心算法应该是最好的选择之一。

贪心算法的性质:

贪心选择性质是指所求问题的整体最优解可以通过一系列局部最优的选择,即贪心选择来达到。

i.这是贪心算法可行的第一个基本要素,也是贪心算法与动态规划算法的主要区别。

(1)在动态规划算法中,每步所做的选择往往依赖于相关子问题的解,因而只有在解出相关子问题后,才能做出选择。

(2)在贪心算法中,仅在当前状态下做出最好选择,即局部最优选择,然后再去解出这个选择后产生的相应的子问题。

最优子结构性质:

①当一个问题的最优解包含其子问题的最优解时,称此问题具有最优子结构性质。

i.运用贪心策略在每一次转化时都取得了最优解。问题的最优子结构性质是该问题可用贪心算法或动态规划算法求解的关键特征。

②贪心算法的每一次操作都对结果产生直接影响,而动态规划则不是。

i.贪心算法对每个子问题的解决方案都做出选择,不能回退;动态规划则会根据以前的选择结果对当前进行选择,有回退功能。

ii.动态规划主要运用于二维或三维问题,而贪心一般是一维问题。

贪心算法的求解过程:

使用贪心算法求解问题应该考虑如下几个方面:

(1)候选集合A:为了构造问题的解决方案,有一个候选集合A作为问题的可能解,即问题的最终解均取自于候选集合A。

(2)解集合S:随着贪心选择的进行,解集合S不断扩展,直到构成满足问题的完整解。

(3)解决函数solution:检查解集合S是否构成问题的完整解。

(4)选择函数select:即贪心策略,这是贪心法的关键,它指出哪个候选对象最有希望构成问题的解,选择函数通常和目标函数有关。

(5)可行函数feasible:检查解集合中加入一个候选对象是否可行,即解集合扩展后是否满足约束条件。

贪心算法的一般流程:

//A是问题的输入集合即候选集合
Greedy(A)
{
  S={ };           //初始解集合为空集
  while (not solution(S))  //集合S没有构成问题的一个解
  {
    x = select(A);     //在候选集合A中做贪心选择
    if feasible(S, x)    //判断集合S中加入x后的解是否可行
      S = S+{x};
      A = A-{x};
  }
  return S;
}

           

活动安排问题

活动安排问题就是要在所给的活动集合中选出最大的相容活动子集合,是可以用贪心算法有效求解的很好例子。

该问题要求高效地安排一系列争用某一公共资源的活动。

贪心算法提供了一个简单、漂亮的方法使得尽可能多的活动能兼容地使用公共资源

问题描述:

①设有n个活动的集合E={1,2,…,n},其中每个活动都要求使用同一资源,如演讲会场等,而在同一时间内只有一个活动能使用这一资源。

②每个活动i都有一个要求使用该资源的起始时间si和一个结束时间fi,且si<fi。如果选择了活动i,则它在半开时间区间[si ,fi )内占用资源。若区间[si ,fi )与区间[sj,fj )不相交,则称活动i与活动j是相容的。当 si ≥ fj 或 sj ≥ fi 时,活动i与活动j相容。

③活动安排问题就是在所给的活动集合中选出最大的相容活动子集合。

算法设计与分析_贪心算法①
//数据结构
struct action{
	int s;			//起始时间
	int f;			//结束时间
	int index;		//活动的编号
};
//活动的集合E记为数组:
action a[1000];

           

按活动的结束时间升序排序

排序比较因子:

bool cmp(const action &a, const action &b)
{
	if (a.f<=b.f) return true;
	return false;
}
           

使用标准模板库函数排序(下标0未用):

sort(a, a+n+1, cmp);

算法:

//形参数组b用来记录被选中的活动
void GreedySelector(int n, action a[], bool b[])
{
  b[1] = true;     //第1个活动是必选的
  //记录最近一次加入到集合b中的活动
  int preEnd = 1;
  for(int i=2; i<=n; i++)
    if (a[i].s>=a[preEnd].f)
    {
      b[i] = true;
      preEnd = i;
    }
}
           

背包问题

1.给定一个载重量为M的背包,考虑n个物品,其中第i个物品的重量 ,价值wi (1≤i≤n),要求把物品装满背包,且使背包内的物品价值最大。

2.有两类背包问题(根据物品是否可以分割),如果物品不可以分割,称为0—1背包问题(动态规划);如果物品可以分割,则称为背包问题(贪心算法)。

算法设计与分析_贪心算法①
算法设计与分析_贪心算法①

数据结构

struct bag{

int w; //物品的重量

int v; //物品的价值

double c; //性价比

}a[1001]; //存放物品的数组

排序因子(按性价比降序):

bool cmp(bag a, bag b){

return a.c >= b.c;

}

使用标准模板库函数排序(最好使用stable_sort()函数,在性价比相同时保持输入的顺序):

sort(a, a+n, cmp);

计算背包问题的贪心算法

//形参n是物品的数量,c是背包的容量M,数组a是按物品的性价比降序排序
double knapsack(int n, bag a[], double c)
{
  double cleft = c;        //背包的剩余容量
  int i = 0;
  double b = 0;          //获得的价值
  //当背包还能完全装入物品i
  while(i<n && a[i].w<cleft)
  {
    cleft -= a[i].w;
    b += a[i].v;
    i++;
  }
  //装满背包的剩余空间
  if (i<n) b += 1.0*a[i].v*cleft/a[i].w;
  return b;
}

           

如果要获得解向量 X={x1, x2…xn},则需要在数据结构中加入物品编号:

如果要获得解向量                         ,则需要在数据结构中加入物品编号:
struct bag{
	int w;
	int v;
	double x;		//装入背包的量,0≤x≤1
	int index;		//物品编号
	double c;
}a[1001];
           

计算背包问题的贪心算法,同时得到解向量

double knapsack(int n, bag a[], double c)
{
  double cleft = c;
  int i = 0;
  double b = 0;
  while(i<n && a[i].w<=cleft)
  {
    cleft -= a[i].w;
    b += a[i].v;
    //物品原先的序号是a[i].index,全部装入背包
    a[a[i].index].x = 1.0;
    i++;
  }
  if (i<n) {
    a[a[i].index].x = 1.0*cleft/a[i].w;
    b += a[a[i].index].x*a[i].v;
  }
  return b;
}
           

继续阅读