欧几里得法求最大公约数

欧几里得法求最大公约数:
求a和b的最大公约数
记 a mod b=c ,即a=kb+c 
设a b的最大公约数为d,则a=m*d b=n*d,m和n互质。
c=a-kb=md-knd=(m-kn)d,m和n互质,则n和m-kn互质,c和b的最大公约数也是d
所以: "a和b(a>b)的最大公约数等于b和a Mode b的最大公约数",递归或迭代计算,直到余数为0,此时除数为最大公约数


也可以这样理解:
a=m*d b=n*d,m和n互质,c=a-b=(m-n)*d ,m-n 和 n互质,

即,"a 、b(a>b)的最大公约数等于 b 和 a-b 的最大公约数。"递归或迭代计算,直到两数相等,此时的值为最大公约数。


/*  欧几里得法 求最大公约数
 *
 */
#include <iostream>
#include <stdlib.h>
using namespace std;
//m>n
int gcd1(int m,int n)
{
	if(m<0)
		m=0-m;
	if(n<0)
		n=0-n;
	for(int c=m%n;c!=0;)
	{
		m=n;
		n=c;
		c=m%n;
	}
	return n;
}
// m>2
int gcd2(int m,int n)
{
	if(m<0)
		m=0-m;
	if(n<0)
		n=0-n;
	while(m!=n)
	{
		m=m-n;
		if(m<n)
		{
			int c=m;
			m=n;
			n=c;
		}
	}	
	return m;
}

int main()
{
	for (int i=0; i<10;i++)
	{
		int n=rand()%10+1;
		int m=rand()%10+10+1;
		cout<<m<<" gcd1 "<<n<<"="<<gcd1(m,n)<<endl;
		cout<<m<<" gcd2 "<<n<<"="<<gcd2(m,n)<<endl;

	}	
}


评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值