扩展中国剩余定理

扩展中国剩余定理

x = a[1] (b[1])
x = a2
x = a[1] + k1[b1];
x = a[2] + k2[b2]
k1[b1] + (-k2[b2]) = a[2] - a[1] ;
c = a[2] - a[1]
c = a[2] - x ;

lt excrt()
{
    lt x,y,k;
    lt M=bi[1],ans=ai[1];//第一个方程的解特判
    for(int i=2; i<=n; i++)
    {
        lt a=M,b=bi[i],c=(ai[i]-ans%b+b)%b;//
        lt gcd=exgcd(a,b,x,y),bg=b/gcd;
        if(c%gcd!=0) return -1; //判断是否无解，然而这题其实不用
        x=mul(x,c/gcd,bg);
        ans+=x*M;//更新前k个方程组的答案
        M*=bg;//M为前k个m的lcm
        ans=(ans%M+M)%M;
    }
    return (ans%M+M)%M;
}