乘法逆元详解（费马小定理与拓展欧几里得）

什么是乘法逆元？

对于整数a,m，如果满足gcd(a,m)=1，怎必定存在x使得,那么x称为a对m的乘法逆元。

乘法逆元的用途？

对于除法取模(b/a)%m，由于不存在公式(b/a)%m=(b%m)/(a%m)，所以必须先计算b/a再取模，但是当a很大是会造成b/a爆精度。乘法逆元的作用就是将这个除法取模变成乘法取模，有(b/a)%m=b*x%m=b%m*x%m，其中x为a对m的乘法逆元。

乘法逆元的计算方式

1.费马小定理（m必须为质数）

费马小定理：a^（p-1）≡1（mod p）其中p为质数。

对于a的逆元x有,同时根据飞马小定理又有a^(m-1)≡1(mod m），因此有ax=a^(m-1)，故x=a^(m-2)。利用快速幂可在logn时间内求出a对m的逆元。

2.扩展欧几里得

费马小定理只能用于m为质数时，所以必须寻求一般的方法。对于,可以等价为ax=km+1，那么有ax-km=1，我们令-m=y,则有ax+my=1。因此我们只要寻找到这个方程的解x,y其中x即为a对m的乘法逆元。该方程可用拓展欧几里得算法求解。拓展欧几里得算法如下:

基本算法：对于不完全为 0 的非负整数 a，m，gcd（a，m）表示 a，m 的最大公约数，必然存在整数对 x，y ，使得 gcd（a，m）=ax+my。

我们观察到：欧几里德算法停止的状态是： a= gcd ， b = 0 ，那么，这是否能给我们求解 x y 提供一种思路呢？因为，这时候，只要 a = gcd 的系数是 1 ，那么只要 b 的系数是 0 或者其他值（无所谓是多少，反正任何数乘以 0 都等于 0 但是a 的系数一定要是 1），这时，我们就会有： a*1 + m*0 = gcd

当然这是最终状态，但是我们是否可以从最终状态反推到最初的状态呢？

假设当前我们要处理的是求出 a 和 b的最大公约数，并求出 x 和 y 使得 a*x + m*y= gcd ，而我们已经求出了下一个状态：b 和 a%m 的最大公约数，并且求出了一组x1 和y1 使得： m*x1 + (a%m)*y1 = gcd ， 那么这两个相邻的状态之间是否存在一种关系呢？

我们知道： a%m= a - (a/m)*m（这里的 “/” 指的是整除，例如 5/2=2 , 1/3=0），那么，我们可以进一步得到：

        gcd = m*x1 + (a-(a/m)*m)*y1

               = m*x1 + a*y1 – (a/m)*m*y1

               = a*y1 + m*(x1 – a/m*y1)

    对比之前我们的状态：求一组 x 和 y 使得：a*x + b*y = gcd ，是否发现了什么？

    这里：

        x = y1

        y = x1 – a/m*y1

    以上就是扩展欧几里德算法的全部过程。算法实现如下

int exgcd(int a,int m,int &x,int &y)
{
    if(m==0)
    {
        x=1;
        y=0;
        return a;
    }
    int ans=exgcd(m,a%m,x,y);
    int t=x;
    x=y;
    y=t-a/m*y;
    return ans;
}