青少年编程与数学 02-018 C++ 数据结构与算法 17课题、数论算法-CSDN博客

本文链接：https://blog.csdn.net/qq_40071585/article/details/147617734

青少年编程与数学 02-018 C++ 数据结构与算法 17课题、数论算法

课题摘要:
数论是数学的一个分支，它研究整数的性质和关系。在计算机科学中，数论算法被广泛应用于密码学、编码理论、计算机安全等领域。

一、最大公约数（GCD）算法

最大公约数（GCD）是两个或多个整数共有的最大正整数因子。求解最大公约数的常用算法是欧几里得算法。

欧几里得算法

欧几里得算法基于以下原理：两个整数 (a) 和 (b) 的最大公约数等于 (b) 和 (a \mod b) 的最大公约数。

示例代码：

int gcd(int a, int b) {
    while (b != 0) {
        int temp = b;
        b = a % b;
        a = temp;
    }
    return a;
}

二、最小公倍数（LCM）算法

最小公倍数（LCM）是两个或多个整数共有的最小正整数倍数。求解最小公倍数可以通过最大公约数来计算。

最小公倍数算法

两个整数 (a) 和 (b) 的最小公倍数可以通过以下公式计算：(\text{LCM}(a, b) = \frac{|a \times b|}{\text{GCD}(a, b)})。

示例代码：

int lcm(int a, int b) {
    return abs(a * b) / gcd(a, b);
}

三、素数判断算法

素数是大于1的自然数，除了1和它本身外没有其他正因数。判断一个数是否为素数的常用方法是试除法。

试除法

试除法通过检查一个数 (n) 是否能被从2到 (\sqrt{n}) 的任何整数整除来判断其是否为素数。

示例代码：

#include <cmath>

bool is_prime(int n) {
    if (n <= 1) {
        return false;
    }
    if (n <= 3) {
        return true;
    }
    if (n % 2 == 0 || n % 3 == 0) {
        return false;
    }
    int i = 5;
    while (i * i <= n) {
        if (n % i == 0 || n % (i + 2) == 0) {
            return false;
        }
        i += 6;
    }
    return true;
}

四、素数生成算法

生成一定范围内的所有素数的常用方法是埃拉托斯特尼筛法。

埃拉托斯特尼筛法

埃拉托斯特尼筛法通过逐步标记合数来生成素数。具体步骤如下：

创建一个从2到 (n) 的整数列表。
从列表中的第一个数开始，标记其所有倍数为合数。
移动到下一个未标记的数，重复步骤2。
未标记的数即为素数。

示例代码：

#include <cmath>
#include <vector>

std::vector<int> sieve_of_eratosthenes(int n) {
    std::vector<bool> primes(n + 1, true);
    primes[0] = primes[1] = false;
    for (int i = 2; i <= std::sqrt(n); ++i) {
        if (primes[i]) {
            for (int j = i * i; j <= n; j += i) {
                primes[j] = false;
            }
        }
    }
    std::vector<int> result;
    for (int i = 2; i <= n; ++i) {
        if (primes[i]) {
            result.push_back(i);
        }
    }
    return result;
}

五、模运算相关算法

模运算在数论中非常重要，它涉及到整数除法的余数。以下是一些模运算相关的算法。

模幂运算

模幂运算计算 (a^b \mod m)。可以通过快速幂算法来高效计算。

示例代码：

int modular_exponentiation(int a, int b, int m) {
    int result = 1;
    a = a % m;
    while (b > 0) {
        if (b % 2 == 1) {
            result = (result * a) % m;
        }
        a = (a * a) % m;
        b = b / 2;
    }
    return result;
}

模逆运算

模逆运算找到一个数 (a) 在模 (m) 下的逆元，即找到一个数 (x) 使得 (a \times x \equiv 1 \mod m)。可以通过扩展欧几里得算法来计算。

示例代码：

int extended_gcd(int a, int b, int& x, int& y) {
    if (b == 0) {
        x = 1;
        y = 0;
        return a;
    }
    int x1, y1;
    int gcd = extended_gcd(b, a % b, x1, y1);
    x = y1;
    y = x1 - (a / b) * y1;
    return gcd;
}

int modular_inverse(int a, int m) {
    int x, y;
    int gcd = extended_gcd(a, m, x, y);
    if (gcd != 1) {
        return -1; // 不存在逆元
    }
    return (x % m + m) % m; // 确保结果为正数
}