【算法学习笔记】34:扩展欧几里得算法

裴蜀定理

描述

对于任意正整数 a a a、 b b b,一定存在整数系数 x x x, y y y,使得:
a x + b y = g c d ( a , b ) ax + by = gcd(a, b) ax+by=gcd(a,b)

并且 g c d ( a , b ) gcd(a, b) gcd(a,b)是对于任意的系数 x x x和 y y y放在 a a a和 b b b上能凑出的最小正整数。

证明

如下如果有整数系数 x x x、 y y y,那么 a x + b y ax + by ax+by一定是 g c d ( a , b ) gcd(a, b) gcd(a,b)的倍数,因为 a a a和 b b b都分别是 g c d ( a , b ) gcd(a, b) gcd(a,b)的倍数。因此能凑出来的数字最小就是 g c d ( a , b ) gcd(a, b) gcd(a,b)。

接下来证明一定存在 x x x、 y y y能凑出 g c d ( a , b ) gcd(a, b) gcd(a,b)。证明存在性的东西可以用构造法,只要把这个东西构造出来了,那么就一定存在了。这个构造的方法就是扩展欧几里得算法,使得对于任意的 a a a、 b b b,都能构造出来 x x x、 y y y使得 a x + b y = g c d ( a , b ) ax + by = gcd(a, b) ax+by=gcd(a,b)。

扩展欧几里得算法

欧几里得算法是 g c d ( a , b ) = g c d ( b , a m o d b ) gcd(a, b) = gcd(b, a~mod~b) gcd(a,b)=gcd(b,a mod b),递归出口是当 b b b为 0 0 0时返回 a a a,因为 0 0 0和任何数的最大公约数就是那个数。

扩展欧几里得算法则是在欧几里得算法的基础上,把系数 x x x和 y y y构造出来。

考虑到欧几里得算法的递归特性,可以用之前递归出的结果计算出前面的结果。

如,欧几里得算法计算gcd的递归出口是 b b b为 0 0 0时返回 a a a,此时 a a a和 b b b的最大公约数就是 a a a,要使得 a x + b y = a ax + by = a ax+by=a,显然有一组整数解 x = 1 x = 1 x=1且 y = 0 y = 0 y=0。

考虑除了递归出口之外的递归分支,需要计算 g c d ( b , a m o d b ) gcd(b, a~mod~b) gcd(b,a mod b),假设此时已经求出一组解 y y y和 x x x(这里将 y y y和 x x x调换便于计算,不调换也是可以的,结论是另一种写法公式),即使得:
b y + ( a m o d b ) x = g c d ( b , a m o d b ) = g c d ( a , b ) by + (a~mod~b)x = gcd(b, a~mod~b) = gcd(a, b) by+(a mod b)x=gcd(b,a mod b)=gcd(a,b)

考虑到 a m o d b = a − ⌊ a b ⌋ ⋅ b a~mod~b = a - \lfloor \frac{a}{b} \rfloor \cdot b a mod b=a−⌊ba⌋⋅b,整理一下 a a a和 b b b的系数,得到:
a x + ( y − ⌊ a b ⌋ ⋅ x ) b = g c d ( a , b ) ax + (y - \lfloor \frac{a}{b} \rfloor \cdot x) b = gcd(a, b) ax+(y−⌊ba⌋⋅x)b=gcd(a,b)

因此,从上一组解的 x x x不变, y y y减去 ⌊ a b ⌋ ⋅ x \lfloor \frac{a}{b} \rfloor \cdot x ⌊ba⌋⋅x就是为 a a a和 b b b构造的系数解。

cpp 复制代码
#include <iostream>

using namespace std;

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() {
    int t; cin >> t;
    
    while (t -- ) {
        int a, b; cin >> a >> b;
        int x, y;
        exgcd(a, b, x, y);
        cout << x << ' ' << y << endl;
    }
    
    return 0;
}
相关推荐
阳洞洞3 分钟前
leetcode 二分查找应用
算法·leetcode·二分查找
猎猎长风4 分钟前
【数据结构和算法】1. 数据结构和算法简介、二分搜索
数据结构·算法
Pasregret7 分钟前
模板方法模式:定义算法骨架的设计模式
算法·设计模式·模板方法模式
左灯右行的爱情41 分钟前
JVM-卡表
java·jvm·算法
奋斗者1号2 小时前
逻辑回归:使用 S 型函数进行概率预测
算法·机器学习·逻辑回归
CodeJourney.2 小时前
基于DeepSeek与Excel的动态图表构建:技术融合与实践应用
数据库·人工智能·算法·excel
Gerry_Liang2 小时前
LeetCode热题100——283. 移动零
数据结构·算法·leetcode
刘大猫262 小时前
Arthas sc(查看JVM已加载的类信息 )
人工智能·后端·算法
web安全工具库2 小时前
Python内存管理之隔代回收机制详解
java·jvm·算法
大学生亨亨3 小时前
蓝桥杯之递归
java·笔记·算法·蓝桥杯