「营业日志 2021.3.14」代数整数上线性筛的两个实例

在两年前我加强出了 这么 一道题目，当时得到的复杂度是 $\Theta (k^2\log n)$ 的，因为最后要做 $n^2$ 个快速幂。

经过了最近一段时间的摸鱼，突然发现这个复杂度可以再稍微降一点。

限于本人水平有限，接下来会列出一些我不会证明的结果，但这些结论应当是容易在教科书或网上查询确认的。

$d=4$（Gaussian 整数）

对于 $d=4$ 的情况，我们需要对于 $|a|+|b|\le k$，求出全体 $(a+b\mathrm{i}{)}^n$。

我们的想法和 $\mathbb{Z}$ 上的情形是类似的，我们如果有 $z=xy$，那么 $z^n=x^n{y}^n$，也就不用对 $z$ 再额外算一次快速幂了。

由于 $\mathbb{Z}[\mathrm{i}]$ 是欧几里得整环，所以它是 UFD。我们只需要对所有的不可约元来算一下快速幂。而 $\mathbb{Z}[\mathrm{i}]$ 上的不可约元只有以下三种：

1. $1+\mathrm{i}$
2. $a+b\mathrm{i}$，其中 $a^2+b^2$ 为 $4k+1$ 形式的素数。
3. $4k+3$ 形式的素数。

那么在我们筛的范围内，1 中只有 $1$ 个，3 中只有 $\Theta (k/\log k)$ 个，2 的比较多，考虑所有 $p\le k^2$ 的满足该形式的素数都会统计在内，但是每个 $p$ 只对应于一组 $a,b$ 和 $b,a$，因此大概只有 $\Theta (k^2/\log k)$ 个不可约元。

那么我们的算法大概是 $\Theta (k^2{\log }_{k}n)$ 的。

至于怎么线性筛，我们照猫画虎，将全体 $a+b\mathrm{i}$ 按照 $a^2+b^2$ 从小到大排列，然后再按照这个顺序从小到大枚举不可约元就行了。

理论上需要取全体圆内整点 $a^2+b^2\le k^2$。但是因为这里的图形非常特殊，如果 $z$ 在圆内但不在正方形内的话，它是没有贡献的，因为它只能乘以 $\pm 1,\pm \mathrm{i}$（乘以更大的数就在圆外了），但是这四个点都在正方形外。进一步地，为了减少计算量可以只筛 $a>0,b\ge 0$的部分。

$d=6$（Eisenstein 整数）

对于 $d=6$ 的情况，我们需要对于 $|z|\le k$，求出全体 $z^n$，其中 $z=a+b\omega \in \mathbb{Z}[\omega ]$，其中 $\omega$ 是 $3$ 次单位根。

其实 $\mathbb{Z}[\omega ]$ 就是 $\mathbb{Q}[\sqrt{-3}]$ 上的代数整数，可以打网格证明是欧几里得整环，也就是 UFD。

$\mathbb{Z}[\omega ]$ 上的不可约元的分类也表明其中 $|p|\le k$ 的不可约元数量是 $\Theta (k^2/\log k)$ 的。因此复杂度也是 $\Theta (k^2{\log }_{k}n)$。

人到中年，不太想写代码了，有兴趣再补，哈哈.jpg