组合数

根据这个公式求，不仅慢，而且中间值有阶乘会溢出，不采用

因此我们可以取对数来解决数太大的问题：

#include<cmath>
#include<iostream>
using namespace std;
double C(int n,int m){
    double lg=0;
    for(int i=m+1;i<=n;i++)
        lg+=log10(i);
    for(int i=1;i<=n-m;i++)
        lg-=log10(i);
    return pow(10,lg);
}
int main(){
    cout<<C(1000,500);
}

这个方法能处理大数速度也很快，但是又取对数，又是还原，难免存在精度问题，代码也比较长，来看下面简洁又不涉及精度的递推思路：

根据递推式写代码就很简单了：

int C(int n,int m){
    if(n==m||m==0)return 1;
    else return C(n-1,m)+C(n-1,m-1);
}

这个虽然代码短但是很慢，计算 C( 30 , 15 ) 居然都花了1秒，别提求大数的组合数了

来看下面这个递推式：

这个就非常快了，可以求大数的组合数，但大数的组合数会超过整型数据的范围，所以使用double类型

#include<iostream>
using namespace std;
double C(int n,int m){
    if(n==m||m==0)return 1;
    else return C(n,m-1)*(n-m+1)/m;
}
int main(){
    cout<<C(1000,500);
}

算得 C( 1000 , 500 ) = 2.70288e+299 速度非常快，时间小于0.1ms

组合数打表，就是打出杨辉三角，这个可以取模

#define N 1000
#define MOD 1000000007
long long C[N+10][N+10];
void Init(){
    for(int n=0;n<=N;n++){
        for(int m=0;m<=n;m++){
            if(n==m||m==0)C[n][m]=1;
            else C[n][m]=C[n-1][m]%MOD+C[n-1][m-1]%MOD;
        }
    }
}

还可以对组合数公式取对数

int n,m;
cin>>n>>m;
double lg=0;
for(int i=m+1;i<=n;i++)
    lg+=log10(i);
for(int i=1;i<=n-m;i++)
    lg-=log10(i);
int a=(int)lg;//整数部分，决定小数点位置
double b=lg-(int)lg;//小数部分，决定每一位上的数字是多少
cout<<pow(10,b)<<"e+"<<a<<endl;//这样可以就不会溢出了