为了账号安全,请及时绑定邮箱和手机立即绑定

求最大公约数的4种算法

标签:
C
for(z=0; z<10000000; z++) 循环只是为了增加程序的运行时间,
让我们体会算法的时间复杂度。

算法一:短除法

想法,采用短除法找出2个数的所有公约数,将这些公因子相乘,结果就是2个数的最大公约数。【找公因子,只能使用蛮力法】

#include<stdio.h>
#include<time.h>
void main()
{ 
int m=28,n=72;

int i,f=1;
 int z;

clock_t start,finish;

double duration;

start= clock();

for(z=0; z<10000000; z++)
{

  for(i=2;i<=m&&i<=n;)
  {

    while(m%i==0&&n%i==0)
    {

     f*=i;
     m/=i;
     n/=i;
    }
    i++;    
    }

}

finish=clock();

duration=(double)(finish-start)/CLOCKS_PER_SEC;

printf("time=%lf seconds\n",duration);

printf("result=%d\n",f);

}

image.png

算法二:辗转相除法

辗转相除法, 又名欧几里德算法(Euclidean algorithm),是求最大公约数的一种方法。它的具体做法是:用较小数除较大数,再用出现的余数(第一余数)去除除数,再用出现的余数(第二余数)去除第一余数,如此反复,直到最后余数是0为止。如果是求两个数的最大公约数,那么最后的除数就是这两个数的最大公约数。


#include<stdio.h>
#include<time.h>
void main()
{ 
int m=28,n=72;

 int i,f=1;
 int z;

clock_t start,finish;

double duration;

start= clock();

for(z=0; z<10000000; z++)
{

   if(m<n)
   {
    i=m;
    m=n;
    n=i;
   }

 i=m%n;
 while(i)
 {

  m=n;
  n=i;
  i=m%n;
}

f=n;

}

finish=clock();

duration=(double)(finish-start)/CLOCKS_PER_SEC;

printf("time=%lf seconds\n",duration);

printf("result=%d\n",f);

}
time=0.037000 seconds
result=4
Press any key to continue

算法三:蛮力法,从2个公约数中较小的数开始递减,二个公约数除以它,可以同时除尽,变是最大公约数,我想的,很笨的一种。


#include<stdio.h>
#include<time.h>
void main()
{ 
int m=28,n=72;

int i,f=1;
 int z;

clock_t start,finish;

double duration;

start= clock();

for(z=0; z<10000000; z++)
{

f=m>n?n:m;

for(i=f;i>0;i--)
{
  if(m%i==0&&n%i==0)
  {
     f=i;
   break;
  }

}

}

finish=clock();

duration=(double)(finish-start)/CLOCKS_PER_SEC;

printf("time=%lf seconds\n",duration);

printf("result=%d\n",f);

}
time=0.992000 seconds
result=4
Press any key to continue

算法四:辗转相减法

辗转相减法是一种简便的求出两数最大公约数的方法。(更相减损术)辗转相减法(求最大公约数),即尼考曼彻斯法,其特色是做一系列减法,从而求得最大公约数。例如 :两个自然数35和14,用大数减去小数,(35,14)->(21,14)->(7,14),此时,7小于14,要做一次交换,把14作为被减数,即(14,7)->(7,7),再做一次相减,结果为0,这样也就求出了最大公约数7

#include<stdio.h>
#include<time.h>
void exchange(int *m,int *n);
void main()
{ 
int m=28,n=72;

int i,f=1;
 int z;

clock_t start,finish;

double duration;

start= clock();

if(m<n)
exchange(&m,&n);

for(z=0; z<10000000; z++)
{

    i=m-n;
    while(i)
    {

     if(i>n)
     {
         m=i;
     }else
     {
       m=n;
       n=i;
     }

     i=m-n;
    }

     f=m;
}

finish=clock();

duration=(double)(finish-start)/CLOCKS_PER_SEC;

printf("time=%lf seconds\n",duration);

printf("result=%d\n",f);

}

void exchange(int *m,int *n)
{
 int temp;
  temp=*m;
  *m=*n;
  *n=temp;

}
time=0.020000 seconds
result=4
Press any key to continue

看看4个算法的运行时间,还是我自己想的算法时间最久,来一个
大神拯救我吧。

点击查看更多内容
8人点赞

若觉得本文不错,就分享一下吧!

评论

作者其他优质文章

正在加载中
全栈工程师
手记
粉丝
73
获赞与收藏
369

关注作者,订阅最新文章

阅读免费教程

感谢您的支持,我会继续努力的~
扫码打赏,你说多少就多少
赞赏金额会直接到老师账户
支付方式
打开微信扫一扫,即可进行扫码打赏哦
今天注册有机会得

100积分直接送

付费专栏免费学

大额优惠券免费领

立即参与 放弃机会
意见反馈 帮助中心 APP下载
官方微信

举报

0/150
提交
取消