第三次课贪心算法精选文档.ppt

《第三次课贪心算法精选文档.ppt》由会员分享，可在线阅读，更多相关《第三次课贪心算法精选文档.ppt（73页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、第三次课贪心算法本讲稿第一页，共七十三页第3章 内容提纲o3.1 贪心算法的基本思想o3.2 删数问题o3.3 装箱问题本讲稿第二页，共七十三页 顾名思义，贪心算法总是作出在当前看来最好的选择。也就是说贪心算法并不从整体最优考虑，它所作出的选择只是在某种意义上的局部最优局部最优选择。当然，希望贪心算法得到的最终结果也是整体最优的。虽然贪心算法不能对所有问题都得到整体最优解，但对许多问题它能产生整体最优解。如单源最短路经问题，最小生成树问题等。在一些情况下，即使贪心算法不能得到整体最优解，其最终结果却是最优解的很好近似。本讲稿第三页，共七十三页活动安排问题 活动安排问题就是要在所给的活动集合中选

2、出最大的相容活动子集合，是可以用贪心算法有效求解的很好例子。该问题要求高效地安排一系列争用某一公共资源的活动。贪心算法提供了一个简单、漂亮的方法使得尽可能多的活动能兼容地使用公共资源。本讲稿第四页，共七十三页活动安排问题 设有n个活动的集合E=1,2,n，其中每个活动都要求使用同一资源，如演讲会场等，而在同一时间内只有一个活动能使用这一资源。每个活动i都有一个要求使用该资源的起始时间si和一个结束时间fi,且si fi。如果选择了活动i，则它在半开时间区间si,fi)内占用资源。若区间si,fi)与区间sj,fj)不相交，则称活动i与活动j是相容的。也就是说，当sifj或sjfi时，活动i与活

3、动j相容。本讲稿第五页，共七十三页活动安排问题templatevoid GreedySelector(int n,Type s,Type f,bool A)A1=true;int j=1;for(int i=2;i=fj)Ai=true;j=i;else Ai=false;下面给出解活动安排问题的贪心算法GreedySelectorGreedySelector:各活动的起始时间和结束时各活动的起始时间和结束时间存储于数组间存储于数组s s和和f f中且按结束中且按结束时间的非减序排列时间的非减序排列 本讲稿第六页，共七十三页活动安排问题 由于输入的活动以其完成时间的非非减序减序排列，所以算法g

4、reedySelectorgreedySelector每次总是选择具有最早完成时间具有最早完成时间的相容活动加入集合A中。直观上，按这种方法选择相容活动为未安排活动留下尽可能多的时间。也就是说，该算法的贪心选择的意义是使剩余的可安排使剩余的可安排时间段极大化时间段极大化，以便安排尽可能多的相容活动。本讲稿第七页，共七十三页 算法greedySelectorgreedySelector的效率极高。当输入的活动已按结束时间的非减序排列，算法只需O(n)O(n)的时间安排n个活动，使最多的活动能相容地使用公共资源。如果所给出的活动未按非减序排列，可以用O(nlogn)O(nlogn)的时间重排。本讲

5、稿第八页，共七十三页活动安排问题 例：例：设待安排的11个活动的开始时间和结束时间按结束时间的非减序排列如下：i1234567891011Si130535688212fi4567891011121314本讲稿第九页，共七十三页活动安排问题 算法算法greedySelector greedySelector 的计算过程的计算过程如左图所示。图中每行相应于算法的一次迭代。阴影长条表示的活动是已选入集合A的活动，而空白长条表示的活动是当前正在检查相容性的活动。本讲稿第十页，共七十三页活动安排问题 若被检查的活动i的开始时间Si小于最近选择的活动j的结束时间fi，则不选择活动i，否则选择活动i加入集合

6、A中。贪心算法并不总能求得问题的整体最优解整体最优解。但对于活动安排问题，贪心算法greedySelector却总能求得的整体最优解，即它最终所确定的相容活动集合A的规模最大。这个结论可以用数学归纳法证明。本讲稿第十一页，共七十三页贪心算法的基本要素 本节着重讨论可以用贪心算法求解的问题的一般特征。对于一个具体的问题，怎么知道是否可用贪心算法解此问题，以及能否得到问题的最优解呢?这个问题很难给予肯定的回答。但是，从许多可以用贪心算法求解的问题中看到这类问题一般具有2个重要的性质：贪心选择性质贪心选择性质和最优子结构性质最优子结构性质。本讲稿第十二页，共七十三页贪心算法的基本要素1 1、贪心选择

