[AcWing]204. 表达整数的奇怪方式(C++实现)求解线性同余方程组的最小非负整数x(扩展中国剩余定理)模板题
1. 题目
2. 读题(需要重点注意的东西)
思路:
中国剩余定理
中国剩余定理是为了解如下线性同余方程组:
以上只为了解释中国剩余定理,本题的主要思路见下方
。
本题的思路(扩展中国剩余定理)
核心思想:将线性同余方程组按如下方式两两合并即可!
3. 解法
---------------------------------------------------解法---------------------------------------------------
#include <iostream>
#include <algorithm>
using namespace std;
typedef long long LL;
// 扩展欧几里得算法求最大公约数和k1、k2
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;
cin >> n;
// 现有方程m1、a1
LL x = 0, a1, m1;
cin >> a1 >> m1;
// -----核心:用现有的m1和a1合并m和a--------
for (int i = 0; i < n - 1; i ++ )
{
LL a2, m2;
cin >> a2 >> m2;
// 算出一组k1、k2
LL k1, k2;
// d是a1和a2的最大公约数
LL d = exgcd(a1, a2, k1, k2);
// 如果无解
if ((m2 - m1) % d)
{
x = -1; // 无解,输出-1
break;
}
k1 *= (m2 - m1) / d;
// 避免溢出,将k1变成最小的正整数解
k1 = (k1 % (a2/d) + a2/d) % (a2/d);
//
x = k1 * a1 + m1;
// 求a1和a2的最小公倍数
LL a = abs(a1 / d * a2);
m1 = k1 * a1 + m1; // 将k1 * a1 + m1赋值给m1,用于下一次合并
a1 = a; // 将a1和a2的最小公倍数再赋值给a1,用于下一次合并
}
// 求出最小正整数解
if (x != -1) x = (m1 % a1 + a1) % a1;
cout << x << endl;
return 0;
}
4. 可能有帮助的前置习题
5. 所用到的数据结构与算法思想
- 中国剩余定理
- 扩展欧几里得算法
6. 总结
求解线性同余方程组的最小非负整数x(扩展中国剩余定理)模板题,理解公式并背下代码。