中国剩余定理（不互质版）

最新推荐文章于 2022-11-27 16:43:12 发布

您的好友急速上线

最新推荐文章于 2022-11-27 16:43:12 发布

阅读量1.3k

点赞数 2

分类专栏：数论

本文链接：https://blog.csdn.net/qq_36979930/article/details/81532811

版权

数论专栏收录该内容

8 篇文章 0 订阅

订阅专栏

博客介绍了中国剩余定理在不互质情况下的应用，通过记录一道HDU 3579的模板题来帮助理解。目前博主主要侧重于记忆模板，待深入理解原理后再进行详细阐述。

摘要由CSDN通过智能技术生成

中国剩余定理...有点难懂呀，还是记一下模板好了刚好接触到一道模板题，是在hdu上的 3579

中国剩余定理（不互质）模板题

暂时先把模板放在这里，等什么时候理解了中国剩余定理再回来写理解吧，目前还是就学着用一下模板吧

///不互质的中国剩余定理
///HDU3579亲测AC
#include <iostream>
#include <cstdio>
#include <algorithm>
#include <cstring>

using namespace std;
const int maxn=1000;
int mod[maxn];
int sty[maxn];

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

int Chinese(int sty[],int mod[],int len){
    bool flag=false;
    int n1=mod[0],b1=sty[0],x,y;
    for(int i=1;i<len;i++){
        int n2=mod[i],b2=sty[i];
        int bb=b2-b1;
        int d=ex_gcd(n1,n2,x,y);
        if(bb%d){
            flag=true;
            break;
        }
        int k=bb/d*x;
        int t=n2/d;
        if(t<0)
            t=-t;
        k=(k%t+t)%t;
        b1=b1+n1*k;
        n1=n1/d*n2;
    }
    if(flag)
        return -1;
    return b1;
}
int main()
{
    int T,n,cas;
    cas=1;
    scanf("%d",&T);
    while(T--){
        scanf("%d",&n);
        for(int i=0;i<n;i++)
            scanf("%d",&mod[i]);
        for(int i=0;i<n;i++)
            scanf("%d",&sty[i]);
        if(n==1&&sty[0]==0){
            printf("Case %d: %d\n",cas++,mod[0]);
            continue;
        }
        printf("Case %d: %d\n",cas++,Chinese(sty,mod,n));
    }
    return 0;
}