题目链接
题意:给定2n个整数a1,a2,…,an和m1,m2,…,mn,求一个最小的非负整数x,满足∀i∈[1,n],x≡mi(mod ai)。
输入格式
第1行包含整数n。
第2…n行:每i+1行包含两个整数ai和mi,数之间用空格隔开。
输出格式
输出最小非负整数x,如果x不存在,则输出-1。
如果存在x,则数据保证x一定在64位整数范围内。
数据范围
1≤ai≤2^31−1,
0≤mi<ai
1≤n≤25
输入样例:
2
8 7
11 9
输出样例:
31
思路:要求找到最小的一个x满足∀i∈[1,n],x % ai =mi,直接就有个中国剩余定理(前提条件模数ai必需两两互质)
- A = a1 * a2 * ……*an
- Ai = A / ai
- Ai^-1 表示 Mi模mi的逆元
- x = m1 * A1 *A1^-1 + m2 * A2 * A2^-1……mn * An * An^-1
比如把x带入x ≡ m1(mod a1)中证明:因为A1*A1^-1相当于1,而后面的A2,A3,An模上a1都是0,所有这个x对于每个表达式成立。
显然这个题并不满足互质的条件,就得从定理的根本原理来下手了。
转换一下:
x = k1 * a1 + m1
x = k2 * a2 + m2
k1 * a1 - k2 * a2 = m2 - m1
( a1,a2,m1,m2都已知,相当于求a1 * x + a2 * y = c;就是扩展欧几里得了啊。)
如果gcd(a1, a2) | m2 - m1,便有解
整理下就有:
x = a1 * k1 + m1 + k * gcd(a1, a2)
x = x0 + k * a (用x0和a来替代),显然可以把n个式子两两结合最后变成一个 这个时候要求的x % a = x0.
代码实现:
#include<iostream>
#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 -= a / b * x;
return d;
}
int main()
{
int n;
scanf("%d", &n);
ll a1, m1;
cin >> a1 >> m1;
bool has_answer = 1;
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 = 0;
break;
}
k1 *= (m2 - m1) / d;
ll t = a2 / d;
//避免k1很大,让它变成最小的一个正解
k1 = (k1 % t + t) % t;
m1 = a1 * k1 + m1;
a1 = abs(a1 / d * a2);
}
if(has_answer){
//取一个正的余数
cout << (m1 % a1 + a1) % a1 << endl;
}
else cout << -1 << endl;
return 0;
}