中国剩余定理

最新推荐文章于 2020-03-26 14:19:43 发布

*笑三笑

最新推荐文章于 2020-03-26 14:19:43 发布

阅读量121

点赞数

分类专栏：中国剩余定理

原文链接：https://blog.csdn.net/niiick/article/details/80229217

版权

中国剩余定理专栏收录该内容

3 篇文章 0 订阅

订阅专栏

求解同余方程组
在这里插入图片描述
设M=m1m2m3m4mi…;
则Mi=M/mi;
x=a1M1pow(M1,-1)+a2M2pow(M2,-1)+a3M3pow(M3,-1)…
又因为求最小正整数解 x=(x%M+M)%M;

void exgcd(int a,int b,int &x,int &y)
{
    if(b==0){ x=1; y=0; return;}
    exgcd(b,a%b,x,y);
    int tp=x;
    x=y; y=tp-a/b*y;
}

int china()
{
    int ans=0,lcm=1,x,y;
    for(int i=1;i<=k;++i) lcm*=b[i];
    for(int i=1;i<=k;++i)
    {
        int tp=lcm/b[i];
        exgcd(tp,b[i],x,y);
        x=(x%b[i]+b[i])%b[i];//x要为最小非负整数解
        ans=(ans+tp*x*a[i])%lcm;
    }
    return (ans+lcm)%lcm;
}