中国剩余定理

最新推荐文章于 2019-08-04 12:34:05 发布

nowting_csdn

最新推荐文章于 2019-08-04 12:34:05 发布

阅读量123

点赞数 2

分类专栏：程序设计与分析（c++实现）文章标签：数论

个人笔记，仅供复习

本文链接：https://blog.csdn.net/weixin_42373330/article/details/93234697

版权

程序设计与分析（c++实现）专栏收录该内容

16 篇文章 0 订阅

订阅专栏

定理：

设 $m_{1},m_{2},···,m_{n}$ ,是两两互素的n个正整数，记 $m=m_{1}m_{2}···m_{n},M_{i}=\frac{m}{m_{i}},1\leq i \leq n$ ,那么下列模线性方程组
$(*)=\begin{cases}\\ &x\equiv b_{1} (mod\ \ m_{1}) \\ &x\equiv b_{2} (mod\ \ m_{2}) \\ &\ \ \ \ \ \ \ \ \ · \\ &\ \ \ \ \ \ \ \ \ · \\ &\ \ \ \ \ \ \ \ \ · \\ &x\equiv b_{n} (mod\ \ m_{n}) \\ \end{cases}$
的解为 $x_0\equiv (b_1M_1y_1+b_2M_2y_2+b_3M_3y_3+···+b_nM_ny_n)mod \ m$ ,其中， $y_i$ 是 $M_iy_i\equiv 1 mod(\ m_i\ )$ 的解， $i = 1, 2, 3, \cdot \cdot \cdot, n$ 。
证明略。（可参考特殊情况：https://blog.csdn.net/weixin_42373330/article/details/82114105）
程序实现：


#define ll long long 

ll extend_gcd(ll a, ll b, ll& x, ll& y)
{
	if(b==0)
	{
		x=1;
		y=0;
		return a;
	}
	ll t=extend_gcd(b,a%b,y, x);
	y-=(a/b)*x;
	return t;
}

ll ChinaRemainder(int m0[], int b[])
{
	ll d, x, y, n, m=1, a = 0;
	for(int i=1;i<= n;i++) m=m*m0[i];
	for(int i=1;i<=n;i++)
	{
		ll MM=m/m0[i];
		extend_gcd(MM,m0[i], x, y);//求MM关于m0[i]的逆元yi
		a=(a+b[i]*m0[i]*x)%m;
	}
	return a;
}