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

---

解题思路

可以先参考看看》》中国剩余定理（CRT ）(跳过前面的铺垫直接看解法）

其实思路就是上面那个文章中的思路。首先，依题意得：
$x\equiv b_i(mod{a}_i)$
而这样的方程有很多，自然变成一个同余方差组。

我们先举个解方程的例子，再带回题目中看：
再将问题进行分解，看 $x_1$ 如何求：
我们可以先求出 同余方程$x{x}_1\equiv 1(3),x{x}_1\equiv 0(5),x{x}_1\equiv 0(7)$ 的解$x{x}_1$
那么 $x_1=2\ast x{x}_1$

求 $x{x}_1$：
因 $x{x}_1\equiv 1(3),x{x}_1\equiv 0(5),x{x}_1\equiv 0(7)$ 所以xx一定是 5和 7的公倍数，我们设 $y_1=35t$ 那么现在就变成了 $35t\equiv 1(3)$，可以转化成 $35t-1=3k$ 再转化一下： $35t-3k=1$。

这样就变成一个$exgcd$，扩欧求出t。

照这样求出所有的 $xx，x=2\ast y_1+3\ast y_2+2\ast y_3$
但是此时 x 还不是最小解， x 要不断减去 $3,5,7$ 的最小公倍数直到不能减(==直接取模)

回到题目，把$x\equiv b_i(mod{a}_i)$转化为：$b_i+a_{i}y=x$ ——> $x-a_{i}y=b_i$
WOW，和上面的方程如此相似，扩欧走起。。
（因为题目保证 $a_i$ 与 $a_j$ 互质，所以我们可以扩欧）

---

代码

#include<iostream>
#include<cmath>
#include<cstring>
#include<cstdio>
#include<algorithm>
#define ll long long

using namespace std;

int n;
ll a[100010],b[100010],maxn,ans;

void exgcd(ll a,ll b,ll &x,ll &y)
{
	if(b==0)
	{
		x=1,y=0;
		return;
	}
	exgcd(b,a%b,x,y);
	ll z=x;
	x=y;
	y=z-(a/b)*y;
}

int main() {
	scanf("%d",&n);
	maxn=1;
	for(int i=1;i<=n;i++)
	{
		scanf("%lld%lld",&a[i],&b[i]);
		maxn=maxn*a[i];
	}
	for(int i=1;i<=n;i++)
	{
		ll k=maxn/a[i];
		ll x=0,y=0;
		exgcd(k,a[i],x,y);
		if(x<=0)x+=a[i];
		ans=ans+x*k*b[i];
	}
	printf("%lld",ans%maxn);
}