扩展欧几里得算法及其应用

前言

由于数论的板子真的很抽象,也很难背,所以特此记录扩展欧几里得算法的板子和它的用途 本篇文章只涉及应用,不涉及证明,如需理解证明还请各位移步其他优秀的讲解!

扩展欧几里得算法 

先粘一下板子的代码 typedef long long LL ; LL exgcd(LL a, LL b, LL &x, LL &y) { if (!b) { x = 1, y = 0 ; return a ; } LL d = exgcd(b, a % b, y, x) ; y -= a / b * x ; return d ; }
变量解释 对于方程:ax + by = d 其中 a 和 b 都是常数 (已知量),d 是 a 和 b 的最大公约数 x 和 y 是我们希望求得的一组满足方程的解

应用例题 

题目链接🔗:

题目分析 

AC代码

#include <iostream>
#include <algorithm>
#include <cstring>

using namespace std ;

typedef long long LL ; 

LL exgcd(LL a, LL b, LL &x, LL &y) 
{
    if (!b) 
    {
        x = 1, y = 0 ; 
        return a ; 
    }

    LL d = exgcd(b, a % b, y, x) ; 
    y -= a / b * x ; 

    return d ; 
}

int main() 
{
    ios::sync_with_stdio(false) ; 
    
    LL a, b, m, n, L ; 
    cin >> a >> b >> m >> n >> L ;

    LL x, y ; 
    LL d = exgcd(m - n, L, x, y) ; 

    if ((b - a) % d) cout << "Impossible" << endl ;
    else 
    {
        x *= (b - a) / d ; 
        LL t = abs(L / d) ; 
        cout << (x % t + t) % t << endl ; // 求最小正整数解
    }

    return 0 ; 
}

难点解释

为什么要计算 t ? 解释:

再来一道题目巩固一下

同余方程模版题

题目描述

题目分析 

a * x % b = 1 等价于找到两个数 x 和 y 使得 a * x + b * y = 1 这恰好是我们扩展欧几里得算法的基本解决对象,直接套板子就行了,由于题目保证输入一定有解,所以我们可以认为 a 和 b 是互质的,因此可以使用扩展欧几里得算法。 最后记得对b取模保证答案为最小正数。

AC代码 

#include <iostream>
#include <algorithm>
#include <cstring>

using namespace std ;

typedef long long LL ; 

int exgcd(int a, int b, int &x, int &y) 
{
    if (!b) 
    {
        x = 1, y = 0 ; 
        return a ; 
    }
    
    int d = exgcd(b, a % b, y, x) ; 
    y -= a / b * x ; 
    
    return d ; 
}

int main() 
{
    ios::sync_with_stdio(false) ; 
    int a, b ; 
    cin >> a >> b ; 
    
    int x, y ; 
    exgcd(a, b, x, y) ; 
    
    cout << (x % b + (LL)b) % b << endl ; 

    return 0 ; 
}

END

经验分享 程序员 微信小程序 职场和发展