[论文解读] Method of Hidden Parameters and Pell's Equation

S. N. Arteha|ArXiv.org|Mar 14, 2003

Neural Networks and Applications参考文献 8被引用 3

一句话总结

本文提出了一种名为'隐参数法'的方法，用于推导佩尔方程 $x^2 - Ay^2 = 1$ 的最小解的显式参数公式，通过固定参数 $l$ 和 $m$ 将不同 $A$ 值的解关联起来，显著减少了计算时间。该方法能够利用小的辅助参数直接计算无限多组解，与标准连分数方法相比，计算时间最多可减少 1.83 倍。

ABSTRACT

Using the representation of numbers, the methods of decreasing the number of calculation steps for Pell's equation are developed. The formulae relating the natural solution x0 of Pell's equation for A0 with the infinite number of minimal natural solutions xi for other Ai are obtained for various cases. The parametric representation of Pell's equation solutions are obtained with the help of new "method of hidden parameters".

研究动机与目标

• 为任意非完全平方数 $A$ 减少求解佩尔方程 $x^2 - Ay^2 = 1$ 所需的计算步骤数量。
• 推导出将给定 $A_0$ 的最小解 $x_0$ 与其它 $A_i$ 的无限多组最小解 $x_i$ 相关联的显式参数公式。
• 识别并利用解中的结构模式，特别是针对解的大小达到最大值的素数和'准素数' $A$ 值。
• 开发一种系统化方法——'隐参数法'，通过引入和消除辅助参数，将解题过程转化为参数恒等式。

提出的方法

• 引入两类关系：'水平'关系中，基本参数 $l$ 和 $m$ 固定，$x_i$ 随索引 $i$ 线性变化；'垂直'关系中，$l$ 和 $m$ 非线性变化。
• 通过 $a_i = a_0 + i|l^2 - 3m^2|$，$b_i = b_0 + 2ilm$，$A_i = a_i^2 + 3b_i^2$ 的参数表示，生成无限多组解。
• 应用恒等式 $(y^2 - Ax^2)(Y^2 - AX^2) = (yY \pm AxX)^2 - A(yX \pm xY)^2$ 推导并验证参数恒等式。
• 使用递推序列 $B_n = d_n B_{n-1} + B_{n-2}$，其中 $B_{-2} = 1$，$B_{-1} = 0$，将连分数渐近分数与解参数关联。
• 利用条件 $dl - 3gm = \pm 1$，将初始值 $a_0$ 和 $b_0$ 参数化为 $d$，$g$，$l$，$m$ 的函数，以构建解族。
• 通过引入隐参数，推导其相互关系，并消去这些参数，最终以初始 $A_0$ 和 $x_0$ 表示解。

实验结果

研究问题

• RQ1 在不同 $A$ 值下，佩尔方程最小解的参数结构是什么，特别是针对素数和'准素数' $A$？
• RQ2 如何通过利用解的结构模式来减少求解佩尔方程的计算步骤？
• RQ3 能否从一个已知解出发，通过显式参数公式生成佩尔方程的无限多组解？
• RQ4 不定线性方程 $ax - by = 1$ 在将佩尔方程与隐参数法联系起来的过程中起什么作用？
• RQ5 解的大小 $x_i$ 的局部最大值和绝对最大值如何与 $A_i$ 的算术结构相关联？

主要发现

• 隐参数法能够将佩尔方程的解以显式参数恒等式形式表示，将单个 $A_0$ 的解 $x_0$ 与其它 $A_i$ 的无限多组解 $x_i$ 关联起来。
• 对于 $A < 1,000,000$，使用该方法的改进算法将计算时间从标准连分数法的 3 小时 2 分钟减少到 1 小时 43 分钟，提速 1.827 倍。
• 对于 $A < 2,000,000$，提速提升至 1.828 倍，表明该方法在最困难情况（如绝对最大值）下尤为有效。
• 所有 $A > 46$ 的绝对最大值均出现在 $A = 4n+1$ 的素数上，局部最大值则出现在素数 $p_i$ 或双倍素数 $2p_j$ 上，表明解的复杂性存在结构模式。
• 该方法识别出'准素数' $A$ 值（如 $A = 1729$），其解结构类似于素数，但在逐差法中不呈现标准形式。
• 当 $A = A_0^2 + 3B_0^2 = 4N+1$ 且 $B_0$ 为偶数时，$x = 2(l^2 + 3m^2)$ 的参数解成立，并通过代换 $a = A_0 \pm 3B_0$，$b = A_0 \mp B_0$ 生成解。