精读笔记
Problem Setting
论文标题：Feasibility Determination for Subjective Probability Constraints（arXiv preprint / 2026）。
这篇论文实际处理的是 simulation-based R&S 里的 feasibility determination，但约束对象是概率：对每个有限备选系统 i、每个 constraint ℓ，需要判断 p_{iℓ} 是否不超过一个或多个主观阈值 h。观测不是一般连续输出，而是事件是否发生的 Bernoulli 输出。
困难点有三层。第一，Bernoulli 数据在小概率或大概率区域高度偏态，直接套正态均值约束程序需要 batch mean，而 batch size 一大，fully sequential 的样本节省基本被抵消。第二，subjective constraints 意味着同一 constraint 上有多个阈值，naive 做法会把每个阈值当成独立 feasibility problem，错误概率和样本量都会随阈值数膨胀。第三，决策者可能不是一次性给出所有阈值，而是看到结果后继续 tighten 或 relax，这和传统固定阈值 feasibility guarantee 有冲突。
关键矛盾是：统计有效性要求对所有 system × constraint × threshold 同时控制错误，而实际效率要求尽早停、复用样本，并且不能依赖粗糙的正态近似。
Motivation
已有路线不够的根源在于建模层级错了：它们把 Bernoulli probability constraint 当成普通 mean constraint 来处理，然后用 batching 去制造近似正态性。这在 p 接近 0 或 1、阈值也接近极端时尤其不自然，因为 Bernoulli 的离散结构和 odds 结构被抹掉了。
作者的核心观察是：既然约束本质上是 probability threshold comparison，就应该直接构造 Bernoulli 与 Bernoulli threshold 的比较，而不是估计 variance 后走 normal boundary。更进一步，subjective thresholds 并不是一堆无结构的标准，它们在同一 constraint 上有天然排序；如果 dummy threshold samples 被单调耦合，多阈值判定可以共享信息。
这篇填的缺口是：在 Bernoulli probability constraints + multiple subjective thresholds 下，给一个不依赖 batching 的 statistically valid first-pass procedure，并提供一个更贴近决策流程的 sequential-threshold heuristic。
Core Idea
核心想法很干净：把阈值 h 也随机化成一个 Bernoulli(h) dummy system，然后比较真实输出 Y 和 dummy 输出 I。若真实系统事件概率 p 小于阈值 h，随机游走 S(r)=∑(Y-I) 应有向下漂移；若 p 大于 h，则有向上漂移。于是 feasibility decision 就变成 hitting -H 或 +H 的 gambler’s ruin 问题。
这改变了建模方式：prior 是“估计 p 的均值并围绕 h 构造置信边界”，本文是“构造一个概率为 h 的参考过程，与系统过程做序贯胜负比较”。这个变化带来的 inductive bias 是 Bernoulli/odds-ratio 结构内生化：错误概率由 odds-ratio 和 hitting boundary 控制，而不是由 batch-normal approximation 控制。
和 prior 的本质区别不在 recycling 这个词，而在 reference object 从 deterministic threshold 变成 stochastic Bernoulli comparator；这让理论分析从 Brownian/normal approximation 回到 finite-state random walk，且多阈值可以通过 uniform coupling 形成单调信息流。
Method
1. Odds-ratio IZ：它解决的是 difference IZ 在概率边界附近不合适的问题。绝对差值在 p 接近 0 或 1 时尺度不对，而 odds-ratio 会自动把极端概率区的相同绝对差看得更重要。更关键的是，它正好给出 random walk 上下跳概率比，因此能直接进入 gambler’s ruin 公式。
2. Dummy Bernoulli threshold process：它解决的是“如何不用正态近似比较 p 和 h”。每个阈值 h 生成 I~Bernoulli(h)，然后用 ∑(Y-I) 的 hitting time 判定。这样不需要初始样本估计方差，也不需要 batch mean。核心变化是把未知 p 与固定 h 的估计问题变成两个 Bernoulli 流的序贯比较。
3. Boundary H 的选择：H 由 β_ℓ ≥ 1/(1+θ_ℓ^{H_ℓ}) 决定，本质是在 gambler’s ruin 的最不利配置下控制单 constraint/threshold 的错误概率。它是该方法统计保证的锚点。
4. 多阈值 coupling：同一 constraint 下所有 threshold 的 dummy Bernoulli 由同一个 uniform 生成，因此 h_1<h_2 时 I(h_1)≤I(h_2)。这使 correctness events 有嵌套关系，证明中只需管两侧最困难阈值，而不是对所有阈值做线性 union bound。这是 subjective constraints 可扩展的关键。
5. Multi-pass recycling heuristic：后续新增阈值时，算法保存历史 envelope vLB/vUB 和 LAST，试图复用过去样本判断新阈值。BRF_B 用历史 dummy Bernoulli 平均比较，BRF_N 用阈值常数比较，BRF_BN 两者取先决策者。这部分解决的是交互式阈值调整的样本效率，但理论上明显弱于 first-pass。
Key Insight / Why It Works
最核心的贡献是把 feasibility against threshold 还原成 Bernoulli-vs-Bernoulli 的 sequential comparison，并用 odds-ratio IZ 精确匹配 random walk 的 drift ratio。方法有效不是因为更复杂的 algorithm engineering，而是因为统计对象选对了：Bernoulli 数据不再被正态化，阈值也被放进同一 Bernoulli family 中。
多阈值效率来自 latent monotone structure，而不是单纯 recycling。对于同一 constraint，阈值越大 dummy success 越多；因此如果在较严格/较宽松边界上已经正确，其他阈值的正确性可以被包含。这解释了为什么 threshold 数量增加不必等比例增加错误预算或样本量。这里的关键 insight 可以迁移：当多个 standards/thresholds 有天然序关系时，应通过 coupling 保留单调性，而不是独立构造 confidence intervals。
BRF 相比 RF 的增益主要不是 scaling，而是 better inductive bias：RF 的正态假设在 Bernoulli 小概率处错位，batching 只是用样本量修补近似；BRF 直接利用离散随机游走结构，避免了这种错配。实验中 RF 在大 batch 下才能恢复 PCD，说明 prior 的瓶颈确实是 approximation layer，而不是实现细节。
Multi-pass 的收益则更像 memory reuse + test-time adaptivity。它根据当前可行集决定是否加阈值，本质上把“不必要的阈值组合”从计算图里剪掉。这部分很实用，但不是同等强度的理论贡献；它的收益很大程度依赖阈值策略是否刚好能快速 prune。若所有阈值最终都必须密集检查，增益会明显下降。
Relation To Prior Work
最接近的谱系有三条：Bernoulli selection 的 odds-ratio IZ / gambler’s ruin 方法，constrained R&S 的 feasibility determination，以及 Zhou et al. 的 subjective constraints / recycled observations。
与 Bernoulli selection 的关系：本文借用了 odds-ratio IZ 和 random walk absorption analysis，但目标从 selecting best Bernoulli population 变成 checking feasibility against thresholds。这个迁移是实质性的，因为 comparator 不再是另一个真实系统，而是一个由阈值生成的 dummy Bernoulli process。
与 constrained R&S / chance constrained selection 的关系：已有方法多关注固定 constraint set 或正态均值约束；即便处理 probability constraint，也常绕回 batching 或两阶段 selection。本文的新增信息是针对 Bernoulli probability constraints 直接给 feasibility 判定 procedure，并把 subjective thresholds 纳入同一保证。
与 Zhou et al. RF 的关系最关键。RF 已经有 subjective thresholds 和 recycling，但假设 normal observations，需要 tolerance ε 和 variance estimation。本文把同样的问题换成 Bernoulli-native procedure。看似只是替换分布假设，实质上改变了误差控制机制：从 Brownian/normal-style boundary 转为 finite-state random-walk boundary。
多 pass 部分更像 Zhou et al. sequentially added thresholds 思路的 Bernoulli 版本重组，创新性弱于 first-pass BRF；真正的新东西是 BRF 的 dummy Bernoulli comparator 与多阈值单调耦合证明。
Dataset / Evaluation
Evaluation 覆盖了三个层面：单 Bernoulli constraint 的 stopping time 与 RF 对比、多 pass heuristic 的有效性与样本效率、以及 inventory simulation 作为较真实的多系统多约束案例。任务覆盖对论文 claim 基本是对口的，因为 claim 本来就是 simulation procedure 的统计有效性和样本效率，而不是泛化到真实数据分布。
实验较好地验证了两个核心点：第一，RF 对 batch size 和 Bernoulli skewness 敏感，尤其小概率下需要很大 batch 才能达到 nominal PCD；第二，BRF/MPB 在 subjective thresholds 下能显著减少 observations，inventory 例子也说明在实际 simulation optimization setting 中可用。
但 evaluation 也有明显边界。多 pass heuristic 没有理论保证，实验只说明“这些设置下没破 PCD”，不能替代 validity proof。实验场景仍是有限、可控、terminating simulation；没有展示 rare-event probability 极小到工程上常见的 10^{-4}/10^{-6} 时是否仍高效。RF baseline 的适配需要 tolerance conversion 和 batch size 选择，这里虽然讨论充分，但公平性仍依赖这些设定。总体上，实验支持核心机制，但不能证明 multi-pass 的普适可靠性。
Limitation
第一，统计保证的边界很清楚：BRF first-pass / 所有阈值预先给定时有 PCD 保证；后续 sequentially added thresholds 的 MPB_B、MPB_N、MPB_BN 没有证明。论文建议若需要 validity 就把所有阈值放 first pass，这实际上削弱了 multi-pass 的理论地位。
第二，方法依赖 Bernoulli IID terminating replications。跨 constraint 可相关，但 replication 之间要 IID；对 steady-state simulation、dependent output、adaptive sampling with changing systems 等场景不直接成立。
第三，dummy Bernoulli 变量带来额外随机性。作者认为它换来了 simple random walk 和免方差估计，但这也可能在接近边界时增加判定波动。与 deterministic threshold + martingale concentration boundary 的比较文中未充分说明。
第四，odds-ratio IZ 是合理但强假设的 decision geometry。它默认 decision maker 对概率差异的敏感性按 odds 尺度变化；很多实际 chance constraints 可能更关心绝对风险差或非对称 loss。本文没有解决 IZ 选择本身的决策建模问题。
第五，scalability 上限来自最困难系统/阈值。若大量 p 接近 threshold 的 acceptable band 边界，或者阈值网格密集覆盖极端小概率区域，stopping time 仍可能很大。所谓多阈值扩展不随阈值数线性恶化，依赖单调耦合和“最多两个有效困难阈值”的证明结构；样本时间并非完全与阈值数无关。
第六，multi-pass 的大幅收益可能主要来自 threshold policy / pruning structure，而不是算法本身普遍更强。若决策者最终需要所有阈值的完整可行性表，MPB 相对 BRF 的优势会变小。
Takeaway
- 1. 对 Bernoulli probability constraints，不要先 batch 成 normal mean 再套传统 feasibility procedure；直接在 Bernoulli 层建模能同时获得更自然的误差控制和更好的样本效率。
- 2. Subjective thresholds 的关键结构是单调性。
- 通过 common-uniform coupling，把多个阈值组织成嵌套随机过程，是本文最值得迁移的 insight。
- 3. Odds-ratio IZ 不只是一个 tolerance choice，而是让 Bernoulli decision problem 与 gambler’s ruin 分析对齐的建模选择；它把概率尺度上的 difficulty 重新参数化了。
一句话总结
这篇论文把 subjective probability feasibility determination 从“Bernoulli 数据的正态近似问题”改造成“带 odds-ratio IZ 的 Bernoulli 随机游走比较问题”，其核心贡献是 Bernoulli-native 的非渐近有效性与多阈值单调复用。