求最大公约数的四种算法

最新推荐文章于 2024-06-09 11:58:49 发布
最新推荐文章于 2024-06-09 11:58:49 发布
阅读量1.1w 收藏 31
点赞数 12
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/qq_41123884/article/details/88363584
版权

一． 题目分析

如果有一个自然数a能被自然数b整除，则称a为b的倍数，b为a的约数。几个自然数公有的约数，叫做这几个自然数的公约数。公约数中最大的一个公约数，称为这几个自然数的最大公约数。

根据约数的定义可知，某个数的所有约数一定不大于这个数本身，几个自然数的最大公约数也一定不大于其中任何一个数。所以当要求任意两个正整数的最大公约数，即是求出一个不大于其中两数中的任何一个，但又能同时整除两个整数的最大自然数。

二．算法构造（流程图）

1．辗转相除法


2.穷举法
在这里插入图片描述
3.更相减损法
在这里插入图片描述
4.stein算法

在这里插入图片描述

三．算法实现

#include <stdio.h>
#include <time.h>
#include <math.h>

int  gcd(int a,int b)//辗转相除法
{  
	if(a%b==0)//若a能整除b则b为最大公约数
       return b;   
	else  
       return gcd(b,a%b);//否则进行递归
}

void gcd_input(int m,int n)
{ 
	
int t1;
t1=gcd(m,n);
printf("最大公约数是 %d\n",t1);//得出最大公约数
}


void multiple (int a,int b)//穷举法
{
  int p,q,temp;
  p=(a>b)?a:b;   //求两个数中的最大值
  q=(a>b)?b:a;  //求两个数中的最小值
  temp=p;      //将最大值赋给p为变量
  while(1)   //利用循环语句来求满足条件的数值
  {
    if(p%q==0)
      break;  //只要找到变量的和数能被a或b所整除，则中止循
p+=temp;   //如果条件不满足则变量自身相加
  }
 printf("最大公约数是 %d\n",a*b/p);//得出最大公约数a*b/p
}

void gcd2(int m,int n)//更损相减法
{
	int i=0,temp,x;
	while(m%2==0 && n%2==0)  //判断m和n能被多少个2整除
	{
		m/=2;
		n/=2;
		i+=1;//用i来储存2的个数
	}
	if(m<n)     //m保存大的值
	{
		temp=m;
		m=n;
		n=temp;
	}
	while(x)
	{
		x=m-n;//计算两个数的差
		m=(n>x)?n:x;//比较m和n，和差的打消
		n=(n<x)?n:x;
		if(n==(m-n))
			break;
	}
	if(i==0)
		printf("最大公约数是%d\n",n);
	else 
		printf("最大公约数是%d\n",(int )pow(2,i)*n);


}

 int Stein( unsigned int x, unsigned int y )//Stein算法 

  {
int factor = 0;// 当两数均为偶数时将其同时除以2至至少一数为奇数为止，
//记录除掉的所有公因数2的个数factor     
  		int temp;
        if ( x < y )//让 x为两者中大的数
        {
                temp = x;
                x = y;
                y = temp;
        }
        if ( 0 == y )
        {
                return 0;
        }
        while ( x != y )
        {
                if ( x & 0x1 )//判断x是否为奇数
                {
                        if ( y & 0x1 )//再判断y是否为奇数
                        {                           
    y = ( x - y ) >> 1;//将两数之差除于2
                                x -= y;
                        }
                        else//当x为奇数，y为偶数是when x is odd and y is even
                        {
                                y >>= 1;
                        }
                }
                else//当x为偶数
                {
                        if ( y & 0x1 )// 当x为偶数，y为奇数
                        {
                                x >>= 1;
                                if ( x < y )
                                {
                                        temp = x;
                                        x = y;
                                        y = temp;
                                }
                        }
                        else //当x和y都为偶数
                        {
                                x >>= 1;
                                y >>= 1;
                                ++factor;
                        }
                }
        }

		printf("最大公约数是%d\n",x << factor);
        
}


