组合数计算代码

最新推荐文章于 2021-12-04 18:12:32 发布

原创最新推荐文章于 2021-12-04 18:12:32 发布 · 867 阅读

0 ·

CC 4.0 BY-SA版权

杂谈专栏收录该内容

86 篇文章

订阅专栏

double lnchoose(int n, int m)
{

    if (m > n)

    {

        return 0;

    }
    if (m < n/2.0)
    {
        m = n-m;
    }

    double s1 = 0;
    for (int i=m+1; i<=n; i++)
    {
        s1 += log((double)i);
    }

    double s2 = 0;
    int ub = n-m;
    for (int i=2; i<=ub; i++)
    {
        s2 += log((double)i);
    }

    return s1-s2;
}

double choose(int n, int m)
{

    if (m > n)

    {

        return 0;

    }
    return exp(lnchoose(n, m));
}