解线性同余方程组-模板

int n;
LL a[K],b[K];
LL gcd(LL a,LL b)
{
return b==0?a:gcd(b,a%b);
}
void exgcd(LL a,LL b,LL &d,LL &x,LL &y)
{
if(b==0)
{
x=1;y=0;d=a;
return;
}
exgcd(b,a%b,d,y,x);
y-=x*(a/b);
}
LL slove(void)//x=b[i](mod a[i])
{
LL ta=a[0],tb=b[0];
bool flag=true;
for(int i=1;i<n;i++)
{
LL xa=ta,xb=a[i],c=b[i]-tb,d,x,y;
exgcd(xa,xb,d,x,y);
if(c%d)
{
flag=false;break;
}
LL tm=xb/d;
x=(x*(c/d)%tm+tm)%tm;
tb=ta*x+tb;
ta=ta/d*a[i];
}
if(!flag)
return -1;
return tb;
}

解出的是最小正整数解