曹冲养猪（中国剩余定理(CRT)）-CSDN博客

本文链接：https://blog.csdn.net/2301_80065123/article/details/138732137

曹冲养猪（中国剩余定理(CRT)）

题目描述

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

输入格式

第一行包含一个整数 $n$ —— 建立猪圈的次数，接下来 $n$ 行，每行两个整数 $a_i, b_i$ ，表示建立了 $a_i$ 个猪圈，有 $b_i$ 头猪没有去处。你可以假定 $a_1 \sim a_n$ 互质。

输出格式

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

样例 #1

样例输入 #1

样例输出 #1

提示说明

在这里插入图片描述

$\leq n\le10$ ， $\leq b_i\lt a_i\le100000$ ， $\leq \prod a_i \leq 10^{18}$

代码内容

// #include <iostream>
// #include <algorithm>
// #include <cstring>
// #include <stack>//栈
// #include <deque>//队列
// #include <queue>//堆/优先队列
// #include <map>//映射
// #include <unordered_map>//哈希表
// #include <vector>//容器，存数组的数，表数组的长度
#include <bits/stdc++.h>

using namespace std;

typedef long long ll;

//扩展欧几里得算法(裴蜀定理)
//根据费马定理，任意正整数a, b都存在整数x, y
//使得ax + by = gcd(a, b)
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()
{
    ll n;
    cin>>n;
    ll a1,b1;
    cin>>a1>>b1;

    ll op=0;
    for(ll i=0;i<n-1;i++)
    {
        ll a2,b2;
        cin>>a2>>b2;
        ll k1,k2;
        //扩展欧几里得算法
        //k1*a1+k2*a2=gcd(a,b)=d
        /*((b2-b1)/d)*(a*x+b*y)=d*((b2-b1)/d)*/
        //x≡b(mod a)
        //x=k1*a1+b1;
        //x=k2*a2+b2;
        //k1*a1+k2*a2=b2-b1
        ll d=exgcd(a1,a2,k1,k2);
        if((b2-b1)%d)
        {
            op=-1;
            break;
        }

        //更新状态
        k1*=(b2-b1)/d;
        ll t=a2/d;
        //将解变成一个最小的正整数解
        k1=(k1%t+t)%t;

        op=k1*a1+b1;
        //更新a和m(k为变量,取余时会自动消失)
        b1=k1*a1+b1;
        a1=abs(a1/d*a2);
    }

    ll res;
    if(op!=-1) res=(b1%a1+a1)%a1;

    cout<<res<<endl;
    return 0;
}