exgcd&&中国剩余定理

拓展欧几里得模板

$求解gcd(a,b)=a\ast x+b\ast y，x与y的值$
$由gcd(a,b)=$ $gcd(b,a$%$b$$)$
$a\ast x+b\ast y=b\ast x+(a-\lfloor \frac{a}{b}\rfloor \ast b)\ast y$
$a\ast x+b\ast y=a\ast y+b\ast (x-a/b\ast y)$
$发现a\ast x对应a\ast y只有一项不相等，直接交换即可$
$a\ast x+b\ast y=a\ast x+b\ast (y-\lfloor \frac{a}{b}\rfloor \ast x)$
a*x+b*y=c存在解的充分必要条件c%__gcd(a,b)==0
$a\ast x+b\ast y=c$等价于$a\ast x\equiv c(\mathrm{m}\mathrm{o}\mathrm{d}b)$有解的充分必要条件c%gcd(a,b)=0
例：求解$a\ast x\equiv 1(\mathrm{m}\mathrm{o}\mathrm{d}b)$即为：$a\ast x+b\ast y\equiv 1(\mathrm{m}\mathrm{o}\mathrm{d}b)$的解
当1%__gcd(a,b)==0时可直接用模板解。

int exgcd(int a,int b,int &x,int &y){
	if(b==0){
		x=1,y=0;
		return a;
	}
	int d=exgcd(b,a%b,y,x);
	y-=a/b*x;
	return d;
}

中国剩余定理

$$\{\begin{array}{r}x\equiv a_1(\mathrm{m}\mathrm{o}\mathrm{d}{m}_1)\\ x\equiv a_2(\mathrm{m}\mathrm{o}\mathrm{d}{m}_2)\\ ...\\ x\equiv a_n(\mathrm{m}\mathrm{o}\mathrm{d}{m}_n)\end{array}(m_1...m_n两两互质)$$
设$M={\prod }_{i=1}^n{m}_i，t_k$为$\frac{M}{m_k}\ast t_k\equiv 1(\mathrm{m}\mathrm{o}\mathrm{d}{m}_k)$的最小非负整数解
其中$\forall {}_{i!=k}\frac{M}{m_i}\ast t_i\equiv 0(\mathrm{m}\mathrm{o}\mathrm{d}{m}_i)$
两边$\ast a[i]$有
$$\{\begin{array}{r}{a}_k\ast \frac{M}{m_k}\ast t_i\equiv a_k(\mathrm{m}\mathrm{o}\mathrm{d}{m}_k)\\ \forall {}_{i!=k}{a}_i\ast \frac{M}{m_k}\ast t_i\equiv 0(\mathrm{m}\mathrm{o}\mathrm{d}{m}_i)\end{array}$$
则有一个解为$x={\sum }_{i=1}^n{a}_i\ast \frac{M}{m_i}\ast t_i$
通解为$x+k\ast M(k\in Z)$
P3868 [TJOI2009] 猜数字

#include <iostream>
using namespace std;
using ll = long long ;
int n;
ll read(){
    char ch=getchar();
    int f=1;
    while(!isdigit(ch)){
        if(ch=='-')f=-1;
        ch=getchar();
    }
    ll x=0;
    while(isdigit(ch)){
        x=(x<<3)+(x<<1)+(ch^48);
        ch=getchar();
    }
    return f*x;
}
ll a[20],b[20];
ll mul(ll a,ll b,ll p){///龟速乘
    ll r=0;
    while(b){
        if(1&b)r=(r+a)%p;
        a=(a+a)%p;
        b>>=1;
    }
    return r;
}
/**
    gcd(a,b)=gcd(b,a%b)
    a*x+b*y=b*x+(a-a/b*b)*y
    a*x+b*y=a*y+b*(x-a/b*y)
*/
void exgcd(ll a,ll b,ll &x,ll &y){
    if(b==0){
        x=1,y=0;
        return ;
    }
    exgcd(b,a%b,y,x);
    y-=a/b*x;
}
ll china(){
    ll M=1;
    for(int i=1;i<=n;i++) M*=b[i];
    ll ans=0;
    ll x,y;
    for(int i=1;i<=n;i++)
    {
        ll tp=M/b[i];
        exgcd(tp,b[i],x,y);
        x=(x%b[i]+b[i])%b[i];
        ans=(ans+mul(mul(tp,x,M),a[i],M))%M;
    }
    return (ans+M)%M;
}
int main()
{
    n=read();
    for(int i=1;i<=n;i++) a[i]=read();
    for(int i=1;i<=n;i++) b[i]=read(),a[i]=(a[i]%b[i]+b[i])%b[i];
    printf("%lld",china());
    return 0;
}