CRT合并

对于两个同余方程组：

x mod a1 = b1 (1)

x mod a2 = b2 (2)

我们令x = a1 * x0 + b1;

带入(2)式，即(a1* x0 + b1) mod a2 = b2，将b1移项，即：

a1 * x0 mod a2 = b2 – b1 (3)

显然，我们可以使用ex_gcd求出它的一个可行解x1，而实际上x0满足通项：

x0 = x1 + a2 / gcd(a1, a2) * k

带回(3)式可得

a1 * x1 + lcm(a1, a2) * k + P = 0

其中a1 * x1 + P是个常数，也就是说，我们将两个同余方程组合并成了一个方程组：

x mod lcm(a1, a2) = a1 * x1 + P

所以，我们只需不断合并，就可以将原方程组化为一个方程，就可以直接求出一组可行解了。

void CRT(int64 a1, int64 b1, int64 a2, int64 b2, int64 &a, int64 &b) {
	if (b1 > b2) swap(b1, b2), swap(a1, a2);
	int64 x, y, g;
	g = gcd(a1, a2); a1 /= g, a2 /= g;
	ex_gcd(a1, a2, x, y);

	a = a1 * a2 / g;
	x = x * g % a2 * (b2 - b1) / g % a;
	b = ((b1 + a1 * x) % a + a) % a;
}

有一些地方比较容易爆long long，需要小心。