中国剩余定理 扩展中国剩余定理

一、中国剩余定理

对于 x ≡ a1(mod m1)
x ≡ a2(mod m2)
x ≡ a3(mod m3)
…………
且 m1 m2 m3…两两互质 此时可运用中国剩余定理
首先求出
s=m1m2m3
M[i]=s/m[i]
t[i]为M[i]关于m[i]的逆元
则最终求得的x =∑(a[i]*t[i]*M[i])%s
代码如下

#include<bits/stdc++.h>
#define ll long long
using namespace std;
ll ksc(ll a,ll b,ll p)  //此处用快速乘，防止爆
{
	ll ret=0;
	while(b)
	{
		if(b&1)
		ret=(ret+a)%p;
		a=(a+a)%p;
		b>>=1;
	}
	return ret;
}
ll ojld(ll a,ll b,ll &x,ll&y) //此处为扩展欧几里得的方法求逆元
{
	if(!b)
	{
		x=1;
		y=0;
		return a;
	}
	ll d=ojld(b,a%b,x,y);
	ll t=x;
	x=y;
	y=t-a/b*y;
	return d;
}
ll inv(ll b,ll p)
{
	ll x,y;
	ll w=ojld(b,p,x,y);
	if(w!=1)
	return -1;
	return (x%p+p)%p;
}
int main()
{
	ll k,a[10],b[10],s=1,m[10],sum=0;
	scanf("%lld",&k);
	for(ll i=0;i<k;i++)
	scanf("%lld",&a[i]);
	for(ll i=0;i<k;i++)
	{
	scanf("%lld",&b[i]);
	s*=b[i]; //求出乘积s
	}
	for(ll i=0;i<k;i++)
	{
	m[i]=s/b[i];
	ll k=inv(m[i],b[i]); //每次的k就相当于上述中的t[i]，为每次求得的逆元
	sum=(sum+ksc(ksc(a[i],m[i],s),k,s))%s;//此处乘法运算时运用了快速乘
	}
	printf("%lld\n",sum);
	return 0;
} 

二、扩展中国剩余定理

扩展中国剩余定理跟中国剩余定理是毛关系都没有的
同样是对于 x ≡ a1(mod m1)
x ≡ a2(mod m2)
x ≡ a3(mod m3)
…………
在没有规定 m1 m2 m3 两两互质的情况下求x
扩展中国剩余定理就是利用拓展欧几里得，每次将两个式子合并为一个式子，依次合并，最终求解

如 对于前两个式子 x=a1+k1 * m1
x=a2+k2 * m2
即可得 m1 * k1 - m2 * k2=a2-a1 就相当于 ax+by=c 其中 k1 k2 相当于x，y
为确保c为正 可 int c=((a2-a1)%m2+m2)%m2
通过扩展欧几里得求出gcd(m1,m2) 为d，且求出了k1，然后k1再乘 (a2-a1)/d 便得到了一个解，
（此处相乘运用了快速乘防爆）
然后将此时得到的k1带入 x=a1+k1 * m1 便得到x
这两个式子合并后的方程设为 X=A+k*N 其中这个A就为刚才求得的x，N就为m1 * m2/d
代码如下

#include <cstdio>
const int maxn = 100010;
typedef long long ll;
int n;
ll a[maxn],b[maxn],a1,b1,x,y;
ll exgcd(ll a,ll b,ll &x,ll &y)
{
    if(!b)
    { 
    x=1; 
    y=0;
    return a;
    }
    ll d=exgcd(b,a%b,x,y);
    ll z=x; 
       x=y;
       y=z-(a/b)*y;
       return d;
}
ll ksc(ll n,ll k,ll mod){
    ll ans=0;
    while(k)
	{
      if(k&1) ans=(ans+n)%mod;
      n=(n+n)%mod;
	  k>>=1;
    }
    return ans;
}
int main(){
    scanf("%d", &n);
    for(int i=1;i<=n;++i)
    scanf("%lld%lld", &b[i], &a[i]);
    a1=a[1];
    b1=b[1];//此处对第一个式子的a[1] b[1]赋值，在后面的for循环中每次更新a1 b1 的值为合并后两式子的A N
    for(int i=2;i<=n;++i)
	{
       	ll c=((a[i]-a1)%b[i]+b[i])%b[i];
       	ll g=exgcd(b1,b[i],x,y);
       	x=ksc(x,c/g,b[i]);
       	a1=a1+b1*x;//将求得的x带入原式子得到值作为合并式子的a1值
       	b1=b1*(b[i]/g);//合并式子的N值为两个式子的b的最小公倍数
       	a1=(a1+b1)%b1;//此处取模
    }
    	printf("%lld\n",a1);
   	    return 0;
}