中国剩余定理模板 Strange Way to Express Integers POJ - 2891

讲解：https://www.cnblogs.com/zylAK/p/9581919.html
取模的数不互素的版本

const int maxn=1e5+5;
void gcd(long long  a,long long b,long long &d,long long &x,long long &y )
{
    if(!b)
    {
        d=a;
        x=1;
        y=0;
    }
    else
    {
        gcd(b,a%b,d,y,x);
        y-=x*(a/b);
    }
}
long long rest[maxn],a[maxn];
long long china(long long n)
{
    long long m=1,d,y,x=0;
    for(int i=0; i<n; i++)
    {
        m*=rest[i];
    }
    for(int i=0; i<n; i++)
    {
        long long w=m/rest[i];
        gcd(rest[i],w,d,d,y);
        x=(x+y*w*a[i])%m;
    }
    return (x+m)%m;
}

题目链接：https://vjudge.net/problem/POJ-2891
题意：给出多个a和r，求一个最小的数x满足所有x%a=r。
注意取模的数不互素

#include<cstdio>
using namespace std;
const int maxn=1e5+5;
void gcd(long long  a,long long b,long long &d,long long &x,long long &y )
{
    if(!b)
    {
        d=a;
        x=1;
        y=0;
    }
    else
    {
        gcd(b,a%b,d,y,x);
        y-=x*(a/b);
    }
}
long long rest[maxn],a[maxn];
long long china(long long n)
{
    long long m=a[1],z=rest[1];
    long long x,y,d;
    for(int i=2; i<=n; i++)
    {
        gcd(m,a[i],d,x,y);
        if((rest[i]-z)%d!=0)
            return -1;
        x=(rest[i]-z)/d*x%a[i];
        z+=(x*m);
        m=m/d*a[i];
        z%=m;
    }
    return (z%m+m)%m;
}
int main()
{
    long long k;
    while(~scanf("%lld",&k))
    {
        for(int i=1; i<=k; i++)
        {
            scanf("%lld%lld",&a[i],&rest[i]);
        }
        printf("%lld\n",china(k));
    }
    return 0;
}