fix(gcd): cont 单项情形 sign 不一致 (#17) — 对齐 Maple/SymPy/FLINT 约定

[lihaokun]

Summary

修复 issue #17：CLPoly cont() 对单项与多项多项式返回不一致的 sign：

• 单项：cont(-x) = -1（返 leading coef 带符号）
• 多项：cont(-x+1) = 1（gcd 折正）

而 Maple / SymPy / FLINT 统一约定 content always non-negative。这也是 issue #14 (squarefreebasis unit leak) 的根因之一（PR #16 仅在 sqf 层兜底修了表象）。

Changes

3 处修复（详细见 docs/fixes/2026-05-12-cont-single-term-sign.md，含半形式化正确性证明）：

1. clpoly/polynomial_gcd.hh:1067 (univariate ZZ cont 特化版)：init ZZ c = 0，gcd(0, c_i) 走 GMP mpz_gcd 取绝对值，自动归正
2. clpoly/upolynomial.hh:220 (template Tc cont 版本)：同样模式
3. clpoly/polynomial_gcd.hh:631+ (多变量 cont)：return 前加 4 行 sign normalize，cover single-group 路径

Before / After

                  Before              After (matches Maple/SymPy/FLINT)
cont(x)         = 1     ✓             1     ✓
cont(-x)        = -1    ✗             1     ✓
cont(-2x)       = -2    ✗             2     ✓
cont(-x+1)      = 1     ✓             1     ✓
cont(-2x+2)     = 2     ✓             2     ✓
cont(常数 -6)   = -6    ✗             6     ✓
cont(-2xy mv)   = -2y   ✗             2y    ✓
cont(0 empty)   = 1     ✓             1     ✓

Test plan

• 新增 test/test_cont_sign_consistency.cc (16 用例 / 16 PASS)
  • Univariate cont 7 case（含 issue cont() 单项 vs 多项情形 sign 不一致（违反 Maple/SymPy/FLINT 约定） #17 复现）
  • Multivariate cont (lex) 3 case
  • cont 非负 invariant 6 case
• 全量回归 27 test files / 0 failed（含 issue squarefreebasis produces constant -1 in output basis #14 test_squarefreebasis_bugfix）
• 参考实现 SymPy 实测对照通过

Semi-Formal Correctness Proof

详 docs/fixes/2026-05-12-cont-single-term-sign.md §3：

• Spec 数学定义
• Lemma 1（loop invariant，归纳证明）
• Main Theorem（实现 = Spec，分 case）
• 等价性 Claim（修后多项情形与原实现等价）
• 边缘情形枚举表

Lean side

Lean Model.lean 端有自己的 MvPolyZZ.contImpl / SparsePolyZZ.contImpl（在 feature/formal-proofs 分支），同样问题。本 PR 后会在 feature/formal-proofs 分支独立同步修（不在本 PR 范围）。

Closes #17

🤖 Generated with Claude Code