P1495 中国剩余定理(CRT)/ 曹冲养猪 超详细题解
前言
因为洛谷的题目本身就有很优秀的题解,所以本系列的主要是用于自身的学习打卡记录,如果在此基础上也能帮助到其他同学,那就再好不过了。顺带一提,可能有一定的内容是洛谷题解基础上增加了自己的理解构成的,侵,联,删。
一、题面
1. 原题(贴图)
二、解析及AC代码
1.解析
这个题是个很标准的中国剩余定理类的题目,为了避免有同学不知道中国剩余定理,我们贴一张图片
那么不难发现定理的a即是题目中的b,定理中的m即使题目中的a,那么题目就很简单了,就是套公式而已,这里提醒一下
M
i
M_i
Mi^-1不是-1次方,而是乘法逆元(【不懂点这】)
关于其计算,常用扩展欧几里得算法、快速幂算法计算,这里采用扩展欧几里得算法,关于其计算原理,这里贴一张 _ WZT _ 大佬的图
那么,此时同学们就可以自己去尝试一下了,如果还是不太理解,可以参考下博主下面的AC代码,会带有一定的注释哦
2.AC代码
#include<iostream>
#include<algorithm>
#include<cstring>
#include<map>
#include<vector>
using namespace std;
typedef long long LL;
const int N = 1e6 + 10, MAXI = 0x3f3f3f3f;
LL t, n, m, k, z;
LL a[N], b[N];
LL M = 1;
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()
{
cin >> n;
for (int i = 1; i <= n; i++)
{
cin >> a[i] >> b[i];
M *= a[i];
}
LL res = 0;
for (int i = 1; i <= n; i++)
{
LL mi = M / a[i];
LL x = 0, y = 0;
//扩展欧几里得算法计算逆元
exgcd(mi, a[i], x, y);
LL MI = x;
res = ((res + b[i] * mi * MI) % M + M) % M;
}
cout << res;
}
如果觉得有用还请点个赞吧,拜托拜托