针对线性方程组
中国剩余定理给出了以下的一元线性同余方程组:
![中国剩余定理1](https://i-blog.csdnimg.cn/blog_migrate/52cb94f8a9f873d733de233eda6c7a8d.jpeg)
中国剩余定理说明:假设整数
m1,
m2, ... ,
mn两两互质,则对任意的整数:
a1,
a2, ... ,
an,
方程组(S)
有解,并且通解可以用如下方式构造得到:
设
是整数
m1,
m2, ... ,
mn的乘积,并设
是除了
mi以外的
n- 1个整数的乘积。
![中国剩余定理2](https://i-blog.csdnimg.cn/blog_migrate/1c06faac8d71acfb76a014c29293e49f.jpeg)
![中国剩余定理3](https://i-blog.csdnimg.cn/blog_migrate/567227f09ad901fd20afd36f5e0e5534.jpeg)
设
这个就是逆元了
![中国剩余定理4](https://i-blog.csdnimg.cn/blog_migrate/60116bdde261afd187a57fdfdce064a0.jpeg)
![中国剩余定理5](https://i-blog.csdnimg.cn/blog_migrate/c011095ae89875c2f73b305542e339ee.jpeg)
通解形式为
在模M的意义下,方程组(S)只有一个解:
![中国剩余定理6](https://i-blog.csdnimg.cn/blog_migrate/11f25abb9c0beb73b07215426036c380.jpeg)
![中国剩余定理7](https://i-blog.csdnimg.cn/blog_migrate/bfaaecfce3a214f5477bf5b825d4b48a.jpeg)
ll a[N],m[N];
void ex_gcd(ll a,ll b,ll &x,ll &y,ll &d)
{
if(!b)
{
d=a;x=1;y=0;
}
else
{
ex_gcd(b,a%b,y,x,d);
y-=x*(a/b);
}
}
ll inv(ll a,ll n)
{
ll x,y,d;
ex_gcd(a,n,x,y,d);
return d==1?(x+n)%n:-1;//检查gcd(a,n)是否为1
}
ll china(int n,ll *a,ll *m)//X=a[i] (mod m[i]) 这个线性方程组一共有n个方程 0~n-1
{
ll mul = 1, res = 0;
for (int i = 0; i < n; i++)
{
mul *= m[i];
}
for (int i = 0; i < n; i++)
{
ll now = mul / m[i];//now表示是除了当前的这个m[i]外,其他m的公倍数
res = (res + now * inv(now, m[i])*a[i]) % mul;
}
return (res + mul) % mul;//线性方程组在0~mul区间上有唯一解
}