c++ - 超大数的运算问题

Question

在一个问题中遇到的需要求这样的一个表达式：

应该如何计算可以得到结果？如果用龙贝格积分计算是算不出来的，伽马函数用定义算也求不出来。请问应该使用什么算法？

Answer 1

前面倒是好说，后面的积分不太好算呀

Answer 2

怀疑题目有笔误。左下角数字应该是251528。这样分式部分就正好是Beta函数的定义：

$$frac{Gamma(m+n)}{Gamma(m) Gamma(n)}=frac{1}{B(n,m)}$$

积分部分，是不完全Beta函数的定义：

$$int_0^{frac{1}{2}} theta ^{n-1} (1-theta )^{m-1} , dtheta=B_{frac{1}{2}}(n,m)$$

两者乘积叫做正则化不完全Beta函数（Regularized Beta Function）。

$$frac{B_{frac{1}{2}}(n,m)}{B(n,m)}=I_{frac{1}{2}}(n,m)$$

这个函数正好是Beta分布的累积分布函数（CDF），所以许多软件包都能计算。比如Excel：

结果并没有超出计算机浮点数的表示范围。