7、性质、贪心选择性质 所谓贪心选择性质贪心选择性质是指所求问题的整体最优解整体最优解可以通过一系列局部最优局部最优的选择，即贪心选择来达到。这是贪心算法可行的第一个基本要素，也是贪心算法与动态规划算法的主要区别。动态规划算法通常以自底向上自底向上的方式解各子问题，而贪心算法则通常以自顶向下自顶向下的方式进行，以迭代的方式作出相继的贪心选择，每作一次贪心选择就将所求问题简化为规模更小的子问题。对于一个具体问题，要确定它是否具有贪心选择性质，必须证明每一步所作的贪心选择最终导致问题的整体最优解。本讲稿第十三页，共七十三页贪心算法的基本要素 当一个问题的最优解包含其子问题的最优解时，称此问题具有最优子

8、结构性质最优子结构性质。问题的最优子结构性质是该问题可用动态规划算法或贪心算法求解的关键特征。2 2、最优子结构性质、最优子结构性质本讲稿第十四页，共七十三页贪心算法的基本要素 贪心算法和动态规划算法都要求问题具有最优子结构性质，这是2类算法的一个共同点。但是，对于具有最优子结构最优子结构的问题应该选用贪心算法还是动态规划算法求解?是否能用动态规划算法求解的问题也能用贪心算法求解?下面研究2个经典的组合优化问题组合优化问题，并以此说明贪心算法与动态规划算法的主要差别。3、贪心算法与动态规划算法的差异本讲稿第十五页，共七十三页贪心算法的基本要素o0-10-1背包问题：背包问题：给定n种物品和一个

9、背包。物品i的重量是Wi，其价值为Vi，背包的容量为C。应如何选择装入背包的物品，使得装入背包中物品的总价值最大?在选择装入背包的物品时，对每种物品在选择装入背包的物品时，对每种物品i i只有只有2 2种种选择，即装入背包或不装入背包。不能将物品选择，即装入背包或不装入背包。不能将物品i i装入背包多次，也不能只装入部分的物品装入背包多次，也不能只装入部分的物品i i。本讲稿第十六页，共七十三页贪心算法的基本要素o背包问题：背包问题：与0-1背包问题类似，所不同的是在选择物品i装入背包时，可以选择物品可以选择物品i i的一的一部分部分，而不一定要全部装入背包，1in。这2类问题都具有最优子结构

10、最优子结构性质，极为相似，但背包问题可以用贪心算法求解，而0-1背包问题却不能用贪心算法求解。本讲稿第十七页，共七十三页贪心算法的基本要素 首先计算每种物品单位重量的价值Vi/Wi，然后，依贪心选择策略，将尽可能多的单位重量价值最高单位重量价值最高的物品装入背包。若将这种物品全部装入背包后，背包内的物品总重量未超过C，则选择单位重量价值次高的物品并尽可能多地装入背包。依此策略一直地进行下去，直到背包装满为止。具体算法可描述如下页：用贪心算法解背包问题的基本步骤：本讲稿第十八页，共七十三页贪心算法的基本要素void Knapsack(int n,float M,float v,float w,f

11、loat x)Sort(n,v,w);int i;for(i=1;i=n;i+)xi=0;float c=M;for(i=1;ic)break;xi=1;c-=wi;if(i=n)xi=c/wi;算法算法knapsackknapsack的主要的主要计算时间在于将计算时间在于将各种物品依其单各种物品依其单位重量的价值从位重量的价值从大到小排序。因大到小排序。因此，算法的计算此，算法的计算时间上界为时间上界为O O（nlognnlogn）。）。为了证明算法的为了证明算法的正确性，还必须正确性，还必须证明背包问题具证明背包问题具有贪心选择性质有贪心选择性质。本讲稿第十九页，共七十三页贪心算法的基本要

12、素 对于0-10-1背包问题背包问题，贪心选择之所以不能得到最优解是因为在这种情况下，它无法保证最终能将背包装满，部分闲置的背包空间使每公斤背包空间的价值降低了。事实上，在考虑0-1背包问题时，应比较选择该物品和不选择该物品所导致的最终方案，然后再作出最好选择。由此就导出许多互相重叠的子问题。这正是该问题可用动态规划算法动态规划算法求解的另一重要特征。实际上也是如此，动态规划算法的确可以有效地解0-1背包问题。本讲稿第二十页，共七十三页3.2 删数问题给定一个高精度的正整数N（不超过240位），去掉任意S个数字后剩下的数字按原左右次序组成一个新的正整数。编程对于给定的N和S，寻找一种方案使得剩

