扩展中国剩余定理（excrt）
本来应该先学中国剩余定理的。但是有了扩展中国剩余定理，朴素的 CRT 就没用了。
扩展中国剩余定理用来求解如下形式的同余方程组：
\[\begin{cases} x \equiv a_1\ ({\rm mod}\ b_1) \\ x\equiv a_2\ ({\rm mod}\ b_2) \\ ... \\ x \equiv a_n\ ({\rm mod}\ b_n)\end{cases} \]扩展中国剩余定理的基本思想是合并，通过 \(n - 1\) 次合并，将一个大的同余方程组合并成一个同余方程。
假设现在有两个同余方程：
\[\begin{cases} x \equiv a\ ({\rm mod}\ A) \\ x\equiv b\ ({\rm mod}\ B) \end{cases} \]现在要将他们合并。首先转化成不定方程：
\[\Rightarrow \begin{cases} A | (x - a) \\ B | (x - b) \end{cases} \]\[\Rightarrow \begin{cases} x - a = Ak_1 \\ x - b = Bk_2 \end{cases} \]\[\Rightarrow \begin{cases} x = a + Ak_1 \\ x = b + Bk_2 \end{cases} \]\[\Rightarrow a + Ak_1 = b + Bk_2 \]\[\Rightarrow Ak_1 - Bk_2 = b - a \]成功转化成了系数为 \((A, -B, b - a)\) 的不定方程，使用 exgcd 求出他的一个根。因此转化成了一个同余方程：\(x \equiv Ak_1 + a (\bmod \ \text{lcm}(A, B))\)。合并完成。
// 合并 x = a(mod A), x = b(mod B) 两个方程
// 返回的是新的 a', b'，满足 x = a'(mod b')
PII merge(int a, int A, int b, int B) {
int k1, k2; int d = exgcd(k1, k2, A, B);
k1 = k1 * (b - a) / d;
int p = A / d * B; return {(A * k1 + a) % p, p};
}
bonus：
- 如果 \(x\) 的系数不为 \(1\)。
也就是 P4774 [NOI2018] 屠龙勇士。
求解形如：
\[\begin{cases} p_1 x \equiv a_1\ ({\rm mod}\ b_1) \\ p_2 x\equiv a_2\ ({\rm mod}\ b_2) \\ ... \\ p_n x \equiv a_n\ ({\rm mod}\ b_n)\end{cases} \]Excrt 因为 \(x\) 的系数是一，因此可以直接联立两个不定方程。也尝试将这个东西转化成不定方程的形式。假设现在需要合并的两个同余方程是：
\[\begin{cases} px \equiv a (\bmod \ b) \\ Px \equiv A(\bmod \ B)\end{cases} \]\[\Rightarrow \begin{cases} b | (px - a) \\ B | (Px - A) \end{cases} \]\[\Rightarrow \begin{cases} px - a = k_1 b \\ Px - A = k_2 B \end{cases} \]\[\Rightarrow \begin{cases} px - k_1 b = a \\ Px - k_2 B =A\end{cases} \]然后发现两个 \(x\) 的系数不同，不能直接合并了。而这两个柿子两边又不能同时除以 \(p\) 或者 \(P\)，因为不保证逆元存在。这就非常难搞。
一个神奇的思路是直接解出两个方程。以第一个方程为例，方程中只有两个未知数 \(x\) 和 \(-k_1\)，可以解出一个特解 \(x_0\)。那么所有 \(x\) 就可以表示成：
\[x = x_0 + \dfrac{b}{(p, b)} \times \alpha \]同理解第二个方程，可以得到
\[x = x_1 + \dfrac{B}{(P, B)} \times \beta \]我们惊奇的发现这两个 \(x\) 的系数相同了。所以可以合并一下：
\[x_0 + \dfrac{b}{(p, b)} \times \alpha = x_1 + \dfrac{B}{(P, B)} \times \beta \]里面只有 \(\alpha, \beta\) 两个未知数，解出他们两个就可以得到 \(x\)。
- 扩展中国定理进行模数非质数的合并
即 古代猪文。
求 \(\dbinom{n}{m} \bmod \ 999911658\) 的值。
将 \(999911658\) 质因数分解得到：\(999911658 = 2 \times 3 \times 4679 \times 35617\)。
因此可以对 \(2, 3, 4679, 35617\) 分别做一遍 \(\text{Lucas}\)，得到下面的同余方程：
\[\begin{cases} x \equiv a_1(\bmod \ 2) \\ x \equiv a_2(\bmod \ 3) \\ x \equiv a_3(\bmod \ 4679) \\ x \equiv a_4(\bmod \ 35617) \\ \end{cases}\]可以直接用 excrt 合并一下。
另外一个应用是扩展卢卡斯。其基本思路也是将模数拆成若干质因数的次方，计算后 excrt 合并。
标签：begin,end,定理,方程,笔记,Excrt,cases,equiv,mod From： https://www.cnblogs.com/LcyRegister/p/17922691.html