算法笔记5——贪心算法

一、定义

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

二、贪心选择性质

三、贪心求解过程

一般流程

（1）候选集合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;
｝

五、例题分析

例1（活动安排问题）：

设有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;		//活动的编号
};

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

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

//形参数组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;
　　　　}
}

例2（背包问题）：

代码如下：

//形参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;
}

例3（最优装载问题）：

有一批集装箱要装上一艘载重量为c的轮船，其中集装箱i的重量为wi。最优装载问题要求确定在装载体积不受限制的情况下，将尽可能多的集装箱装上轮船。

代码如下：

struct load {
	int index;		//集装箱编号
	int w;			//集装箱重量
}box[1001];

bool cmp (load a, load b) {
	if (a.w<b.w) return true;
	else return false;
}
stable_sort(box, box+n+1, cmp);

while (scanf("%d%d", &c, &n)!=EOF)
{
　　memset(box, 0, sizeof(box));
　　memset(x, 0, sizeof(x));
　　for (int i=1; i<=n; i++) 
　　{
　　　　scanf("%d", &box[i].w);
　　　　box[i].index = i;
　　}
　　//按集装箱的重量升序排序
　　stable_sort(box, box+n+1, cmp);
　　if (box[1].w>c) {
　　　　printf("No answer!\n");
　　　　continue;
　　}
//贪心算法的实现，重量最轻者先装载
　　int i;
　　for (i=1; i<=n && box[i].w<=c; i++)
　　{
　　　　x[box[i].index] = 1;
　　　　c -= box[i].w;
　　}
　　//输出装载的集装箱数量
　　printf("%d\n", i-1);
　　//输出装载的集装箱编号
　　for (i=1; i<=n; i++)
　　　　if (x[i]) printf("%d ", i);
　　printf("\n");
｝