算法最大公约数.docVIP

实验项目——求最大公约数 1. 实验题目 求两个自然数m和n的最大公约数。 2. 实验目的 ⑴ 复习数据结构课程的相关知识，实现课程间的平滑过渡； ⑵ 掌握并应用算法的数学分析和后验分析方法； ⑶ 理解这样一个观点：不同的算法能够解决相同的问题，这些算法的解题思路不同，复杂程度不同，解题效率也不同。 3. 实验要求 ⑴ 至少设计出三个版本的求最大公约数算法； ⑵ 对所设计的算法采用Θ符号进行时间复杂性分析； ⑶上机实现算法，并用计数法和计时法分别测算算法的运行时间； ⑷ 通过分析对比，得出自己的结论。 一．连续整数检测： #includeiostream.h int f1(int m,int n) { int t; if(mn)t=n; else t=m; while(t) { if(m%t==0n%t==0)break; else t=t-1; } return t; } void main() { int a,b; cout请输入两个整数： ; cinab; cout最大公约数是： ; coutf1(a,b)endl; } 根据代码考虑最坏情况他们的最大公约数是1，循环做了t-1次，最好情况是只做了1次，可以得出O（n）=n/2; 二．欧几里得算法#includestdio.h void main() { int r,n,m,t; printf(请输入两个正整数:\n); scanf(%d,%d,n, m); if(nm)/*选出较大数*/ { t=n; n=m; m=t; } while(m!=0)/*辗转取余算法*/ { r=n%m; n=m; m=r; } printf(它们的最大公约数为:%d\n,n); } 根据代码辗转相除得到欧几里得的O(n)= log n 三．递归算法 #includestdio.hint f(long a,long b){long s,t=a%b; if(t==0)s=b;else s=f(b,t); return s;}int main(){long a,b,t; scanf(%ld%ld,a,b); if(ab){t=a;a=b;b=t;} printf(%ld\n,f(a,b));}} 时间复杂度：T(n)=O(n); 基本语句执行次数C(n)=n#includeiostream.h #includestdio.h #includestdlib.h #includetime.h #define N 100 int w,w2,w3;//用于计数 int f1(int m,int n) { int t; if(mn)t=n; else t=m; while(t) { if(m%t==0n%t==0)break; else t=t-1; w++; } return t; } int f2(int m,int n) { int r; r=m%n;w2=1; while(r!=0) { m=n; n=r; r=m%n; w2++; } return n; } int f3(int m,int n) { int i=2,j=0,h=0; int a[N],b[N],c[N]; while(in) { if(n%i==0) { j++; a[j]=i; n=n/i; w3++; } else { i++; w3++; } } j++; a[j]=n; i=1; int u; u=j; while(i=j) { //printf(%d ,a[i]); i++; w3++; } //printf(\n); i=2; j=0; while(im) { if(m%i==0) { j++; b[j]=i; m=m/i; w3++; } else { i++; w3++; } } j++; b[j]=m; i=1; while(i=j) { //printf(%d ,b[i]); i++; w3++; } int k=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(h1) { k=k*c[h]*c[h