分治算法(金块问题)

分治算法(金块问题)
分治算法(金块问题)

问题描述:

有一个老板有一袋金块。每个月将有两名雇员会因其优异的表现分别被奖励一个金块。按规矩,排名第一的雇员将得到袋中最重的金块,排名第二的雇员将得到袋中最轻的金块。根据这种方式,除非有新的金块加入袋中,否则第一名雇员所得到的金块总是比第二名雇员所得到的金块重。如果有新的金块周期性的加入袋中,则每个月都必须找出最轻和最重的金块。假设有一台比较重量的仪器,我们希望用最少的比较次数找出最轻和最重的金块

算法思想:

分而治之方法与软件设计的模块化方法非常相似。为了解决一个大的问题,可以:1) 把它分成两个或多个更小的问题; 2) 分别解决每个小问题; 3) 把各小问题的解答组合起来,即可得到原问题的解答。小问题通常与原问题相似,可以递归地使用分而治之策略来解决。

问题分析:

一般思路:假设袋中有n 个金块。可以用函数M a x(程序1 - 3 1)通过n-1次比较找到最重的金块。找到最重的金块后,可以从余下的n-1个金块中用类似的方法通过n-2次比较找出最轻的金块。这样,比较的总次数为2n-3。

分治法:当n很小时,比如说,n≤2,识别出最重和最轻的金块,一次比较就足够了。当n 较大时(n>2),第一步,把这袋金块平分成两个小袋A和B。第二步,分别找出在A和B中最重和最轻的金块。设A中最重和最轻的金块分别为HA 与LA,以此类推,B中最重和最轻的金块分别为HB 和LB。第三步,通过比较HA 和HB,可以找到所有金块中最重的;通过比较LA 和LB,可以找到所有金块中最轻的。在第二步中,若n>2,则递归地应用分而治之方法

程序设计

据上述步骤,可以得出程序1 4 - 1的非递归代码。该程序用于寻找到数组w [ 0 : n - 1 ]中的最小数和最大数,若n < 1,则程序返回f a l s e,否则返回t r u e。

当n≥1时,程序1 4 - 1给M i n和M a x置初值以使w [ M i n ]是最小的重量,w [ M a x ]为最大的重量。

首先处理n≤1的情况。若n>1且为奇数,第一个重量w [ 0 ]将成为最小值和最大值的候选值,因此将有偶数个重量值w [ 1 : n - 1 ]参与f o r循环。当n 是偶数时,首先将两个重量值放在for 循环外进行比较,较小和较大的重量值分别置为Min和Max,因此也有偶数个重量值w[2:n-1]参与for循环。

在for 循环中,外层if 通过比较确定( w [ i ] , w [ i + 1 ] )中的较大和较小者。此工作与前面提到的分而治之算法步骤中的2) 相对应,而内层的i

f负责找出较小重量值和较大重量值中的最小值和最大值,这个工作对应于3 )。for 循环将每一对重量值中较小值和较大值分别与当前的最小值w [ M i n ]和最大值w [ M a x ]进行比较,根据比较结果来修改M i n和M a x(如果必要)。

复杂性分析

注意到当n为偶数时,在for 循环外部将执行一次比较而在f o r循环内部执行3 ( n / 2 - 1 )次比较,比较的总次数为3 n / 2 - 2。当n 为奇数时,f o r循环外部没有执行比较,而内部执行了3(n-1)/2次比较。因此无论n 为奇数或偶数,当n>0时,比较的总次数为「3n/2ù-2次。

关键步骤:

程序14-1 找出最小值和最大值的非递归程序

template

bool MinMax(T w[], int n, T& Min, T& Max)

