1.防溢出:
如果直接用c(n,m)=n!/(n-m)!*m! 来编程很可能会在算n!时就爆了long long,所以每一步最好把除分母也算上。
所以对于c(n,m)来说取m=min(m,n-m)来算c(n,m)=n*(n-1)*(n-2)*...*(n-m+1)/m*(m-1)*(m-2)*...*1 。显然分子分母都是m项相乘,从后往前去算:
先算(n-m+1)/1 ,then * (n-m+2)/2 ,then....,then * n/m 这个时候也满足每一步的相除都是整除,因为任何连续的m个整数一定会出现一个因数m。
long long com(long long n,long long m)
{
long long ans=1;
for(int i=min(m,n-m);i>0;i--)
ans*=(n-i+1)/(1+min(m,n-m)-i);
return ans;
}
同时也满足m=0,或者n=0的情况。
如果要约定m>n或m<0 时c(n,m)=0,那上面代码就要添上几句。
2.高效:
如果多次用到连续的组合数,不用一个一个计算,可以用迭代计算:[i][j]=(c[i-1][j]+c[i-1][j-1])
有组合数的关系:c(i,j)=c(i-1,j)+c(i-1,j-1);所以可以时空权衡开个二维数组储存结果:
c[0][0]=1;
for(int i=1;i<=limax;i++)
{
c[i][0]=1;
for(int j=1;j<=i;j++)
c[i][j]=(c[i-1][j]+c[i-1][j-1])%mod;
}
未完待续...