[Luogu P1495]曹冲养猪

aideng1990

于 2019-04-16 17:05:00 发布

阅读量115

点赞数

CC 4.0 BY-SA版权

原文链接：http://www.cnblogs.com/clockwhite/p/10718311.html

本文介绍了一个基于中国剩余定理（CRT）的算法实现，通过使用C++编程语言，展示了如何解决数学中的一类问题。该算法利用了扩展欧几里得算法求解逆元，并通过具体实例演示了CRT的计算过程。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

题目链接

中国剩余定理（孙子定理）的裸题。直接放代码。

#include<iostream>
#include<cstdio>
#include<algorithm>
using namespace std;
#define ll long long
ll read(){
    ll res=0,f=1;
    char ch=getchar();
    while(ch<'0'||ch>'9'){
        if(ch=='-')f=-1;
        ch=getchar();
    }
    while(ch>='0'&&ch<='9'){
        res=res*10+(ch-'0');
        ch=getchar();
    }
    return res*f;
}
ll n,a[15],b[15],ans;
void exgcd(ll a,ll b,ll &x,ll &y){
    if(b==0)x=1,y=0;
    else{
        exgcd(b,a%b,x,y);
        ll xt=x;
        x=y;
        y=xt-a/b*y;
    }
}
void intchina(ll *m,ll *u){
    ll Mk=1,M[15];
    for(ll i=1;i<=n;++i)Mk*=m[i];
    for(ll i=1;i<=n;++i){
        M[i]=Mk/m[i];
        ll x,y;
        exgcd(M[i],m[i],x,y);
        ans=(ans%Mk+M[i]*x*u[i]%Mk+Mk)%Mk;
    }
    while(ans<0)ans+=Mk;
    ans%=Mk;
}
int main(){
n=read();
for(int i=1;i<=n;++i)a[i]=read(),b[i]=read();
intchina(a,b);
cout<<ans;
return 0;
}