贪心算法

贪心算法

在求最优解问题的过程中，依据某种贪心标准，从问题的初始状态出发，直接去求每一步的最优解，通过若干次的贪心选择，最终得出整个问题的最优解，这种求解方法就是贪心算法。
贪心算法不是从整体上考虑问题，它所做出的选择只是在某种意义上的局部最优解，而由问题自身的特性决定了该题运用贪心算法可以得到最优解。

1、背包问题
给定一个载重量为M的背包，考虑n个物品，其中第i个物品的重量 ，价值wi （1≤i≤n），要求把物品装满背包，且使背包内的物品价值最大。
有两类背包问题（根据物品是否可以分割），如果物品不可以分割，称为0—1背包问题（动态规划）；如果物品可以分割，则称为背包问题（贪心算法）。

用贪心算法，按性价比从高到底的顺序选取物品，获得最优值240。由于物品可以分割，剩下的空间装入物品3的一部分，而获得了更好的性能。

#include<iostream>
#include<algorithm>
using namespace std;
struct bag
{
	int value, weight,index;
	double c,x;
}a[1000];
bool cmp(bag a, bag b) {
	return a.c >= b.c;
}
void TxBag(int n,bag a[],double c)
{
	double cleft = c;
	int i = 0;
	double b = 0;
	while (i < n && a[i].weight <= cleft)
	{
		cleft -= a[i].weight;
		b += a[i].value;
		//物品原先的序号是a[i].index，全部装入背包
		a[a[i].index].x = 1.0;
		i++;
	}
	if (i < n) 
	{
		a[a[i].index].x = 1.0 * cleft / a[i].weight;
		b += a[a[i].index].x * a[i].value;
	}
	cout << "最优值为：" << b;
}
int main()
{
	int n=0;
	while (cin >> a[n].weight>>a[n].value)
	{
		a[n].index = n;
		a[n].c = a[n].value / a[n].weight;
		n++;
	}
	sort(a, a + n,cmp);
	TxBag(n, a, 50);
}

2、木棒加工问题
现有n根木棒，已知它们的长度和重量。要用一部木工机一根一根地加工这些木棒。该机器在加工过程中需要一定的准备时间，是用于清洗机器，调整工具和模版的。
木工机需要的准备时间如下：
（1）第一根木棒需要1min的准备时间；
（2）在加工了一根长为l ，重为w的木棒之后，接着加工一根长为l ’ (l ≤ l’ )，重为 w’ ( w≤w’)的木棒是不需要任何准备时间的，否则需要一分钟的准备时间。
给定n根木棒，找到最少的准备时间。
举例：
现在有长和重分别为（4，9），（5，2），（2，1），（3，5）和（1，4）的五根木棒，那么所需准备时间最少为2min，顺序为（1，4），（3，5），（4，9），（2，1），（5，2）。
输入
输入有多组测试例。输入数据的第一行是测试例的个数T。
每个测试例两行：第一行是一个整数n（1≤n≤5000），表示有多少根木棒；第二行包括n×2个整数，表示l1，w1，l2，w2，l3，w3，…，ln，wn，其中li和wi表示第i根木棒的长度和重量。数据由一个或多个空格分隔。
输出
输出是以分钟为单位的最少准备时间，一行一个

本题仅仅使用贪心算法是不够的，排序之后还要使用动态规划的算法。

#include<iostream>
#include<algorithm>
using namespace std;
#define maxN 100
struct stick
{
	int l;//长度
	int w;//重量
};
int cmp(stick a, stick b)
{
	//长度相等时，按重量排序
	if (a.l == b.l) 	return a.w < b.w;
	//优先按长度排序
	else if (a.l < b.l) return true;
	return false;
}
int LIS(int n, stick a[])
{
	//数组b表示木棒分组的序号
	int b[maxN];
	memset(b, 0, sizeof(b));//清空数组b（将数组初始化为0）
	int i, j, k;
	b[0] = 1;
	for (i = 1; i < n; i++)
	{
		//计算第i个木棒的的分组序号
		k = 0;
		for (j = 0; j < i; j++)
			if (a[i].w < a[j].w && k < b[j]) k = b[j];
		b[i] = k + 1;
	}
	//查找最大的分组序号（数组b中的最大值）
	int max = 0;
	for (i = 0; i < n; i++)
		if (b[i] > max) max = b[i];
	return max;
}
int main()
{
	stick data[maxN];
	int T, n,b[maxN];
	cin >> T;
	for (int i = 0; i < T; i++)
	{
		cin >> n;
		for (int j = 0; j < n; j++)
			cin >> data[j].l >> data[j].w;
		sort(data, data + n, cmp);
		b[i] = LIS(n,data);
	}
	for (int i = 0; i < T; i++)
		cout << b[i] << endl;
	return 0;
}

