【算法学习笔记】35:扩展欧几里得算法求解线性同余方程

线性同余方程问题

线程同余方程问题是指 a x ≡ b ( m o d m ) ax \equiv b~(mod~m) ax≡b (mod m),给定 a a a、 b b b和 m m m,找到一个整数 x x x使得该方程成立,即使得 a x m o d m = b ax~mod~m=b ax mod m=b,随便返回任何一个解都可以。

例如 4 x ≡ 3 ( m o d 5 ) 4x \equiv 3~(mod~5) 4x≡3 (mod 5),那么 x x x的一个可能的解可以是 2 2 2。

接下来用扩展欧几里得算法尝试构造这个解。从 a x ≡ b ( m o d m ) ax \equiv b~(mod~m) ax≡b (mod m)可知,一定存在一个 y y y使得:
a ⋅ x = m ⋅ y + b a \cdot x = m \cdot y + b a⋅x=m⋅y+b

也就是说,因为 a x ax ax模 m m m的余数是 b b b,所以 a x ax ax一定可以表示成 m m m的整数 y y y倍再加上一个 b b b。也就是:
a x − m y = b ax - my = b ax−my=b

令 y ′ = y y' = y y′=y,那么就是:
a x + m y ′ = b ax + my' = b ax+my′=b

因此原线性同余方程问题求 x x x有解,等价于这个方程求 x x x和 y ′ y' y′有解。而根据扩展欧几里得算法里所讨论的, a a a是 g c d ( a , m ) gcd(a,~m) gcd(a, m)的倍数, m m m也是 g c d ( a , m ) gcd(a,~m) gcd(a, m)的倍数,所以它们拼到一起也必须是 g c d ( a , m ) gcd(a,~m) gcd(a, m)的倍数。

因此,这个方程有解的充要条件 是 b b b必须是 g c d ( a , m ) gcd(a,~m) gcd(a, m)的倍数,也即 g c d ( a , m ) ∣ b gcd(a,~m)~|~b gcd(a, m) ∣ b 。

例题:AcWing 878. 线性同余方程

这题最终结果要限制在int范围内,因为 m m m也是在int范围内的,并且:
a x + m y = b ⇔ a ( k m + r ) + m y = b ⇔ a r + m ( a k + y ) = b ax + my =b \\ \Leftrightarrow a(km + r) + my = b \\ \Leftrightarrow ar + m(ak + y) = b ax+my=b⇔a(km+r)+my=b⇔ar+m(ak+y)=b

也就是说,把系数 x x x变成 r = x m o d m r = x~mod~m r=x mod m时,另一个系数只要从 y y y变成 a k + y ak+y ak+y就可以了,其中 k = ⌊ x m ⌋ k = \lfloor \frac{x}{m} \rfloor k=⌊mx⌋。

所以可以直接把结果 x x x模 m m m,一定也是一个合法的解,并且满足在int范围内的要求。

cpp 复制代码
#include <iostream>

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);
    // d = b * y + (a % b) * x = b * y + (a - a / b * b) * x
    //   = a * x + b * (y - a / b * x)
    y -= a / b * x;
    return d;
}

int main() {
    int t; cin >> t;
    
    while (t -- ) {
        int a, b, m; cin >> a >> b >> m;
        // ax % m = b, ax + my' = b, iff gcd(a, m) = d | b
        int x, y;
        int d = exgcd(a, m, x, y);
        if (b % d) puts("impossible");
        else cout << (LL)x * (b / d) % m << endl;
    }
    
    return 0;
}
相关推荐
胖咕噜的稞达鸭1 小时前
算法入门:专题攻克一---双指针(3)有效三角形的个数 查找总价格为目标值的两个商品(剑指offer题目)
算法
逻辑留白陈6 小时前
Adaboost进阶:与主流集成算法对比+工业级案例+未来方向
算法
Learn Beyond Limits6 小时前
Mean Normalization|均值归一化
人工智能·神经网络·算法·机器学习·均值算法·ai·吴恩达
天选之女wow6 小时前
【代码随想录算法训练营——Day28】贪心算法——134.加油站、135.分发糖果、860.柠檬水找零、406.根据身高重建队列
算法·leetcode·贪心算法
Gohldg6 小时前
C++算法·贪心例题讲解
c++·数学·算法·贪心算法
远远远远子6 小时前
类与对象 --1
开发语言·c++·算法
Aaplloo6 小时前
【无标题】
人工智能·算法·机器学习
西望云天7 小时前
The 2024 ICPC Asia Nanjing Regional Contest(2024南京区域赛EJKBG)
数据结构·算法·icpc
10岁的博客7 小时前
容器化安装新玩法
算法
不会算法的小灰7 小时前
HTML简单入门—— 基础标签与路径解析
前端·算法·html