{// 寻找w [ 0 : n - 1 ]中的最小和最大值

// 如果少于一个元素,则返回f a l s e

// 特殊情形: n <= 1

if (n < 1) return false;

if (n == 1) {Min = Max = 0;

return true;}

/ /对Min 和M a x进行初始化

int s; // 循环起点

if (n % 2) {// n 为奇数

Min = Max = 0;

s = 1;}

else {// n为偶数,比较第一对

if (w[0] > w[1]) {

Min = 1;

Max = 0;}

else {Min = 0;

Max = 1;}

s = 2;}

// 比较余下的数对

for (int i = s; i < n; i += 2) {

// 寻找w[i] 和w [ i + 1 ]中的较大者 // 然后将较大者与w [ M a x ]进行比较 // 将较小者与w [ M i n ]进行比较

if (w[i] > w[i+1]) {

if (w[i] > w[Max]) Max = i;

if (w[i+1] < w[Min]) Min = i + 1;}

else {

if (w[i+1] > w[Max]) Max = i + 1;

if (w[i] < w[Min]) Min = i;}

}

return true; }

算法分析与设计习题集整理

算法分析与设计习题集整理 第一章算法引论 一、填空题: 1、算法运行所需要的计算机资源的量,称为算法复杂性,主要包括时间复杂度和空间复杂度。 2、多项式10()m m A n a n a n a =+++L 的上界为O(n m )。 3、算法的基本特征:输入、输出、确定性、有限性 、可行性 。 4、如何从两个方面评价一个算法的优劣:时间复杂度、空间复杂度。 5、计算下面算法的时间复杂度记为: O(n 3) 。 for(i=1;i<=n;i++) for(j=1;j<=n;j++) {c[i][j]=0; for(k=1;k<=n;k++) c[i][j]= c[i][j]+a[i][k]*b[k][j]; } 6、描述算法常用的方法:自然语言、伪代码、程序设计语言、流程图、盒图、PAD 图。 7、算法设计的基本要求:正确性 和 可读性。 8、计算下面算法的时间复杂度记为: O(n 2) 。 for (i =1;i

分治算法例题

目录 1031 输油管道问题 (2) 解题思路 (2) 程序代码 (2) 1032 邮局选址 (5) 解题思路 (5) 程序源代码 (5) 1034 集合划分2 (7) 解题思路: (7) 程序源代码: (7) 1033 集合划分 (9) 解题思路 (9) 程序源代码 (9)

1031 输油管道问题 解题思路 本题目可以分为两个步骤: 1、找出主管道的位置; 2、根据主管道的位置,计算各个油井到主管道的长度之和。 根据题意,设主管道贯穿东西,与y 轴平行。而各个子油井则分布在主输油管道的上下两侧。如下图: 由上图,其实只需要确定主管道的y 坐标,而与各个子油井的x 坐标无关!根据猜测,易知:主管道的y 坐标就是所有子油井y 坐标的中位数。(可以用平面几何知识证明,略) 求中位数的方法可以用排序后取a[(left+right)/2],当然更推荐用书上的线性时间选择算法解决。记求得的主管道为m y ,最后要输出的结果只需要计算:1||n i m i y y =-∑,输出即可。 另外要提醒的是本题多Case 。 程序代码 #include #include void swap (int &a ,int &b ) { int tmp = a ; a = b ; b = tmp ; }

//本函数求arr[p:q]的一个划分i,使arr[p:i-1]都小于arr[i],arr[i+1,q]都大于arr[i] int partition(int *arr,int p,int q) { int index = p-1,start = p,base = arr[q]; for(;start

大学算法分析与设计复习总结

大学算法分析与设计复习总结 第1章绪论 考点: 1、算法的5个重要特性。(P3) 答:输入、输出、有穷性、确定性、可行性 2、描述算法的四种方法分别是什么,有什么优缺点。(P4) 答: 1.自然语言优点:容易理解;缺点:容易出现二义性,并且算法都很冗长。 2.流程图优点:直观易懂;缺点:严密性不如程序语言,灵活性不如自然语言。 3.程序设计语言优点:用程序语言描述的算法能由计算机直接执行;缺点:抽象性差,是算法设计者拘泥于描述算法的具体细节,忽略了“好”算法和正确逻辑的重要性,此外,还要求算法设计者掌握程序设计语言及其编程技巧。 4.伪代码优点:表达能力强,抽象性强,容易理解 3、了解非递归算法的时间复杂性分析。(P13) 要点:对非递归算法时间复杂性的分析,关键是建立一个代表算法运行时间的求和表达式,然后用渐进符号表示这个求和表达式。 非递归算法分析的一般步骤是: (1)决定用哪个(或哪些)参数作为算法问题规模的度量。 (2)找出算法的基本语句。 (3)检查基本语句的执行次数是否只依赖问题规模。 (4)建立基本语句执行次数的求和表达式。 (5)用渐进符号表示这个求和表达式。

[例1.4]:求数组最小值算法 int ArrayMin(int a[ ], int n) { min=a[0]; for (i=1; i

通用分支递归式: 使用扩展递归技术求解下列递推关系式(1) (2)

分治法--凸包

//纯代码 #include #define PPmax 30 typedef struct node { float x,y; }Point; Point DingDian[PPmax];//用于存放凸边形的顶点 int DingDnum=0; typedef struct Pointss { Point p1,p2; }SDian; SDian BDD[PPmax];//存放凸边形边的双点 int BDDnum=0;//记录已经存入的双点数 int abc=0;//全局变量,用来记录划边顺序 void dianxu(Point P[],int n) //对已经输入的点数组进行排序 {//存放点的数组P[], 点的个数n // ——点排序 // 具体实现:(选择排序实现) // 1.按x坐标从小到大排 // 2.扫描最前端点元素 // 2.1若不存在相同x坐标的点,则执行3 // 2.2若存在相同x坐标的点,按|y|排列大->小 // 3.扫描最后端点元素 // 3.1若不存在相同x坐标的点,则执行4 // 3.2若存在相同x坐标的点,按|y|排列小->大 // 4.返回。 Point temp; int i,j; int min,max; float px; int pxnum; //按x坐标排序 for(i=0;i<=n-2;i++) { min=i; for(j=i+1;j<=n-1;j++) { if(P[j].x

{ min=j; } } temp=P[i]; P[i]=P[min]; P[min]=temp; } //数组两端相同x按|y|排序 //前端排序 px=P[0].x; pxnum=1; while(P[pxnum].x==px) { pxnum++; } for(i=0;i<=pxnum-2;i++) { //max=i; min=i; for(j=i+1;j<=pxnum-1;j++) { if( P[j].y < P[min].y ) { min=j; } } temp=P[i]; P[i]=P[min]; P[min]=temp; } //后端排序 px=P[n-1].x; pxnum=1; while(P[n-1-pxnum].x==px) { pxnum++; } for(i=n-pxnum;i<=n-2;i++) { min=i;

分治法解决凸包问题

#include typedef struct { int x,y; }Point; int arrayLength=0; int main() { printf("**************软件工程3班**************"); printf("*****************成冠辉*****************"); Point array[15],result[15]; InitialData(array); printf("请输入数据:\n\n"); Put(array,arrayLength); Qsort(array,0,arrayLength-1); printf("排序后:\n"); Put(array,arrayLength); GetResult(array,result); printf("结果集为:\n"); Put(result,resultCount); getchar(); return 0; } //返回分裂位置 int Partition(Point *array,int l,int r) { int p=l,i=l,j=r+1; do { do { i++; }while(array[i].x

j--; }while(array[j].x>array[p].x); swap(array,i,j); }while(iarray[14]射线左边的点DealWithLeft(14,0,array,result); //处理array[14]->array[0]射线右边的点} void InitialData(Point *array) { FILE *fp=freopen("input.txt","r",stdin); char ch;

算法分析——分治法

分治算法 小组的汇报内容: 一、分治算法的基本概念 (2) 二、分治算法的基本思想及策略 (2) 三、分治法适用的情况 (3) 四、分治法的基本步骤 (3) 五、分治法的复杂性分析 (4) 六、快速傅里叶变换 (5) 七、可使用分治法求解的一些经典问题 (9) 八、依据分治法设计程序时的思维过程 (9) 九、分治法与其它常见算法的比较 (9) 小组的分工情况: 彭勇讲前五个部分,天西山讲第六个部分,胡化腾讲最后三个部分,吕璐负责汇报的资料收集,小组分工以及最后的资料整理。

三、分治法适用的情况 分治法所能解决的问题一般具有以下几个特征: 1) 该问题的规模缩小到一定的程度就可以容易地解决 2) 该问题可以分解为若干个规模较小的相同问题,即该问题具有最优子结构性质。 3) 利用该问题分解出的子问题的解可以合并为该问题的解; 4) 该问题所分解出的各个子问题是相互独立的,即子问题之间不包含公共的子子问题。 第一条特征是绝大多数问题都可以满足的,因为问题的计算复杂性一般是随着问题规模的增加而增加; 第二条特征是应用分治法的前提它也是大多数问题可以满足的,此特征反映了递归思想的应用;、 第三条特征是关键,能否利用分治法完全取决于问题是否具有第三条特征,如果具备了第一条和第二条特征,而不具备第三条特征,则可以考虑用贪心法或动态规划法。 第四条特征涉及到分治法的效率,如果各子问题是不独立的则分治法要做许多不必要的工作,重复地解公共的子问题,此时虽然可用分治法,但一般用动态规划法较好。 四、分治法的基本步骤 分治法在每一层递归上都有三个步骤: step1 分解:将原问题分解为若干个规模较小,相互独立,与原问题形式相同的子问题; step2 解决:若子问题规模较小而容易被解决则直接解,否则递归地解各个子问题

