中国剩余定理

最新推荐文章于 2022-02-14 16:47:39 发布

DAYH

最新推荐文章于 2022-02-14 16:47:39 发布

阅读量346

点赞数

文章标签：算法线性代数机器学习

本文链接：https://blog.csdn.net/qq_56427460/article/details/122334580

版权

中国剩余定理用于求解同余方程

$\begin{cases} x\equiv a_1 (\mod m_1)\\ x\equiv a_2(\mod m_2)\\ \cdots\cdots\\ x\equiv a_n(\mod m_n) \end{cases}$

公式 $x = (a_1\times M_1\times M_1^{-1} + \cdots+ a_2 \times M_n\times M_i^{-1}) \% M$

其中 $M=m_1\times m_2 \times \cdots \times m_n\\$

$M_i=M/bi\\$

$M_i^{-1}$ 为 $M_i$ 模 $m_i$ 下的逆元

例题：生物节律阿九大战朱最学

//阿九大战朱最学
#include<iostream>

using namespace std;

typedef long long LL;
//扩展欧几里得求逆元
LL exgcd(LL a, LL b, int &x, int &y)
{
  if (!b){x = 1, y = 0; return a;}
  LL d = exgcd(b, a % b, y, x);
  y -= (a / b) * x;
  
  return d;
}
//中国剩余定理
LL rct(LL *a, LL *m,int n)
{
    LL ans = 0;
    LL M = 1;
    //m为除数的积
    for (int i = 0; i < n; i++) M *= m[i];
    
    for (int i = 0; i < n; i++) 
    {
        int x , y;
        LL mi = M / m[i];
        exgcd(mi, m[i], x, y);
        ans = (ans + a[i] * mi % M * x) % M;
    }
    //cout << M << endl;
    return (ans + M) % M;
}


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