费马大定律代码化和定理《计算机科学中的数学》外扩学习1

前言

• 《计算机科学中的数学信息与智能时代的必修课》
  • 第1章 什么是证明
    • 1.1 命题
      • 1.15 费马大定律-(又称费马最后定理)
      • xn+yn=znx^{n}+y^{n}=z^{n}xn+yn=zn没有正整数解.

基础算式展开

• 当n=1时,方程就是x+y=z,显然有无数正整数解(1+2=3)当n=1时,方程就是x+y=z,显然有无数正整数解(1+2=3)当n=1时,方程就是x+y=z,显然有无数正整数解(1+2=3)
• 当n=2时,方程就是x+y=z,显然有无数正整数解(32+42=52)当n=2时,方程就是x+y=z,显然有无数正整数解(3^{2}+4^{2}=5^{2})当n=2时,方程就是x+y=z,显然有无数正整数解(32+42=52)
• 费马断言,当n≥3时,这样的正整数解不存在费马断言,当n≥3时,这样的正整数解不存在费马断言,当n≥3时,这样的正整数解不存在
  • 我们用python代码简单试一下

n=3x=1y=2print(f"x={x},{n}+ y={y},{n}",pow(x,n)+pow(y,n))print(3**n)

从图中可以看到当n≥3的时候,13+23=331^{3}+2^{3}=3^{3}13+23=33是不对的.或者说是无解的

尝试设计复杂代码

我尝试利用while循环来帮助我实现对费马大定律的实现

n=1x=1y=2z=0whilepow(x,n)+pow(y,n)==pow(z,n):print(f"({x}^{n})+({y}^{n})=({z}^{n})")n+=1else:print(f"n={n},费马断言成立")

这个并不是我想要的结果.

用ai分析.问题如下

1. 核心逻辑错误while循环条件初始不成立.

• 左边的x=1 y=1结果是pow(1,1) + pow(2,1) = 1 + 2 = 3
• 右边的z=0, 结果是pow(0,1)
• 条件判断为3==0即False
• 而 while 循环的执行规则是：只有条件为True时，才会进入循环体执行内部代码，因此该循环体（print 和 n += 1）永远不会被执行，直接跳转到else分支输出结果。

我从这个点上想想办法
有点变化

n=1x=1y=2z=x+ywhilepow(x,n)+pow(y,n)==pow(z,n):print(f"({x}^{n})+({y}^{n})=({z}^{n})")n+=1else:print(f"n={n},费马断言成立")

2「概念错误」：费马大定理的条件与变量取值完全错误

• 费马大定理的核心条件是 n > 2，而代码从 n=1 开始验证，且混淆了 “有解” 与 “无解” 的逻辑：
• n=1 时，方程 x + y = z 是普通加法，存在无数正整数解（如 1+2=3），与费马大定理无关；
• n=2 时，方程是勾股定理，存在无数正整数解（如 3²+4²=5²），同样与费马大定理无关；
• 费马大定理关注的是 n ≥ 3 时，方程无正整数解。
• 变量 z 取 0 违反 “正整数解” 的要求：费马大定理讨论的是 x、y、z 均为正整数（大于 0 的整数）， 而代码中 z=0 不符合定理的前提条件。

n=1#核心条件是大于2,而为1则是普通加法,为2时勾股定理.我们需要的是大于2x=1y=2z=x+y#这段重新修改了,算是可以了

3「逻辑缺陷」：while-else 结构的使用场景错误，逻辑闭环缺失

• Python 中 while-else 的特性是：else 分支仅在 while 循环条件自然变为 False（未被 break 强制终止）时执行，此处用于验证费马大定理完全不适用。
• 代码没有对 x、y、z 进行任何更新逻辑：即使循环条件初始成立，循环体内也只更新了 n，x、y、z 始终保持初始值，会导致无限循环（若条件为 True），无法完成多组数值的验证。

n=1x=1y=2z=0#这个初始就是有问题.whilepow(x,n)+pow(y,n)==pow(z,n):print(f"({x}^{n})+({y}^{n})=({z}^{n})")n+=1#这里缺乏x,y的更新逻辑 +=1,会导致无限循环else:#不适用费马定律,print(f"n={n},费马断言成立")

4「语法细节缺陷」：pow() 函数在 n=0 时存在隐含风险

虽未触发，但需注意

• 若后续代码导致 n 变为 0，pow(0, 0) 在 Python 中会返回 1（数学上无明确定义），而 pow(x, 0) 对任意 x≠0 均返回 1，会导致验证结果失真，但此问题在当前代码中未触发。

利用ai写的代码,学习用,非本人写

# 验证费马大定理：n>2时，x^n + y^n = z^n 无正整数解defverify_fermat_simple(n_start,max_n,max_xyz):""" 简单验证费马大定理（有限范围） :param n_start: 起始指数（建议≥3） :param max_n: 最大指数 :param max_xyz: x、y、z的最大取值（正整数） """# 遍历指定范围的指数nforninrange(n_start,max_n+1):has_solution=False# 遍历所有可能的正整数x、y、zforxinrange(1,max_xyz+1):foryinrange(1,max_xyz+1):forzinrange(1,max_xyz+1):ifpow(x,n)+pow(y,n)==pow(z,n):print(f"找到解：{x}^{n}+{y}^{n}={z}^{n}")has_solution=Truebreakifhas_solution:breakifhas_solution:break# 输出当前n的验证结果ifnothas_solution:print(f"n={n}，在x、y、z≤{max_xyz}范围内，费马大定理断言成立")else:print(f"n={n}，找到有效解，费马大定理断言不成立（有限范围）")# 调用函数：验证n=3到n=5，x、y、z最大为100if__name__=="__main__":verify_fermat_simple(n_start=3,max_n=5,max_xyz=100)