C++实现（同余方程组）中国剩余定理

最新推荐文章于 2025-01-27 21:08:14 发布

lgq2016

最新推荐文章于 2025-01-27 21:08:14 发布

阅读量3.6k

点赞数 2

分类专栏：网络安全理论文章标签：中国剩余定理 C++实现

网络安全理论专栏收录该内容

1 篇文章

订阅专栏

本文详细介绍了一种解决中国剩余数定理问题的算法实现，该算法能够处理多个模数互质或不互质的情况，通过扩展欧几里得算法求解未知数x，满足一系列同余方程。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

#include<iostream>
#include<stdio.h>
using namespace std;
#define LL long long

/***********************************
中国剩余数定理
ai=x(mod mi)    求x
eg:1=x(mod 11) 2=x(mod 7)
最后求得为满足条件的最小值
************************************/
LL ex_gcd(LL a,LL b,LL &x,LL &y)
{
    if(b==0)
    {
        x=1;y=0;
        return a;
    }
    LL r=ex_gcd(b,a%b,x,y);
    LL t=x;
    x=y;
    y=t-a/b*y;
    return r;
}

///n个mi互质
//const LL maxn = 20;
//LL a[maxn], m[maxn], n;
//LL CRT(LL a[], LL m[], LL n)
//{
//    LL M = 1;
//    for (int i = 0; i < n; i++) M *= m[i];
//    LL ret = 0;
//    for (int i = 0; i < n; i++)
//    {
//        LL x, y;
//        LL tm = M / m[i];
//        ex_gcd(tm, m[i], x, y);
//        ret = (ret + tm * x * a[i]) % M;
//    }
//    return (ret + M) % M;
//}

///n个mi不互质
const LL maxn = 1000;
LL a[maxn], m[maxn], n;
LL CRT(LL a[], LL m[], LL n) {
    if (n == 1) {
        if (m[0] > a[0]) return a[0];
        else return -1;
    }
    LL x, y, d;
    for (int i = 1; i < n; i++) {
        if (m[i] <= a[i]) return -1;
        d = ex_gcd(m[0], m[i], x, y);
        if ((a[i] - a[0]) % d != 0) return -1;  //不能整除
        LL t = m[i] / d;
        x = ((a[i] - a[0]) / d * x % t + t) % t;
        a[0] = x * m[0] + a[0];
        m[0] = m[0] * m[i] / d;
        a[0] = (a[0] % m[0] + m[0]) % m[0];
    }
    return a[0];
}

int main ()
{
   LL T;
   while( ~scanf("%lld",&T) )
   {
       for(int i=0; i<T; i++) scanf("%lld%lld",&m[i],&a[i]);
       cout<<CRT(a,m,T)<<endl;
   }

}

转自：https://blog.youkuaiyun.com/u010468553/article/details/38346195