同余方程&孙子定理&模线性同余方程

参考博客

Acdreamer

---

一次同余方程

$$a * x \equiv b (mod\ n)$$

分情况讨论

---

1. 若有$(a,n) = 1$,根据 欧几里得扩展，存在 l,m使得$a*l+n*m = 1$，等式两边同时模n，得 $a*l \equiv 1 (mod \ n )$,对于 $$a * x \equiv b (mod\ n)$$ $$a*l*x \equiv b*l (mod \ n)$$ $$x \equiv b*l (mod\ n)$$

---

1. 若有 $(a,n) = d$, 对于 $$a * x \equiv b (mod\ n)$$

移项

$$a * x-b \equiv 0 (mod\ n)$$
即 $$n|(a*x-b)$$
已知 $$d|a,d|n$$
于是 $$d|b$$
所以同余方程有解的条件是 $d|b$
令 $a_1 = a/d,b_1 = b/d,n_1 = n/d$
则有 $$n_1 |(a_1*x-b_1)$$
即为 $$a_1*x \equiv b_1(mod \ n_1)\ \ \ \ \ 且 (a_1,n_1) = 1$$
就变成情况1，可以解出此时的 $$x \equiv b_1*l (mod\ n_1)$$
即 $$x = b_1*l + n_1*t$$ t的取值范围 为(0,d-1)

---

中国剩余定理

解线性同余方程，运用拉格朗日插值法，得孙子定理

模板

int m[maxn];
int  China(int a[],int n,int b[])//其中a[]是模数，b[] 是余数
{
    int M = 1;
    for(int i = 0;i < n;++i)
        M *= a[i];
    int ans = 0;
    for(int i = 0;i < n;++i)
    {
        int Mi = M/a[i];
        LL x,y;
        ex_gcd(Mi,a[i],x,y);
        ans += b[i] * x*Mi;
    }
    ans = (ans%M + M) %M;
    return ans;
}

扩展

合并两个模数
假设
         $y= a_1 (mod \ m_1)$
         $y =a_2 (mod \ m_2)$
令
         $x \equiv a_1 + t_1*m_1; （1）$
         $x \equiv a_2 + t_2*m_2; （ 2）$
两式相减得         $t_1*m_1- t_2 *m_2 = a_1 - a_2（3）$
        利用欧几里得扩展求出 $t_1$带入（1） 求出一个特解$x$

$$y = x (mod \ lcm(m_1,m_2))$$
证明 x与y同余于$lcm（m_1,m_2)$
$m_1|x-y$
$m_2|x-y$
所以 $lcm(m_1,m_2)|x-y$,所以x和y关于$lcm(m_1,m_2)$ 同余