初等数论总结(填坑)

1.欧几里得算法(GCD)

---

1.1 普通数的GCD

• 表达式

  $$gcd\left ( a,b\right ) = gcd\left (b, a\%b\right )$$

• 代码实现

  //非递归
  int gcd(int a, int b) {
      for (int c = a%b; c != 0; c = a%b) a = b, b = c;
      return b;
  }
  

• 复杂度：$O(log\left (n\right ))$

---

1.2 高精度的GCD

• 表达式

  $$gcd\left (a, b\right ) = gcd\left (a-b, b\right )\\ gcd\left (2\cdot a, b\right ) = gcd\left (a, b\right )\\ gcd\left (2\cdot a, 2\cdot b\right ) = gcd\left (a, b\right )\\ $$

• 复杂度：$O(log\left (a\right ) + log\left (b\right ))$

---

2.扩展欧几里得算法(exGCD)

---

• 方程形式：

  $$aX + bY = gcd\left (a, b\right )$$

• 求解方法：

  设
  

  $$a\cdot X_{1} + b\cdot Y_{1} = gcd\left ( a,b\right )\\ b\cdot X_2 + (a \% b)\cdot Y_2 = gcd\left(b, a \% b \right )$$

  其中
  

  $$gcd\left (a, b\right ) = gcd\left (b, a \% b\right)\\ a \% b = a-(a/b)·b$$

  整理可得
  

  $$a·X_1 + b·Y_1 = b·X_2 + (a-(a / b)·b)·Y_2\\ a·X_1 + b·Y_1 = a·Y_2 + b·[X_2-(a / b)]·Y_2$$

  所以
  

  $$X_1 = Y_2\\ Y_1 = X_2-(a / b)·Y_2$$

• 代码实现：

  //递归
  void exGcd(int a, int b, int &x, int &y) {
      if (b == 0) { x = 1; y = 0; return; }
      exGcd(b, a%b, x, y);
      int t = x; x = y; y = t-a/b*y;
  }

---

3.裴蜀定理

---

• 推论：

  设 $gcd\left (a_{1}, a_{2}, a_{3}...a_{n}\right) = d$

  存在 $X_{1}, X_{2}, X_{3}...X_{n}$

  使得 $\sum_{i = 1}^{n} a_{i}\cdot X_{i} = d$

• 求解方法：

  先求 $a_{1}\cdot X_{1} + a_{2}\cdot X_{2} = gcd\left (a_{1}, a_{2}\right )$

  再将 $gcd\left (a_{1}, a_{2}\right )$ 带回原式，反复迭代就可得解

---

4.数论筛法

---

4.1 Eraosthens筛法

• 代码实现

  void getPrime(const int n, int *p) {
      int tot = 0;
      memset(vis, 0, sizeof(vis));
      for (int i = 2; i <= n; i++) {
          if (vis[i]) continue;
          p[++tot] = i;
          for (int j = 2; i*j <= n; j++) vis[i*j] = 1;
      }
  }
  

• 复杂度：$O(n\ log\ n)$

---

4.2 Euler筛法：

• 详情请百度 贾志鹏_线性筛法与积性函数.

• 代码实现

  void getPrime(const int n, int *p) {
      int tot = 0;
      memset(vis, 0, sizeof(vis));
      for (int i = 2; i <= n; i++) {
          if (!vis[i]) p[++tot] = i;
          for (int j = 1; j <= tot; j++) {
              if (i*p[j] > n) break;
              vis[i*p[j]] = 1;
              if (i%p[j] == 0) break;
          }
      }
  }
  

• 复杂度：$O(n)$

---

5.素数测试

---

5.1 Miller-Rabin算法

• 填坑

---

5.2 二次探测

• 我都不会

---

6.质因数分解

---

6.1 暴力分解法

• 代码实现自行脑补

• 复杂度：$O(sqrt\left (n\right ))$

---

6.2 Pollard-Rho 椭圆曲线分解法

• 连怎么读我都不会

---

7.快速幂

---

• 代码实现

  int pow(int a, int b, int p) {
      int ret = 1;
      while (b) {
          if (b & 1) ret = (ret*a)%p;
          a = (a*a)%p; b >>= 1;
      } return ret;
  }
  

• 复杂度：$O(log\left (b\right ))$

---

8.积性函数

---

8.1 定义(引用)

在非数论的领域，积性函数指所有对于任何a,b都有性质f(ab)=f(a)f(b)的函数。
在数论中的积性函数：对于正整数n的一个算术函数 f(n)，若f(1)=1，且当a,b互质时f(ab)=f(a)f(b)，在数论上就称它为积性函数。
若对于某积性函数 f(n) ，就算a, b不互质，也有f(ab)=f(a)f(b)，则称它为完全积性的。

8.2 重要的我知道的积性函数

• $\varphi\left (n\right )$_欧拉函数

  • 定义

    $\varphi\left (n\right )=$ 比 $n$ 小且与 $n$ 互质的正整数的数目.

  • 计算

    详情请见 wikipedia

    $$\varphi \left ( n \right ) = \prod_{i = 1}^{r}p_{i}^{k_{i}-1}(p_{i}-1) = n\prod_{p|n} (1-\frac{1}{p})$$

• $\mu\left (n\right )$_莫比乌斯函数

  • 定义 智商捉急只能打成这样

    $$\mu\left (n\right ) = \left\{\begin{matrix} 1 & n=1\\ (-1)^{k} & n=p_{1}p_{2}...p_{k}\\ 0 & other \end{matrix}\right.$$

  • 简单来说就是如果 $n=1$ 则为 $1$ ，如果有平方因子就为 $0$ ，否则就是 $(-1)^{k}$ ， $k$ 是因数的个数。

8.3 几个等式只知道一个

$$\sum_{d|n} \varphi\left (d\right )=n$$

---