从恒等式到算法：五步与复杂度

• 命题框：把积分恒等式变成 LCU 电路，核心指标 ‖c‖₁ = 1+O(ε)——LCU 开销不爆炸，成功概率只含物理衰减
• 五步流水线（每步带公式与数值）：

1. 截断 [−K, K]：误差 = 1 − (2/π)arctan K；演算 K = 1 → 误差 = 1 − (2/π)(π/4) = 0.5；由 arctan(1/K) ≈ 1/K，误差 ≈ 2/(πK) ≤ 0.01 → K ≥ 200/π ≈ 64

2. 梯形离散化：cⱼ = wⱼ/(π(1+kⱼ²)) > 0 全正（免相位修正）；‖c‖₁ ≈ ∫_{−K}^{K} dk/(π(1+k²)) = (2/π)arctan K = 1 − O(1/K)

3. 实现酉因子：等差结构 e^{−iLkⱼs} = e^{iLKs}·(e^{−iL(2K/M)s})ʲ——二进制分解后受控选择仅需 O(log M) 次受控幂次

4. 组装 LCU：O_coef → 受控演化序列 → O_coef†

5. 后选择 + 振幅放大：P ≈ ‖u(T)‖²/‖u₀‖²（纯耗散例 ‖u(T)‖ = e^{−γT}‖u₀‖ → P ≈ e^{−2γT}）；放大轮数 O(q)，q = ‖u₀‖/‖u(T)‖

• 演算（教材例一，完整数值）：A = diag(1, 2, 4)、T = 1、u₀ = (1,1,1)ᵀ/√3：

1. 逐分量 Poisson 核：e^{−At} = diag(e^{−t}, e^{−2t}, e^{−4t})

2. ‖u(T)‖² = (e^{−2} + e^{−4} + e^{−8})/3 ≈ (0.1353 + 0.0183 + 0.0003)/3 ≈ 0.0513

3. ‖u(T)‖ ≈ 0.227 → q ≈ 4.4 → 振幅放大约 5 轮

• 结果框：总查询 Õ(Γₚ^{1+1/p}·q^{2+2/p}·T^{1+1/p}/ε^{1+2/p})；态制备 O(q) 匹配下界 Ω(q)（亚线性则蕴含 BQP = PP）