拓展欧几里得算法及其应用 二元一次不定方程 线性同余方程 线性同余方程组 拓展中国剩余定理 模板

拓展欧几里得算法

欧几里得算法直接使用g++中的<algorithm>库中__gcd()函数即可。

$(a,b)=(b,a\mathrm{m}\mathrm{o}\mathrm{d}b)$.

拓展欧几里得算法用于求出不定方程$ax+by=(a,b)$的一个特解$x_0,y_0$,
顺带求出$(a,b)$,[通解$x=x_0+\frac{b}{(a,b)}t,y=y_0-\frac{a}{(a,b)}t(t\in Z)$]{}.

解不定方程

不定方程$ax+by=c$有解等价于$(a,b)\mid c$.据此判断是否有解，若有解，假设有一组特解$x_0^{{}^{\prime }},y_0^{{}^{\prime }}$,则它们的$\frac{c}{(a,b)}$倍显然是原不定方程的一组特解。$x_0=\frac{c}{(a,b)}{x}_0^{{}^{\prime }},y_0=\frac{c}{(a,b)}{y}_0^{{}^{\prime }}$,而通解依旧是$x=x_0+\frac{b}{(a,b)}t,y=y_0-\frac{a}{(a,b)}t(t\in Z)$.

求解模线性方程（线性同余方程）

$ax\equiv c(\mathrm{m}\mathrm{o}\mathrm{d}m)⟺ax+my=c$.

求乘法逆元

$ab\equiv 1(\mathrm{m}\mathrm{o}\mathrm{d}m)$,
则a关于模m的乘法逆元是b，b关于模m的乘法逆元是a。或者说$Z_m$群中a和b互为乘法逆元。

用乘法逆元有$\frac{A}{b}\equiv A\times b^{-1}(\mathrm{m}\mathrm{o}\mathrm{d}c)$.当左边的式子A是很大的数，而b是小规模数，且除出来的数一定是整数的时候，可以用右式边算边模。

求解
$ax\equiv 1(\mathrm{m}\mathrm{o}\mathrm{d}m)⟺ax+my=1$.解出的x即为解，只是注意需要用通解公式将$x$调整到$Z_m$范围内。

线性同余方程组

方程组 a i x ≡ c i ( m o d m i ) ( i = 1 , 2 , 3 , … , n ) a_ix \equiv c_i \pmod{m_i} \quad (i=1, 2, 3, \ldots, n)