本博客主要讲解如何求一个数的乘法逆元。
定义
对于任意整数 a , m , b a,m,b a,m,b,若 a , m a,m a,m互质,且 a ∣ b a|b a∣b ,则存在一个整数 x x x 使得 b / a ≡ b × x ( m o d m ) b/a\equiv b\times x(mod\ m) b/a≡b×x(mod m),则称 x x x 是 a a a 的模 m m m 乘法逆元,记为 a − 1 ( m o d m ) a^{-1}(mod\ m) a−1(mod m)。
求解
方法1:费马小定理
因为 b / a ≡ b × a − 1 ≡ b / a × b × b − 1 ( m o d m ) b/a\equiv b\times a^{-1}\equiv b/a \times b\times b^{-1}(mod\ m) b/a≡b×a−1≡b/a×b×b−1(mod m),所以 a × a − 1 ≡ 1 ( m o d m ) a\times a^{-1}\equiv 1(mod\ m) a×a−1≡1(mod m)。
当 m m m 为质数(此时用 p p p 代替 m m m)且 a < p a<p a<p,依据费马小定理, a p − 1 ≡ 1 ( m o d p ) a^{p-1}\equiv 1(mod\ p) ap−1≡1(mod p),即 a × a p − 2 ≡ 1 ( m o d p ) a\times a^{p - 2} \equiv 1(mod\ p) a×ap−2≡1(mod p)。因此,当模数 p p p 为质数时, a p − 2 a^{p - 2} ap−2 即为 a a a 的乘法逆元。
方法2:扩展欧几里得算法
当仅保证 a , m a,m a,m互质时,我们可以通过求线性同余方程 a × x ≡ 1 ( m o d m ) a\times x\equiv1(mod\ m) a×x≡1(mod m) 来得到所求的逆元。这一步很容易用扩展欧几里得算法来写。
ps
关于组合计数中所用的逆元,当其模数为质数时可以由另一种方法推出来。
首先,我们可以递推求出 k ! ( k ∈ [ 1 , n ] ) k!(k\in [1,n]) k!(k∈[1,n]),对于 n ! n! n! 我们可以用费马小定理来求出它的逆元。然后,我们反过来递推:由于 n ! n! n!的逆元相当模上 m m m 后的一个分母,所以我们可以由 ( n − 1 ) ! − 1 ( m o d p ) = n − 1 ( m o d p ) ] × n (n-1)!^{-1}(mod\ p) = n^{-1}(mod\ p)]\times n (n−1)!−1(mod p)=n−1(mod p)]×n来求出 ( n − 1 ) − 1 ( m o d p ) (n-1)^{-1}(mod\ p) (n−1)−1(mod p)从而求出所求的组合数。