同余与二元一次方程求解(下)

逆

同余方程 $ax\equiv 1(\mathrm{mod}m)$ 的解，为$a$模$m$的逆，记为$a^{-1}$。上一篇说到这个同余方程有解的条件是$gcd(a,m)=1$,以下说到的逆都是在这个条件下。
求逆可以用扩展欧几里得或者费马小定理(条件:$m$是素数)

逆元与除法取模

我们知道取模操作满足以下性质

(a+b)mod m = ((a mod m)+(b mod m))mod m
(a-b)mod m = ((a mod m)-(b mod m))mod m
(a*b)mod m = ((a mod m)*(b mod m))mod m

但是对除法取模进行类似的操作是错误的，例如
$(60/4)(\mathrm{mod}25)=15$
$((60(\mathrm{mod}25))/(4(\mathrm{mod}25)))(\mathrm{mod}25)=2$
这里就需要用到逆元，$4$ 模 $25$ 的逆是$19$
$(60\ast 19)(\mathrm{mod}25)=1140(\mathrm{mod}25)=15$

设$b$模$m$的逆是$d$，所以$b\ast d(\mathrm{mod}m)=1$
$(\frac{a}{b})(\mathrm{mod}m)=(\frac{a}{b})(\mathrm{mod}m)\ast (bd)(\mathrm{mod}m)=\left(\frac{a}{b}bd\right)(\mathrm{mod}m)=(ad)(\mathrm{mod}m)$
把除法取模转化成乘法取模

逆元与求解二元一次线性方程$ax+my=b(ax\equiv b(\mathrm{mod}m))$

得到逆之后就可以求形如$ax\equiv b(\mathrm{mod}m)$的任意同余方程
令$a^{\prime }$为$a$模$m$的逆,$a^{\prime }a\equiv 1(\mathrm{mod}m)$
在$ax\equiv b(\mathrm{mod}m)$的左右同乘以$a^{\prime }$,得到$a{a}^{\prime }x=a^{\prime }b(\mathrm{mod}m),x\equiv a^{\prime }b(\mathrm{mod}m)$

例如求$6x\equiv 40(\mathrm{mod}23)$$,6$模$23$的逆是$4$,所以$x=4\ast 40(\mathrm{mod}23)=22$

转一个代码注释

/*
--------------------佛祖保佑，永无bug--------------------
	                   _ooOoo_
	                  o8888888o
	                  88" . "88
	                  (| -_- |)
	                  O\  =  /O
	               ____/`---'\____
	             .'  \\|     |//  `.
	            /  \\|||  :  |||//  \
	           /  _||||| -:- |||||-  \
	           |   | \\\  -  /// |   |
	           | \_|  ''\-/''  |   |
	           \  .-\__  `-`  ___/-. /
	         ___`. .'  /-.-\  `. . __
	      ."" '<  `.___\_<|>_/___.'  >'"".
	     | | :  `- \`.;`\ _ /`;.`/ - ` : | |
	     \  \ `-.   \_ __\ /__ _/   .-` /  /
	======`-.____`-.___\_____/___.-`____.-'======
	                   `=-='
*/