为何scipy.special.comb与math.comb的组合数计算结果存在差异？

• math.comb 是Python标准库为精确整数计算设计的实现，全程采用任意精度整数运算，返回结果为原生int类型，计算过程无近似、无截断，和数学定义的组合数结果完全对齐。
• scipy.special.comb 是为大规模科学计算场景优化的实现，默认配置下走浮点加速路径，返回结果为双精度浮点数（float64），不会对所有整数结果做精确还原，计算过程天然受浮点精度限制。

• 当n=25时：所有组合数的总和为2^25 = 33554432，单组最大组合数C(25,12)=5200300，所有数值都远小于双精度浮点数可无误差表示整数的上限（2^53 ≈ 9e15）。这种场景下scipy计算出的每一项浮点值，和math返回的整数转同精度浮点的结果完全相等，逐次累加后两边结果自然匹配，等值判断返回True。
• 当n=250时：所有组合数的总和为2^250，单组最大组合数C(250,125)的数值规模已经远超2^53，超出了双精度浮点数的精确整数表示范围。此时scipy在计算大数值组合数时，每一项都会产生不同程度的浮点舍入误差，误差逐次叠加后，最终的求和结果和math库全程精确整数计算的结果出现可见偏差，等值判断自然返回False。

scipy.special.comb 本身提供了exact参数，只有显式指定exact=True时，才会切换到整数精确计算路径返回int类型结果；默认的exact=False配置是面向大规模批量计算的性能优化，用极小的精度损失换取更高的运算速度。如果调用时加上该参数：

sum(scipy.special.comb(250, i, exact=True) for i in range(0, 251))

得到的结果会和math库的求和值完全一致。

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