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

>Link

luogu P4777

---

>Description

求解一元线性同余方程组$\{\begin{array}{c}x\equiv b_1(mod{a}_1)\\ ......\\ x\equiv b_n(mod{a}_n)\end{array}$
不保证$a$两两互质

---

>解题思路

由于模数不保证两两互质，所以我们就不能用中国剩余定理的方法exgcd了
要用扩展中国剩余定理

参考自：建议看看这个巨爷的笔记QAQ，这里做个总结
我们先考虑两个式子$\{\begin{array}{c}x\equiv b_1(mod{a}_1)\\ x\equiv b_2(mod{a}_2)\end{array}$
可以转换成$\{\begin{array}{c}x=b_1+a_1\ast k_1\\ x=b_2+a_2\ast k_2\end{array}$
联立得 $b_1+a_1\ast k_1=b_2+a_2\ast k_2$
移项得 $a_1\ast k_1-a_2\ast k_2=b_2-b_1$
这个式子就可以用exgcd来解决了，但是得满足$gcd(a_1,a_2)|b_2-b_1$这个前提（回顾下exgcd：$a_1\ast k_1-a_2\ast k_2=gcd(a_1,a_2)$），最后求出来两边再乘上$(b_2-b_1)/gcd(a_1,a_2)$就行了

我们设用exgcd求出来的其中两个解为 $K_1$和$K_2$，即 $a_1\ast K_1-a_2\ast K_2=gcd(a_1,a_2)$
按照前面说的，两边乘上$\frac{(b_2-b_1)}{gcd(a_1,a_2)}$，得到 $k_1=K_1\ast \frac{(b_2-b_1)}{gcd(a_1,a_2)}$，$k_2=K_2\ast \frac{(b_2-b_1)}{gcd(a_1,a_2)}$
现在我们就可以求出一个$x$了

可是显然不只有一个$x$
设求出来的两个解为$x_0$和$x_1$，满足$x_0<x_1$
两个解都满足同余方程，所以$\{\begin{array}{c}{x}_0\equiv b_1(mod{a}_1)\\ x_1\equiv b_1(mod{a}_1)\end{array}$$\{\begin{array}{c}{x}_0\equiv b_2(mod{a}_2)\\ x_1\equiv b_2(mod{a}_2)\end{array}$
设存在 $s_1$和$s_2$，使得 $x_1=x_0+a_1\ast s_1$， $x_1=x_0+a_2\ast s_2$
移项得 $x_1-x_0=a_1\ast s_1$， $x_1-x_0=a_2\ast s_2$，所以 $a_1\ast s_1=a_2\ast s_2$
设 $d_1=\frac{a_1}{gcd(a_1,a_2)}$，$d_2=\frac{a_2}{gcd(a_1,a_2)}$（我们知道$gcd(d_1,d_2)=1$），得 $d_1\ast gcd(a_1,a_2)\ast s_1=d_2\ast gcd(a_1,a_2)\ast s_2$，即 $d_1\ast s_1=d_2\ast s_2$
所以我们得到 $d_1|s_2$，
$d_1\ast d_2|s_2\ast d_2$
$d_1\ast d_2\ast gcd(a_1,a_2)|s_2\ast d_2\ast gcd(a_1,a_2)$
$lcm(a_1,a_2)|x_1-x_0$
所以$x_1-x_0$一定是$lcm(a_1,a_2)$的倍数！！！求出一个$x$，其它加上任意倍$lcm(a_1,a_2)$就行了

so

$M$表示当前所有$a$的$lcm$，$X$表示当前解
枚举到当前数，解表示为 $X+t\ast M$
满足 $X+t\ast M\equiv b_i(mod{a}_i)$，移项得 $t\ast M\equiv b_i-X(mod{a}_i)$

---

>代码

//还是不是很懂 需要再看一下QAQ

#include <iostream>
#include <cstdio>
#include <cstring>
#include <algorithm>
#define N 100010
#define LL long long
using namespace std;

int n;
LL a[N], b[N], M, ans;

LL mul (LL aa, LL bb, LL Mod)
{
	LL ret = 0;
	aa %= Mod;
	for (; bb; bb >>= 1, aa = (aa + aa) % Mod)
	  if (bb & 1) ret = (ret + aa) % Mod;
	return ret;
}
LL ex_gcd (LL aa, LL bb, LL &x, LL &y)
{
	if (!bb) {x = 1, y = 0; return aa;}
	LL X = x, Y = y;
	LL ret = ex_gcd (bb, aa % bb, X, Y);
	x = Y;
	y = X - (aa / bb) * Y;
	return ret;
}

int main()
{
	scanf ("%d", &n);
	for (int i = 1; i <= n; i++)
	  scanf ("%lld%lld", &a[i], &b[i]);
	M = a[1], ans = b[1];
	for (int i = 2; i <= n; i++)
	{
		LL x = 0, y = 0;
		LL g = ex_gcd (M, a[i], x, y);
		//if ((b[i] - ans) % g != 0) nosolution~
		/*x = x * ((b[i] - ans) / g);
		LL t = a[i] / g;
		x = (x % t + t) % t;*/
		LL t = (b[i] - ans % a[i] + a[i]) % a[i];
		x = mul (x, t / g, a[i] / g);
		ans = ans + x * M;
		M = M * (a[i] / g);
		ans = (ans % M + M) % M;
	}
	ans = (ans % M + M) % M;
	printf ("%lld", ans);
	return 0;
}