中国剩余定理-表达整数的奇怪方式

表达整数的奇怪方式 

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

中国剩余定理：

举例：

推导：

合并方程组：

Example

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

输入格式

第1行包含整数n。

第2..n行：每i+1行包含两个整数aiai和mimi,数之间用空格隔开。

输出格式

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

数据范围

1≤ai≤231−11≤ai≤231−1,
0≤mi<ai0≤mi<ai
1≤n≤251≤n≤25

输入样例：

2
8 7
11 9

输出样例：

31

#include<iostream>
#include <algorithm>

using namespace std;

typedef long long LL;

//扩展欧几里得算法模版
LL exgcd(LL a, LL b, LL &x, LL &y)
{
    if (!b)//如果b = 0,返回一组最小解
    {
        x = 1, y = 0;
        return a;
    }
    by + (a mod b)x = (a, b)最大公约数d
    LL d = exgcd(b, a % b, y, x);//得到最大公约数
    //a mod b = a - (a / b)下取整 * b带入
    //by + (a - (a / b)下取整 * b)x = d;
    //展开整理：ax + b (y - (a/b)*x)
    //得到a的系数不变为a,y的系数变成：y - (a/b)*x
    y -= a / b * x;
    return d;
}


int main()
{
    //第1行包含整数n。
    int n;
    cin >> n;
   // scanf("%d", &n);
    //第2..n行：每i+1行包含两个整数ai和mi,数之间用空格隔开。
    //合并思路：每次把新的方程合并成现有方程中
    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)//如果余数不是0则无解
        {
            has_answer = false;
            break;
        }
        //否则有解
        //k1a1 - k2a2 = d,d是m1 - m2的倍数，所以得翻若干倍
        //翻成m2 - m1
        k1 *= (m2 - m1) / d;
        //中间过程，k1变成方程最小的整数解
        LL t = a2 / d;//首先把a2 / d存下来
        //然后用k1变成最小正整数解
        k1 = (k1 % t + t) % t;
        
        //最后变成x  = ka + m
        m1 = a1 * k1 + m1;
        
        //新的方程求a1，a2最小公倍数
        a1 = abs(a1 / d * a2); 
    }
    if (has_answer)//如果有解
    {
        //C++上取余可能去负的，比如-5 mod 3 余 -2
        //要想得到正余数，需要-5 mod 3加上3然后在mod 3就可以了
        cout << (m1 % a1 + a1) % a1 << endl;//求的是m1 % a1正的余数，非负的余数
    }//如果无解
    else puts("-1");
    
    return 0;
}