/****************************************************************************
File name:求最大公约数
Description: 1.辗转相除法(又名欧几里德法))C语言中用于计算两个正整数a,b的最大公约数和最小公倍数,实质它依赖
于下面的定理:
a b=0
gcd(a,b)={
gcd(b,a mod b) b!=0
根据这一定理可以采用函数嵌套调用形式进行求两个数的最大公约数和最小公倍数。
①函数嵌套调用
其算法过程为:?前提:设两数为a,b设其中a 做被除数,b做除数,temp为余数
1、大数放a中、小数放b中;
2、求a/b的余数;
3、若temp=0则b为最大公约数;
4、如果temp!=0则把b的值给a、temp的值给a;
5、返回第二步;
②函数递归调用
3.穷举法(也叫枚举法),穷举法求两个正整数的最大公约数的解题步骤:从两个数中较小数开始由大到小列
举,直到找到公约数立即中断列举,得到的公约数便是最大公约数。
4.更相减损法:
第一步:任意给定两个正整数;判断它们是否都是偶数。若是,则用2约简;若不是则执行第二步。
第二步:以较大的数减较小的数,接着把所得的差与较小的数比较,并以大数减小数。继续这个操作,直到所
得的减数和差相等为止。
则第一步中约掉的若干个2与第二步中等数的乘积就是所求的最大公约数。
其中所说的“等数”,就是最大公约数。求“等数”的办法是“更相减损”法。所以更相减损法也叫等值算法。
5.Stein算法:对两个正整数 x>y :
①均为偶数 gcd( x,y ) = 2gcd( x/2,y/2 );
②均为奇数 gcd( x,y ) = gcd( (x+y)/2,(x-y)/2 );
③x奇y偶 gcd( x,y ) = gcd( x,y/2 );
④x偶y奇 gcd( x,y ) = gcd( x/2,y )? 或 gcd( x,y )=gcd( y,x/2 );
Author: 软件工程1704班 赵茹平
Version: 用C语言实现
Date: 2019年3月5日
The last revision time:2019年3月9日
*****************************************************************************/
#include <stdio.h>
#include <stdlib.h>
#include<time.h>
#include<math.h>
clock_t start,stop;/*clock_t是clock()函数返回的变量类型*/
double duration1,duration2,duration3,duration4,duration5;/*记录被测函数运行时间,以秒为单位*/
int divisor1 (int a,int b) /*自定义函数求两数的最大公约数*/
{
int temp; /*定义整型变量*/
if(a<b) /*通过比较求出两个数中的最大值和最小值*/
{ temp=a;a=b;b=temp;} /*设置中间变量进行两数交换*/
while(b!=0) /*通过循环求两数的余数,直到余数为0*/
{
temp=a%b;
a=b; /*变量数值交换*/
b=temp;
}
return (a); /*返回最大公约数到调用函数处*/
}
int divisor2(int a,int b) /*自定义函数求两数的最大公约数*/
{
int temp; /*定义整型临时变量*/
if(a<b) /*通过比较求出两个数中的最大值和最小值*/
{
temp = a;/*a为大数*/
a = b;
b = temp;
}
while(b!=0) /*通过循环求两数的余数,直到余数为0*/
{
temp = a%b;
a = b;
b = temp;
}
return a; /*返回最大公约数到调用函数处*/
}
int divisor3 (int a,int b) /*自定义函数求两数的最大公约数*/
{
int temp; /*定义义整型变量*/
temp = (a>b)?b:a; /*采种条件运算表达式求出两个数中的最小值*/
while(temp>0)
{
if (a%temp == 0&&b%temp == 0) /*只要找到一个数能同时被a,b所整除,则中止循环*/
break;
temp--; /*如不满足if条件则变量自减,直到能被a,b所整除*/
}
return temp; /*返回满足条件的数到主调函数处*/
}
int gcd1(int a,int b)/*自定义函数求两数的最大公约数*/
{
int i,temp,x;
i = 0;/*i用来记约掉2的个数*/
while(a%2 == 0&&b%2 == 0) /*判断m和n能被多少个2整除*/
{
a /= 2;
b /= 2;
i += 1;
}
if(a<b) /*m保存大的值*/
{
temp = a;
a = b;
b = temp;
}
while(x)/*通过循环求两数的差,直到差与减数相等*/
{
x = a-b;
a = (b>x)?b:x;/*差值与n作比较,m存大的数*/
b = (b<x)?b:x;/*n存较小的数*/
if(b == (a-b))
break;
}
if(i == 0)
return b;
else
return (int)pow(2,i)*b;/*返回最大公约数*/
}
int gcd2(int u,int v)
{
if(u == 0) return v;/*若u为0,则v为最大公约数*/
if(v == 0) return u;/**//*若v为0,则u为最大公约数*/
if(~u&1) /*u是偶数*/
{
if(v&1) /*v是奇数*/
return gcd2(u >> 1,v);
else /*u和v都是偶数*/
return gcd2(u >> 1,v >> 1) << 1;
}
else if(~v&1) /*u是奇数,v是偶数*/
return gcd2(u,v >> 1);
else /*u和v都是奇数*/
if (u > v)
return gcd2((u - v) >> 1, v);
return gcd2((v - u) >> 1, u);/*返回最大公约数*/
}
int main()
{
int i,t1,m[20],n[20];
srand((unsigned)time(NULL));/*初始化随机数*/
printf("第一组数为:");
for(i=0;i<20;i++)/*输出20组数据*/
{//scanf("%d",&m[i]);
m[i]=rand()%100;/*随机数范围是1到100*/
printf(" %d ",m[i]);/*输出随机数*/
}
printf("\n");
printf("第二组数为:");
for(i=0;i<20;i++)/*输出20组数据*/
{//scanf("%d",&n[i]);
n[i]=rand()%100;/*随机数范围是1到100*/
printf(" %d ",n[i]);/*输出随机数*/
}
printf("\n函数嵌套调用求最大公约数:\n");
start=clock();
for(i=0;i<20;i++)/*输出20组数据*/
{
t1 = divisor1(m[i],n[i]);/*求两数的最大公约数*/
printf("第%d组数为:%d %d\n",i+1,m[i],n[i]);
printf("The higest common divisor is :%d\n",t1); /*输出最大公约数*/
}
stop=clock();
duration1=((double)(stop-start))/CLK_TCK;
printf("\nThe First Running Time Is:%fs\n",duration1);
printf("函数递归调用求最大公约数:\n");
start=clock();
for(i=0;i<20;i++)/*输出20组数据*/
{
t1 = divisor2(m[i],n[i]); /*求两数的最大公约数*/
printf("第%d组数为:%d %d\n",i+1,m[i],n[i]);
printf("The higest common divisor is :%d\n",t1); /*输出最大公约数*/
}
stop=clock();
duration2=((double)(stop-start))/CLK_TCK;
printf("\nThe Second Running Time Is:%fs\n",duration2);
printf("穷举法求最大公约数:\n");
start=clock();
for(i=0;i<20;i++)/*输出20组数据*/
{
t1 = divisor3(m[i],n[i]); /*求两数的最大公约数*/
printf("第%d组数为:%d %d\n",i+1,m[i],n[i]);
printf("The higest common divisor is: %d\n",t1); /*输出最大公约数*/
}
stop=clock();
duration3=((double)(stop-start))/CLK_TCK;
printf("\nThe Third Running Time Is:%fs\n",duration3);
printf("更相减损法求最大公约数:\n");
start=clock();
for(i=0;i<20;i++)/*输出20组数据*/
{
t1=gcd1(m[i],n[i]); /*求两数的最大公约数*/
printf("第%d组数为:%d %d\n",i+1,m[i],n[i]);
printf("The higest common divisor is:%d\n",t1); /*输出最大公约数*/
}
stop=clock();
duration4=((double)(stop-start))/CLK_TCK;
printf("\nThe Fourth Running Time Is:%fs\n",duration4);
printf("Stein算法求最大公约数:\n");
start=clock();
for(i=0;i<20;i++)/*输出20组数据*/
{
t1=gcd2(m[i],n[i]); /*求两数的最大公约数*/
printf("第%d组数为:%d %d\n",i+1,m[i],n[i]);
printf("The higest common divisor is:%d\n",t1); /*输出最大公约数*/
}
stop=clock();
duration5=((double)(stop-start))/CLK_TCK;
printf("\nThe Fifth Running Time Is:%fs\n",duration5);
return 0;
}
流程图如下:
1.函数嵌套调用
2.函数递归调用
3.穷举法
4.更相减损法
经验总结如下:
在本次程序算法设计中我学会了:
1.如何调用时间函数
#include<time.h>是头文件
clock_t start,stop;/*clock_t是clock()函数返回的变量类型*/
double duration;/*记录被测函数运行时间,以秒为单位*/
int main()
{
start=clock();/*开始时间*/
MyFunction();/*把被测函数加在这里*/
stop=clock();/*停止时间*/
duration=((double)(stop-start))/CLK_TCK;/*计算运行时间*/
printf("\nThe First Running Time Is:%fs\n",duration);
return 0;
}
2.如何调用随机数组
#include<math.h>
int main()
{
int i,m[20];
srand((unsigned)time(NULL));/*初始化随机数*/
for(i=0;i<50;i++)/*输出50组数据*/
{
m[i]=rand()%100;/*随机数范围是1到100*/
printf(" %d ",m[i]);/*输出随机数*/
}
}
这次程序算法设计我收获很大,在设计过程中虽然遇到了很多问题,但我通过向同学求助与上网查找资料逐一击破了难题,也许我的收获是别人早都会的,但我会尽自己最大的努力完成每一次编程,找到自己的问题,弥补自己的差距。