3、鱼塘钓鱼问题
约翰有h（1≤h≤16）个小时的时间，在该地区有n（2≤n≤25）个湖，这些湖刚好分布在一条路线上，该路线是单向的。约翰从湖1出发，他可以在任一个湖结束钓鱼。但他只能从一个湖到达另一个与之相邻的湖，而且不必每个湖都停留。
假设湖i（i＝1～n—1），以5分钟为单位，从湖i到湖i＋1需要的时间用ti（0＜ti≤192）表示。例如t3＝4，是指从湖3到湖4需要花20分钟时间。
已知在最初5分钟，湖i预计钓到鱼的数量为fi（fi≥0）。以后每隔5分钟，预计钓到鱼的数量将以常数di（di≥0）递减。如果某个时段预计钓到鱼的数量小于或等于di，那么在下一时段将钓不到鱼。为简单起见，假设没有其它的钓鱼者影响约翰的钓鱼数量。
输入
对每组测试例，第一行是n，接下来一行是h。下面一行是n个整数fi（1≤i≤n），然后是一行n个整数di（1≤i≤n），最后一行是n—1个整数ti（1≤i≤n—1）。
输入样例

输出
对每个测试例，输出在每个湖上花费的时间，这是约翰要实现钓到最多的鱼的计划（必须使整个计划在同一行输出）。
接下来一行是钓到的鱼的数量：
如果存在很多方案，尽可能选择在湖1钓鱼所耗费的时间，即使有些时段没有钓到鱼；如果还是无法区分，那就尽可能选择在湖2钓鱼所耗费的时间，以此类推。

#include<iostream>
using namespace std;
#define NUM 30
int f[NUM];//每个湖钓到鱼的数量
int d[NUM];//鱼的递减数量
int t[NUM];//湖与湖之间的时间
int plan[NUM];//
int best=0;
void greedy(int pos, int time)
{
	if (time <= 0) return;//时间已经用完
	int i, j;
	int fish[NUM];
	int p[NUM];
	int t = 0;
	for (i = 0; i < pos; ++i)
		fish[i] = f[i];
	memset(p, 0, sizeof(p));
	for (i = 0; i < time; ++i)
	{
		int max = 0;//鱼最多的湖中，鱼的数量
		int id = -1;//鱼最多的湖的编号
		//查找鱼最多的湖中，鱼的数量和湖的编号
		for (j = 0; j < pos; ++j)
			if (fish[j] > max) {
				max = fish[j];
				id = j;
			}
		if (id != -1)//找到了，进行钓鱼处理
		{
			++p[id];
			fish[id] -= d[id];
			t += max;
		}
		//没有找到（从湖1起到湖pos全部钓完了），就把时间放在湖1上
		else ++p[0];
	}
	//处理最优方案
	if (t > best)
	{
		best = t;//最优值
		memset(plan, 0, sizeof(plan));
		for (i = 0; i < pos; ++i)//最优解
			plan[i] = p[i];
	}
}

int main()
{
	int n, h,time=0;
	cin >> n >> h;
	h = h * 60 / 5;
	for (int i = 0; i < n; i++)
		cin >> f[i];
	for (int i = 0; i < n; i++)
		cin >> d[i];
	for (int i = 0; i < n - 1; i++)
		cin >> t[i];
	for (int i = 1; i <= n&& h-time; i++)
	{
		greedy(i, h - time);
		time += t[i];
	}
	for (int i = 0; i < n-1; i++)
		cout << plan[i] * 5 << ",";
	cout << endl;
	cout << plan[n-1] * 5 << endl;
	cout << "最优值为：" << best << endl;
}