乘法逆元的四种求法（拓展欧几里得、费马小定理、递归、递推）

前言

逆元：如果$a\ast x\equiv 1(modp)$，且a与p互质，则称x是a关于p的逆元。
对于这个概念和倒数有本质的区别，因为除法不能将mod数化进去。引用一个例子：

(a +  b) % p = (a%p +  b%p) %p  （对）

(a  -  b) % p = (a%p  -  b%p) %p  （对）

(a  *  b) % p = (a%p *  b%p) %p  （对）

(a  /  b) % p = (a%p  /  b%p) %p  （错）

为什么除法错的?
证明是对的难，证明错的只要举一个反例
(100/50)%20 = 2 ≠ (100%20) / (50%20) %20 = 0

对于一些题目，我们必须在中间过程中进行求余，否则数字太大，电脑存不下，那如果这个算式中出现除法，我们是不是对这个算式就无法计算了呢？

答案当然是 NO (>o<)
这时就需要逆元了

我们知道：
如果
ax = 1
那么x是a的倒数，x = 1/a
但是a如果不是1，那么x就是小数
那数论中，大部分情况都有求余，所以现在问题变了
ax = 1 (mod p)
那么x一定等于1/a吗？

不一定
所以这时候，我们就把x看成a的倒数，只不过加了一个求余条件，所以x叫做 a关于p的逆元
比如2 * 3 % 5 = 1，那么3就是2关于5的逆元，或者说2和3关于5互为逆元
这里3的效果是不是跟1/2的效果一样，所以才叫数论倒数

a的逆元，我们用inv(a)来表示

那么(a / b) % p = (a * inv(b) ) % p = (a % p * inv(b) % p) % p
这样就把除法，完全转换为乘法了。

以上的例子摘自这篇博客，强推，浅显易懂！

---

对于法三用到如下公式变形：

证明：
设x = p % a,y = p / a
于是有 x + y * a = p
(x + y * a) % p = 0
移项得 x % p = (-y) * a % p
x * inv(a) % p = (-y) % p
inv(a) = (p - y) * inv(x) % p
于是 inv(a) = (p - p / a) * inv(p % a) % p

法四就是法三的记忆化版或者说是递推版

Code：

#include <bits/stdc++.h>

using namespace std;
typedef long long LL;

//拓展欧几里得
LL exgcd(LL A, LL B, LL &x, LL &y) {
   
  if (B == 0) {
   
    x = 1;
    y = 0;
    return A;
  }
  LL gcdnum = exgcd(B, A % B, x, y);
  LL tmp = y;
  y = x - A/B*y; //算上来的x,y计算出新的x,y
  x = tmp;
  return gcdnum;
}

//中速乘  : 对于长整型的快速乘要利用到80为的long double的特性O(1)  本算法复杂度O(logN)
LL qmul(LL a, LL b, LL mod) {
   
  LL ans = 0;
  a %= mod;
  while (b) {
   
    if (b&1) {
   
      ans = (ans + a) % mod;
    }
    a = (a << 1) % mod;
    b >>= 1