matlab最大公约数 三种算法讲解学习

matlab最大公约数 三种算法
算法设计与分析
11信本余启盛118632011004
一、上机目的及内容
1.上机内容
求两个自然数m和n的最大公约数。
2.上机目的
（1）复习数据结构课程的相关知识，实现课程间的平滑过渡；
（2）掌握并应用算法的数学分析和后验分析方法；
（3）理解这样一个观点：不同的算法能够解决相同的问题，这些算法的解题思路不同，复杂程度不同，解题效率也不同。
k=f3(m,n);
printf("方法三最大公约数为：%d\n",k);
printf("\n--------------------\n");
printf("\n计数器显示结果：\n\n\n");
end
r=n;
toc;
算法三：
function factor=CommFactor3(m,n)
tic;
factor=1;
themax=max(m,n);
for i=2:1:themax
while (mod(m,i)==0)&&(mod(n,i)==0)
factor=factor*i;
m=m/i;
n=n/i;
{
int t;
if(m>n)t=n;
else t=m;
while(t)
{
if(m%t==0&&n%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++;
2.1如果i是m和n的公因子，则执行下述操作：
2.1.1 factor=factor*i；
2.1.2 m = m / i; n = n / i;
2.2如果i不是m和n的公因子，则i=i+1；
3.输出factor；
根据代码考虑最坏情况他们的最大公约数，循环做了i-1次;最好情况是只做了1次，可以得出:O（n）=n/2;
2.循环直到r等于0
2.1 m = n;
2.2 n = r;
2.3 r = m % n;
3.输出n ;
根据代码辗转相除得到欧几里得的:O(n)=logn
算法三：蛮力法（短除法）
CommFactor3
输入：两个自然数m和n
输出：m和n的最大公约数
1.n），执行下述操作：
2.如果m%t==0&&n%t==0，结束
2.1如果t不是m和n的公因子，则t=t-1;
3.输出t ;
根据代码考虑最坏情况他们的最大公约数是1，循环做了t-1次，最好情况是只做了1次，可以得出O（n）=n/2;
算法二：欧几里德算法
CommFactor2
输入：两个自然数m和n
输出：m和n的最大公约数
1. r = m % n;
tic;
ifm>n)
t=n;
else t=m;
while(t)
if(m%t==0&&n%t==0)break;
else t=t-1;
end
end
r=t
toc;
算法二：
function r=CommFactor2(m,n)
tic;
r=mod(m,n);
while r~=0
m=n;
n=r;
r=mod(m,n);
}
return n;
}
int f3(int m,int n)
{
int i, factor = 1;
for (i = 2; i <= m && i <= n; i++)
{
while (m % i == 0 && n % i == 0) //此处不能用if语句
{
factor = factor * i;
m = m / i; n = n / i;
二、实验原理及基本技术路线图
（1）至少设计出三个版本的求最大公约数算法；
（2）对所设计的算法采用大O符号进行时间复杂性分析；
（3）上机实现算法，并用计数法和计时法分别测算算法的运行时间；
（4）通过分析对比，得出自己的结论。
三、所用仪器、材料（设备名称、型号、规格等或使用软件）
1台PC及VISUAL C++6.0软件matlab .2008
end
end
toc;
三种算法时间复杂度比较：（c++语言）
#include"iostream.h"
#include"stdio.h"
#include"stdlib.h"
#include"time.h"
#define N 100
int w,w2,w3;//用于计数
int f1(int m,int n)
w3++;
}
}
return factor;
}
int main(void)
{
int m,n;
printf("请输入m,n :\n");
scanf("%d%d",&m,&n);
int k;
k=f1(m,n);
printf("方法一最大公约数为：%d\n",k);
k=f2(m,n);
printf("方法二最大公约数为：%d\n",k);
MATLAB程序代码：
main.m
x=fix(rand(1,1000)*1000);
y=fix(rand(1,1000)*1000);
fori=1:1000
A(i)=CommFactor2(x(i),y(i));
end
x=x';
y=y';
算法一：
function r=CommFactor1(m,n)
四、实验方法、步骤（或：程序代码或操作过程）
实验采用三种方法求最大公约数
1、连续整数检测法。
2、欧几里得算法
3、蛮力法（短除法）
根据实现提示写代码并分析代码的时间复杂度：
算法一：连续整数检测法。
CommFactor1
输入：两个自然数m和n
输出：m和n的最大公约数
1.判断m和n哪个数小，t=min（m，n）