[lnsyoj517/luoguP4777]扩展中国剩余定理
题意
原题链接
求线性同余方程组
的最小非负整数解。
sol
与[lnsyoj163/luoguP1495]曹冲养猪不同的是,本题无法保证互质,这就导致中国剩余定理无法使用,需要一种新的方式来解决问题。
我们思考:如果我们可以将两个方程合并为一个方程,那么整个方程组都可以合成为一个方程,问题也就解决。
下面进行推导:
对于两个方程 、,,我们将其分别转换为 和 ,整理得 ,利用扩展欧几里得算法,我们可以求出一组可行解 。此时,新的方程的 即为原来的 ,新的方程的 即为原来的 的最大公倍数。
最终,合并成的一个方程,其 的值即为最终答案.
注意:勤取模
代码
#include <iostream>
#include <algorithm>
#include <cstring>
using namespace std;
typedef long long LL;
const int N = 100005;
int n;
LL a[N], m[N];
LL exgcd(LL a, LL b, __int128 &x, __int128 &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(){
while (scanf("%d", &n) != EOF){
for (int i = 1; i <= n; i ++ ) scanf("%lld%lld", &m[i], &a[i]);
LL b = a[1];
__int128 M = m[1];
bool flag = false;
for (int i = 2; i <= n; i ++ ){
__int128 p, q;
LL d = exgcd(M, m[i], p, q);
if (((a[i] - b % m[i] + m[i]) % m[i]) % d) {
flag = true;
puts("-1");
break;
}
__int128 u = 1;
p = u * p * ((a[i] - b % m[i] + m[i]) % m[i]) / d;
LL t = abs(m[i] / d);
p = (p % t + t) % t;
LL m1 = M;
M = M / d * m[i];
b = ((u * m1 * p % M + b % M) % M + M) % M;
}
if (flag) continue;
printf("%lld\n", b);
}
return 0;
}
分类:
题解 / 2024训练
标签:
数学
【推荐】国内首个AI IDE,深度理解中文开发场景,立即下载体验Trae
【推荐】编程新体验,更懂你的AI,立即体验豆包MarsCode编程助手
【推荐】抖音旗下AI助手豆包,你的智能百科全书,全免费不限次数
【推荐】轻量又高性能的 SSH 工具 IShell:AI 加持,快人一步
· 地球OL攻略 —— 某应届生求职总结
· 周边上新:园子的第一款马克杯温暖上架
· Open-Sora 2.0 重磅开源!
· 提示词工程——AI应用必不可少的技术
· .NET周刊【3月第1期 2025-03-02】