【模板】中国剩余定理求法

#include<cstdio>
#define ll long long
using namespace std;
const int maxn=15;
int n;
long long a[maxn],b[maxn];
long long exgcd(long long a,long long b,long long &x,long long &y) {
    if(!b) {
        x=1;y=0;
        return a;
    }
    long long q=exgcd(b,a%b,y,x);
    y-=a/b*x;
    return q;
}
long long china() {
    long long M,res; 
    M=1;res=0;
    for(int i=1;i<=n;i++) M*=a[i];
    for(int i=1;i<=n;i++) {
        long long m,p,d,y;
        m=M/a[i];d=exgcd(m,a[i],p,y);
        res=(res+m*p*b[i])%M;
    }
    if(res<0) res+=M;
    return res;
} 
int main() {
    scanf("%d",&n);
    for(int i=1;i<=n;i++) scanf("%lld%lld",&a[i],&b[i]);
    printf("%lld",china());
}
posted @ 2020-03-20 19:58  仰望星空的蚂蚁  阅读(8)  评论(0)    收藏  举报  来源