扩展欧几里得算法

同余定理

(a+b)%c=(a%c+b%c)%c
(a*b)%c=(a%c*b%c)%c

同余定理的加法证明，乘法与之相似

a = k1 * m + r1
b = k2 * m + r2
( a + b )% m =(( k1 * m + r1 )+( k2 * m + r2 ))% m
= (( k1 + k2 )* m +( r1 + r2 ))% m
= ( r1 + r2 )% m
= ( a % m + b % m )% m
( a + b )% m = ( a % m + b % m )% m

什么是拓展欧几里得？简单的说，就是求关于x,y的方程 ax + by = gcd(a,b) 的所有整数解，注意这里是等号右边是gcd(a,b),所以做题的时候必须等号右边是gcd（a,b）或者gcd(a,b)的整数倍才能求解。（如果是整数倍n的话，先两边除以整数n求解方程，最后再乘回n，不能直接不除n求解方程）。

扩展欧几里得算法的作用是：求解a和b的最大公约数r和a*x+b*y的一个解。exgcd会返回a和b的最大公约数r, x和y是在递归的过程中求出来的，一开始x和y是不用赋初值的，结束递归的时候，x=1，y=0，在回溯的过程中x、y的值逐步改变，从而得到x、y的一组解。

如何用欧几里得算法求一个数的逆元呢？

由a*x ≡ 1(mod m)这个方程可以知道x是a的逆元。

a*x≡1(mod m) 等价于 a*x+m*y=1 可以用扩展欧几里得求得一组解，（x+m)mod m即为a的逆元。其实x就是a的逆元了，但是在完全剩余系m里，a的逆元有无数个，所以要使a的逆元x尽量小，让它模上m就行了。但是有可能它本来就比m小，一模m可能就等于0了，所以要提前加上m再模m，这就是a的逆元等于（x+m)mod m的原因。

至于y，随便一个整数就行了叭。

用扩展欧几里得算法求逆元的原因是a*x≡1(mod m) 可以轻易写成a*x+m*y=1 (mod m)的形式，如果把模m作为b代入，求得的值就是a的逆元，如果b的值不是m的话，求解的值只是方程a*x+b*y=gcd(a,b)的值了，而如果b不是m（或者m的倍数），那么这个值一定不是a的逆元。

参考文章：https://blog.csdn.net/sdutstudent/article/details/78795643