======二次域及有理逼近相关问题======

=====引言=====

二次域是初等数论中篇幅很大的一部分，主要集中在不定方程篇章。

这里采用了较高一点的观点统一了这些简单的二次不定方程。希望读者在阅读之前，已经提前了解了简单的代数学。

=====绪论=====

====定义====

首项系数为1的整系数二次多项式$x^2+px+q=0$的零点是：

$$\frac{-p±\sqrt{p^2-4q}}{2}$$

称为“含有根号d的二次整数”，全体记作二次整环$Z(\sqrt{d})$，对于加减乘封闭。不同的d对应于不同的整环。普通的整数环是每一个二次整环的理想。

第一种情况：对于所有的d，$a+b\sqrt{d}$一定是二次整数。

第二种情况：当d模4余1，a与b是奇数的时候，$\frac{a+b\sqrt{d}}{2}$也是二次整数。因为这种情况也是首系数为1的整系数多项式的零点：

$$x^2-ax+\frac{a^2-db^2}{4}=0$$

奇数的一半称半整数。两个半整数配上除以4余1的d开二次根号，也是二次整数。

以上的d全部可正可负。当d为正时就是普通意义的二次根号，当d为负的时候可以理解成对绝对值开根号，并乘以虚数单位i。

二次整数有两个线性无关的分量，因此二次整数是二维的。

同类二次整数的比构成二次有理数，全体构成二次域$Q(\sqrt{d})$。

若d为正，集合中所有数均为实数，称为实二次整环或者实二次域。若d为负，集合中除了一般的有理数以外全部不是实数，称为虚二次整环或虚二次域。

====范数与单位数====

对于同一个整系数二次方程，有两个根。如果它们不是一般的有理数，那么它们在形式上只在二次根号前相差一个正负号。

如果两个二次有理数只在二次根号之前相差正负号，称它们互为**共轭**关系。因为一般的有理数在二次根号前面的系数是0，因此一般的有理数与它自身为共轭关系。

显然，在虚二次域中，某数的共轭的概念，与复数共轭的概念一致。但是在实二次域中这两个概念不一致。

在二次域中，由加减乘除（非0）四则运算产生的等式，无法区分共轭关系。也就是说，在等式中将每一个数换成它的共轭数，即将每一个二次根号的符号改变，等式仍然成立。

二次有理数与它的共轭的和称为**迹**。某数的迹就是它的有理数部分的2倍，形式简单，因此很少研究迹。

二次有理数与它的共轭的积称为**范数**：

$$N(a+b\sqrt{d})=a^2-db^2$$

显然，在虚二次域中，范数的概念，与复数的模的平方的概念一致。但是在实二次域中这两个概念不一致。由于d不含平方因子，不可能是平方数，因此只有0的范数是0。

范数具有保持乘法和除法（非0）的良好性质。

$$N(a_1+b_1\sqrt{d})N(a_2+b_2\sqrt{d})=N((a_1+b_1\sqrt{d})(a_2+b_2\sqrt{d}))$$

$$\frac{N(a_1+b_1\sqrt{d})}{N(a_2+b_2\sqrt{d})}=N\left(\frac{a_1+b_1\sqrt{d}}{a_2+b_2\sqrt{d}}\right)$$

一个二次有理数与它的共轭相乘为这个数的范数，因此它的倒数就是它的共轭与范数之比。

$$a+b\sqrt{d}=\frac{a-b\sqrt{d}}{N(a+b\sqrt{d})}$$ 

如果一个二次整数的倒数还是二次整数，称这个二次整数为**单位数**。二次整数是单位数的充要条件是它的范数为1或-1。

单位数对于乘法封闭，构成单位群。有一个核心位置的定理（证明极难）：

**狄利克雷单位定理：数域的单位群是有限生成阿贝尔群。**

狄利克雷单位定理表明：单位群维数有限，存在一组基。所有的单位数可以由基的乘积表示。这组基（不含1和-1）称为**基本单位数**。

=====虚二次域=====

====简介====

在虚二次域中，仅当d为-1和-3的时候，存在除了1和-1以外的单位数。当d为负数且不为-1或-3的时候，单位数只有1和-1。

当d为-1的时候，单位数有4个：1，-1，i，-i。当d为-3的时候，单位数有6个：1，-1，$\frac{1+\sqrt{3}i}{2}$，$\frac{1-\sqrt{3}i}{2}$，$\frac{-1+\sqrt{3}i}{2}$，$\frac{-1-\sqrt{3}i}{2}$。

