【学习笔记】快速 GCD 算法

目标：在 $\mathcal{O}(V)$ 预处理的前提下，做到 $\mathcal{O}(1)$ 求算 $V$ 以内任意数对的 $\gcd$ 。

Lemma. 对于任意正整数 $n$，可以找到三个正整数 $a,b,c(a⩽b⩽c)$ 使得 $abc=n$，满足 $a,b⩽\sqrt{n}\wedge (c⩽\sqrt{n}\vee c\in \text{Prime})$ 。

Proof. 当 $n$ 存在超过 $\sqrt{n}$ 的质因子时，令 $c$ 为该质因子即成立。因此只需证 $n$ 只含小质因数的情况。

构造其分配方案：取 $p$ 为 $n$ 最小质因数，由归纳法，将 $\frac{n}{p}$ 划分为 $a,b,c(a⩽b⩽c)$ 。显然 $a⩽\sqrt[3]{\frac{n}{p}}$ 。由归纳假设 $b⩽\sqrt{\frac{n}{p}}$，且 $c⩽\sqrt{\frac{n}{p}}\vee c\in \text{Prime}$ 。但是 $n$ 的质因数都不超过 $\sqrt{n}$，且 $c\mid n$，所以总是有 $c⩽\sqrt{n}$ 。对 $p$ 的大小进行分类讨论：

• 若 $p⩽\sqrt[4]{n}$，则 $ap⩽\sqrt[3]{\frac{n}{p}}\cdot p=\sqrt[3]{n{p}^2}⩽\sqrt{n}$，于是 $ap,b,c$ 就是合法的划分。
• 否则，若 $a⩾p$ 则 $n=pabc⩾p^4>n$，矛盾，故 $a<p$ 。又因 $p$ 是最小质因数 $\Rightarrow a=1$ 。显然最小质因子 $p⩽\sqrt{n}$，所以 $p,b,c$ 就是合法的划分。

Q.E.D.

顺便给出了求解方法：线性筛的时候知道了 $p$，直接叠加到 $\frac{n}{p}$ 划分中的最小数上就行。

知道了这个之后呢？那么原来的 $\gcd (n,v)$ 可以变为 $\gcd (abc,v)$ 。预处理 $\mathcal{O}(\sqrt{n}\times \sqrt{n})=\mathcal{O}(n)$ 的 $\gcd$ 表，可以直接查询 $\gcd (a,v)=\gcd (a,v\mathrm{mod}a)$，因为与质数求 $\gcd$ 是平凡的。

而 $\gcd$ 是质因子指数取 $\min$，相当于两边同时除以 $d$ 使得二者不再能同时除以任意一个数。所以每次 $v$ 跟 $a$ 或 $b$ 求 $\gcd$ 后，直接将 $v$ 除以这么多就可以了。具体可见代码。

所以这是 $\mathcal{O}(V)$ 预处理之后 $\mathcal{O}(3)$ 查询的方法。

---

上面的 $\text{Lemma}$ 似乎很需要些猜测才能得到。不妨试试 $\text{HALF-GCD}$ 的方法：将每个数视为 $x=2$ 的关于 $x$ 的多项式即可。有复杂度
$$\mathcal{T}(n)=2\mathcal{T}\left(\frac{n}{2}\right)+\mathcal{O}(1)$$

其中 $n$ 是多项式的长度。不难发现，若我们预处理 $V^{\alpha }$ 以内的数对的 $\gcd$ 其中 $\alpha$ 是 $(0,1)$ 内常数，则递归层数是常数 ${\log }_2({\alpha }^{-1})$ ，故这也是 $\mathcal{O}(1)$ 的查询。

据 $\text{skip2004}$ 所说，处理 $V^{1/3}$ 内的数对的 $\gcd$ 是跑得很快的。