组合数学--逆元求组合数

逆元求解

公式大概就是这么个公式。具体的实现代码如下：

#include<iostream>
#include<cstdio>
#include <stdio.h>
#include<algorithm>
#include<cstring>
#include<cmath>
#include<cstdlib>
#include<queue>
#include<map>
#include<vector>
#include <set>
#define ll long long
using namespace std;
const ll mod=1e9+7;
const ll N=1e5+5;
ll powm(ll a,ll b,ll mod)//快速幂
{
    ll ans=1;
    while(b)
    {
        if(b&1)
        {
            ans=ans*a%mod;
        }
        b=b>>1;
        a=a*a%mod;
    }
    return ans;
}
ll inv(ll a,ll mod)//求逆元
{
   return powm(a,mod-2,mod);
}
ll A(int n,int m)//全排列公式
{
	return fac[m]*inv(m-n)%mod;
}
 
ll C(int n,int m)//组合数公式
{
	return fac[m]*inv(n)%mod*inv(m-n)%mod;
}


ll fac[N];
int main()
{
    ll n,m,k;
    cin>>n>>m;
    fac[0]=1;
    for(int i=1;i<=n;i++)
    {
        fac[i]=fac[i-1]*i%mod;//计算阶乘
    }
    //C(n,m) = n!*(m!%mod的逆元)*((n-m)!%mod的逆元)%mod
    ll res=fac[n]*inv(fac[m],mod)%mod*inv(fac[n-m],mod)%mod;
    cout<<res<<endl;
    return 0;
}