【基础算法】数论的基本概念

简介

• 本文主要介绍数论的一些基本概念以及一些基本求解算法（如试除法求约数、欧几里得算法及其扩展等）。

整除

定义

• 设a，b属于整数，且a不等于0。如果存在整数c，使得b=ca，那么我们就称b被a整除，记为a|b。
• 上面的a是b的约数，b是a的倍数。

性质

试除法求约数

• 显然，对于整数n，它的约数必然小于根号n。
• 那么我们可以从1枚举到根号n，设当前枚举的为i，如果n%i=0则i是n的约数，特别地，如果i≠n/i，那么n/i也是n的约数。

code

    for (int i = 1; i <= x / i; i++) {
        if (x % i == 0) {
            approximate_number.push_back(i);
            if (i != x / i) {
                approximate_number.push_back(x / i);
            }
        }
    }

带余数除法

定义

• 设a、b为整数，且a≠0，那么必然存在等式b=ca+r。
• 如果没有特别说明，一般这个余数指的是最小非负余数。

性质

1. 任意整数被a除所得到的余数r必然落在区间[0, a - 1]上；
2. 任意连续a个整数被a除，必然会得到上述所说的a个数。

同余定义

• 设a、b为整数，如果存在整数p（p≠0），使得a mod p = b mod p，那么我们称a和b同余，也可以记为a 三 b (mod p)。

素数和合数

素数定义

• 素数指那些大于1的，且除了1和它自身之外再没有其它约数的自然数。

合数定义

• 合数是指除了1和它自身之外还有其它约数的自然数。

其他

• 特别地，1既不是素数又不是合数。

最大公约数和最小公倍数

最大公约数定义

最大公因数，也称最大公约数、最大公因子，指两个或多个整数共有约数中最大的一个。a，b的最大公约数记为（a，b），同样的，a，b，c的最大公约数记为（a，b，c），多个整数的最大公约数也有同样的记号。求最大公约数有多种方法，常见的有质因数分解法、短除法、辗转相除法、更相减损法。与最大公约数相对应的概念是最小公倍数，a，b的最小公倍数记为[a，b]。

最小公倍数定义

两个或多个整数公有的倍数叫做它们的公倍数，其中除0以外最小的一个公倍数就叫做这几个整数的最小公倍数。整数a，b的最小公倍数记为[a，b]，同样的，a，b，c的最小公倍数记为[a，b，c]，多个整数的最小公倍数也有同样的记号。

求解最大公约数

• 我们可以使用欧几里得算法求解这个问题

欧几里得算法

• Q：假设有两个整数a、b，如何求两个数的最大公约数呢？
• A：不妨假设a > b，此时有两种情况：
  1. b整除a，显然此时b就是a和b的最大公约数。
  2. b不整除a，即：a = q * b + r，r >= 0，我们通过证明可以得到gcd(a, b) = gcd(b, a % b)

  

  此时我们可以通过递归来求解a和b的最大公约数，递归停止的条件为b=0。
  根据上面的分析，我们可以写出以下求解代码。

  int gcd(int a,int b)
  {
      if(!b){
          return a;
      }
      return gcd(b,a%b);
  }
  

  也可以使用迭代的方法写代码

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

• 以上就是欧几里得算法，时间复杂度为O(n)。

求解多个数的最大公约数

• 对于序列a，我们每次求解t = gcd(a[i], a[i+1])，用t代替序列a中的a[i]和a[i+1]，重复多次，直到序列a中剩下一个数，这个数就是序列的最大公约数。

求解最小公倍数

扩展欧几里得算法

code

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;
}