【洛谷P4777】拓展中国剩余定理(Excrt)

$link$

分析：

前置知识：$exgcd$, $crt$
$exgcd$用来解同余方程 $crt$用来解同余方程组 其中$p_i$互质
$excrt$也可以解同余方程组 但$p_i$可以不互质

先不用管合并所有方程组 考虑合并两个
$$\{\begin{array}{l}x\equiv A(modP)\\ x\equiv a_i(mod{p}_i)\end{array}$$
化成正常形式：
$$\{\begin{array}{l}x=A+P{k}_1\\ x=a_i+p_i{k}_2\end{array}$$
上下相减 得
$$P{k}_1-p_i{k}_2=a_i-A$$
用$exgcd$可求出$k_1=p_{i}l+b$ $(l$为任意整数$)$
代入回去就可以得出 $x=A+P{p}_{i}l+Pb$
然后这是$p_i$和$P$互质的时候 也就是普通$crt$

再考虑$p_i$和$P$不互质 先除去$gcd(p_i,P)$
即
$$\frac{P}{gcd(p_i,P)}{k}_1-\frac{p_i}{gcd(p_i,P)}{k}_2=\frac{a_i-A}{gcd(p_i,P)}$$
同理 $exgcd$得出
$$k_1=\frac{p_i}{gcd(p_i,P)}l+b$$
然后同样代入回去
$$x=A+lcm(p_i,P)l+Pb$$
复杂度$O(nlogn)$
如果防溢出 还要用龟速乘 求最小非负数解每次要处理下

CODE：

#include<iostream>
#include<cmath>
#include<cstring>
#include<cstdio>
#include<algorithm>
using namespace std;
typedef long long ll;
const int N=1e6+5;
int n;
ll a[N],p[N],x,y;
ll mul(ll a,ll k,ll Mod)
{
	ll res=0;
	while(k)
	{
		if(k&1) (res+=a)%=Mod;
		k>>=1;
		(a+=a)%=Mod;
	}
	return res;
} 
ll exgcd(ll a,ll b,ll &x,ll &y)
{
	if(b==0)
	{
		x=1;
		y=0;
		return a;
	}
	ll res=exgcd(b,a%b,x,y);
	ll k=x;
	x=y;
	y=k-a/b*y;
	return res;
}
ll excrt()
{
	ll Mod=p[1],ans=a[1];
	for(int i=2;i<=n;i++)
	{
		ll k=((a[i]-ans)%p[i]+p[i])%p[i];
		ll gcd=exgcd(Mod,p[i],x,y);
		ll qwq=p[i]/gcd;
		x=mul(x,k/gcd,qwq);
		ans+=x*Mod;
		Mod*=qwq;
		ans=(ans%Mod+Mod)%Mod;
	}
	return ans;
}
int main(){
	scanf("%d",&n);
	for(int i=1;i<=n;i++)
		scanf("%lld%lld",&p[i],&a[i]);
	printf("%lld",excrt());
	
    return 0;
}