解方程 · 精确的根，与一堵数学的墙 · 符号计算的艺术

solve 给出方程的精确根：二次是求根公式，三次四次也有（繁琐的）根式解。但到了五次，一堵墙立在那里——Abel-Ruffini 定理证明：五次以上方程没有通用的根式解。这一章看符号求解的威力，以及计算机如何「知道」并诚实面对这堵墙。

解方程 · 精确的根，与一堵数学的墙

进入精确求解。解方程是符号计算最「解气」的能力之一——不是逼近一个根，而是给出精确的根。但这条路走到五次方程时，会撞上数学史上一堵著名的墙。这一章，既看它的威力，也看它如何诚实地面对自己的极限。

solve：精确的根

solve 求方程的精确解：

from sympy import *
x = symbols('x')

print(solve(x**2 - 2, x))     # [-sqrt(2), sqrt(2)]   —— 精确,不是 ±1.414
print(solve(x**2 + 1, x))     # [-I, I]               —— 复数根也精确

$x^2-2=0$ 的根是精确的 $\pm\sqrt2$ ，而不是数值解给的 $\pm1.41421356$ ； $x^2+1=0$ 干脆给出虚数根 $\pm i$ 。更漂亮的是，它能解带符号系数的方程，直接吐出你中学背过的求根公式：

a, b, c = symbols('a b c')
print(solve(a*x**2 + b*x + c, x))
# [(-b - sqrt(b**2 - 4*a*c))/(2*a), (-b + sqrt(b**2 - 4*a*c))/(2*a)]

$\dfrac{-b\pm\sqrt{b^2-4ac}}{2a}$ ——一字不差。这是符号计算的典型场景：不是解某一个具体方程，而是推导出一整个公式。

三次、四次：繁琐但仍有根式解

二次有求根公式，三次、四次也有（卡尔达诺公式、费拉里公式），只是繁琐得多。SymPy 照样能给：

sols = solve(x**3 - 2, x)      # x³=2 的三个立方根(1 实 2 复),都用根式精确表示
print(sols[0])                  # 2**(1/3)

$x^3=2$ 的三个根——一个实的 $\sqrt[3]2$ 和两个复的——都能用根式（ $\sqrt[3]{\ }$ 、 $i$ 、 $\sqrt3$ ）精确写出。四次也一样，虽然公式长到没人愿意手写，但存在，CAS 能给。

五次的墙：Abel-Ruffini 定理

到五次方程，天变了：

print(solve(x**5 - x + 1, x))
# [CRootOf(x**5 - x + 1, 0), CRootOf(x**5 - x + 1, 1), ...]  —— 不是根式!

SymPy 没有给出根式解，而是返回一串 CRootOf。这不是它偷懒——是一条深刻的数学定理：

Abel-Ruffini 定理（1824）：五次及以上的一般多项式方程，没有通用的「根式解」。 也就是说，不存在一个像求根公式那样、只用加减乘除和开方就能表达任意五次方程根的公式——这是被证明的不可能，不是「还没找到」。这堵墙背后，是 19 世纪 Évariste Galois 创立的群论 / Galois 理论：一个方程能否用根式解，取决于它的「对称性群」的结构。五次方程的对称群「太复杂」（ $S_5$ 不可解），于是根式无能为力。

所以计算机知道这堵墙的存在。面对五次，它不硬凑，而是诚实地用 CRootOf(多项式, k) 表示「这个多项式的第 k 个根」——这仍是一个精确的符号对象，能数值地算到任意精度，只是无法用根式写出。这和积分那章的 erf 异曲同工：当已有的符号不够用，就诚实地造一个新符号精确承载它。

r = CRootOf(x**5 - x + 1, 0)
print(N(r, 20))    # -1.1673039782614186843   —— RootOf 能精确求值到任意位数

超越方程：无穷多解怎么办

再看含三角、指数的超越方程，它可能有无穷多解。这时 solve 会力不从心，要用更严谨的 solveset——它返回一个解集，能表示「无穷多」：

print(solveset(sin(x), x, domain=S.Reals))
# Union(ImageSet(Lambda(_n, 2*_n*pi), Integers), ImageSet(...2*_n*pi + pi...))

$\sin x=0$ 的解是所有 $n\pi$ （ $n$ 取任意整数）。solveset 用 ImageSet 精确表达这个无穷集合——这是 solve（只返回有限个解）做不到的。方程的解可能是无穷集，符号系统用「集合」这个精确语言来承载它。

[1] Niels Henrik Abel (1824) / Paolo Ruffini (1799) —— 五次方程根式不可解。
[2] Ian Stewart.《Galois Theory》—— 方程可解性与群论，极佳的普及。
[3] SymPy 文档.「solve, solveset, CRootOf.」https://docs.sympy.org/latest/modules/solvers/solvers.html

解方程 · 精确的根，与一堵数学的墙