凸包问题

凸包问题 摘要:凸包问题是计算机几何中的一个经典问题,它要求将平面内的点集用最少的凸点将所有的顶点封闭。凸包问题的应用十分广泛,很多最优化问题经过抽象后可以发现它最终是凸包问题模型;它还可以用于人际关系网络的最小化搜索,通过人际关系,可以合理推断出某人的身份,职位等个人特征。目前求取凸包的常用算法有:穷举法,格雷厄姆扫描法(Graham),分治法,蛮力法和Jarris 步进法。其中穷举法与蛮力法都是建立在穷举的算法思想上,它们的时间复杂度比较大,格雷厄姆扫描法采用几何方面的知识,降低了求解过程的时间复杂度。关键词:凸包问题;计算机几何;格雷厄姆扫描法 一、引言 凸包问题的完整描述:令S 是平面上的一个点集,封闭S 中所有顶点的最小凸多边形,称为S 的凸包,表示为CH(S)。如下图一所示,由红色线段表示的多边形就是点集Q={p0,p1,...p12}的凸包。 图一 凸包问题是计算机几何的一个经典问题,它可以解决很多优化模型,目前目前求取凸包的常用算法有:穷举法,格雷厄姆扫描法(Graham),分治法,蛮力法和Jarris步进法。本文主要讨论穷举法,蛮力法,以及格雷厄姆扫描法,通过算法思想的描述,计算出相应的时间复杂度,比较这几种算法的优劣。

