扩展中国剩余定理

概述

  扩展中国剩余定理(exCRT)是中国剩余定理(CRT)更普适的版本。

   exCRT 没有 CRT 中的模数互质的限制，这一变化虽然导致可能出现无解的情况，但却令该算法能解决更为普通的线性同余方程组。此外，不一定互质也说明可能无法用 CRT 的方法解 exCRT 问题，所以我们需要考虑新的算法。

求解

给定方程组
$$\{\begin{array}{l}x\equiv a_1(mod{m}_1)\\ x\equiv a_2(mod{m}_2)\\ x\equiv a_3(mod{m}_3)\\ ⋮\\ x\equiv a_n(mod{m}_n)\end{array}$$
求 $x$ 。

  本着难题化简的想法，我们先考虑如何解仅含两个同余式的同余方程组：

  对于方程组

$$\{\begin{array}{l}x\equiv a_1(mod{m}_1)\\ x\equiv a_2(mod{m}_2)\end{array}$$

  首先，依据同余的定义，我们可以将原本的两个线性同余方程式转换为普通的线性方程式：

$$\{\begin{array}{l}x=a_1+k_1\times m_1\\ x=a_2+k_2\times m_2\end{array}(k_1\in \mathbb{Z},k_2\in \mathbb{Z})$$

  联立可得 $a_1+k_1\times m_1=a_2+k_2\times m_2$ ，再此基础上移项得 $k_1\times m_1-k_2\times m_2=a_2-a_1$ 。现在，设 $A=m_1,B=m_2,C=a_2-a_1$ ，则该式可转换为 $A\times k_1+B\times k_2=C$ 。

  可以发现这个式子略感熟悉。回忆一个知识点：扩展欧几里得算法求线性方程。

  可以发现，该式为扩展欧几里得的标准求解方程。所以我们可以利用扩展欧几里得算法求解 $k_1$ ，然后代入求解，得出一组特解 $x_0$ 。可以发现对于这两个同余式来说存在通解式 $x=x_0+lcm(m_1,m_2)\times k(k\in \mathbb{Z})$ 。同时，该通解式可以转换为 $x\equiv x_0(modlcm(m_1,m_2))$ 。而这个同余式与之前的 $x\equiv a_i(mod{m}_i)$ 没有本质区别。所以可以用这个式子代替

$\{\begin{array}{l}x\equiv a_1(mod{m}_1)\\ x\equiv a_2(mod{m}_2)\end{array}$

  这样的话就将两个同余式合并为了一个。可以直接用扩展欧几里得算法解得线性同余方程式。

  由于总共有 $n$ 项，所以如果我们运行 $n-1$ 次该过程，就可以将原同余方程组转换为一个线性同余方程式。同时，该线性同余方程式的 $a$ 项为所求的一个特解 $x_0$ ，通解式为 $x=x_0+k\times lcm(m_1,m_2,m_3\cdots ,m_n)(k\in \mathbb{Z})$

code

long long x,y,m[Size],A[Size],M1,A1,n;

long long mul(long long X,long long Y,long long mod)//龟速乘，防止数据过大
{
	X%=mod;
	Y%=mod;
	long long v=(long double)X*Y/mod;
	long long Ans=X*Y-v*mod;
	if(Ans<0)
		Ans+=mod;
	else
	{
		if(Ans>mod)
		{
			Ans-=mod;
		}
	}
	return Ans;
	
}

long long exgcd(long long a,long long b)
{
	if(b==0)
	{
		x=1;
		y=0;
		return a;
	}
	long long Tem=exgcd(b,a%b);
	long long z=x;
	x=y;
	y=z-y*(a/b);
	return Tem;
}

bool exCRT()
{
	M1=m[0];
	A1=A[0];
	for(int i=1;i<n;i++)
	{
		long long Tem=A[i]-A1;
		long long temp=exgcd(M1 , m[i]);
		if(Tem%temp)
		{
			return false;
		}
		long long TEm=m[i]/temp;
		x=((x*(Tem/temp))%TEM+TEm)%TEm;
		A1+=x*M1;
		M1=M1/temp*m[i];
	}
	return true;
}