|
|
已按你的要求移除表格,将数值示例改为条目式叙述,保持A4单页紧凑排版。
---
有效三元组临界指数理论及数值计算(A4单页终版)
作者:朱火华
---
一、有效三元组定义
设正整数 a,b,c 满足:
a \le b < c,\qquad a+b > c
此类三元组称为有效三元组。此时方程 a^n+b^n=c^n 在 n=2 时左边大于右边,n\to\infty 时左边小于右边,由介值定理存在唯一实数解 n_{\text{crit}}>1,称为该三元组的临界指数。
若 a+b \le c,则临界指数 \le 1,属次费马区域,本文不讨论。
---
二、等腰三元组闭式公式
当 a=b 时,方程化为 2a^n=c^n,直接解得:
\boxed{n_{\text{crit}} = \frac{\ln 2}{\ln(c/a)}}
常用数值:\ln 2 \approx 0.6931472。
特例:模三元组 (K+1,\,K+1,\,K+2) 的统一临界指数:
n_{\text{mod}} = \frac{\ln 2}{\ln\!\bigl(\frac{K+2}{K+1}\bigr)}
---
三、一般三元组数值求解
当 a<b 时无闭式解,需数值求解方程:
a^n + b^n = c^n
\quad\text{或}\quad
\left(\frac{a}{c}\right)^n + \left(\frac{b}{c}\right)^n = 1
二分法:函数 f(n)=a^n+b^n-c^n 严格递减,上界可取 n_{\max}= \frac{\ln 2}{\ln(c/a)} 与 \frac{\ln 2}{\ln(c/b)} 中较大者乘以安全系数2,迭代至所需精度。
牛顿法(快速精化):
n_{k+1} = n_k - \frac{a^{n_k}+b^{n_k}-c^{n_k}}{a^{n_k}\ln a + b^{n_k}\ln b - c^{n_k}\ln c}
初始近似:
n_0 = \frac{\ln(c/a) + \ln(c/b)}{2\ln c - \ln a - \ln b}
推荐先用二分法收缩区间,再以牛顿法精化,2~3步可得高精度解。
快速近似公式(a,b 接近时):
n_{\text{crit}} \approx \frac{\ln 2}{\ln\!\bigl(\frac{2c}{a+b}\bigr)}
---
四、数值示例
· (6,6,7)(模三元组):闭式解得 n \approx 0.6931472 / 0.1541506 \approx 4.4966,该值大于 a=6 且非整数。
· (6,6,8)(等腰):闭式解得 n \approx 0.6931472 / 0.2876821 \approx 2.409,小于3。
· (5,6,7)(非等腰):二分+牛顿得 n \approx 2.973,小于3。
· (4,5,6)(非等腰):精确值 \approx 2.48705;近似公式给出 2.409,相对偏差约3%。
所有有效三元组的临界指数均为无理数。当 n_{\text{crit}}<3 时,方程在整数指数 n\ge 3 下无解;当 n_{\text{crit}}\ge 3 时,其值恒非整数。此数值现象与费马大定理完全吻合。
---
五、逻辑骨架(三元组分类与计算路径)
```
全部有效三元组 (a≤b<c, a+b>c)
│
├─ 等腰 (a=b) → 闭式 n = ln2/ln(c/a)
│ ├─ 模三元组: (K+1,K+1,K+2)
│ └─ 一般等腰
│
└─ 非等腰 (a<b) → 数值求解
├─ 二分法定界
└─ 牛顿法精化
```
---
六、附录:关键常数与公式速查
· \ln 2 \approx 0.6931471805599453
· 等腰闭式:n = \dfrac{\ln 2}{\ln c - \ln a}
· 模三元组:n = \dfrac{\ln 2}{\ln(K+2) - \ln(K+1)}
· 牛顿增量:\Delta n = -\dfrac{a^n+b^n-c^n}{a^n\ln a + b^n\ln b - c^n\ln c}
· 有效条件:a \le b < c 且 a+b > c
---
全文毕
排版适配 A4 单页,纯文本无表格,可直接打印或作为数值计算参考卡使用。 |
|