P4777 【模板】扩展中国剩余定理（EXCRT）

题目链接

题面：

代码：

#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<string>
#include<map>
#define ll long long
using namespace std;
const int maxn=100010;
ll a[maxn],m[maxn];
int k;
ll q_pow(ll a,ll b,ll p)
{
    ll ans=0;
    while(b)
    {
        if(b&1) ans=(ans+a)%p;
        a=(a+a)%p;
        b>>=1;

    }
    return ans%p;
}

ll mypow(ll a,ll b,ll p)
{
    ll ans=1;
    while(b)
    {
        if(b&1) ans=q_pow(ans,a,p);
        a=q_pow(a,a,p);
        b>>=1;
    }
    return ans;
}

ll exgcd(ll a,ll b,ll &x,ll &y)
{
    if(b==0)
    {
        x=1;
        y=0;
        return a;
    }

    ll d=exgcd(b,a%b,x,y);

    ll z=x;
    x=y;
    y=z-a/b*y;
    return d;
}

ll china(void)
{
    ll ans=a[1];
    ll x,y;
    ll M=m[1];
    for(int i=2;i<=k;i++)
    {

        ll c=(a[i]-ans%m[i]+m[i])%m[i];

        ll d=exgcd(M,m[i],x,y);
        if(c%d!=0) return -1;

        x=q_pow(x,c/d,m[i]/d);

        ans=(ans+x*M);
        M=M*(m[i]/d);
        ans=(ans%M+M)%M;
    }
    return (ans%M+M)%M;
}




int main(void)
{

    scanf("%d",&k);

    for(int i=1;i<=k;i++)
    {
        scanf("%lld%lld",&m[i],&a[i]);
    }

    printf("%lld\n",china());
    return 0;
}