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

题目描述

自从曹冲搞定了大象以后，曹操就开始捉摸让儿子干些事业，于是派他到中原养猪场养猪，可是曹冲满不高兴，于是在工作中马马虎虎，有一次曹操想知道母猪的数量，于是曹冲想狠狠耍曹操一把。举个例子，假如有 16 头母猪，如果建了 3 个猪圈，剩下 1 头猪就没有地方安家了。如果建造了 5 个猪圈，但是仍然有 1 头猪没有地方去，然后如果建造了 7 个猪圈，还有 2 头没有地方去。你作为曹总的私人秘书理所当然要将准确的猪数报给曹总，你该怎么办？

输入格式
第一行包含一个整数 n——建立猪圈的次数，接下来 n 行，每行两个整数 $a_i{b}_i$
表示建立了 $a_i$ 个猪圈，有 $b_i$ 头猪没有去处。你可以假定 $a_i$，$a_j$互质。

输出格式

输出包含一个正整数，即为曹冲至少养母猪的数目。

输入输出样例

输入

3
3 1
5 1
7 2

输出

16

说明/提示

$1\le n\le 10，0\le b_i<a_i\le 1000000，1\le \prod a_i\le 10^{18}$

题解

不是吧不是吧不会有人看不出来 这是模板CRT，
即做法为套公式：
$m=\prod a_i(即最小公倍数)t=m/a[i]x=t^{-1}(mod{a}_i)$
$ans=\sum x\ast t\ast b[i]$
其中逆元用扩展欧几里得求
当然这只是CRT的其中一个解，要求最小解还得模 $m$。
具体的CRT证明这里不多讲，本人是看一篇巨佬的博客懂的

代码

#include<iostream>
#include<cstdio>
#define ll long long //不开longlong见祖宗
using namespace std;
ll n,m=1,ans,t,x,y,a[20],b[20];
void exgcd(ll a,ll b)//扩展欧几里得算法求逆元
{
	if(b==0)
	{
		x=1;
		y=0;
		return;
	}
	exgcd(b,a%b);
	int t=x;
	x=y;
	y=t-a/b*y;
}
int main()
{
	scanf("%lld",&n);
	for(int i=1; i<=n; i++)
	{
		scanf("%lld%lld",&a[i],&b[i]);
		m*=a[i];//求最小公倍数
	}
	for(int i=1; i<=n; i++)
	{
		t=m/a[i];
		x=y=0;
		exgcd(t,a[i]);
		x=(x%a[i]+a[i])%a[i];
		ans=(ans+x*t*b[i])%m;//CRT 记得模数
	}cout<<ans;
}