贪心算法与动态规划----基于HDOJ2059龟兔赛跑问题分析。

题目

据说在很久很久以前，可怜的兔子经历了人生中最大的打击——赛跑输给乌龟后，心中郁闷，发誓要报仇雪恨，于是躲进了杭州下沙某农业园卧薪尝胆潜心修炼，终于练成了绝技，能够毫不休息得以恒定的速度(VR m/s)一直跑。兔子一直想找机会好好得教训一下乌龟，以雪前耻。
最近正值HDU举办50周年校庆，社会各大名流齐聚下沙，兔子也趁此机会向乌龟发起挑战。虽然乌龟深知获胜希望不大，不过迫于舆论压力，只能接受挑战。
比赛是设在一条笔直的道路上，长度为L米，规则很简单，谁先到达终点谁就算获胜。
无奈乌龟自从上次获胜以后，成了名龟，被一些八卦杂志称为“动物界的刘翔”，广告不断，手头也有了不少积蓄。为了能够再赢兔子，乌龟不惜花下血本买了最先进的武器——“"小飞鸽"牌电动车。这辆车在有电的情况下能够以VT1 m/s的速度“飞驰”，可惜电池容量有限，每次充满电最多只能行驶C米的距离，以后就只能用脚来蹬了，乌龟用脚蹬时的速度为VT2 m/s。更过分的是，乌龟竟然在跑道上修建了很多很多（N个)的供电站，供自己给电动车充电。其中，每次充电需要花费T秒钟的时间。当然，乌龟经过一个充电站的时候可以选择去或不去充电。
比赛马上开始了，兔子和带着充满电的电动车的乌龟并列站在起跑线上。你的任务就是写个程序，判断乌龟用最佳的方案进军时，能不能赢了一直以恒定速度奔跑的兔子。

Input
本题目包含多组测试，请处理到文件结束。每个测试包括四行：
第一行是一个整数L代表跑道的总长度
第二行包含三个整数N，C，T，分别表示充电站的个数，电动车冲满电以后能行驶的距离以及每次充电所需要的时间
第三行也是三个整数VR，VT1，VT2，分别表示兔子跑步的速度，乌龟开电动车的速度，乌龟脚蹬电动车的速度
第四行包含了N(N<=100)个整数p1,p2…pn,分别表示各个充电站离跑道起点的距离，其中0<p1<p2<…<pn<L
其中每个数都在32位整型范围之内。

Output
当乌龟有可能赢的时候输出一行 “What a pity rabbit!"。否则输出一行"Good job,rabbit!";
题目数据保证不会出现乌龟和兔子同时到达的情况。

注：设起点为第0个充电站，终点为第N+1个充电站。

动态规划与贪心算法

贪心算法

1.将问题Q(n)（n是问题的规模），分解成问题Q（n-1）和递推问题W(0）。
2.Q（n）的最优解必然是建立在Q（n-1）的最优解上。
3.解出Q（1）的最优解，递推得到Q（n）的最优解。

本题分析：

Q（n）:到达第n个供电站的最短时间。（将起点设置为第0个供电站，终点设置为最后一个供电站）
W(0）：从上一个供电站到下一个供电站所需的最短时间。

这种分解是错误的。因为问题Q（n-1）的最优解中并不一定包含Q（n）的最优解，即到达供电站（n-1）花费的时间最短并不意味着到达供电站(n)所花费的时间最短。因此用贪心算法难以解决。

动态规划

1.对于原问题Q（n）,列出其各阶子问题Q（1），Q（2）…Q(n-1)以及递推问题W(0);
2.Q（n）的最优解必然是建立在其各阶子问题中的某一个的最优解上。
3.解出Q（1）的最优解，递推得到Q（n）的最优解。

本题分析

Q（n）:到达第n个供电站的最短时间。（将起点设置为第0个供电站，终点设置为最后一个供电站）
W(0):从第i个充电站充满电之后，一直行驶到第n个充电站，中间不再进站充电所需的时间。（0<=i<n）

for(i=0;i<n;i++){
	Q(n)=min(Q(n),Q(i)+W(0))
}

到达站点n的所花时间最短，必然是先到达某个站点i进行充电，然后不再充电，行驶到站点n。（注意起点也是站点0）

总结

在贪心算法中，原问题Q（n）的最优解，必然建立在其上一阶子问题Q（n-1）的最优解上。
而在动态规划中，原问题Q（n）的最优解，不一定建立在其上一阶子问题上，而是建立在其各阶子问题Q（1）~Q（n-1）的某一个子问题的最优解上。故，需要数组来存储其各阶子问题的最优解，然后遍历子问题最优解，寻求原问题最优解。
可以说，贪心算法，是动态规划的一种特例，只需要遍历子问题Q（n-1）.

代码

#include<stdio.h>
#include <math.h>
#define MAX 0xfffff  

//NT函数：返回从pa到pb不充电所花的时间。
double NT(double i,double pa,double pb,double c,double V1,double V2,double T){
	double t2;
	if(pb-pa<c){
		t2=(pb-pa)/V1+T;
	}else{
		t2=T+c/V1+(pb-pa-c)/V2;
	}
	if(!i){
		return t2-T;
	}
	return t2;
	
}
int main()
{
	int i,j,N;
	double L,C,T,VR,VT1,VT2;
	double p[102],dp[102];//p[]记录每个充电站的位置，dp记录到达该充电站的最少时间

	while(scanf("%lf",&L)!=EOF){
		scanf("%d %lf %lf",&N,&C,&T);
		scanf("%lf %lf %lf",&VR,&VT1,&VT2);
		for(i=1;i<=N;i++){
			scanf("%lf",&p[i]);
			dp[i]=MAX;
		}
		p[0]=0;p[i]=L;
		dp[0]=0;dp[i]=MAX;
		for(i=1;i<=N+1;i++){
			
			for(j=0;j<i;j++){
				if(dp[j]+NT(j,p[j],p[i],C,VT1,VT2,T)<dp[i]){
					dp[i]=dp[j]+NT(j,p[j],p[i],C,VT1,VT2,T);
				}
			}
		}
		if(L/VR<dp[N+1]){
			printf("Good job,rabbit!\n");
		}else{
			printf("What a pity rabbit!\n");
		}

	}

	getchar();
	getchar();
	return 0;
}