Skip to content

算法Algorithm

区间根覆盖证书

Interval root cover certificate · Validated branch-and-prune root coverage

以互不相交的局部根证书和完整二叉盒树认证紧区域内的全部根,保留未决叶,证明覆盖不变量并明确有限停止与计算成本。

形式陈述 ​

输出是一张完整覆盖,而非一串近似根 ​

设目标区域 D⊂Rn 为紧盒,F 在其开邻域连续可微。先给出有限个(允许 m=0)互不相交、包含于 D 的正半径闭盒 B1,…,Bm,每个都带一张成功的严格Krawczyk证书,因而各含恰一个根且根位于盒内部。候选盒可以由普通近似求根产生,但通过认证以前不计入根数。

全域证书是一棵以 D 为根的有限二叉树。每个内部节点保存分割坐标,子盒沿该坐标中点精确分成两半,闭边界允许重叠。每个叶子只能有以下两种已证明的标签:

  1. 无根叶: 对某个分量 i,可靠区间值 Fi(X) 不含零。
  2. 已覆盖叶: X⊂intBj,并记录相应编号 j。

若全部叶子都有这两种标签,便证明 D 中恰有 m 个不同根。只要还有未决叶,结果就必须保留“全域根数未认证”的状态,即使已经找到了若干正确根。

构造算法与包含不变量 ​

维护待处理盒队列,初始只有 D。取出一个盒 X 后,先尝试已覆盖判据,再计算 F(X) 尝试无根判据;两者都失败时,沿最长边二分,将两个子盒放回队列并保存父子关系。若达到节点、位数或时间预算,将尚未处理部分保留为未决叶,不能把它们默认当成空区域。

每一步维护如下不变量:原区域 D 等于全部已经结束的叶盒与当前待处理盒的并;任何尚未归入候选根盒的真根都留在待处理盒内。局部证书证明存在性,覆盖不变量负责没有遗漏,互不相交负责没有重复计数。

这个算法把“猜候选”和“检查完整性”分开。没有承诺从任意黑箱函数自动发现一切根;它给出了能独立复核一份候选清单是否完整的计算任务。

直觉

找到两个交点,回答的是“至少有两个”。要回答“恰有两个”,还需检查所有没有画出的地方。覆盖树把整个矩形逐块交代:这块某条方程不可能等于零;那块完全落在已证明只含一个根的盒里。每一次分割都记住两个孩子,检查者不必相信搜索程序没有悄悄跳过一片区域。

同一根附近可能出现好几个已覆盖叶。它们共用一个根盒编号,根数仍只计一次。相反,两张相交的局部根证书可能指向同一个根,所以本页从输入就要求候选根盒互不相交。

完备证书为何推出精确根数 ​

每个局部证书给出一个根,盒子两两不交,故至少有 m 个不同根。对任意 x∈D,从树根沿包含它的子盒下行,有限树最终到达某个叶子;若 x 在分割边界,可以任选包含它的孩子。无根叶不可能含零点,已覆盖叶中的根只能是对应 Bj 中的唯一根。因此没有第 m+1 个根。

验证器可以从父盒和分割坐标重新生成子盒,而不接受子盒自行声称的端点。这让“无缝覆盖”变成局部可检查规则;所有端点为有理数时,中点分割也是精确运算。

例子与边界

圆与指数曲线的全域计数 ​

取

D=[−1,1]×[0,2],F(x,y)=(x2+y2−1y−ex/2).

使用Krawczyk页中的 B−,B+ 及两个有理预条件矩阵。公开证书生成了93个分割节点、90个无根叶和4个已覆盖叶,最大深度24;总计187个节点,无未决叶。两张局部证书加上这棵树,证明目标矩形内恰有两个根。

四个已覆盖叶并不表示四个根。它们只分配到两个互不相交的已认证盒。完整树和验证脚本可在终点页下载,检查者可以重新计算每个叶子的方程包围。

若只允许处理5个节点,同一算法留下6个未决叶。此时两张局部证书仍各自有效,但不能报告全域恰有两个根。扩大预算与已证正确性是两件不同的事。

小盒也可能仍然未决 ​

候选清单允许为空。对 f(x)=x2+1、D=[−1,1],函数区间为 [1,2],树根本身就是无根叶;这是一张 m=0 的完整证书。

对 f(x)=x2,没有根盒覆盖零时,每次细分都会留下含零的盒,函数值区间也始终包含零。盒宽低于 10−20 并不能改变这个逻辑状态。重根需要别的局部证明,不能把“很窄”当作存在且唯一。

复多项式的重根簇可以改交圆盘总重数证书:例如 (z−c)3 在任意正半径圆盘中可认证总重数三,无须假装导数可逆。这个输出不宣称三个不同根,也不替代本页一般实向量黑箱的存在且唯一合同;它利用的是全纯结构和多项式系数。

本页的根盒必须在 D 内,且其根在盒内部。因此 f(x)=x、D=[0,1] 的边界根不能由当前格式的严格内包含盒认证。可以扩展目标区域后另外处理原边界,或引入适合边界的精确证书;不能未经说明把域外的一半盒当作域内覆盖。

推论与应用

有限停止的清楚版本 ​

假设候选盒确实覆盖了 D 的全部根,且区间求值满足一致收敛性:当 X⊂D 的直径趋零时,wid(Fi(X)) 对全部这样的小盒一致趋零。端点精度或特殊函数尾界应同步提高,以免固定误差底破坏这个条件。最长边二分、每个有限深度节点最终都会处理的搜索,将在有限步结束。

证明可由有限维紧性给出。反设有任意小的未决盒 Xk,选 xk∈Xk,再取收敛子列 xk→x∈D。未决表示每个分量区间都含零;包含性和宽度趋零迫使 F(x)=0。但所有根都落在某个 intBj,所以充分小的 Xk 必完全落入该盒,应该被标为已覆盖,矛盾。于是存在统一尺度,所有更小盒都能结束。最长边二分达到该尺度所需的深度有统一上界,有限深度的二叉树只有有限节点。

这一定理的“候选完整”是假设;实际有限成功证书反过来证明它。若候选漏根,搜索可能一直保留未决区域;若求值平台阻止区间收窄,即使候选完整也可能无法自动停止。公开算例使用固定18阶指数包围并确实完成有限树,因此其正确性由实际叶证据建立,不依赖声称该固定阶扩张在无限细分下必然收敛到零宽。

证书大小与检查成本 ​

若树有 V 个节点、m 张根盒证书,逐叶扫描全部候选盒的直接实现连同分割与坐标验证需 O(V(1+m)n) 次坐标操作,加 O(V) 次函数区间求值。验证 m 个稠密Krawczyk证书通常另需 O(mn3) 算术操作。树的显式存储为 O(V) 个记录,深度优先的工作栈为 O(H),其中 H 为最大深度;有理端点位长和特殊函数求值次数都须另计。

这些成本以实际证书大小为参数,没有给任意 F 关于输入长度的多项式上界。根很近、导数接近奇异或区间表达式过宽时,V 可以很大。结构良好的证书允许独立验证器换用不同搜索顺序,而不必复制作者寻找候选根的过程。

参考资料
关系图谱7 个相邻概念 · 2 类关系

拖动节点调整位置。

显示关系

显示:依赖

  1. 前置三跳
  2. 前置二跳
  3. 前置一跳
  4. 当前条目
  5. 后续一跳
  6. 后续二跳
  7. 后续三跳
文字版关系按与当前条目的最短距离分组
类型化关系