1、贪心算法
(1)原理:在对问题求解时,总是做出在当前看来是最好的选择。也就是说,不从整体最优上加以考虑,他所做出的仅是在某种意义上的局部最优解。贪心算法不是对所有问题都能得到整体最优解,但对范围相当广泛的许多问题他能产生整体最优解或者是整体最优解的近似解。
(2)特性:贪心算法采用自顶向下,以迭代的方法做出相继的贪心选择,每做一次贪心选择就将所求问题简化为一个规模更小的子问题,通过每一步贪心选择,可得到问题的一个最优解,虽然每一步上都要保证能获得局部最优解,但由此产生的全局解有时不一定是最优的,所以贪婪法不要回溯。能够用贪心算法求解的问题一般具有两个重要特性:贪心选择性质和最优子结构性质。
1)贪心选择性质
所谓贪心选择性质是指所求问题的整体最优解可以通过一系列局部最优的选择,即贪心选择来达到。这是贪心算法可行的第一个基本要素。贪心算法则通常以自顶向下的方式进行,以迭代的方式作出相继的贪心选择,每作一次贪心选择就将所求问题简化为规模更小的子问题。
对于一个具体问题,要确定它是否具有贪心选择性质,必须证明每一步所作的贪心选择最终导致问题的整体最优解。证明的大致过程为:首先考察问题的一个整体最优解,并证明可修改这个最优解,使其以贪心选择开始。做了贪心选择后,原问题简化为规模更小的类似子问题。然后用数学归纳法证明通过每一步做贪心选择,最终可得到问题的整体最优解。其中,证明贪心选择后的问题简化为规模更小的类似子问题的关键在于利用该问题的最优子结构性质。
2)最优子结构性质
当一个问题的最优解包含其子问题的最优解时,称此问题具有最优子结构性质。
(3)贪心算法与动态规划算法的差异:
动态规划和贪心算法都是一种递推算法,均有最优子结构性质,通过局部最优解来推导全局最优解。两者之间的区别在于:贪心算法中作出的每步贪心决策都无法改变,因为贪心策略是由上一步的最优解推导下一步的最优解,而上一部之前的最优解则不作保留,贪心算法每一步的最优解一定包含上一步的最优解。动态规划算法中全局最优解中一定包含某个局部最优解,但不一定包含前一个局部最优解,因此需要记录之前的所有最优解。
(4)基本思路:
1)建立数学模型来描述问题。
2)把求解的问题分成若干个子问题。
3)对每一子问题求解,得到子问题的局部最优解。
4)把子问题的解局部最优解合成原来解问题的一个解。
2、活动安排问题
活动安排问题就是要在所给的活动集合中选出最大的相容活动子集合,是可以用贪心算法有效求解的很好例子。该问题要求高效地安排一系列争用某一公共资源的活动。贪心算法提供了一个简单、漂亮的方法使得尽可能多的活动能兼容地使用公共资源。
问题描述:
设有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相容。活动安排问题就是要在所给的活动集合中选出最大的相容活动子集合。
求解思路:
将活动按照结束时间进行从小到大排序。然后用i代表第i个活动,s[i]代表第i个活动开始时间,f[i]代表第i个活动的结束时间。按照从小到大排序,挑选出结束时间尽量早的活动,并且满足后一个活动的起始时间晚于前一个活动的结束时间,全部找出这些活动就是最大的相容活动子集合。事实上系统一次检查活动i是否与当前已选择的所有活动相容。若相容活动i加入已选择活动的集合中,否则,不选择活动i,而继续下一活动与集合A中活动的相容性。若活动i与之相容,则i成为最近加入集合A的活动,并取代活动j的位置。
下面给出求解活动安排问题的贪心算法,各活动的起始时间和结束时间存储于数组s和f中,且按结束时间的非减序排列。如果所给的活动未按此序排列,可以用O(nlogn)的时间重排。具体代码如下:
//4d1 活动安排问题 贪心算法
#include "stdafx.h"
#include <iostream>
using namespace std;
template<class Type>
void GreedySelector(int n, Type s[], Type f[], bool A[]);
const int N = 11;
int main()
{
//下标从1开始,存储活动开始时间
int s[] = {0,1,3,0,5,3,5,6,8,8,2,12};
//下标从1开始,存储活动结束时间
int f[] = {0,4,5,6,7,8,9,10,11,12,13,14};
bool A[N+1];
cout<<"各活动的开始时间,结束时间分别为:"<<endl;
for(int i=1;i<=N;i++)
{
cout<<"["<<i<<"]:"<<"("<<s[i]<<","<<f[i]<<")"<<endl;
}
GreedySelector(N,s,f,A);
cout<<"最大相容活动子集为:"<<endl;
for(int i=1;i<=N;i++)
{
if(A[i]){
cout<<"["<<i<<"]:"<<"("<<s[i]<<","<<f[i]<<")"<<endl;
}
}
return 0;
}
template<class Type>
void GreedySelector(int n, Type s[], Type f[], bool A[])
{
A[1]=true;
int j=1;//记录最近一次加入A中的活动
for (int i=2;i<=n;i++)//依次检查活动i是否与当前已选择的活动相容
{
if (s[i]>=f[j])
{
A[i]=true;
j=i;
}
else
{
A[i]=false;
}
}
}
3.背包问题
注意:这是背包问题,而不是0-1背包问题,背包问题可以用贪心算法进行求解,但0-1无法用贪心算法求解,需要用动态规划算法求解;
首先对贪心算法做一下总结,以及它与动态规划算法的区别:
贪心算法两个最重要的性质:
(1)贪心选择性质;
(2)最优子结构性质;
其中,贪心选择性质:自顶向下进行决策,每次做出的决策都是局部最优解,且每次做出决策后问题规模都变小了;最优子结构性质:即问题的最优解结构中包含子问题的最优解;
动态规划算法的两个最重要的性质:
(1)重叠子问题性质;
(2)最优子结构性质;
其中最优解子结构性质和贪心算法相似,唯一不同的是重叠子问题性质,因为动态规划算法是自底向上的算法,它需要首先将原始问题分解为若干个相互有联系的子问题,在计算的时候有的子问题可能会被计算很多次,所以动态规划算法会将这些子问题的解存在一个表格中,使得最终对于这些子问题只需要求解一次(可以使原来需要再指数时间内解决的问题可以在多项式问题中得到解决)
背包问题求解代码如下:
//背包问题,使用贪心算法进行求解
//======================================================
#include <iostream>
#include "merge_sort.cpp"
using namespace std;
void init_data(float *v, float *w, float *x, int n) //初始化数据
{
cout << "请输入每类物体的价值:" << endl;
for (int i = 0; i < n; i++)
{
cin >> v[i];
}
cout << "请输入每类物体的重量: " << endl;
for (int i = 0; i < n; i++)
{
cin >> w[i];
}
for (int i = 0; i < n; i++)
x[i] = 0;
}
void Knapsack(int n, float M, float *v, float *w, float *x)
//n是物体的种类数,M是背包容量,v是每类物体的价值,w是每类物体的重量,x是每类物体装入的份额,属于[0,1]
{
int i = 0;
float c = M;
MergeSort(w , v , n); //v[i]/w[i]是每一类物体的单位重量价值,然后对它们进行降序排序
for (i = 0; i < n; i++)
{
if (w[i] > c)
break;
x[i] = 1;
c -= w[i];
}
if (i < n)
x[i] = c / w[i];
}
int main(void)
{
float M = 0.0; //背包容量
cout << "请输入背包容量: ";
cin >> M;
int n = 0; //物体数量
cout << "\n请输入物体数量 n : ";
cin >> n;
float *v = new float[n];
float *w = new float[n];
float *x = new float[n];
init_data(v, w, x, n); //初始化数据
Knapsack(n, M, v, w, x);
cout << "排好序的w[i]: " << endl;
for (int i = 0; i < n; i++)
{
cout << w[i] << " ";
}
cout << "\n\n输出最后的决策x[i] : " << endl;
for (int i = 0; i < n; i++)
{
cout << x[i] << " ";
}
/*MergeSort(w, v, n);
cout << "输出排好序的w[i] : " << endl;
for (int i = 0; i < n; i++)
{
cout << w[i] << " ";
}*/
system("pause");
delete v , w , x;
return 0;
}
(其中使用的排序算法基于合并排序算法)