在虚二次域中，仅当d为-1和-3时，存在基本单位数i和$\frac{1+\sqrt{3}i}{2}$。其他情况不存在1和-1以外的其他单位数，也就不存在基本单位数。

因此，两个整环$Z(i)$和$Z(\sqrt{3}i)$是特殊的整环，称为高斯整环和艾森斯坦整环。它们直观上分别构成复平面上正方形点阵和正六边形点阵（正三角形格点），研究虚二次域的时候最经常用到这两个整环。

虚二次域中对范数的研究可以转化为**椭圆上整点问题**，有名的“圆上整点问题”可以转化为对d为-1的虚二次域的研究。

====三种整环====

虚二次域中，我们经常研究的整环占一小部分，原因是虚二次域中很多优秀的性质无法满足，使得研究有很大困难。

有三种整环的概念：

**Euclid整环**（欧几里得）：满足**辗转相除法**的整环。

**主理想整环**：每一个理想都是主理想的整环，即满足**Bezout定理**（裴蜀定理）。

**唯一分解整环**：每个元素的非相伴分解都唯一的整环，满足**唯一分解定理**。

三个概念是层层嵌套包含的关系，唯一分解整环在最外面，欧几里得整环在最里面。欧几里得整环一定是主理想整环，主理想整环一定是唯一分解整环，而反之则不然。因此三个定理也有层层递推的关系。

虽然唯一分解整环不一定是主理想整环，例如在取模多项式整环中可以找到反例，但是在二次域中，这两个概念是重合的，即二次域的主理想整环与唯一分解整环范畴重合。因此，二次域只分为辗转相除和唯一分解两种特殊情形。

在虚二次域中，只有-1、-2、-3、-7和-11对应的虚二次整环是Euclid整环，其余均不满足辗转相除法。

在实二次域中，只有2、3、5、6、7、11、13、17、19、21、29、33、37、41、57和73，共16个整环是Euclid整环。

对于二次域，有很重要的概念叫类数。理想的全体除以理想构成的群，得到商群的大小就称为类数。类数为1，说明相应的整环是主理想整环。

Gauss猜想有无穷个类数为1的实二次域，这个问题至今没有得到解决——关于实二次域的大多数此类研究进展都很慢。

已经得到解决的是，虚二次域中，加上上面的5个，只有-19、-43、-67和-163也是主理想整环。

好在之前的嵌套关系成立。我们只需知道高斯整环（-1）和艾森斯坦整环（-3）都是Euclid整环，满足辗转相除法和唯一分解定理就够了。

参见OEIS：

