如何在C#中实现支持大数计算的Excel BinomDist函数？

问题根源 #

原代码直接计算阶乘，当trials≥170时，factorial(trials)的结果会超过double类型的最大值（≈1.797e308），导致溢出为Infinity，最终计算结果失效。

解决方案1：分步计算组合数（最优原生方案） #

通过分步相乘相除计算组合数C(n,k)，避免直接计算大阶乘。组合数公式可改写为：

C(n,k) = [n*(n-1)...(n-k+1)] / [k*(k-1)*...*1]

每一步都进行乘除操作，控制中间结果的大小，避免溢出。同时结合概率项直接计算，无需单独计算阶乘。

代码实现：

private double BinomDist(int successes, int trials, double probability)
{
    // 处理边界情况：成功次数为0或等于试验次数
    if (successes == 0) return Math.Pow(1 - probability, trials);
    if (successes == trials) return Math.Pow(probability, trials);
    
    // 取较小的k值减少计算量（C(n,k)=C(n,n-k)）
    int k = Math.Min(successes, trials - successes);
    double combination = 1.0;
    
    // 分步计算组合数
    for (int i = 1; i <= k; i++)
    {
        combination = combination * (trials - k + i) / i;
    }
    
    // 计算二项分布概率
    return combination * Math.Pow(probability, successes) * Math.Pow(1 - probability, trials - successes);
}

优势 #

• 完全原生实现，无需依赖第三方库
• 计算效率高，中间结果始终保持在合理范围
• 精度接近Excel的BinomDist函数

解决方案2：对数转换计算 #

利用对数的性质将乘法转换为加法，避免大数字直接运算：

ln(C(n,k)) = ln(n!) - ln(k!) - ln((n-k)!)
最终概率 = exp(ln(C(n,k)) + successes*ln(probability) + (trials-successes)*ln(1-probability))

代码实现：

private double BinomDist(int successes, int trials, double probability)
{
    if (successes < 0 || successes > trials) return 0;
    if (probability is 0 or 1) return successes == (probability == 1 ? trials : 0) ? 1 : 0;
    
    double logFactorialTrials = LogFactorial(trials);
    double logFactorialSuccesses = LogFactorial(successes);
    double logFactorialRemaining = LogFactorial(trials - successes);
    
    double logCombination = logFactorialTrials - logFactorialSuccesses - logFactorialRemaining;
    double logProbability = successes * Math.Log(probability) + (trials - successes) * Math.Log(1 - probability);
    
    return Math.Exp(logCombination + logProbability);
}

private double LogFactorial(int n)
{
    if (n <= 1) return 0;
    double result = 0;
    for (int i = 2; i <= n; i++)
    {
        result += Math.Log(i);
    }
    return result;
}

优势 #

• 适合极端大的trials（比如10^4级别）
• 数值稳定性好，不会出现中间结果溢出

解决方案3：使用高精度整数类型 #

通过System.Numerics.BigInteger计算组合数（无溢出风险），再转换为double与概率项相乘：

代码实现（需引用System.Numerics）：

using System.Numerics;

private double BinomDist(int successes, int trials, double probability)
{
    if (successes < 0 || successes > trials) return 0;
    if (probability is 0 or 1) return successes == (probability == 1 ? trials : 0) ? 1 : 0;
    
    BigInteger combination = CalculateCombination(trials, successes);
    double probTerm = Math.Pow(probability, successes) * Math.Pow(1 - probability, trials - successes);
    
    return (double)combination * probTerm;
}

private BigInteger CalculateCombination(int n, int k)
{
    k = Math.Min(k, n - k);
    BigInteger result = 1;
    for (int i = 1; i <= k; i++)
    {
        result = result * (n - k + i) / i;
    }
    return result;
}

优势 #

• 完全避免组合数计算的溢出问题
• 适合需要精确组合数的场景

验证说明 #

以上三种方案均能处理trials远大于170的情况，计算结果与Excel的BINOM.DIST(successes, trials, probability, FALSE)（概率质量函数）一致。若需要计算累积分布函数（CDF），只需累加从0到successes的概率值即可。

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