乘法逆元

求解逆元有如下几种方式：

1. 如果 mod 为 素数，那么可以使用费马小定理+快速幂求解

2. 如果 gcd(a,mod)==1 其中，a为求逆元素，并且满足前面的等式保证a存在逆元。 那么可以使用 扩展欧几里得

3. 递推法快速打表求解多个逆元， 要求 mod 必须为奇质数

4. 如果 mod 不为素数，那么可以使用通用解法，这里不详细讲（因为我也不会

下面讲解1，2，3方法

方法一：（如果 mod 为 素数，那么可以使用费马小定理+快速幂求解）

1. 费马小定理

如果 gcd(a,p) == 1， 并且 p 为素数，那么就有 

公式变形：

所以：求解 a 的逆元，就是求解 。因此需要使用到快速幂

2. 代码：求 a 在 mod 下的乘法逆元

//求 m^n % k 
ll quick_power(ll m, ll n, ll k){
    ll ans = 1;
    m %= k;
    while(n){
        if(n & 1)
            ans = ans * m % k;
        n >>= 1;
        m = m * m % k;
    }
    return ans;
}

ll inverse_quick_pow(ll a, ll mod){
    return quick_power(a,mod-2,mod);
}

 

方法二：如果 gcd(a,mod)==1 其中，a为求逆元素，并且满足前面的等式保证a存在逆元。 那么可以使用 扩展欧几里得

1. 扩展欧几里得：传送门

2. 由于求解逆元，其实就是解方程 ax+py ≡ 1 (mod p)

3. 代码: 求 a 在 mod 下的乘法逆元

//ax + by = gcd(a,b)
//gcd 保存了 a,b 的最大公约数
void ext_gcd(ll a, ll b, ll &gcd, ll &x, ll &y) {
    if (!b) { gcd=a; x=1; y=0; }
    else { ext_gcd(b,a%b,gcd,y,x); y-=(a/b)*x; }
    // 理解上面这个这句话就用上面的 extend_gcd 来理解，主要是使用了引用的概念
    // 并且 ext_gcd 在 y-=(a/b)*x 上面，也就是说 ext_gcd 参数中的 x,y 都是
    // 下一层函数返回的 x2, y2
    // 然后使用 y-=(a/b)*x 处理一下本层函数 y1 即可
}


ll inverse_ext_gcd(ll a, ll mod) {
    ll x, y, gcd;
    ext_gcd(a, mod, gcd, x, y);
    // gcd(x,mod) ==1 是有逆元的充要条件
    // 当满足有逆元的时候，需要调整 x 到 0-mod 之间
    // 对于负数 x+mod 即可调整正确，
    // 为了编码统一将正数页进行调整， (x+mod)%mod
    return gcd == 1 ? (x + mod) % mod : -1;
}

 

 方法三：递推法快速打表求解多个逆元， 要求 mod 必须为奇质数

1. 原理：基于一个公式  

2. 初始条件： inv[1]=1;

3. 代码：

// 递推法快速求 逆元
const int maxn = 3e3 + 5;
ll inv[maxn];

//要求： mod 为奇质数
//公式: inv[i]=(mod - mod /i) * inv[mod%i]%mod
//初始条件: inv[1]=1;
void inverse_recursion(ll mod) {
    inv[1] = 1;
    for (int i = 2; i < maxn; i++)
        inv[i] = (mod - mod / i) * inv[mod % i] % mod;
}