[[http://oeis.org/A048981|A048981 Squarefree values of n for which the quadratic field Q(sqrt(n)) is norm-Euclidean]]

[[http://oeis.org/A003172|A003172 Q(sqrt(n)) is a unique factorization domain (or simple quadratic field)]]

====相伴与唯一分解====

如果一个二次整数乘一个单位数得到另一个二次整数，那么这两个二次整数是**相伴**关系。

唯一分解定理一定要考虑相伴关系才有可能成立。例如，若不考虑相伴关系，由于-1是单位数，整数不满足唯一分解：

$$10=2*5=(-2)*(-5)$$

我们必须在相伴这个等价关系构成的诸多等价类中，为每个类指定一个数作为这个类的代表，即定义**本原数**，才可能有唯一分解。

例如在上面的例子中，如果指定2和5为本原数，那么-2和-5就不是本原数，此时10的分解才变得唯一了。

本原数的规定是人为的，即如果定义-2和5、2和-5或者-2和-5为本原数，在唯一分解的角度不会引起矛盾。一般会根据实际问题的研究方便定义本原数。例如，如果我们习惯于在正整数范畴研究问题，那么将正整数定义为本原数即可。

我们看到，事实上只需为所有的素数（在唯一分解前提下与不可约数等价）定义本原数就够了，其他的非素数的本原数定义必然由素数的本原数定义合成。

====Gauss整数====

一般将$Q(i)$称为高斯域，相应的$Z(i)$为高斯整环，高斯整环中的每个元素为高斯整数，即复平面上正方形格点。

高斯域恰好是四次分圆域，因此常用来解决**四次互反律**问题。

高斯整数中，一个数有四个相伴数（含本身）。

高斯整数中的全体素数分为三类：

**分歧**数：1+i，为原来的2的因子。分歧数的共轭是它的相伴数，因此可以指定任一分歧数为本原数代表。

**惯性**数：所有正整数中4k+3形式的素数，在高斯整数中仍旧为素数。在整环扩张中保持了素数的特性，因此称为“惯性”。

**分裂**数：所有正整数中4k+1形式的素数，在高斯整数中可以拆成一对共轭的两个素数，这两个素数不相伴。这样的新素数是分裂的。

当然，这两个共轭的素数是不同的，即共轭的两个分裂数是互素的。

对于素数中的分裂数和惯性数，本原数的指定往往有着严格的规定，这是为了解决四次剩余问题的方便。

规定：高斯整数中的本原素数$\pi$有：

$$\pi\equiv 1 \mod 2(1+i)$$

在2(1+i)的缩系中有4个剩余类，除了1+i的每个素数的每个相伴数恰好落入其中一类。

对于1+i与它的相伴数，一般指定1+i是本原素数。

高斯整数最简单的应用，是解决勾股方程的解。勾股方程，就是满足下面形式的方程：

$$x^2+y^2=z^2$$

左边恰好构成高斯整数的范数，即：

$$N(x+yi)=z^2$$

通过模4的分析，我们知道右边模4必然余1，即如果含模4余3的惯性数因子，必然含偶数个。

由于分歧数和分裂数的范数都是一般整数中的素数，将左边唯一分解后必然也只能成对出现（在共轭与相伴的意义下）。即：

$${(u+vi)}^2=x+yi$$

$$N(u+vi)=z$$

用一般的整数写出来就是：

$$u^2-v^2=x$$

$$2uv=y$$

$$u^2+v^2=z$$

勾股方程的几何意义，就是单位圆上的圆周角定理，或者半正切的外能代换公式。

高斯整数同样能证明相应的四次形式无解。即：

$$x^4+y^4=z^4$$

利用高斯整数的唯一分解，可以解决圆上整点问题。即给定范数为n的条件下，有多少个高斯整数满足这个范数n：

$$N(x+yi)=n$$

仍旧将左边和右边唯一分解。左边在高斯整数意义下唯一分解，右边在正整数范畴唯一分解。

$$N({(1+i)}^a)N(u_1+v_1i)N(u_2+v_2i)=2^apq$$

对于分歧和分裂的素数，范数是原整数中的素数，而4k+3形式惯性的素数，范数是原素数的平方。因此n中4k+3形式的素数必须成对出现，否则无解。

然后利用简单的计数法就知道，在n中4k+3形式的素数成对出现前提下，整点个数与含多少个2（或1+i）无关，只与4k+1形式的素数个数有关，每一个4k+1形式的素数提供2中选择方法，在计数中扩大2倍。最后由对称性，整点个数乘4即可。

====Eisenstein整数====

注：Eisenstein（艾森斯坦）是Gauss的得意门生。

一般将$Q(\sqrt{3}i)$称为艾森斯坦域，相应的$Z(\sqrt{3}i)$为艾森斯坦整环，艾森斯坦整环中的每个元素为艾森斯坦整数，即复平面上正六边形格点。

艾森斯坦域恰好是三次分圆域，也是六次分圆域，因此常用来解决**三次互反律**问题。结合已经解决的二次互反律，就能给出六次剩余的手动计算。同样，如果结合高斯域中的四次互反律，就能解决十二次剩余的手动计算。

艾森斯坦整数中，一个数有六个相伴数（含本身）。

同样，艾森斯坦整数中的全体素数分为三类：

**分歧**数：$\frac{3+\sqrt{3}i}{2}$，为原来的3的因子。

**惯性**数：所有正整数中3k+2形式（2和6k+5形式）的素数，在艾森斯坦整数中仍旧为素数。

**分裂**数：所有正整数中3k+1形式（6k+1形式）的素数，在高斯整数中可以拆成一对共轭的两个素数，这两个素数不相伴。这样的新素数是分裂的。同样，这两个共轭的素数是不同的，即共轭的两个分裂数是互素的。

对于素数中的分裂数和惯性数，本原数的指定也有着严格的规定，这是为了解决三次剩余问题的方便。

规定：艾森斯坦整数中的本原素数$\pi$有：

$$\pi\equiv 1 \mod 3$$

在3的缩系中有6个剩余类，除了$\frac{3+\sqrt{3}i}{2}$的每个素数的每个相伴数恰好落入其中一类。注意，这与通常的3的剩余类不同。艾森斯坦整数中3的全部剩余类有9个，而缩系中有6个。

对于$\frac{3+\sqrt{3}i}{2}$与它的相伴数，可以指定$\frac{3+\sqrt{3}i}{2}$是本原素数。

艾森斯坦整数可以解决下面形式的方程的解：

$$x^2+3y^2=z^2$$

或者：

$$x^2-xy+y^2=z^2$$

或者：

$$x^2+xy+y^2=z^2$$

后两个在整数范畴是等价的。这里的求解完全仿照勾股方程即可，不再赘述。

艾森斯坦整数同样能证明相应的三次形式无解。即：

$$x^3+y^3=z^3$$

这是因为立方和公式拆分后，左边会多出艾森斯坦整数范数的一项。

利用艾森斯坦整数的唯一分解，可以解决一种椭圆上整点问题。即给定范数为n的条件下，有多少个艾森斯坦整数满足这个范数n。三种形式为：

$$x^2+3y^2=n$$

或者：

$$x^2-xy+y^2=n$$

或者：

$$x^2+xy+y^2=n$$

方法仍旧完全一样，不再赘述。

=====实二次域=====

====简介====

对于每一个实二次域，单位数有无限个。

对于实二次域，基（不含1和-1）只有一个。因此有如下结论：取ε为正的单位数中除了1以外两分量（任一）最小的，那么全体单位数可以表示成：

$$±\varepsilon ^n \quad n\in Z$$

研究实二次域的结构，首先就要求解基本单位数，因此就有了Pell方程（佩尔方程）。佩尔方程的解法已经被研究透彻，可以采用连分数算法。

====连分数====

连分数是一种记号。例如，长为4的连分数：

$$<a_0,a_1,a_2,a_3>=a_0+\frac{1}{a_1+\frac{1}{a_2+\frac{1}{a_3}}}$$

只是为了形式上简洁，才记成等号左边的样子。这里的四个变元可以任意取值。

连分数中，从前面一直到第k个变元构成的连分数，是它的**第k渐进分数**。连分数各变元的下标从0开始。

渐进分数有递推关系：

$$\frac{p_n}{q_n}=\frac{x_np_{n-1}+p_{n-2}}{x_nq_{n-1}+q_{n-2}}$$

这个式子和下面的Farey数列的递推关系很像。

上面的例子是有限连分数。如果分式无限地写下去，有无限个变元，就得到无限连分数。无限连分数收敛等价于渐进分数收敛。

有定理：

**无限连分数，如果各变元均大于等于1，那么一定收敛。**

因为只要各变元为正，无限连分数的偶渐进分数单调递增（都比它小），奇渐进分数单调递减（都比它大）。而在均大于等于1时，相邻（奇偶间）两个渐进分数之间距离可以给出估计式，趋于0，因此收敛。

显然可以看到，连分数关于下标为偶数的变元单调递增，关于下标为奇数的变元单调递减。这无论它有限或无限都成立。

还有一个事实，对于有限连分数，全体结尾为1的有限连分数和全体结尾不为1的有限连分数一一对应，即同一个连分数有两种表示：

$$<a_0,a_1,a_2,a_3>=<a_0,a_1,a_2,a_3-1,1>$$

**简单连分数：连分数从第1项开始全都是正整数。如果有限，要求最后一项不为1。**（第0项可以任意）

简单连分数的值，一定大于偶数的渐进分数，一定小于奇数的渐进分数。无限简单连分数一定收敛。

仿照一般分数的概念，第0项是0的连分数称为“真分数”。显然如果这之后的所有变元都大于等于1，那么得到的真分数一定落在0到1之间。

如果“规定第0项是该数的取整”，那么全体实数都有“唯一的简单连分数表示”。其中：

如果两个无限简单连分数的值相等，必然逐项相等。

如果两个有限简单连分数的值相等，不仅要逐项相等，而且必然项数也相同。

无限简单连分数不能与有限简单连分数值相等。有理数与有限简单连分数具有一一对应关系，因此无限简单连分数全都是无理数。

如果要求0、1区间内某个数的简单连分数表示（第0项为0），只需：

第一步：取倒数，得到的数大于1。

第二步：取整，得到的数在0到1之间。这个数称为余数。

第三步：对余数重复上述操作。

这样就得到了相应的表示。

====循环连分数====

仿照循环小数的概念，如果在连分数后面形成了循环，则形成**循环连分数**。如果循环节从第0项开始，称为**纯循环连分数**，否则称为**混循环连分数**。例如纯循环连分数：

$$<a_0,a_1,a_2,a_3,a_0,a_1,a_2,a_3,…>=<\overline{a_0,a_1,a_2,a_3}>$$

混循环连分数:

$$<a_0,a_1,a_2,a_3,a_1,a_2,a_3,…>=<a_0,\overline{a_1,a_2,a_3}>$$

混循环连分数后面循环的部分，一般要找最早循环的部分，称为它的“纯循环部分”。循环节一般取最小正整数。

对二次有理数执行“无限简单连分数”计算，即取倒数、取整交替，得到的余数还是二次有理数。

这里考虑一类特殊形式的二次有理数。

设$\xi$是如下形式的二次有理数，d只需是正数：

$$\xi_k=\frac{1}{q_k}\left(c_k+\sqrt{d}\right)=<a_k,a_{k+1},a_{k+2},a_{k+3},…>\quad q_k|N(c_k-\sqrt{d})$$

称$\xi_k$为**第k余数**。

只要让分母整除分子的范数，取倒数、取整的交替过程会变得简单很多。事实上如果这条件不成立，只要分子分母同时乘以分母的绝对值，并强行压入根号，条件就成立了。因此，任何二次有理数都能写成这形式。

如果这条件成立，那么每个余数都满足这条件。

根据上文的简单连分数算法：对余数取整可以得到a，进而得到新的c。

$$a_k=\frac{c_{k+1}+c_k}{q_k}$$
 
取整后得到的新的c为负数，且绝对值一定比根号d小，因此范数为负。取倒数，得到新的分母q，q总是正的。

$$q_k q_{k+1}=-N(c_{k+1}-\sqrt{d})$$

由于范数为负，取倒数之后根号d前面的符号不变，而c的符号由负变正（负数前面加负号变为正数）。

余数循环显然是循环连分数的循环条件。由于计算方法固定，只要余数重复必循环。

余数$\xi$里面，每个c都为负数，且绝对值一定比根号d小，因此c的个数有限。每个q都整除对应c构成二次整数的范数，因此q的个数有限。

余数有限必重复。因此**所有二次有理数都可以写成循环连分数**。反过来容易证明**循环连分数一定是二次有理数**。

由于无限简单连分数写法唯一，因此就找到了二次有理数的一种唯一表示。

====Stern-Brocot树====

Stern-Borcot树从两个简单的分数开始：0和无穷。

$$\frac{0}{1}\quad\frac{1}{0}$$

每次我们在相邻的两个分数$\frac{a}{b}$、$\frac{c}{d}$中间插入一个分数$\frac{a+c}{b+d}$，这样就完成了一次迭代，得到下一个序列。

$$\frac{0}{1}\quad\frac{1}{1}\quad\frac{1}{0}$$

$$\frac{0}{1}\quad\frac{1}{2}\quad\frac{1}{1}\quad\frac{2}{1}\quad\frac{1}{0}$$

将每一层新加入的数写在上一层下面，最终形成树状结构：

{{ :2020-2021:teams:namespace:正有理数树.png?400 |}}

每一个序列都是 Stern-Brocot 树的中序遍历。在每一层的序列中，真分数是单调递增的，并且序列中的分数都是最简分数。因此，Stern-Borcot树可以当成一棵平衡树，建立和查询就向平衡树一样做就行了。

====Farey数列====

第k层Farey数列，是分母小于等于k的所有最简真分数按大小顺序排列形成的序列，即满足单调性。可以理解为上面的树的左半分支的变形。

$$\frac{0}{1}\quad\frac{1}{2}\quad\frac{1}{1}$$

$$\frac{0}{1}\quad\frac{1}{3}\quad\frac{1}{2}\quad\frac{2}{3}\quad\frac{1}{1}$$

$$\frac{0}{1}\quad\frac{1}{4}\quad\frac{1}{3}\quad\frac{1}{2}\quad\frac{2}{3}\quad\frac{3}{4}\quad\frac{1}{1}$$

$$\frac{0}{1}\quad\frac{1}{5}\quad\frac{1}{4}\quad\frac{1}{3}\quad\frac{2}{5}\quad\frac{1}{2}\quad\frac{3}{5}\quad\frac{2}{3}\quad\frac{3}{4}\quad\frac{4}{5}\quad\frac{1}{1}$$

区别在于，仅当满足分母条件，需要插入特定的分数的时候，才在这个位置插入这个特定的分数。由于第n层插入的均为分母为n的最简分数，第n层插入的分数总共有$\varphi(n)$个。

====最佳有理逼近====

讨论如何用有理数“最佳地”逼近无理数，不妨假设无理数落入0、1区间。

“最佳”一词的概念：选定的有理数必须保证，比它与无理数的距离更近的有理数，分母都比它大。不存在分母比它小的有理数，到给定无理数的距离更近。

**最佳有理数：在法雷数列的某一行中，与给定无理数距离最近的那个有理数。**

这个有理数可能在下面几行依旧与无理数“距离最近”，但一定有某一行，会找到一个新的有理数，与无理数距离更近。因此去重复后可以得到**最佳逼近有理数列**，分母严格递增，距离递减。

比无理数大的称为上逼近，否则为下逼近。由于无理数和有理数之间一定有有理数，最佳逼近有理数列必然为若干个上逼近，之后若干个下逼近，交替进行的形式。

有结论：

**渐进分数必然为最佳逼近。偶项渐进分数全都是下逼近，奇项渐进分数全都是上逼近。渐进分数列是下上交错的逼近。**

在最佳逼近有理数列中，渐进分数是下上关系改变之前的倒数第一个数。如果将最佳逼近有理数列都写成有限简单连分数，那么在渐进分数之后（下上关系改变之后），连分数长度加一。

例如，$<0,2,4,2>$是根号6减去2的一个渐进分数。那么它的渐进分数列是：

$$<0,2>=\frac{1}{2},<0,2,4>=\frac{4}{9},<0,2,4,2>=\frac{9}{20}…$$
它的最佳逼近有理数列是：
$$<0,1>=1,<0,2>=\frac{1}{2},<0,2,1>=\frac{1}{3},<0,2,2>=\frac{2}{5},<0,2,3>=\frac{3}{7},<0,2,4>=\frac{4}{9},<0,2,4,1>=\frac{5}{11},<0,2,4,2>=\frac{9}{20}…$$

从每个渐进分数（不包含）开始，到下一个渐进分数（包含）为止，同为上逼近，或同为下逼近。

在最佳逼近列中，每一个最佳分数是上一个最佳分数与再往前一个渐进分数的分子分母对应求和。

====Dirichlet逼近定理====

Dirichlet（狄利克雷）逼近定理是说，对于任意的一个无理数$\theta$，均能找到无穷个有理数逼近它，满足不等式：

$$\left|\frac{p}{q}-\theta\right|\leqslant\frac{1}{q^2}$$

更强的结论是，在右边的分母上能放一个$\sqrt{5}$。这个$\sqrt{5}$是最优的，在无理数为黄金分割的时候取等。因此，在上式中的小于等于号事实上也可以是小于号。

更进一步的定理是Kronecker（克罗内克）的逼近定理。

如果$\theta$为无理数，$\alpha$为实数，则对任意正数$\varepsilon$, 存在整数n和整数m，使得:

$$\left|n\theta−m−\alpha\right|<\varepsilon$$

这两个定理都可以用抽屉原理来解决。事实上，历史上第一次正式提出抽屉原理，就是Dirichlet为了解决Pell方程而研究这个有理数逼近条件才正式提出来抽屉原理的。

当然，采用抽屉原理的证明可以发现，上文中提到的最佳逼近有理数列，每项满足定理中右边改成分母为一次式的不等式。

进一步有结论，渐进有理数列中，每一项均满足Dirichlet逼近定理的不等式。

我们用Dirichlet逼近定理来逼近二次根式$\sqrt{d}$。即有无穷个有理数（显然为正有理数）满足：

$$\left|\frac{x}{y}-\sqrt{d}\right|\leqslant\frac{1}{y^2}$$

于是，下面的范数就有：

$$\left|N(x+y\sqrt{d})\right|=|x-y\sqrt{d}||x+y\sqrt{d}|\leqslant\frac{1}{y}(\frac{1}{y}+2y\sqrt{d})\leqslant2\sqrt{d}+1$$

这是对范数拆出的两项进行估值。这也直观地说明只要有理数与根号d越接近，范数越小。

因此，范数较小的二次整数有无限个，进而采用一些手段，就可以推出范数为正负1的单位数存在，也存在无限个。

进而可以发现，对于所有$\sqrt{d}$的渐进分数，配上系数之后得到的二次整数的范数都落在非常小的区间。由于$\sqrt{d}$的渐进分数是余数循环的，只要其中出现使得范数为正负1的渐进分数，经过一个循环之后新的渐进分数凑成的二次整数也应当满足范数为正负1，即这个新渐进分数也是单位数。由于第-1个渐进分数规定为$\frac{1}{0}$，对应的二次整数范数为1，那么只要计算每个循环节处前一个渐进分数即可。

根据上逼近与下逼近的结论，第奇数个渐进分数得到的范数为负，偶数个为正。即是否存在范数为负1的二次整数取决于循环连分数的循环节长度是否为奇数。

最后还有一个结论，每经过一个循环，相当于旧的二次整数乘上了一个单位数，得到新的二次整数。因此上面得到的单位数是基本单位数。

====Pell方程的定义====

实二次域中，对范数的研究最终可以转化为如下形式的二元二次不定方程，称为佩尔方程：

$$x^2-dy^2=k$$

其中d和k是给定的整数，并且d是正数，要求x和y的整数解。显然x和y与符号无关，将x与y均为正的解，称为**正解**。

和二次域略有区别的是，这里d并没有要求“不含平方因子”的限制，但是可以将d的平方因子并入y中，变为某种特殊类型的解。

佩尔方程也可以理解成**双曲线上整点问题**，考虑双曲线上有没有整点，有多少个整点。

佩尔方程等价于，在实二次域中，求范数为k的两分量均为整数的特殊形式二次整数。

$$N(x+y\sqrt{d})=k$$

只要这样的二次整数存在，佩尔方程就有解，并且不同类型的二次整数对应不同类型的解。

对于上面的单位数求解问题，可以化为以下的**四个佩尔方程**：

$$x^2-dy^2=1$$

对应于范数为1的一种情况。

$$x^2-dy^2=-1$$

对应于范数为-1的一种情况。

$$x^2-dy^2=4$$

这里d除4余1，x和y的奇偶性必然相同。当x和y均为偶数时，两边同时除以4，方程与第一个方程等价。当x与y均为奇数时，对应于范数为1，并且x与y均为半整数的情况。

$$x^2-dy^2=-4$$

这里d除4余1，x和y的奇偶性必然相同。当x和y均为偶数时，两边同时除以4，方程与第二个方程等价。当x与y均为奇数时，对应于范数为-1，并且x与y均为半整数的情况。

因此在后两个方程中，更注重讨论x与y的奇数解。

====非4k+1形式二次域的Pell方程====

类似实二次域的基本单位数结论，可以给出这四类佩尔方程的“基本解”的概念。

有结论：只要d不是完全平方数，方程

$$x^2-dy^2=1$$

必然有解。并且记$\rho_1=x+y\sqrt{d}$为最小的正解，于是所有解都可以用

$$\{\pm\rho^n  | n\in Z\}$$

表示。

在d不含平方因子，并且ε的范数是1，并且ε不是半整数形式的时候，这里的ρ(1)就是对应的实二次域中基本单位数ε。此时对应的-1和-4的方程均无解，4的方程有解，但是无奇数解。

例如d为3时，将根号3表示为连分数：

$$\sqrt{3}=<1,\overline{1,2}>$$

根据连分数算法，可以得到$Q(\sqrt{3})$的基本单位数是$2+\sqrt{3}$。

方程

$$x^2-dy^2=-1$$

很多时候无解。有解的时候，记$\rho_{-1}=x+y\sqrt{d}$为最小的正解，于是所有解都可以用

$$\{\pm\rho^n  | n\equiv 1 \quad \mod 2\}$$

表示。

在d不含平方因子，并且ε的范数是-1，并且ε不是半整数形式的时候，这里的ρ(-1)就是对应的实二次域中基本单位数ε。并且ε的平方范数为1，ρ(-1)的平方是对应的ρ(1)。此时对应的4和-4的方程均有解，但是均无奇数解。

例如d为2时，将根号2表示为连分数：

$$\sqrt{2}=<1,\overline{2}>$$

根据连分数算法，可以得到$Q(\sqrt{2})$的基本单位数是$1+\sqrt{2}$。

在1到100中，二次域2、10、26、58、74和82存在这种形式的基本单位数，后文4k+1类中也存在退化成为此类的基本单位数。


====4k+1形式二次域的Pell方程====

这一类二次域中的基本单位数，需要将相应的二次根号除以2，再将连分数的第0项改为半整数，最后得到的才是相应的解。这是一种新的连分数算法。

方程

$$x^2-dy^2=4$$

总是有解，但是很多时候无奇数解。有奇数解的时候，记$\rho_4=\frac{x+y\sqrt{d}}{2}$为最小的正解，于是所有奇数解都可以用

$$\{\pm\rho^n  | n\equiv 1 ,2 \quad \mod 3\}$$

的分子表示。

在d不含平方因子，并且ε的范数是1，并且ε是半整数形式的时候，这里的ρ(4)就是对应的实二次域中基本单位数ε。此时对应的-1与-4的方程均无解。可以发现，ε的立方恰好分母为1，范数为1，因此ρ(4)的立方恰好是对应的ρ(1)。

例如d为21的时候。

又例如d为69的时候。

事实上1到100中，只有二次域69、77和93的基本单位数属于这种分母为2的半整数形式，而21、33和57属于非半整数形式。

方程

$$x^2-dy^2=-4$$

有解等价于对应的-1的方程有解，在此基础上，有奇数解等价于对应的4的方程有解。此时，记$\rho_{-4}=\frac{x+y\sqrt{d}}{2}$为最小的正解，于是所有奇数解都可以用

$$\{\pm\rho^n  | n\equiv 1 ,5 \quad \mod 6\}$$

的分子表示。

在d不含平方因子，并且ε的范数是-1，并且ε是半整数形式的时候，这里的ρ(-4)就是对应的实二次域中基本单位数ε。可以发现，ε的立方恰好分母为1，范数为-1，因此ρ(-4)的立方恰好是对应的ρ(-1)。同样有，ρ(-4)的平方是对应的ρ(4)，ρ(-4)的六次方是对应的ρ(1)。

例如当d为5的时候，$\frac{\sqrt{5}}{2}$的半整数连分数表示为：

$$\frac{\sqrt{5}}{2}=<\frac{1}{2},\overline{1}>$$

于是解得$Q(\sqrt{5})$基本单位数为$\frac{1}{2}+\frac{\sqrt{5}}{2}$。

但是d为17的时候，$\frac{\sqrt{17}}{2}$的半整数连分数表示为：

$$\frac{\sqrt{17}}{2}=<\frac{3}{2},\overline{1,1,3}>$$

于是解得$Q(\sqrt{17})$基本单位数为$4+\sqrt{17}$。我们失望地发现，它不属于这种形式。

在1到100中，5、13、29、53、61和85属于这种分母中含2的半整数形式，而17、37、41、65、73、89和97属于非半整数形式。

综上，以上所有情况可以简记为-4基本解的若干次乘方，构成6循环的序列：

-4、4、-1、4、-4、1、-4、4、-1、4、-4、1……

只有1的情况是一定有解的，仅当-1有解并且4有奇数解，二者同时成立的时候，-4才有奇数解。

当d含有平方因子的时候，1的方程仍旧一定有解。这时如果将根号d中的平方因子全部提出来，依旧能对应于一个二次域。由范数为1，可以得到结论：

这时1方程的基本解是对应二次域基本单位数的若干次方。

同样地，-1和4和-4的方程如果有解，基本解也是对应二次域基本单位数的若干次方。因此解的情况和对应二次域的情况有关。

====一般的Pell方程====

对于一定有解的1的方程，任意一个基本解的范数都是1。因此有结论：一般的佩尔方程

$$x^2-dy^2=k$$

一旦有解就一定有无数个解，即这样的双曲线上一旦有整点就一定有无数个整点。这是因为，将范数为k的对应的二次整数乘以范数为1的基本解对应二次整数，得到依然二次整数范数依然为k，而范数为1的基本解有无限个，相乘得到的范数为k的基本解不同。

依靠同一个解，与不同的单位数相乘，得到的全体解，称为**一族解**。同一个方程解的族数，可以大于1，也可以为1，也可以无解为0。

上面的讨论表明：方程-1的族数，为0或1，4的奇数解情形族数，为2或0，-4的奇数解情形族数，也为2或0。