欧几里得（辗转相除）算法与扩展欧几里得算法保姆级讲解

一.欧几里得算法（gcd）

我们先思考一个问题，如果我们现在要求两个数的最大公约数，你应该怎么求呢？我们这里用gcd表示为求两个最大公约数的函数，那么你不妨在草稿纸上写一下，gcd(a,b)是否等于gcd(b,a%b)呢？

这就是我们今天要讲的，其实我们叫这个算法的原理实际上就是辗转相除法。我们对于两个数a和b，如果我们要求a和b的最大公约数，当a>b时，我们用a除以b，得到结果q和余数r,再用除数b除以余数r再得到一个余数，再用余数除以余数，直到我们余数为零时，最后结果便为最大公余数。

看个代码：

#include<bits/stdc++.h>
using namespace std;
//板子
int gcd(int num1,int num2){
    int a=max(num1,num2);
    int b=min(num1,num2);
    if(a%b==0){
        return b;
    }
    return gcd(a%b,b);//如果余数不为零，就用辗转相除法一直求即可
}
int main(){
    int x=0;
    int y=0;
    while(cin>>x>>y){
        if(x==0&&y==0){
            break;
        }else{
            cout<<gcd(x,y)<<endl;//求x和y的最大公约数
        }
    }
    system ("pause");
    return 0;
}
  

二.扩展欧几里得算法（exgcd）

一.算法开场介绍

为什么我们大家用扩展欧几里得算法都命名为exgcd，实际上就是extend-gcd的意思。

我们先介绍一下裴蜀定理，若a，b是整数，且gcd(a,b)=d，那么对于任意的整数x，y，ax+by=m中的m一定是d的倍数。（特别地，如果a、b是整数，那么一定存在整数x、y使ax+by=gcd(a,b)。）

扩展欧几里得算法解决的是这样的一个问题：我们现在给定两个非零的整数a和b，求一组整数解（x，y）使得ax+by = gcd(a,b)成立，其中gcd(a,b)表示a和b的最大公约数。

欧几里得算法总是把gcd(a,b)转化为求解gcd(b,a%b)，而当b变为0时返回a，此时的a就等于gcd。也就是说，欧几里得算法结束时变量a中存放的是gcd，变量b中存放的是0，因此此时显然有a × 1 + b × 0 = g c d 成立，此时有x=1、y=0成立。那么我们不妨利用上面的欧几里得算法的过程来计算x和y。目前已知的是递归边界成立时为x=1、y=0，需要想办法反推出最初始的x和y。

从这里我们可以退出x1和y1，直到我们推出我们想要求的x和y即可。

void exgcd(int a, int b)
{
	//当前目的：求解 ax + by = gcd(a, b) 这么一个方程
	
	if(b == 0) //a, b不断改变的过程中，b最终必然会成为0
	{
		//在 b = 0 时方程还要成立？ 使 x = 1, y = 0 ，必然成立 
		x = 1;
		y = 0; //建议返回0。
		return;
	} 
	
	exgcd(b, a % b);//把下一层系数传进去（先求下一个方程的解 ）
	
	//现在我们已经拿到了下一个方程的解x, y
	int tx = x;//暂时存一下x，别丢了
	x = y;
	y = tx - a / b * y; 
}

二.实战

题目描述

求关于 x 的同余方程 ax≡1(modb) 的最小正整数解。

输入描述

一行，包含两个整数 a,b，用一个空格隔开。

输出描述

一个整数 x0​，即最小正整数解。输入数据保证一定有解。

样例输入 1 

3 10

样例输出 1 

7

#include<bits/stdc++.h>
using namespace std;
long long x,y;
void exgcd(int a,int b){
    if(b==0){
        x=1;
        y=0;
        return ;
    }
    exgcd(b,a%b);//实际上就是递归
    int temp=x;
    x=y;
    y=temp-(a/b)*y;
}
int main(){
    long long a,b;
    cin>>a>>b;
    exgcd(a,b);
    x=(x+b)%b;
    cout<<x<<endl;
    system("pause");
    return 0;
}

完结撒花，祝你所有题必然AC！

轻点关注，愿与你相遇与下一篇文章！