组合数取余[C(n,m) %p]

如果 n， m很小（不超过50），可以用C++
的库函数 double tgamma(double x) ，这是一个欧拉积分.

ll C(ll n,ll m){
    return (ll)round(tgamma(n+1)/tgamma(m+1)/tgamma(n-m+1));
}

如果 n, m 不大，可以开 O(n^2) 的空间，可以利用杨辉恒等式来预处理组合数表

const ll mo=1e9+7;
ll C[1005][1005];
void getC(int n){
    for(int i=0;i<=n;i++){
        for(int j=0;j<=i;j++){
            if(j==0 || j==i)
                C[i][j]=1;
            else
                C[i][j]=(C[i-1][j-1]+C[i-1][j])%mo;
        }
    }
}
。

如果 n， m比较大，可以开 O(n)的空间，可以利用逆元来求解，当然，要保证 p 是素数

const ll mo=1e9+7;
ll C(ll n,ll m){
    static ll M=0,inv[N],mul[N],invMul[N];
    while(M<=n){
        if(M){
            inv[M]=M==1?1:(mo-mo/M)*inv[mo%M]%mo;
            mul[M]=mul[M-1]*M%mo;
            invMul[M]=invMul[M-1]*inv[M]%mo;
        }
        else mul[M]=1,invMul[M]=1;
        M++;
    }
    return mul[n]*invMul[m]%mo*invMul[n-m]%mo;
}