中国剩余定理以及扩展中国剩余定理

最新推荐文章于 2024-06-25 21:10:30 发布

呱太

最新推荐文章于 2024-06-25 21:10:30 发布

阅读量2.9k

点赞数 8

分类专栏：算法与数据结构文章标签：算法 c++ 扩展中国剩余定理算法算法模板

本文链接：https://blog.csdn.net/qq_35650513/article/details/127469912

版权

算法与数据结构专栏收录该内容

19 篇文章 1 订阅

订阅专栏

在这里插入图片描述

中国剩余定理必须有两两互质的条件；而扩展中国剩余定理没有限制（可能互质，也能不互质）。所以只记忆一个扩展中国剩余定理的板子就行.

在这里插入图片描述

代码

#include <iostream>
#include <algorithm>

using namespace std;
typedef long long LL;
int n;

LL exgcd(LL a, LL b, LL &x, LL &y){
    if(b==0){
        x = 1, y = 0;
        return a;
    }
    LL d, x1, y1;
    d = exgcd(b, a%b, x1, y1);
    x = y1, y = x1 - a/b*y1;
    return d;
}

LL excrt(LL m[], LL r[]){
    LL m1, m2, r1, r2, p, q;
    m1 = m[0], r1 = r[0];
    for(int i = 1; i < n; ++i){
        m2 = m[i], r2 = r[i];
        LL d = exgcd(m1, m2, p, q);
        if( (r2-r1)%d !=0 ) return -1; // 不能整除，说明无解
        p = p*(r2-r1)/d;
        p = (p%(m2/d) + (m2/d))%(m2/d);
        r1 = m1*p+r1;
        m1 = m1*m2/d;
    }
    return (r1%m1+m1)%m1;
    
}


int main(){
    cin >> n;
    LL m[30], r[30];
    for(int i = 0;i < n; ++i){
        cin >> m[i] >> r[i];
    }
    cout << excrt(m, r) << endl;
    
    
    return 0;
}