IEEE 754中是否存在亚O(n)方法精确复现浮点重复加法？

具体而言，有没有更快的方式，能精确复现以下伪代码的执行结果：

f64 notQuiteFMA(f64 acc, f64 f, bigint n) {
    for (bigint i = 0; i < n; i++) {
        acc = acc + f;
    }
    return acc;
}

在实数运算中，直接返回acc + f*n即可，但有限精度运算下情况要复杂得多。比如notQuiteFMA(0, 1, 1 << 256)的结果约为2^{53，而非无限精度下预期的2}256，这也凸显了加速该计算的必要性。

• 注1：已知可通过含2^128项的查表法实现O(1)时间复杂度，但此类极端算法不在讨论范围内；
• 注2：需在当前舍入模式下实现位级精确（暂不考虑非信号NaN）；
• 注3：上述代码为伪代码，仅代表双精度浮点加法序列，不考虑编程语言的高精度/SIMD优化等情况；
• 注4：64位浮点数的尾数为52位，当x约为2^53时，x + 1 == x，具体值取决于舍入模式；
• 注5：假设大整数对有界常数取模为O(1)时间复杂度。

内容的提问来源于stack exchange，提问作者TLW