void main()
{  

	clock_t start, finish; 
	double duration; 

	int i,j,x;
    int flag=1;
	int str[20]={345,12350,34363,3345,422,6755,536,57,375,57,37,975,775,4735,673,5633,563,339,596,84};
	
    while(flag)
	{
		
		printf("请选择您所使用的方法,按F键结束:\n");

		printf("1.辗转相除法\n");

		printf("2.穷举法\n");

		printf("3.更损相减法\n");

		printf("4.stein算法\n");

		scanf("%d",&x);

		if (x=='F')
			break;

	switch(x)

	{
		case 1:{
					printf("辗转相除法所得到的结果为:\n");

					start = clock(); 

					for( i=0;i<20;i++)
					{
						if(str[i]<0||str[i+1]<0||str[i] != (int)str[i]||str[i+1] != (int)str[i+1])
							printf("输入非法请重新输入!");//输入非法的判断

						gcd_input(str[i],str[i+1]);
					}
	
						finish = clock(); 

						duration = (double)(finish - start) / CLOCKS_PER_SEC; 
		 
						break;
					
			   }
		case 2:{
					printf("穷举法所得到的结果为:\n");

					start = clock(); 

					for( i=0;i<20;i++)
					{
						if(str[i]<0||str[i+1]<0||str[i] != (int)str[i]||str[i+1] != (int)str[i+1])
							printf("输入非法请重新输入!");
				
						multiple(str[i],str[i+1]);
					}
		
						finish = clock(); 
	
						duration = (double)(finish - start) / CLOCKS_PER_SEC; 
		 
						break;
				
				}

		case 3:{
					printf("更损相减法所得到的结果为:\n");

					start = clock(); 

					for( i=0;i<20;i++)
					{
						if(str[i]<0||str[i+1]<0||str[i] != (int)str[i]||str[i+1] != (int)str[i+1])
							printf("输入非法请重新输入!");
		
						gcd2(str[i],str[i+1]);
					}
		
						finish = clock(); 

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

		case 4:{
					printf("stein法所得到的结果为:\n");
			   			   
				 start = clock(); 

					for( i=0;i<20;i++)
					{
					if(str[i]<0||str[i+1]<0||str[i] != (int)str[i]||str[i+1] != (int)str[i+1])
							printf("输入非法请重新输入!");
						Stein(str[i],str[i+1]);
					}
		
						finish = clock(); 

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

	printf( "所花费时间为%f seconds\n\n", duration);

	}

}

四．调试、测试及运行结果
1.测试
1.1辗转相除法

#include <stdio.h>
#include <time.h>
#include <math.h>
int  gcd(int a,int b)//辗转相除法
{  
	if(a%b==0)
       return b;   
	else  
       return gcd(b,a%b);
}

void gcd_input(int m,int n)
{ 
	
int t1;

t1=gcd(m,n);

printf("最大公约数%d\n",t1);/*最大公约数*/

}

void main()
{
	long k = 10000000; 
	clock_t start, finish; 
	double duration; 
	int i,j,x;
    int flag=1;
	int str[20]={90,177,83,45,42,35,66,54,89,54,64,96,785,64,63,543,35,345,56,54};
	start = clock(); 
	for( i=0;i<20;i++)
		gcd_input(str[i],str[i+1]);
	finish = clock(); 
	duration = (double)(finish - start) / CLOCKS_PER_SEC; 
	printf( "所花费时间为%f seconds\n\n", duration);
}

辗转相除法测试结果
1.2穷举法

#include <stdio.h>
#include <time.h>
#include <math.h>

void multiple (int a,int b)//穷举法
{
  int p,q,temp;
  p=(a>b)?a:b;   /*求两个数中的最大值*/
  q=(a>b)?b:a;  /*求两个数中的最小值*/
  temp=p;      /*最大值赋给p为变量自增作准备*/
  while(1)   /*利用循环语句来求满足条件的数值*/
  {
    if(p%q==0)
      break;  /*只要找到变量的和数能被a或b所整除，则中止循环*/
    p+=temp;   /*如果条件不满足则变量自身相加*/
  }
 printf("最大公约数是 %d\n",a*b/p);
}

void main()
{
		clock_t start, finish; 
	double duration; 
	int i,j,x;
    int flag=1;
	int str[20]={90,177,83,45,42,35,66,54,89,54,64,96,785,64,63,543,35,345,56,54};
	start = clock(); 
	for( i=0;i<20;i++)
		multiple(str[i],str[i+1]);
	finish = clock(); 
	duration = (double)(finish - start) / CLOCKS_PER_SEC; 
	printf( "所花费时间为%f seconds\n\n", duration);
}

穷举法测试结果
1.3更相减损法

#include <stdio.h>
#include <time.h>
#include <math.h>

void gcd2(int m,int n)//更损相减法
{
	int i=0,temp,x;
	while(m%2==0 && n%2==0)  //判断m和n能被多少个2整除
	{
		m/=2;
		n/=2;
		i+=1;
	}
	if(m<n)     //m保存大的值
	{
		temp=m;
		m=n;
		n=temp;
	}
	while(x)
	{
		x=m-n;
		m=(n>x)?n:x;
		n=(n<x)?n:x;
		if(n==(m-n))
			break;
	}
	if(i==0)
	//	return n;
	printf("最大公约数是%d\n",n);
	else 
	//return (int )pow(2,i)*n;
	printf("最大公约数是%d\n",(int )pow(2,i)*n);


}


void main()
{
		clock_t start, finish; 
	double duration; 
	int i,j,x;
    int flag=1;
	int str[20]={90,177,83,45,42,35,66,54,89,54,64,96,785,64,63,543,35,345,56,54};
	start = clock(); 
	for( i=0;i<20;i++)
	gcd2(str[i],str[i+1]);
	finish = clock(); 
	duration = (double)(finish - start) / CLOCKS_PER_SEC; 
	printf( "所花费时间为%f seconds\n\n", duration);
}

更损相减法测试代码结果

1.4.Stein算法

#include <stdio.h>
#include <time.h>
#include <math.h>

 int Stein( unsigned int x, unsigned int y )//Stein算法 

  {
int factor = 0;// 当两数均为偶数时将其同时除以2至至少一数为奇数为止，
//记录除掉的所有公因数2的个数factor     
  		int temp;
        if ( x < y )//让 x为两者中大的数
        {
                temp = x;
                x = y;
                y = temp;
        }
        if ( 0 == y )
        {
                return 0;
        }
        while ( x != y )
        {
                if ( x & 0x1 )//判断x是否为奇数
                {
                        if ( y & 0x1 )//再判断y是否为奇数
                        {                           
    y = ( x - y ) >> 1;//将两数之差除于2
                                x -= y;
                        }
                        else//当x为奇数，y为偶数是when x is odd and y is even
                        {
                                y >>= 1;
                        }
                }
                else//当x为偶数
                {
                        if ( y & 0x1 )// 当x为偶数，y为奇数
                        {
                                x >>= 1;
                                if ( x < y )
                                {
                                        temp = x;
                                        x = y;
                                        y = temp;
                                }
                        }
                        else //当x和y都为偶数
                        {
                                x >>= 1;
                                y >>= 1;
                                ++factor;
                        }
                }
        }

		printf("最大公约数是%d\n",x << factor);
        
}

void main()
{
	l	clock_t start, finish; 
	double duration; 
	int i;
	int str[20]={90,177,83,45,42,35,66,54,89,54,64,96,785,64,63,543,35,345,56,54};
	start = clock(); 
	for( i=0;i<20;i++)
		Stein(str[i],str[i+1]);
	finish = clock(); 
	duration = (double)(finish - start) / CLOCKS_PER_SEC; 
	printf( "所花费时间为%f seconds\n\n", duration);
}

图4.stein算法测试代码截图
2.运行结果
总共用了4组不同的数据来做测试，求其平均时间。
2.1辗转相除法
在这里插入图片描述在这里插入图片描述
在这里插入图片描述
在这里插入图片描述
2.2 穷举法

在这里插入图片描述
在这里插入图片描述
在这里插入图片描述
在这里插入图片描述

平均时间为0.0495s

2.3更损相减法

![在这里插入图片描述![](https://img-blog.csdnimg.cn/20190309110746745.png?x-oss-pro![在这里插入图片描述](https://img-blog在这里插入图片描述
在这里插入图片描述
在这里插入图片描述
在这里插入图片描述
平均时间为0.032s

2.4 Stein算法
在这里插入图片描述在这里插入图片描述
在这里插入图片描述
在这里插入图片描述
平均时间为0.029s

五.总结
1.计算运行时间的方法

计算运行时间可以运用到clock()函数。clock函数定义如下： clock_t clock( void ); 这个函数返回从“开启这个程序进程”到“程序中调用clock()函数”时之间的CPU时钟计时单（clock tick）数。其中clock_t是用来保存时间的数据类型。

在time.h文件中，我们可以找到对它的定义：
#ifndef _CLOCK_T_DEFINED
typedef long clock_t;
#define _CLOCK_T_DEFINED
#endif

为了获取 CPU 所使用的秒数，仍需要除以 CLOCKS_PER_SEC。CLOCKS_PER_SEC 等于 1000000。在开头记录下时间start=clock(),再在末尾记录下时间end=clock(),两个时间相减即再除以CLOCKS_PER_SEC就为时间差。

2．运算符号的使用

2.1按位与(&)
a&x为对数a的二进制形式的取位操作，即去a二进制形式的第x位。在stein算法中a&1可以用于判断数a的奇偶性，即a末位为0即为偶数，末位为1的为奇数。

2.2按位左移(<<)
a<<=x等价于a*2^x，即为a乘以2的x次幂。原理是让a的二进制形式左移x位。这样与2的幂次方相乘使运算更快更方便。

2.3按位右移(>>)
a>>=x等价于a/(2^x),即使a除以2的x次幂。原理是让a的二进制形式右移x位。这样与2的幂次方相除使运算也更快。

Stein算法的优点之一就在于运用了移位运算操作，求最大公约数更加简洁和精确。

最爱吃芝士
关注
  • 12
    点赞
  • 踩
  • 31
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
C++求最大公约数四种方法解析
12-16
C++求最大公约数的四种方法思路，供大家参考，具体内容如下 将最近学的求最大公约数的四种方法总结如下： 第一种：穷举法之一 解释：拿其中一个数出来，用一个临时变量（tem）保存，每次都把那两个数除以这个临时...
浅析求两个数的最大公约数的三种算法
小葱的博客
11-05 4万+
第一种方法：#define _CRT_SECURE_NO_WARNINGS 1 #include<stdio.h> #include<stdlib.h> #include<assert.h>int Max_Com_Divisor(int x, int y)//用辗转相减法求最大公约数 { while (1)//用大数减去小数并将结果保存在大数中 { if (x > y
求最大公约数的4种算法（C++）
热门推荐
Elf.筱焱的博客
03-16 8万+
求最大公约数的4种算法（C++） 一、实验目的 1.计算两个正整数的最大公约数和最小公倍数，并进行程序的调式与测试。 2.理解四种不同的求最大公约数的方法，学习其思维模式。 3.了解算法的概念。对问题的分析时，进行合适的算法设计。 二、算法设计 1.题目分析 （1）首先输入一对正整数，对两个数进行判断，若不是正整数，则给出提示并重新输入，直至两个数均为正整数为止。 （2）将两个数传入计算最大公约数...
第十六讲 最小公倍数
最新发布
tang7mj的博客
06-09 672
最小公倍数（LCM，Least Common Multiple）是指两个或多个整数的所有公倍数中最小的一个。换句话说，最小公倍数是能够被这些整数整除的最小正整数。
辗转相除法、更相减损法、Stein算法
CSDN文章已停止维护，后续文章会在 https://blog.hufeifei.cn 持续更新
07-30 1万+
最大公约数和最小公倍数求解，常用的方法是短除法进行因式分解，然后最大公约数是所有公共因子的乘积，最小公倍数是所有因子的乘积。 本质上求最小公倍数就是求最大公倍数：x=m*a， y=m*b；m是最大公约数，那最小公倍数就是m*a*b。所以可以得到最大公约数与最小公倍数的关系： LCM(A,B)×GCD(A,B)=A×BLCM(A,B)×GCD(A,B)=A×B LCM(A,B)×GCD(A,...
如何求解最大公约数和最小公倍数
公众号（猴子数据分析）
08-24 1万+
1 定义 最大公约数 greatest common divisor，简写为gcd；或highestcommon factor，简写为hcf 如果有一个自然数a能被自然数b整除，则称a为b的倍数，b为a的约数。几个自然数公有的约数，叫做这几个自然数的公约数。公约数中最大的一个公约数，称为这几个自然数的最大公约数。 例： 在2、4、6中，2就是2，4，6的最大公约数。 最小公倍数 最小公
最大公约数的四种方法
Hell_potato777的博客
10-01 1万+
求两数的最大公约数，一共有四种方法：暴力穷举法、更相减损法、辗转相除法、stein 算法，小女不才，花了几天的时间终于把这几种方法全部弄明白，现在就把它们全部分享出来。首先，假设被求的两个数为 x、y，且 x > y。最大公约数 d = gcd (x , y)
欧几里德算法求最大公约数——C++代码
06-08
欧几里得算法，也称为辗转相除法，是计算两个正整数最大公约数（Greatest Common Divisor, GCD）的一种经典方法。该算法基于以下原理：两个正整数a和b（a>b）的最大公约数等于a除以b的余数c和b之间的最大公约数。...
求最大公约数的三种算法
09-21
今天一次无意的思考中想起了最大公约数，想一下最大公约数的算法，第一反映是穷举，再是欧几里得算法，然后是短除。
基础算法-python求最大公约数和最小公倍数
11-05
python求最大公约数和最小公倍数 #辗转相除法 def gcd(a,b): #最大公约数函数，且最小公倍数 = 两个数相乘 / 最大公约数 if b == 0: return a else: return gcd(b,a%b) print("请输入两个数：") j,k = input()....
算法_用欧几里得算法求最大公因数_
09-29
扩展欧几里得算法（Extended Euclidean Algorithm）不仅求出最大公因数，还能同时得到两个数的贝祖等式解，即存在整数x和y使得ax + by = gcd(a, b)。这对于计算模逆、线性同余方程的解等问题非常有用。 扩展...
求最大公约数的几种常见的方法 【详解】
加油！！！
10-19 8244
这里比较推荐是用 辗转相除法（欧几里得算法）和 《九章算术》中的 更相减损法多说一下，因为当时阿明浅浅学过遍辗转相除法，然后不久后就忘干净了，用的时候还要再去反复找，为了方便使用，干脆把 求解最大公约数 的几种常见的方法详细介绍一下，虽然不是最好，但是多少希望对大家有些帮助！希望大家对这些方法，有更加深刻的印象。加油！！！
一文搞定最大公约数（四种方法，赋原理和比较，超详细解答）
Stella-F的博客
06-20 1万+
最大公约数前言1.暴力穷举法代码2.辗转相除法步骤原理代码3.更相减损法步骤原理代码比较4.stein算法比较运算符&移位操作符原理步骤代码 前言 求两数的最大公约数，一共有四种方法：暴力穷举法、更相减损法、辗转相除法、stein 算法，小女不才，花了几天的时间终于把这几种方法全部弄明白（尤其是stein 算法，看了好久才终于看懂……），现在就把它们全部分享出来。 首先，假设被求的两个数为 x、y，且 x > y。最大公约数 d = gcd (x , y) 1.暴力穷举法 正如名字所说，暴击
最大公约数怎么求算法_欧几里得算法（辗转相除算法）2
weixin_39623671的博客
11-24 483
欧几里得95、欧几里得算法(辗转相除算法)2“上个视频，我们学习了如何用分解质因数法求最大公因数和最小公倍数。在使用这个方法时，需要先将每个数分解质因数。例如90=2×32(90=2×3的平方)，105=3×5×7…”网友说。“但并不是每个数都那么好分解…比如这两个数就很难分解：3139=？2117=？…”网友接着说。“遇到这种情况，分解质因数法也无能为力…这怎么办？”网友继续说。“不用担心，这种...
最大公约数的三种求法——（C语言）
m0_69573677的博客
12-06 2万+
最大公约数的三种求法——（C语言）
求最大公约数的几种算法
我的博客
09-12 9109
求两个数的最大公约数是生活中常见的运算。我们普遍知道的有欧几里得法、穷举法，其实还有其它方法：更相减损法。下面面一起来看看： 1、辗转相除法，又叫欧几里得法： int f1(int n1, int m1)//最大公约数子函数 { int n = n1 > m1 ? n1 : m1; int m = n1 < m1 ? n1 : m1; if (n % m == 0)//能直接相除的情况 { return m; } else { while (m) { i
求最大公约数的4种常用算法
AmethystFOB的博客
03-16 5万+
一、问题描述： 运行最大公约数的常用算法 二、问题分析与设计： 1.辗转相除法（又名欧几里德法） ①函数嵌套调用 其算法过程为： 前提：设两数为a,b设其中a 做被除数,b做除数，temp为余数 1、大数放a中、小数放b中； 2、求a/b的余数； 3、若temp=0则b为最大公约数； 4、如果temp!=0则把b的值给a、temp的值给a； 5、返回第二步； 代码： using nam...
四种算法求最大公约数
qq_45915957的博客
04-13 1992
1.辗转相除法 辗转相除法（又名欧几里德法）C语言中用于计算两个正整数a,b的最大公约数和最小公倍数，实质它依赖于下面的定理： 根据这一定理可以采用函数嵌套调用和递归调用形式进行求两个数的最大公约数和最小公倍数,现分别叙述如下： ①函数嵌套调用 其算法过程为： 前提：设两数为a,b设其中a 做被除数,b做除数，temp为余数 1、大数放a中、小数放b中； 2、求a/b的余数； 3、若temp=0...
三种算法求最大公约数c语言
a308285289的博客
03-22 586
#include #include int select=1; /*辗转相除法*/ void ZZXC() { int m,n; printf("请输入两个整数：\n"); scanf("%d,%d",&m,&n); printf("\n"); int t; if(m { m=m+n; n=m-n; m=m-n; } while(n!=0) { t=m
求最大公约数的欧几里德算法
05-28
欧几里德算法，也称辗转相除法，用于求两个正整数的最大公约数（GCD）。算法基于下面的定理： 对于任何整数a和b，它们的最大公约数等于a除以b的余数r和b之间的最大公约数。 算法步骤如下： 1. 如果a小于b，则交换a和b。 2. 用a除以b，得到余数r。 3. 如果r为0，则b是最大公约数。 4. 否则，将b赋值为r，然后返回第二步。 下面是一个Python实现： ```python def gcd(a, b): if a < b: a, b = b, a while b != 0: r = a % b a, b = b, r return a ``` 例如，gcd(60, 48)将返回12，因为60和48的最大公约数是12。

“相关推荐”对你有帮助么？

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值