本系统是一个基于 OpenGL 的 3D 粒子可视化教学平台，专门用于解释和展示 BSD 猜想（Birch–Swinnerton-Dyer 猜想）的核心内容。

BSD 教学系统（3D 可视化）介绍

一、本系统主要做什么

本系统是一个基于 OpenGL 的 3D 粒子可视化教学平台，专门用于解释和展示 BSD 猜想（Birch–Swinnerton-Dyer 猜想） 的核心内容。它将抽象的椭圆曲线、有理点群、局部质数信息、L 函数中心点行为以及 BSD 主公式，通过动态粒子云、柱状图、曲线和标注箭头等直观方式呈现出来。系统包含自动讲解流程（从“什么是 BSD”到“当前研究进展”），并提供多种交互模式（总览、曲线与有理点、局部质数 aₚ、Sato–Tate 分布、中心点 s=1、BSD 主公式、以及“陈恩华证明”演示模式），方便教学和研究表达。

二、BSD 猜想是什么

BSD 猜想是数论中最重要的未解决问题之一（克雷数学研究所千禧年大奖难题之一）。它建立了一条连接算术与分析的深刻桥梁。

核心断言：
设 E 为定义在有理数域 ℚ 上的椭圆曲线，
令 r 为 E(ℚ) 的 代数秩（即 Mordell–Weil 群中自由生成元的个数），
令 ord_s=1 L(E, s) 为 L 函数在 s=1 处的 解析秩（即零点阶数），
则 r = ord_s=1 L(E, s)。

更强版本（强 BSD）还要求 L 函数在 s=1 处的首项系数满足一个精确等式，包含诸如正则化子、周期、Tamagawa 数、Tate–Shafarevich 群大小等算术不变量。

三、陈恩华的研究是什么

在系统中，“陈恩华证明演示” 模式是开发者（陈恩华）个人对 BSD 猜想结构的一个研究表达和教学演示。该模式将 BSD 的代数层（有理点秩）、局部质数层（aₚ 信息）、中心点分析层（L(E,s) 在 s=1 附近的行为）以及主公式同步展示为五幕动态粒子过程：

第一幕：代数秩推进链（左）—— 展示秩的生成结构。
第二幕：局部质数同步柱（中）—— 展示每个质数 p 对应的 aₚ 信息。
第三幕：s=1 塌缩墙（中右）—— 强调中心点附近的变化。
第四幕：中心点陷落云（右）—— 模拟 L 函数在 s=1 的零点行为。
第五幕：顶部同步冠层（上）—— 象征整体一致性锁定。

注意： 该模式是研究视觉表达和教学工具，不代表 BSD 猜想已被严格证明。它的目的是帮助理解 BSD 猜想的各个组成部分及其内在联系。

四、BSD 相关公式排序

以下按从具体到一般、从局部到整体的顺序列出主要公式：

椭圆曲线方程（Weierstrass 形式）：

E: y² = x³ + a·x + b （a, b ∈ ℚ，判别式 ≠ 0）

局部质数因子（欧拉因子）：
对每个质数 p，定义局部因子：
```
Lₚ(E, s) = (1 - aₚ·p⁻ˢ + p·p⁻²ˢ)⁻¹ （若 E 在 p 处有好约化）
```
其中 aₚ = p + 1 - #E(𝔽ₚ)，是约化曲线在有限域上的点数偏差。
完全 L 函数（Euler 乘积）：
```
L(E, s) = ∏ₚ Lₚ(E, s)⁻¹
```
它通过解析延拓成为整个复平面上的亚纯函数，满足函数方程。
Sato–Tate 分布（非 CM 情形）：
归一化 aₚ / (2√p) 对应的角度 θₚ ∈ [0, π] 满足：
```
lim_{X→∞} (1/π(X)) #{p ≤ X : α ≤ θₚ ≤ β} = (2/π) ∫ₐᵦ sin²θ dθ
```
BSD 猜想（弱形式）：
```
rank(E(ℚ)) = ord_{s=1} L(E, s)
```
BSD 猜想（强形式 / 主公式）：
设 r 为代数秩，则 L(E, s) 在 s=1 处的泰勒展开首项系数满足：
```
        L^{(r)}(E, 1) / r! =
        (Ω_E · Reg(E) · ∏ₚ cₚ · |Sha(E)|) / |E(ℚ)_{tors}|²
        
```
其中：
- Ω_E 是实数周期（period）。
- Reg(E) 是正则化子（regulator）。
- cₚ 是 Tamagawa 数（局部因子）。
- Sha(E) 是 Tate–Shafarevich 群（其有限性尚未普遍证明）。
- E(ℚ)_{tors} 是有理挠点群。

总结： 本系统通过交互式 3D 粒子可视化，将上述抽象公式和概念转化为直观的图形，帮助学习者和研究者更深入地理解 BSD 猜想的结构与挑战。开发者陈恩华在此框架内加入了自己对证明路径的视觉化表达，旨在促进数学教育与研究交流。

代码：main.cpp