中国剩余定理

前置知识：求逆元
中国剩余定理 (Chinese Remainder Theorem, CRT) 可求解如下形式的一元线性同余方程组（其中 $n_1,n_2,n_3,...,n_k$ 两两互质）：

$\{\begin{array}{c}x\equiv a_1(mod{n}_1)\\ x\equiv a_2(mod{n}_2)\\ ...\\ x\equiv a_k(mod{n}_k)\end{array}$

算法流程
1、计算所有模数的积M
2、对于第i个方程
$a、计算m_i=\frac{M}{n_i}$
$b、计算m_i在模n_i意义下的乘法逆元m_i^{-1}$
$c、计算c_i=m_i\times m_i^{-1}$
3、方程组的唯一解为$x={\sum }_{i=1}^k{a}_i{c}_i(modM)$

P1495 【模板】中国剩余定理(CRT)/曹冲养猪

完整代码

#include<bits/stdc++.h>

using namespace std;

typedef long long LL;

const int N = 110;

int n;
int A[N], B[N];

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;
	LL M = 1;
	for(int i = 1; i <= n; ++ i)
	{
		cin >> A[i] >> B[i];
		M *= A[i];
	}
	
	LL s = 0;
	for(int i = 1; i <= n; ++ i)
	{
		LL m = M / A[i];
		LL ti, x;
		exgcd(m, A[i], ti, x);
		s += B[i] * m * ti;
	}
	
	cout << (s % M + M) % M << endl;
	
	return 0;
}