问题
生成元 a1,…,ak(gcd 为 1)的 Frobenius 数 F(a1,…,ak) 是不能表示为非负组合的最大整数。k=2 有拉马努金式封闭公式 F = a1a2−a1−a2(01 号命题 2 的证明);k=3 有 Selmer 等的 «算法式» 表达式(通过有限次取整与极值,长度可控);k ≥ 4 时问题变脸:Ramírez Alfonsín(1996)证明一般情形 NP-难。本问题问:NP-难定性之后,是否存在任何 «可用的显式公式»——哪怕超多项式长度但结构清晰(如逐段拟多项式、按参数分层给出)?固定 k 时公式谱的边界在哪里?Curtis(1990)证明:不存在 «有理函数拼接» 型的封闭公式(精确地:任何把 F 表示为生成元的逐段有理函数、段数有限的表示不存在,对 k ≥ 3),这为 «公式» 一词划定了底线,但 «底线之上还有什么» 持续争论。
主要概念
Apéry 集作为统一语言
Apéry 集 Ap(S,a1) = {w0,…,wa1−1},其中 wj 是 S 中模 a1 余 j 的最小元素。它一举决定三个量:F = max Ap − a1、g = (Σ Ap)/a1 − (a1−1)/2(Selmer 公式)、以及表示问题本身。因此 «Frobenius 公式问题» 等价于 «Apéry 集的显式刻画问题»:k=2 时 Apéry 集是等差数列 {0,b,2b,…}(01 号),k=3 时可证明其结构为 «三段分段线性»(故有算法式公式),k ≥ 4 时 Apéry 集可以任意复杂——这就是公式的边界。
NP-难与显式公式的相容性
NP-难是对 «输入长度为 log2a1+…+log2ak 的计算» 的复杂度定性(归约自整数规划,Ramírez Alfonsín 1996;Kannan 对固定 k 给出多项式算法——k 固定时问题在 «维数» 上属参数复杂度阶梯的另一端)。逻辑上 «NP-难» 与 «存在显式公式» 可以相容:公式可以巨长但求值快。真正的障碍是结构性的:Curtis 定理排除了逐段有理函数表示,而现有的一切 «漂亮公式» 都是这种类型。
主要命题与证明
命题 1(Selmer 公式与 Frobenius 公式;完整证明) 设 Ap(S,a1) = {w0,…,wa1−1}。则(i)F(S) = maxj wj − a1;(ii)g(S) = (Σj wj)/a1 − (a1−1)/2。
证明 先证 Apéry 集的基本性质:每个 wj 满足 wj − a1 ∉ S(否则与最小性矛盾),且 S ∩ (模 a1 余 j 类) = wj + a1ℕ。(i)设 wj 最大。wj − a1 不可表示(上面性质),而一切 > wj − a1 的数 x:若 x ≥ wj 则 x ≥ wj 且 x 在某个余类中不小于该类 Apéry 元(x ≥ wj ≥ wj′ 对所有 j′),故 x ∈ S;若 wj−a1 < x < wj,则 x 的余类 j′ 的 Apéry 元 ≤ wj ≤ x+…(x > wj−a1 ⟹ x ≥ wj−a1+1;若 x 的 Apéry 元 wj′ > x,则 wj′ ≥ x+1 > wj−a1+1 ⟹ wj′ − wj ≥ 2−a1+…(更直接:x ≡ j′,x − wj′ ≡ 0(mod a1)且 x < wj′ ⟹ x ≤ wj′−a1 < wj,即 x < wj 时 x ∈ S ⟺ x ≥ wj′…)——完整链条:任何 x > max Ap − a1 时 x 所在类的 Apéry 元 wj′ ≤ x(否则 x ≤ wj′ − a1 ≤ max Ap − a1,矛盾),故 x = wj′ + ta1 可表示。(ii)间隙恰是每个余类中 Apéry 元之下的非元素:第 j 类的间隙是 wj, wj−a1, … 之下不可表示者,即 {wj − a1t : t ≥ 1} ∖ 可表示者——由于 wj−a1t(1 ≤ t ≤ ⌊wj/a1⌋)全部不可表示(减 a1 后仍是同余类中更小者,若可表示则 wj 可表示,矛盾),第 j 类恰有 ⌊wj/a1⌋ 个间隙。求和:g = Σ ⌊wj/a1⌋ = (Σ(wj − rj))/a1,其中 rj = wj mod a1 = j,故 Σ rj = a1(a1−1)/2,得 g = (Σ Ap)/a1 − (a1−1)/2。□
命题 2(Curtis 定理,1990;陈述) 对 k ≥ 3,Frobenius 数不能表示为生成元的 «逐段有理函数»(任意有限段、每段是生成元的有理函数)。特别地,不存在拉马努金公式到 k ≥ 3 的 «同类» 推广。
梗概 证明在参数平面 (3, 2m+1) 或类似切片上追踪 F:在该切片上 F 随参数振荡的 «锯齿» 行为与任何有理函数沿切片的限制(在其定义域上是逐段光滑且单调性受控)矛盾——锯齿次数随参数无界增长,而有限段有理函数只能产生受控的振荡次数。□
命题 3(复杂度的两端,Ramírez Alfonsín 1996 与 Kannan;陈述) 变元个数 k 作为输入一部分时,计算 F 是 NP-难(归约自整数规划的可行性判定);固定 k 时存在多项式算法(Kannan 1992 的 k-次多项式,基于整数规划的固定维数算法)。因此 «公式» 的可得性随 k 逐级坍塌:k=2 封闭公式、k=3 算法式、k 固定多项式算法、一般 NP-难。
现状说明 NP-难归约的机制:把背包式判定(能否恰好凑出 N)嵌入 Frobenius 计算——F 与 «紧邻 F 的可表示性» 携带整数规划的信息。固定 k 的 Kannan 算法则在参数空间的 «格子» 上做极小化,理论上有界但常数巨大;«把 Kannan 算法变成可读公式» 的努力至今没有突破 Curtis 屏障。□
进展与文献
Sylvester 与拉马努金的二生成元公式;Selmer《On the linear Diophantine problem of Frobenius》(1977);Curtis《On formulas for the Frobenius number of a set of two or three integers》(1990);Ramírez Alfonsín《The Diophantine Frobenius Problem》(1996 论文与 2005 专著);Kannan 的固定维数算法(1992)。开放核心:固定 k=4 的 «最简» 公式形态;Apéry 集结构复杂性的精确度量(与 10 号 Betti 复杂度耦合);随机生成元情形的期望公式(08 号)作为 «平均情形公式» 的替代路线。