计算组合数

编写函数，参数是两个非负整数n和m，返回组合数Cnm=n!(m!(n-m)!)，其中m<=n<=25,例如，n=25，m=12时答案为5200300。

 

 

很多人一看到这个题目，第一做法是编写一个计算阶乘的函数，然后返回值。其做法如下

#include <stdio.h>
long long factorial(int n){
    long long m=1;
    for(int i=1;i<n;i++)
    m*=i;
    return m;
}
int main (void)
{
    int m,n;
    scanf("%d %d",&n,&m);
    printf("%d",factorial(n)/(factorial(m)*factorial(m-n)));
    return 0;
}

如果按这样计算，要计算m的阶乘和n的阶乘，这样很容易溢出，我们需要对现有的公式进行化简之后再编写代码。

一个简单的方法是利用n！/m！=(m+1)(m+2)(m+3)......(n-1)n，这样可以很大幅度的减少溢出可能。改良后的完整代码如下

#include <stdio.h>
long long C(int n,int m){
    if(m<n-m) m=n-m;
    long long ans=1;
    for(int i=m+1;i<=n;i++) ans*=i;
    for(int i=1;i<=n-m;i++) ans/=i; 
    return ans;
}
int main (void)
{
    int m,n;
    scanf("%lld %lld",&n,&m);
    printf("%lld",C(n,m));
    return 0;
}

这里的函数C中有一个技巧，if(m<n-m) m=n-m;很多人会想这样做的目的是什么，改变了参数的值难道不会影响到最后的结果吗？

 

我们可以看到Cnm=n!(m!(n-m)!)，如果当m<n-m,那么原式中m就等于n-m，n-m就等于n-(n-m)=m,那么此时的原式Cnm=n!(m!(n-m)!)就变成了Cnm=n!((n-m)!m!)，简单的来说就是把分母中两个数的位置换了一下，把较大的一位换到了前面。举一个极端一点的例子就可以看出这样做的优越性。

我们算一下C25 1的值，在25个东西里面选出一个东西的选法，很简单，有25种选法。

如果我们按照原式执行函数，C25 1=25!/(1!*24!),我们要计算出2*3*4*5*......25的结果，然后再用这个结果除以1，除以2一直除到24为止。

如果我们把分母中两个数换一个位置，把较大的24换到前面来再执行一下函数C，再来看一下运行过程，函数会从m+1一直乘到n，此时m=n-m=24，也就是从25乘到25，结果为25，然后用这个结果除以1，除以2，一直除到n-m为止，但是此时的n-m就等于1，所以25除以1就等于最后的结果。输出为25。

新手自学，纯属手打，不知道能不能帮助到不能理解这个小技巧的童鞋。