扩展中国剩余定理

#define _(x,y) (x%y+y)%y

Exgcd（扩展欧几里得算法）

为了求解方程 $ax+by=c$ 的一组解

我们可以先使用扩展欧几里得算法求解方程 $ax+by=\gcd (a,b)$

不妨令 $a_2=b_1,b_2=a_1\mathrm{m}\mathrm{o}\mathrm{d}{b}_1$

对于方程：
$$\{\begin{array}{l}{a}_1{x}_1+b_1{y}_1=\gcd (a_1,b_1)\\ a_2{x}_2+b_2{y}_2=\gcd (a_2,b_2)\end{array}$$

因为

$$\gcd (a_1,b_1)=\gcd (b_1,a_1\mathrm{m}\mathrm{o}\mathrm{d}{b}_1)=\gcd (a_2,b_2)$$

联立可得：

$$a_1{x}_1+b_1{y}_1=b_1{x}_2+(a_1\mathrm{m}\mathrm{o}\mathrm{d}{b}_1)\cdot y_2$$

$$a_1{x}_1+b_1{y}_1=b_1{x}_2+(a_1-\lfloor \frac{a_1}{b_1}\rfloor \cdot b_1)\cdot y_2$$
$$a_1{x}_1+b_1{y}_1=b_1{x}_2+a_1{y}_2-\lfloor \frac{a_1}{b_1}\rfloor \cdot b_1{y}_2$$
$$a_1{x}_1+b_1{y}_1=a_1{y}_2+b_1\cdot (x_2-\lfloor \frac{a_1}{b_1}\rfloor \cdot y_2)$$

显然 $x_1=y_2,y_1=x_2-\lfloor \frac{a_1}{b_1}\rfloor \cdot y_2$ 时原方程成立

因此，只要求出 $x_2$ 和 $y_2$ ，则可求出 $x_1$ 和 $y_1$

而在函数中，$x_1$ 的返回值为 $y_2$ ，$y_1$ 的返回值为 $x_2$

为了调整当前返回值，$y_1=x_2-\lfloor \frac{a_1}{b_1}\rfloor \cdot y_2=y_1-\lfloor \frac{a_1}{b_1}\rfloor \cdot x_1$

（联系返回值来理解，当然你也可以让 $x_1$ 的返回值为 $x_2$ ，$y_1$ 的返回值为 $y_2$ 但那个亚子就没压行好看了

显然，当 $b=0$ 时，原方程转换为 $ax=a$ ，一定存在一组解
$$\{\begin{array}{l}x=1\\ y=0\end{array}$$

代码如下：

LL exgcd(LL a,LL b,LL &x,LL &y)
{
    if(!b) return x=1,y=0,a;
    LL Q=exgcd(b,a%b,y,x);
    return y-=a/b*x,Q;
}

Excrt （扩展中国剩余定理）

为了求解同余方程组的一个解
$$\{\begin{array}{l}x\equiv a_1(\mathrm{m}\mathrm{o}\mathrm{d}{m}_1)\\ x\equiv a_2(\mathrm{m}\mathrm{o}\mathrm{d}{m}_2)\\ ...\\ x\equiv a_n(\mathrm{m}\mathrm{o}\mathrm{d}{m}_n)\end{array}$$

我们先将前 $i-1$ 个同余方程合并成一个船新的同余方程 $x\equiv A(\mathrm{m}\mathrm{o}\mathrm{d}M)$，且这个同余方程满足前 $i-1$ 个同余方程的约束

现在我们来考虑怎么合并第 $i$ 个同余方程

对于同余方程 $x\equiv A(\mathrm{m}\mathrm{o}\mathrm{d}M)$ ， $x$ 的通解为 $x=A+M\cdot k$

将其代入第 $i$ 个同余方程得 $A+M\cdot k\equiv a_i \pmod{m_i} $

化简得 $x\cdot M+y\cdot m_i=A-a_i$ ，其中 $x,y$ 为未知数

使用扩展欧几里得算法可以解得方程 $x\cdot M+y\cdot m_i=\gcd (M,m_i)$ 的一组特解

记 $x_0=\gcd (M,m_i)$ ，$C=A-a_i$ （注意调整 $C$ 的值 $[0,m_i)$

若 $x_0|C$ 则原方程有解，此时应将 $x$ 调整为原方程的解，即 $x=x\cdot \frac{C}{x_0}$

此时 $x$ 满足前 $i$ 组同余方程

因为 $x\cdot M\equiv A-a_i(\mathrm{m}\mathrm{o}\mathrm{d}{m}_i)$ ，所以 $A-x\cdot M\equiv a_i(\mathrm{m}\mathrm{o}\mathrm{d}{m}_i)$

令 $A^{{}^{\prime }}=A-x\cdot M,M^{{}^{\prime }}=\frac{M\cdot m_i}{\gcd (M,m_i)}$

因为 $A^{{}^{\prime }}\equiv A(\mathrm{m}\mathrm{o}\mathrm{d}M)$ 且 $A^{{}^{\prime }}\equiv a_i(\mathrm{m}\mathrm{o}\mathrm{d}{m}_i)$

则两个同余方程可合并为 $x\equiv A^{{}^{\prime }}(\mathrm{m}\mathrm{o}\mathrm{d}{M}^{{}^{\prime }})$

代码如下：

LL excrt()
{
    LL M=m[1],A=a[1];
    for(int i=2;i<=n;i++)
    {
        LL x,y,x0=exgcd(M,m[i],x,y);
        LL C=_(A-a[i],m[i]);
        if(C%x0) return -1;
        x=ksc(x,C/x0,m[i]);
        x%=m[i];A-=x*M;M*=m[i]/x0;A%=M;
   }
    if(A%M==0) return M;
    return _(A,M);
}