最大公约数的三种算法

希望对大家有所帮助，多谢您的浏览！昆明理工大学信息工程与自动化学院学生实验报告（2011—2012学年第1学期）课程名称：算法设计与分析开课实验室：信自楼机房4442011年10月12日年级、专业、班计科092学号200910405214姓名徐兴繁成绩实验项目名称求最大公约数指导教师吴晟教师评语该同学是否了解实验原理： A.了解口 B.基本了解口C.不了解口该同学的实验能力： A.强口 B.中等 口 C.差 口该同学的实验是否达到要求：A.达到口 B.基本达到口C.未达到口实验报告是否规范： A.规范口 B.基本规范口C.不规范口实验过程是否详细记录： A.详细口 B.一般口 C.没有口教师签名：年 月 日一、上机目的及内容上机内容求两个自然数m和n的最大公约数。上机目的（1） 复习数据结构课程的相关知识，实现课程间的平滑过渡；（2） 掌握并应用算法的数学分析和后验分析方法；（3） 理解这样一个观点：不同的算法能够解决相同的问题，这些算法的解题思路不同，复杂程度不同，解题效率也不同。二、实验原理及基本技术路线图（方框原理图或程序流程图）（1）至少设计出三个版本的求最大公约数算法；（2） 对所设计的算法采用大0符号进行时间复杂性分析；（3） 上机实现算法，并用计数法和计时法分别测算算法的运行时间（4） 通过分析对比，得出自己的结论。三、 所用仪器、材料（设备名称、型号、规格等或使用软件）1台PC及VISUALC++6.0软件四、 实验方法、步骤（或：程序代码或操作过程）希望对大家有所帮助，多谢您的浏览！实验采用三种方法求最大公约数1、连续整数检测法。2、欧几里得算法3、分解质因数算法根据实现提示写代码并分析代码的时间复杂度：方法一：intf1(intm,intn){intt;if(m>n)t=n;elset=m;while(t){if(m%t==0&&n%t==0)break;elset=t-1;}returnt;}根据代码考虑最坏情况他们的最大公约数是1,循环做了t-1次，最好情况是只做了1次，可以得出0(n)=n/2;方法二：intf2(intm,intn){intr;r=m%n;while(r!=0){m=n;n=r;r=m%n;}returnn;}根据代码辗转相除得到欧几里得的0(n)=logn方法三：intf3(intm,intn){inti=2,j=0,h=0;inta[N],b[N],c[N];while(i<n)if(n%i==0){j++;a[j]=i;n=n/i;}elsei++;}j++;a[j]=n;i=1;intu;u=j;while(i<=j){//printf("%d",a[i]);i++;}i=2;j=0;while(i<m){if(m%i==0){j++;b[j]=i;m=m/i;}elsei++;}j++;b[j]=m;i=1;while(i<=j){//printf("%d",b[i]);i++;}intk=1;for(i=1;i<=j;i++){for(k=1;k<=u;k++){if(b[i]==a[k]){h++;c[h]=a[k];//printf("%d",c[h]);a[k]=a[k+1];break;}}}k=1;while(h>1){k=k*c[h]*c[h-1];h=h-2;}if(h==1){k=k*c[1];returnk;}elsereturnk;}根据代码分解质因子算法0(n)=n2+n/2为了计算每种算法运行的次数所用的时间，我将代码稍加改动添加代码如下其中计数器采用的是没做一次循环就加1；计时器是记住开始时间和结束时间，用结束时间减开始时间。#include"iostream.h"#include"stdio.h"#include"stdlib.h"#include"time.h"#defineN100intw,w2,w3;//用于计数intf1(intm,intn){intt;if(m>n)t=n;elset=m;while(t){if(m%t==0&&n%t==0)break;elset=t-1;w++;}returnt;}intf2(intm,intn){intr;r=m%n;w2=1;while(r!=0){m=n;n=r;r=m%n;w2++;returnn;intf3(intm,intn){inti=2,j=0,h=0;inta[N],b[N],c[N];while(i<n){if(n%i==0){j++;a[j]=i;n=n/i;w3++;}else{i++;w3++;}}j++;a[j]=n;i=1;intu;u=j;while(i<=j){//printf("%d",a[i]);i++;w3++;}//printf("\n");i=2;j=0;while(i<m){if(m%i==0){j++;b[j]=i;m=m/i;w3++;elsew3++;}}j++;b[j]=m;i=1;while(i<=j){//printf("%d",b[i]);i++;w3++;}intk=1;for(i=1;i<=j;i++){for(k=1;k<=u;k++){if(b[i]==a[k]){w3++;希望对大家有所帮助，多谢您的浏览！h++;c[h]=a[k];//printf("\n%d",c[h]);a[k]=a[k+1];break;}}}k=1;while(h>1){k=k*c[h]*c[h-1];h=h-2;w3++;}if(h==1){k=k*c[1];returnk;}elsereturnk;intmain(void)intm,n;printf（“ 请输入m,n:\n"）;scanf（"%d%d",&m,&n）;intk;k=f1（m,n）;printf（"方法一最大公约数为：%d\n",k）;k=f2（m,n）;printf（" 方法二最大公约数为：%d\n",k）;k=f3（m,n）;printf（" 方法三最大公约数为：%d\n",k）;printf（"\n \n"）;printf（"\n计数器显示结果：\n\n\n"）;printf（"方法一：%d\n",w2）;printf（"方法二：%d\n",w）;printf（"方法三：%d\n",w3）;printf（"\n \n"）;floata,i;clock_tstart,finish;doubleusetime;i=0;start=clock（）;

while(i<1000000){f1(m,n);i++;}finish=clock();usetime=finish-start;printf(“ 方法一用时％.f*10八(-6)豪秒\n",i=0;start=clock();while(i<1000000){f2(m,n);i++;}finish=clock();usetime=finish-start;printf(“ 方法二用时%.f*10八(-6)豪秒\n",i=0;start=clock();usetime);usetime);while(i<1000000)usetime);usetime);f3(m,n);i++;finish=clock();usetime=finish-start;printf(“ 方法三用时％.f*10八(-6)豪秒\n",usetime);}五、实验过程原始记录(测试数据、图表、计算等)请给出各个操作步骤的截图和说明；三种算法得到结果验证结果：计数器：我想到的是做一次循环就加一计算算法运行时间结果：在计算时间过程中因为计算机的运算速度很快，所以我利用了循环把时间精确得到10-6毫秒六、实验结果、分析和结论(误差分析与数据处理、成果总结等。其中，绘制曲线图时必须用计算纸或程序运行结果、改进、收获)请结合实验的结果分析算法原理；在实验中遇到了些什么问题，如何解决；有什么收获等；在本次实验中代码是独自完成的，一开始我感觉这个代码最多半小时就可以完成，但是第三个算法的时候我分析了好久才写出来，在计算三种方法运行时间的时候，我一开始只精确到毫秒(ms),计算结果都是零，后面我写了一个循环调试才发现是我的精确度还在不够，所以我想到了计算算法执行了1000000次之后所用的时间，然后再求平均每次执行的时间。结果分析：从前面的复杂度0(n)的出欧几里得算法的是最优算法，连续整除法其次，最复杂的是分解质因数算法，再从代码运行的计数器和计算的时间来看结果恰好和前面的复杂度得到的结果一致，所以的出结论：欧几里得算法最