13、下的数字组成的新数最小示例：N=178543，S=4则结果为：13本讲稿第二十一页，共七十三页算法分析oN的存储只能采用数组，可选择字符数组；o可采用贪心算法，即：共删掉S个数，每一步删掉一个数，并且总选择一个使剩下的数最小的数符删去。o为了保证删1个数符后的数最小，可以按照从高位到低位的方向收缩递减区间，如果不存在递减区间则删尾数符；否则删掉递减区间的首字符，这样形成一个新的数串。然后回到串首，重复上面的规则，删下一个数符，。，直至删除S个数符为止。本讲稿第二十二页，共七十三页17854317543154314313本讲稿第二十三页，共七十三页17853417534153413413本讲稿第

14、二十四页，共七十三页 程序#include#include using namespace std;int s;class nodepublic:char c;node*next;class queepublic:node*first;本讲稿第二十五页，共七十三页quee*init()char str250;quee*n;node*tmp;node*tail;n=new quee;n-first=NULL;cinstr;cins;int i=0;while(stri!=0)tmp=new node;tmp-c=stri;tmp-next=NULL;if(n-first=NULL)n-first

15、=tmp;tail=tmp;elsetail-next=tmp;tail=tmp;i+;return n;本讲稿第二十六页，共七十三页void del(quee*start,node*p)node*x;x=start-first;if(x-next)while(x-next!=p)x=x-next;x-next=p-next;delete p;本讲稿第二十七页，共七十三页void print(quee*start)node*tmp;if(start-first=NULL)return;tmp=start-first;while(tmp)printf(%c,tmp-c);tmp=tmp-next

16、;printf(n);本讲稿第二十八页，共七十三页node*find(quee*start)node*tmp;tmp=start-first;while(tmp-next!=NULL&tmp-cnext-c)tmp=tmp-next;return tmp;本讲稿第二十九页，共七十三页int main(int argc,char*argv)quee*n;node*p;n=init();while(s!=0)p=find(n);del(n,p);s-;print(n);return 0;本讲稿第三十页，共七十三页3.3 POJ 1017Packets本讲稿第三十一页，共七十三页Packets题意已

17、知：有6*6 的大箱子和 1*1，2*2，3*3，4*4，5*5，6*6 的木块(平面)问：给定各种木块的数目，求最少需要多少个大箱子来装？例如：输入：0 0 4 0 0 1 -输出 2输入：7 5 1 0 0 0 -输出 1本讲稿第三十二页，共七十三页Packetso解题思想：o贪心准则：先放大的，后放小的1.6*6的木块每个占用一个箱子；2.5*5的木块每个占用一个新箱子，余下11个1*1的空格；3.4*4的木块每个占用一个新箱子，余下5个2*2的空格。4.3*3的木块每4个占用新一个箱子，不足4个也占一个新箱子，分情况余下不同数目的空格；5.2*2的木块先填空格，空格不足开新箱子，每9个

18、2*2的木块占一个新箱子；6.1*1的木块先填空格，空格不足开新箱子，每36个占一个新箱子。本讲稿第三十三页，共七十三页Packetso假设 6*6，5*5，4*4，3*3，2*2，1*1的箱子个数 b6 b5 b4 b3 b2 b1本讲稿第三十四页，共七十三页Packets 4*4，5*5，6*6的块单独开新的箱子。每一个5*5的块还能放下11个1*1的箱子。每一个4*4的块还能放下5个2*2的箱子，或者放下20个1*1的箱子。4*4块5*5块6*6块本讲稿第三十五页，共七十三页 3*3的块1-4块占一个新箱子，一个箱子可放下4个3*3的块。如果一个箱子放下了1个3*3的块，则还可放下5个2

19、*2的块和7个1*1的块。或者可以放下27个1*1的块。如果一个箱子放下了2个3*3的块，则还可放下3个2*2的块和6个1*1的块。或者可以放下18个1*1的块。如果一个箱子放下了3个3*3的块，则还可放下1个2*2的块和5个1*1的块。或者可以放下9个1*1的块。本讲稿第三十六页，共七十三页PacketsnTotal -箱子数1）先放好所有 6*6,5*5,4*4 和3*3 的木块nTotal=b6+b5+b4+(b3+3)/4 4*4,5*5,6*6 单独开新的箱子 3*3 每 1-4 个占一个新箱子 2）再把2*2的塞到放有3*3木块的箱子里 int Contain24=0,5,3,1;

