[洛谷P4777]扩展CRT

题目

传送门 to luogu

思路

板题讲什么思路啊！ 首先重新定义一下变量名：

$$\{\begin{array}{l}x\equiv c_1(\mathrm{m}\mathrm{o}\mathrm{d}{m}_1)\\ x\equiv c_2(\mathrm{m}\mathrm{o}\mathrm{d}{m}_2)\\ x\equiv c_3(\mathrm{m}\mathrm{o}\mathrm{d}{m}_3)\\ \dots \\ x\equiv c_n(\mathrm{m}\mathrm{o}\mathrm{d}{m}_n)\end{array}$$

普通版

就是 $m$ 两两互质的情况。

考虑像拉格朗日插值一样一项一项构造：

第 $i$ 项用来满足第 $i$ 个方程，必然有 ${\omega }_i={\prod }_{j\ne i}{m}_j$ 这一项。考虑将 $c_i$ 乘进去，那必须有一个 $t$ 使得 $t\cdot {\omega }_i\equiv 1(\mathrm{m}\mathrm{o}\mathrm{d}{m}_i)$ 。这不就是求个逆元吗？

记 ${\omega }_i^{-1}$ 为 ${\omega }_i$ 在模 $m_i$ 意义下的逆元，则有：$$x=\sum _{i=1}^n{c}_i{\omega }_i{\omega }_i^{-1}$$

还可以优化，记 $M={\prod }_{i=1}^n{m}_i$ ，则 $x={\sum }_{i=1}^n{c}_i\frac{M}{m_i}{\left(\frac{M}{m_i}\right)}^{-1}$ 。代码就不贴了吧？

扩展版

要是 $m_i$ 不是两两互质怎么办？

要是你想到了通过将 $x$ 扩倍啥的让它们不再互质，恭喜你 你完事儿了

正解是 合并两个方程，直到只有一个方程。假定这两个方程分别是 $x\equiv c(\mathrm{m}\mathrm{o}\mathrm{d}m)$和$x\equiv C(\mathrm{m}\mathrm{o}\mathrm{d}M)$ 。

发现 $x\equiv c(\mathrm{m}\mathrm{o}\mathrm{d}m)$ 等价于 $x=km+c(k\in \mathbb{Z})$ ，尝试对这两个方程进行等价变换：$$x=aM+C=bm+c$$

后面两项相等，等价于：$$bm-aM=C-c$$

因为 $m,c,M,C$ 都是已知，假设得到一个特解 $\{\begin{array}{l}a=a_0\\ b=b_0\end{array}$ ，有通解公式：$$\begin{gathered} a=a_0+k\times \frac{m}{d} \\ \iff aM=a_{0}M+k\times \frac{mM}{d} \end{gathered}$$

毕竟 $M\ne 0$ 嘛。将该式带回 $x=aM+C$ 中得到：$$x=a_{0}M+C+k\times \frac{mM}{d}$$

既然 $k$ 是变量，写成同余形式：$$x\equiv a_{0}M+C(\mathrm{m}\mathrm{o}\mathrm{d}\frac{mM}{d})$$

不要看到有分数线就觉得是小数，$\frac{mM}{d}$ 其实就是 $n$ 和 $m$ 的最小公倍数，肯定是整数。

对了，$a_0$ 怎么算？在 $bm-aM=C-c$ 的方程中，好像只有 $a$ 和 $b$ 是未知数—— $\mathrm{e}\mathrm{x}\mathrm{g}\mathrm{c}\mathrm{d}$ 即可。要是 $a$ 无解，$x$ 也无解。

代码

#include <cstdio>
#include <iostream>
#include <vector>
using namespace std;
typedef __int128 LL; // 防爆措施
inline LL readint(){
	LL a = 0; char f = 1, c = getchar();
	if(c == EOF) return EOF;
	while(c < '0' or c > '9')
		(c == '-' and (f = -1)), c = getchar();
	while('0' <= c and c <= '9')
		a = (a<<3)+(a<<1)+(c^48), c = getchar();
	return a*f;
}

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

LL exCRT(int n, LL c[], LL m[]){
	LL M = m[1], C = c[1];
	for(int i=2; i<=n; ++i){
		LL a, b, r = C-c[i];
		// bm-aM = C-c
		LL d = exgcd(m[i], M, b, a);
		if(r%d != 0)// no solution
			return -1;
		a = r/d*(-a); // 我们得到的是-a
		C = a*M+C; M = M/d*m[i]; C %= M;
	}
	return (C+M)%M;
}

const int MaxN = 100000;
int n; LL c[MaxN+5], m[MaxN+5];
int main(){
	while(~scanf("%d", &n)){
		for(int i=1; i<=n; ++i)
			m[i] = readint(), c[i] = readint();
		printf("%lld\n",exCRT(n, c, m));
	}
    return 0;
}