扩展中国剩余定理模板题-牛客-一个小问题

链接：https://ac.nowcoder.com/acm/problem/15068
来源：牛客网
 

uu遇到了一个小问题，可是他不想答。你能替他解决这个问题吗？
问题：给你k对a和r是否存在一个正整数x使每队a和r都满足：x mod a=r，求最小正解x或无解。

 

输入描述:

第一行是正整数k(k<=100000)
接下来k行，每行有俩个正整数a，r(100000>a>r>=0)

输出描述:

在每个测试用例输出非负整数m，占一行。
如果有多个可能的值，输出最小的值。
如果没有可能的值，则输出-1。

示例1

输入

复制

2
8 7
11 9

输出

复制

31

证明摘抄大佬博客：

https://www.cnblogs.com/zwfymqz/p/8425731.html

 

首先考虑同余方程组只有两个式子的时候

x≡c1(mod m1)

x≡c2(mod m2)

 

将两个式子变形

x=c1+m1k1

x=c2+m2k2

 

联立

c1+m1k1=c2+m2k2

 

移项

m1k1=c2−c1+m2k2

 

我们用(a,b)表示a,b的最大公约数

 

在这里需要注意，这个方程有解的条件是(m1,m2) | (c2−c1)，

因为后面会用到(c2−c1)/(m2,m1)

这一项，如果不整除的话肯定会出现小数。

 

对于上面的方程，两边同除(m1,m2)

 

 

转换一下

 

 

此时我们已经成功把k2消去了。

 

同余式两边同除m1 / (m1,m2)

 

接下来怎么办呢？这个式子已经化到最简了。。

不要忘了，我们刚开始还有两个式子。我们把k1代回去！

 

此时，整个式子中的元素我们都已经知道了

具体一点，这个式子可以看做是

x≡c(modm)

其中

 

推广一下

我们每次把两个同余式合并，求解之后得到一个新的同余式。再把新的同余式和其他的联立，最终就可以求出满足条件的解

 

模板代码：

#include<iostream>
#include<algorithm>
#include<cmath>
#include<stdio.h>
#include<string.h>
#include<vector>
#include<map>
#include<set>
#include<stack>
#define debug cout<<"*"<<endl;
#define input(x) scanf("%d",&x)
#define output(x) printf("%d\n",x);
#define llinput(x) scanf("%lld",&x)
#define lloutput(x) printf("%lld\n",x);
using namespace std;
typedef long long ll;
typedef unsigned long long ull;
typedef pair<int,int> pii;
typedef pair<ll,ll> pll;
typedef pair<string,int> psi;
typedef pair<char,char> pcc;
const int mod=998244353;
const int INF=1e9+7;
const int maxn=1e5+7;
const double PI=acos(-1);
const int N=3007;
void Init()
{
    ios::sync_with_stdio(false);
    cin.tie(0);
    cout.tie(0);
}

ll exgcd(ll a, ll b, ll &x, ll &y)
{
    if(b == 0)
    {
        x = 1;
        y = 0;
        return a;
    }
    ll gcd = exgcd(b, a%b, x, y);
    ll t = x;
    x = y;
    y = t - a/b*y;
    return gcd;
}

ll China(ll B[], ll W[], ll k) //W为除数，B是余数
{
    ll i, a, b, c, d, x, y, t;
    for(i=1; i<k; ++i)
    {
        a = W[i-1], b = W[i];
        c = B[i] - B[i-1];
        d = exgcd(a, b, x, y);
        if(c%d)
            return -1;
        t = b/d;
        x = (x*(c/d)%t+t)%t;
        B[i] = B[i-1] + W[i-1]*x;
        W[i] = W[i]/d*W[i-1];
    }
    return B[k-1];
}
ll m[100101], a[100101];
int main(void) {
    int k;
    while (cin >> k) {
        for (int i = 0; i < k; i++) {
            cin >> m[i] >> a[i];
        }
        printf("%lld\n", China(a,m,k));
    }
    return 0;
}