20、Contain2i 表示当箱子里有i个 3*3的木块时，还能放下多少个2*2的木块3）考虑放有3*3的木块的箱子在塞满2*2的木后还能放多少个1*1的木块：int Contain14=0,7,6,5;Contain1i 表示当箱子里有i个 3*3的木块，并且还填满了2*2的木块后，还能放下多少个1*1的木块。本讲稿第三十七页，共七十三页Packets 如果箱子里放1个3*3木块，那么还能放5个2*2木块,以及7个1*1木块本讲稿第三十八页，共七十三页Packets 如果箱子里放2个3*3木块，那么还能放3个2*2木块,以及6个1*1木块本讲稿第三十九页，共七十三页Packets 如果箱子里放3

21、个3*3木块，那么还能放1个2*2木块,以及5个1*1木块本讲稿第四十页，共七十三页Packets4)计算放好6*6,5*5,4*4,3*3后留下多少空格能放2*2c2=b4*5+Contain2b3%4;留下多少空格能放1*1 c1=b5*11+Contain1b3%4(能放2*2的地方暂不考虑放1*1，因为如果c2b2则多余的b2-c2个位置可以放4*(b2-c2)个1*1的箱子)本讲稿第四十一页，共七十三页Packets5)放 2*2，如果不够放，就开新箱子，放完2*2，计算还剩多少空格能放1*1。if(b2 c1)onTotal+=(b1-c1)/36;oif(b1-c1)%36)on

22、Total+;oocout nTotal endl;Packets本讲稿第四十三页，共七十三页o#include oint Contain14=0,7,6,5;oint Contain24=0,5,3,1;oint nTotal;ovoid main()o int b1,b2,b3,b4,b5,b6;o for(;)o cinb1b2b3b4b5b6;o if(b1=0&b2=0&b3=0&b4=0&b5=0&b6=0)obreak;onTotal=b6+b5+b4+b3/4;oif(b3%4)nTotal+;oint c1;/当前能放 1*1 木块的空格数目oc1=b5*11;/每个放5*5

23、的箱子，还能放11个 1*1 的木块c1+=Contain1b3%4;/加上3*3箱子里能放的数目oint c2;/当前能放 2*2 木块的空格数目（先不考虑将1*1的木块放在 /2*2的空格里）oc2=b4*5;/每个放5*4的箱子，还能放5个 2*2 的木块oc2+=Contain2b3%4;/加上3*3箱子里能放的数目本讲稿第四十四页，共七十三页oif(b2 c1)onTotal+=(b1-c1)/36;oif(b1-c1)%36)onTotal+;oocout nTotal endl;o 本讲稿第四十五页，共七十三页案例3水位 o为使房屋买主可以评估洪水保险的开销，一家真实的房地产公司

24、提供了一种客户端程序，此程序现实了可能被购买得房屋所在地区以100平方米为单位的海拔高度。雨水、雪水和管道破裂流出的水将会汇聚到海拔最低，因为水会从高处流向低处。为简化问题，我们假设水沿着排水沟从高出流向低处，并且水不会被土地所吸收。本讲稿第四十六页，共七十三页o从天气信息的档案中，我们得知了一个区域的典型储水量。作为预期的购房者，我们希望知道低处积水后的水位，以及完全被淹没在水中的面积占此地区的百分比(也就是，10平方米的海拔高度严格低于水平面)。你被要求写一个程序来给出结果。本讲稿第四十七页，共七十三页输入数据o输入包括区域描述所构成的一个序列。每个由一对整数m和n开始，每个都小于30，给

25、定的区域都是100平方米。紧随其后的是m行，每行n个整数，给出地区海平面，以主行序给出。高度的单位是米，分别用正数或者负数表明高于或低于海平面。每个区域最后一个值是一个整数，用以指出有多少立方米的水聚集到此区域。m和n为两个0代表输入结束。本讲稿第四十八页，共七十三页输出描述o对每个区域，显示区域编号(1,2,)，水平面(以米为单位表示高于或低于海平面)以及此区域被水淹没的面积，每个占一行。水位和被水淹没面积的百分率显示时保留2位小数每个区域输出后跟一个空行。本讲稿第四十九页，共七十三页样例o3 3o25 37 45o51 12 34o94 83 27o10000o0 0 oRegion 1o

