8.12.12 ACM-ICPC数学 数论 线性同余方程

8.12.11 ACM-ICPC数学 数论 线性同余方程

引言

线性同余方程在数论中占有重要地位，尤其在密码学、计算机科学和算法竞赛（如ACM-ICPC）中广泛应用。本节将详细讨论线性同余方程的定义、解法以及在实际问题中的应用，并包含代码实现示例。

线性同余方程的定义

线性同余方程的一般形式为： 𝑎𝑥≡𝑏(mod𝑛)ax≡b(modn) 其中，𝑎a、𝑏b 和 𝑛n 为整数，𝑥x 为未知数。需要从区间 [0,𝑛−1][0,n−1] 中求解 𝑥x，当解不唯一时需要求出全体解。

用逆元求解

当 𝑎a 和 𝑛n 互素（gcd⁡(𝑎,𝑛)=1gcd(a,n)=1）时，可以计算 𝑎a 的逆元，并将方程的两边乘以 𝑎a 的逆元，得到唯一解。此时：

𝑥≡𝑏𝑎−1(mod𝑛)x≡ba−1(modn)

扩展欧几里得算法求解

当 gcd⁡(𝑎,𝑛)≠1gcd(a,n)=1 时，方程可能无解或有多个解。具体步骤如下：

1. 计算 gcd⁡(𝑎,𝑛)gcd(a,n)。

2. 如果 gcd⁡(𝑎,𝑛)∣𝑏gcd(a,n)∣b，则方程有解。否则无解。

3. 将方程两边除以 gcd⁡(𝑎,𝑛)gcd(a,n)： 𝑎′𝑥≡𝑏′(mod𝑛′)a′x≡b′(modn′) 其中 𝑎′=𝑎𝑑a′=da​、𝑏′=𝑏𝑑b′=db​ 和 𝑛′=𝑛𝑑n′=dn​。

4. 使用扩展欧几里得算法求解新的方程。

实现代码

扩展欧几里得算法

#include <iostream>
using namespace std;

// 扩展欧几里得算法
int ex_gcd(int a, int b, int& x, int& y) {
    if (b == 0) {
        x = 1;
        y = 0;
        return a;
    }
    int d = ex_gcd(b, a % b, x, y);
    int temp = x;
    x = y;
    y = temp - a / b * y;
    return d;
}

// 线性同余方程求解
bool linearCongruence(int a, int b, int n, int& x, int& y) {
    int d = ex_gcd(a, n, x, y);
    if (b % d != 0) return false;
    int k = b / d;
    x = (x * k % n + n) % n;  // 保证 x 为非负数
    return true;
}

int main() {
    int a = 14, b = 30, n = 100;
    int x, y;
    if (linearCongruence(a, b, n, x, y)) {
        cout << "Solution: x = " << x << endl;
    } else {
        cout << "No solution." << endl;
    }
    return 0;
}

实例讲解

例题

求解线性同余方程： 14𝑥≡30(mod100)14x≡30(mod100)

解法

1. 计算 gcd⁡(14,100)=2gcd(14,100)=2，由于 gcd⁡(14,100)gcd(14,100) 能整除 30，方程有解。

2. 化简方程： 7𝑥≡15(mod50)7x≡15(mod50)

3. 使用扩展欧几里得算法求解 7𝑥≡15(mod50)7x≡15(mod50)： 计算 77 和 5050 的逆元。

4. 根据计算结果，方程的一个解为 𝑥=45x=45。

5. 所有解为： 𝑥≡45+𝑘⋅50x≡45+k⋅50 其中，𝑘k 为任意整数。

应用

线性同余方程在许多实际问题中都有应用，例如：

1. 密码学：线性同余生成器用于伪随机数生成。
2. 计算机科学：哈希函数的设计。
3. 算法竞赛：常用于解决模运算问题，如求逆元、快速幂等。

总结

本节介绍了线性同余方程的基本定义、解的存在性条件、具体求解方法以及实际应用。掌握这些内容，对于理解数论中的许多问题和算法竞赛中的高效解题具有重要意义。并通过代码示例展示了如何实现这些算法，帮助读者更好地理解和应用线性同余方程。