Skip to content

方法Method

输出敏感分析

Output-sensitive analysis · Output-sensitive algorithms

在明确输出接口后,同时用输入规模与实际输出规模刻画算法运行成本的分析方法。

形式陈述 ​

通常的时间复杂度把同样大小的输入合并取最坏情况。输出敏感分析进一步固定计算模型、输入编码与输出表示,令 n(x) 为输入规模、k(x) 为算法必须报告的对象数,以同一组常数对所有合法输入保证

TA(x)≤cF(n(x),k(x))+c0.

这里的 k 是问题规定的答案大小,不能把算法自己的任意中间工作量改名为输出。报告区间实例、报告匹配位置、输出凸包极点,分别定义了不同的 k;同样一个输入改为只返回计数,便改变了接口。若每个输出对象占常数个机器字,显式逐项输出至少需要 Ω(k) 时间。对象为长字符串或大整数时,还应计入写出的总位数。

数据结构查询还要拆开预处理与单次成本:先用 P(n) 时间建索引,再以 Q(n,k) 回答查询。不能把建索引隐藏在一个 O(log⁡n+k) 查询界里,也不能把一次建索引重复算作每次查询的必需工作。

输出敏感界不必是 O(n+k)。O(nh)、O(nlog⁡h) 与 O((n+k)log⁡n) 都可能表达有用的输出依赖;它们是否最优还需要同模型、同接口的下界。枚举复杂性中的“输出多项式时间”允许总时间是输入与输出长度的某个多项式,也不等于这些更精细的保证。

直觉

许多输入对象不会成为答案,却必须被读取或用于排除。分析时给这两类工作分别记账:定位、搜索或预处理支付“找到答案”的成本;每报告一项,再支付与这一项有关的成本。关键不是在公式末尾添一个 k,而是证明所有扫描都能记在少量定位步骤或真实输出上。

这也解释了为什么“小答案”有时仍不便宜。即使凸包只有三个极点,算法通常仍须查看所有输入点,确认没有遗漏一个外部点;输出敏感性会减少超过这份输入检查的额外工作,而不会自动免除读取输入。

例子与边界

区间查询中,停止比较记给谁 ​

在中心区间树的一个节点,中心为 10,驻留的闭区间按左端点排序为

[2,12], [7,11], [9,13], [10,14].

查询 q=8<10 时,所有驻留区间的右端点已经满足 high≥10>8。只需扫描左端点:

所读左端点 比较 动作
2 2≤8 报告 [2,12]
7 7≤8 报告 [7,11]
9 9>8 停止扫描

最后一项未被输出,但并非没有成本。每个访问节点至多有一次这样的停止比较,其余被读项都各自对应一个输出。随后只进入左子树,右子树的区间全在中心右侧,不可能含 8。若树高为 H,整条搜索路径的非输出工作为 O(H);每个区间只驻留一处,输出工作为 O(k),合计 O(H+k)。只有中心骨架确实保持 H=O(log⁡n) 时,才能代入得到 O(log⁡n+k)。

静态构造可用平衡的中心选择保证树高;动态版本还需说明固定坐标宇宙或重建策略。端点表支持有序顺次遍历也是该计费的实现条件:若每取下一项都重新从树根搜索,便平白加入了逐项对数因子。

不知道凸包大小时怎样猜 ​

考虑 n≥3 个互异平面点,先假定无三点共线,输出凸包的 h 个极点。Chan 的分组路线把 Jarvis wrapping 的“一步扫描全部点”换成“一步询问各组凸包”。给定步数预算 H 和组大小 m,执行以下过程:

  1. 分成 ⌈n/m⌉ 组,各组至多 m 点;分别计算只保留极点的局部凸包并按环序存储,总计 O(nlog⁡m)。
  2. 从全局极点出发;每一步在各局部凸包上求支持切线,以二分切线查询取得该组的候选,再从这些候选中选出下一个全局极点。若允许共线点,组内切线选择和组间候选比较都在从当前顶点出发方向相同的候选中选择距离最远者,跳过同一条边上的非极点。每步成本为 O((n/m)log⁡m)。
  3. 回到起点即输出整个环;执行 H 步仍未闭合则返回 incomplete,不能把这段前缀当作最终凸包。

每一步选的是支持整组、进而支持全部输入的下一条边,因此已生成的前缀始终是全局边界的一段。取 m=H≤n,成功或失败的一轮都至多耗费 O(nlog⁡H)。切线查询依赖凸多边形的环序和精确方向测试;这里把它作为几何子程序,不把任意数组二分搜索当作它的替代。

未知 h 时依次取

Ht=min{22t,n},t=1,2,….

设首次成功在第 T 轮。若不是首轮,上一轮失败保证 22T−1<h;所以 2T<2log2⁡h。每轮重新建局部凸包也无妨,因为

∑t=1TO(nlog⁡Ht)≤O(n∑t=1T2t)=O(nlog⁡h).

例如 h=20 时,预算为 4,16,256(受 n 截断),相应对数成本刻度是 2,4,8,总和仍与最后一项同阶。若每轮不复用预处理,却只把预算依次翻倍为 2,4,8,…,得到的账本是 n(1+2+⋯+⌈log2⁡h⌉),只证明 O(n(log⁡h)2);普通翻倍并不能直接替代上述双指数猜测。

上述最远点规则也处理只有部分点共线的情形。例如从 (0,0) 沿同一方向遇到 (1,0) 和 (2,0) 时,应直接选择 (2,0),否则会把边上的中间点误计入极点数 h。空集、单点和全共线情形单独处理;若只报极点,全共线输入输出至多两个端点。若要求全部共线边界输入点,输出参数就变了。把这些情形一起写入,可用 O(nlog⁡max{2,h}) 避免 h≤1 时的对数歧义。

推论与应用

Aho–Corasick 自动机的输出链把每个报告步骤记给一个匹配,因而扫描是 O(n+z);若只需每个模式的出现次数,可以累计状态访问数再沿失败树汇总,未必需要逐个生成全部 z 个位置。这正是报告接口与计数接口的差别。

扫描线报告 K 个交点时,事件处理还可能逐项需要树更新,故出现 O((n+K)log⁡n);输出敏感不保证每个结果只付常数。阅读一个新保证时,可以依次检查“输出是什么、什么工作能记给输出、还有多少未输出工作”。用上面的例子自测:指出 q=8 的第三次比较由树高项支付,再分别求出双指数猜测与普通翻倍的成本和;两个证明回答的是不同的额外开销。

参考资料
  • Timothy M. Chan,Optimal Output-Sensitive Convex Hull Algorithms in Two and Three Dimensions,Discrete & Computational Geometry 16(4),1996,pp. 361–368,§2 的分组、切线与双指数猜测,§4 的退化处理及复用改进。
  • Mark de Berg, Otfried Cheong, Marc van Kreveld, Mark Overmars,Computational Geometry: Algorithms and Applications,第 3 版,Springer,2008,§10.1 的 interval tree 与报告查询。
关系图谱6 个相邻概念 · 2 类关系

拖动节点调整位置。

显示关系

显示:依赖

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