二、穷举法 穷举法的思想很简单直接,它利用凸包性质—如果点集中的两个点的连线属于凸多边形的边,当且仅当点集中其余的点都在这两个点连线的同一侧。 假设点集合S 中有n 个顶点,则这n 个顶点共可以构成(1) 2n n -条边,对于任何一条边来讲,检查其余的n-2 个点的相对于该直线段的正负性。如果它们有相同的正负性,说明该边是凸多边形的边,反之就不属于凸多边形,继续检查。算法流程图如下: 不相同 相同 N Y 图二:算法流程图 上述算法(穷举法)需要执行n(n-1)/2 次,再次检查n-2 个点的正负性,故算法 时间复杂性为2(1)2 n n -∑=3()o n 。显然这并非一个高效的算法凸包问题有许多更加有效的算法可以达到log n n 的渐近时间复杂度。 三、蛮力法 开始 从点集S 中取出两个顶 点u ,v 判断其余的n-2 个点的相对于该直线段的正负性 判断执行次数是否 大于(1)2n n - 把u ,v 加入凸包 结束

凸包

算法分析与设计 课程设计 题目:凸包问题 所属专业:信息与计算科学 年级:1121102 学生姓名:江浩2011213361 学生姓名:张鸿2011213363 学生姓名:杨永涛2011213374 指导教师:于洪 评定成绩:

关于凸包问题的分析与求解 摘要 本文针对凸包问题进行分析和讨论,首先对凸包进行定义性描述并对其做具体分析,它要求将平面内的点集用最少的凸点来封闭,然后给出其具体模型,并对模型进行分析和求解,最后我们讨论几种求解算法(包括穷举法,蛮力法和分治法)的优劣性并给出问题的算法代码。 关键词:凸包问题穷举法分治法 一.引言 凸包(Convex Hull)是一个计算几何(图形学)中的概念。用不严谨的话来讲,给定二维平面上的点集,凸包就是将最外层的点连接起来构成的多边形,它能包含点集中的所有点。 凸包问题的完整描述如下: 令S 是平面上的一个点集,封闭S 中所有顶点的最小凸多边形,称为S 的凸包,表示为CH(S)。如下图一所示,由红色线段表示的多边形就是点集Q={p0,p1,...p12}的凸包。

图一 凸包问题是计算机几何的一个经典问题,它可以解决很多优化模型,目前目前求取凸包的常用算法有:穷举法,格雷厄姆扫描法(Graham),分治法,蛮力法和Jarris步进法。本文主要讨论穷举法,蛮力法,以及格雷厄姆扫描法。 二.问题分析 根据前面的描述我们先分析以下几种算法的特点: 1.穷举法 穷举法的思想很简单直接,它利用凸包性质—如果点集中的两个点的连线属于凸多边形的边,当且仅当点集中其余的点都在这两个点连线的同一侧。 假设点集合S中有n个顶点,则这n个顶点共可以构成 (1) 2 n n- 条边,对于任何一条边来讲, 检查其余的n-2 个点的相对于该直线段的正负性。如果它们有相同的正负性,说明该边是凸多边形的边,反之就不属于凸多边形,继续检查。穷举法需要执行n(n-1)/2次,再次检查n-2 个点的正负性,故算法时间复杂性为 2(1) 2 n n- ∑=3() o n。显然这并非一个高效的算法凸包 问题有许多更加有效的算法可以达到log n n的渐近时间复杂度。 算法流程图如下:

算法分析习题详细答案五

1.最大子段和问题:给定整数序列 n a a a ,,,21 ,求该序列形如 j i k k a 的子段和 的最大值: j i k k n j i a 1max ,0max 1) 已知一个简单算法如下: int Maxsum(int n,int a,int& best i,int& bestj){ int sum = 0; for (int i=1;i<=n;i++){ int suma = 0; for (int j=i;j<=n;j++){ suma + = a[j]; if (suma > sum){ sum = suma; besti = i; bestj = j; } } } return sum; }试分析该算法的时间复杂性。 2) 试用分治算法解最大子段和问题,并分析算法的时间复杂性。 3) 试说明最大子段和问题具有最优子结构性质,并设计一个动态规划算法解最大子段和问题。分析算法的时间复杂度。 (提示:令1()max ,1,2,,j k i j n k i b j a j n L ) 解:1)分析按照第一章,列出步数统计表,计算可得)(2 n O 2)分治算法:将所给的序列a[1:n]分为两段a [1:n/2]、a[n/2+1:n],分别求出这两段的最大子段和,则a[1:n]的最大子段和有三种可能: ①a[1:n]的最大子段和与a[1:n/2]的最大子段和相同; ②a[1:n]的最大子段和与a[n/2+1:n]的最大子段和相同; ③a[1:n]的最大子段和为两部分的字段和组成,即 j n j i l n i j a a a a a 122;

intMaxSubSum ( int *a, int left , int right){ int sum =0; if( left==right) sum = a[left] > 0? a[ left]:0 ; else {int center = ( left + right) /2; int leftsum =MaxSubSum ( a, left , center) ; int rightsum =MaxSubSum ( a, center +1, right) ; int s_1 =0; int left_sum =0; for ( int i = center ; i >= left; i--){ left_sum + = a [ i ]; if( left_sum > s1) s1 = left_sum; } int s2 =0; int right_sum =0; for ( int i = center +1; i <= right ; i++){ right_sum + = a[ i]; if( right_sum > s2) s2 = right_sum; } sum = s1 + s2; if ( sum < leftsum) sum = leftsum; if ( sum < rightsum) sum = rightsum; } return sum; } int MaxSum2 (int n){ int a; returnMaxSubSum ( a, 1, n) ; } 该算法所需的计算时间T(n)满足典型的分治算法递归分式T(n)=2T(n/2)+O(n),分治算法的时间复杂度为O(nlogn)

相关文档
最新文档