当前位置: 代码迷 >> 综合 >> 数学知识:中国剩余定理(孙子定理)
  详细解决方案

数学知识:中国剩余定理(孙子定理)

热度:8   发布时间:2023-12-13 07:39:26.0

ACWing 204. 表达整数的奇怪方式

给定 2n 个整数a1,a2,…,an和m1,m2,…,mn,求一个最小的非负整数 x,满足?i∈[1,n],x≡mi(mod ai)。

输入格式
第1 行包含整数 n。
第 2…n+1行:每 i+1 行包含两个整数ai和mi,数之间用空格隔开。

输出格式
输出最小非负整数 x,如果 x 不存在,则输出 ?1。
如果存在 x,则数据保证 x 一定在64位整数范围内。

数据范围
1≤ai≤231?1,
0≤mi<ai
1≤n≤25
输入样例:
2
8 7
11 9
输出样例:
31

Code:

#include <iostream>
#include <cmath>
#include <algorithm>
using namespace std;typedef long long LL;LL exgcd(LL a, LL b, LL &x, LL &y)
{
    if(!b){
    x = 1;y = 0;return a;}LL d = exgcd(b, a % b, y, x);y = y - a / b *x;return d;
}int main()
{
    int n;cin >> n;bool has_answer = true;LL a1, m1;cin >> a1 >> m1;for(int i = 0; i < n - 1; i ++){
    LL a2, m2;cin >> a2 >> m2;LL k1, k2;LL d = exgcd(a1, a2, k1, k2);if((m2- m1) % d){
    has_answer = false;break;}k1 = k1 * (m2 - m1) / d;LL t = a2 / d;k1 = (k1 % t + t) % t;m1 = a1 * k1 + m1;a1 = abs(a1 * a2 / d);}if(has_answer)cout << (m1 % a1 + a1) % a1 << endl;else puts("-1");return 0;
}
  相关解决方案