26、Water level is 46.67 meters.o66.67 percent of the region is under water.本讲稿第五十页，共七十三页解题思路o由于题目中特别声明水无论如何都会流到当前水面最低的地方，使得问题一下子简化了。很容易想到以下贪心算法：o（1）把每个格子的高度排序；o（2）以低格子到高格子的顺序填水，把水均匀的铺在当前的水面上，并不断更新当前水面面积，和剩余水量；o（3）若剩余水量为0，输出当前水面高度和水面覆盖率。本讲稿第五十一页，共七十三页注意o1水面每达到一个格子的高度，水面的面积也要随之扩大；o2高度有负值，水面的初始高度是最低格子的高度而

27、不是0；o3考虑好若干个格子高度相等的情况；o4每个格子的面积是100；o5没有水时，覆盖率为0；o6水若刚好达到一个格子的高度，此格子不算被水覆盖；o7水面达到最高格子后，要把剩下的水均匀的铺在整个区域上；o8注意浮点数的计算误差，因为这个原因大家在POJ上都没有通过此题:-本讲稿第五十二页，共七十三页#include#include int h1000000;/保存每个格子的高度值，因不知m和n的范围，故设的比较大int cmp(const void*a,const void*b)return*(int*)a)-*(int*)b);本讲稿第五十三页，共七十三页main()int t,m,n

28、,nn;double x,hw;int i;t=0;本讲稿第五十四页，共七十三页 while(1)/输入 scanf(%d%d,&n,&m);if(n=0&m=0)break;t+;nn=m*n;for(i=0;i 0)for(i=1;i nn;i+)/在高度相等的情况下，扩大水面面积 while(i nn&hi=hi-1)i+;if(i=nn|x 0)/若剩下水，就把它铺在当前的范围上 while(i 1且B/A=BA，则C=B/Ao若A1，则C=B,打印1/Co转步骤（2）。本讲稿第六十一页，共七十三页#include#include int num;int a,b;int c;int z

29、c(int a,int b)int xx;xx=b/a;if(xx*a=b)return 1;return 0;本讲稿第六十二页，共七十三页void calab()while(1)c=(b/a)+1;a=a*c-b;b=b*c;printf(%d,c);if(a1&zc(a,b)printf(%dn,b/a);break;else if(a=1)printf(%dn,b);break;本讲稿第六十三页，共七十三页omain()oo int i;o scanf(%dn,&num);o for(i=1;i=num;i+)o o scanf(%d%dn,&a,&b);o calab();o o本讲稿

30、第六十四页，共七十三页最佳浏览路线问题-5017-42-47-19-3-36-34-43-30-4334-23-8-45南北西东本讲稿第六十五页，共七十三页课堂练习 校门外的树o某校大门外长度为L的马路上有一排树，每两棵相邻的树之间的间隔都是1米。我们可以把马路看成一个数轴，马路的一端在数轴0的位置，另一端在L的位置；数轴上的每个整数点，即0，1，2，L，都种有一棵树。o由于马路上有一些区域要用来建地铁。这些区域用它们在数轴上的起始点和终止点表示。已知任一区域的起始点和终止点的坐标都是整数，区域之间可能有重合的部分。现在要把这些区域中的树（包括区域端点处的两棵树）移走。你的任务是计算将这些树都

31、移走后，马路上还有多少棵树。本讲稿第六十六页，共七十三页课堂练习 校门外的树o【输入文件】o输入的第一行有两个整数L（1=L=10000）和 M（1=M=100），L代表马路的长度，M代表区域的数目，L和M之间用一个空格隔开。接下来的M行每行包含两个不同的整数，用一个空格隔开，表示一个区域的起始点和终止点的坐标。o【输出文件】o输出包括一行，这一行只包含一个整数，表示马路上剩余的树的数目。o本讲稿第六十七页，共七十三页课堂练习 校门外的树【样例输入】o500 3o150 300o100 200o470 471o【样例输出】o298本讲稿第六十八页，共七十三页#include#include i

32、nt a10010;int main()int l,n,i,j,s,t,ans;while(scanf(%d%d,&l,&n)!=EOF)memset(a,0,sizeof(a);for(i=0;in;i+)scanf(%d%d,&s,&t);for(j=s;j=t;j+)aj=1;for(ans=0,i=0;i=l;i+)if(!ai)ans+;printf(%dn,ans);return 0;本讲稿第六十九页，共七十三页方法二：区间合并o按区间的开始值对所有区间排序；o扫描合并本讲稿第七十页，共七十三页本讲稿第七十一页，共七十三页程序本讲稿第七十二页，共七十三页方法三：排序扫描本讲稿第七十三页，共七十三页