费马小定理和扩展欧几里得求逆元

什么是逆元：
如果有a*x = 1%p， 则称x是a的逆元

有公式
(a + b)%p = (a%p + n%p)%p
(a - b)%p = (a%p - n%p)%p
(a × b)%p = (a%p × n%p)%p

但没有
(a / b)%p = (a%p / n%p)%p
于是逆元的作用就体现出来了
求出b的逆元inv
就可以有
(a/b)%p = (a*inv)%p
运用乘法的公式，我们就可以在乘的过程中进行取模。

那么对于**(a / b)%p**， 我们又如何求的b的逆元，（这里介绍两种情况， 第一种情况基本可以适用大部分情况）
1.当p为质数
先看费马小定理： b ^ (p - 1) ≡ 1 (modp)。
左边提出一个b，得 b×b ^ (p - 2) ≡ 1 (modp)。
所以b的逆元就是b ^ (p - 2)。
然后我们用取模快速幂求逆元就行了

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

第二种情况：
b, p互质
根据逆元定义有b*inv = 1%p, 整理得binv + kp = gcd(b, p) = 1, 然后又可以用扩展欧几里得算法求inv了

int exgcd(int a, int b, int &x, int &y)
{
    if(!b)
    {
        x = 1;
        y = 0;
        return a;
    }

    int d = exgcd(b, a%b, y, x);
    y = y - a/b*x;
    return d;
}

这里可以稍微对扩展欧几里得的求法做一证明

ax1 + by1 = gcd(a, b)
不妨设gcd(a, b) = d
因为有gcd(b, a%b) == gcd(a, b) = d
于是得， a%bx2 +by2 = d
由a%b = a - a/b*b带入上式
经整理得
ax2 + (y2 - a/bx2) = d = ax1 + b*y1
x1 = x2, y1 = y2 - a/b*x2