数论学习之中国剩余定理 china china!!!

针对线性方程组

中国剩余定理给出了以下的一元线性同余方程组:
中国剩余定理1
 
中国剩余定理说明:假设整数 m1, m2, ... , mn两两互质,则对任意的整数: a1, a2, ... , an,
  方程组(S)
有解,并且通解可以用如下方式构造得到:
  中国剩余定理2
是整数 m1, m2, ... , mn的乘积,并设
  中国剩余定理3
是除了 mi以外的 n- 1个整数的乘积。
  中国剩余定理4
这个就是逆元了
  中国剩余定理5 
通解形式为
  中国剩余定理6 
在模M的意义下,方程组(S) 只有一个解:
  中国剩余定理7
 

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区间上有唯一解
}

你可能感兴趣的:(数学,数论)