组合数取模（卢卡斯定理模版）直接套用即可

#include  
#include  
#include  
#include  
#include  
#include  
#define LL  __int64 
using namespace std;  
LL PowMod(LL a,LL b,LL MOD){  //费马小定理求逆元
    LL ret=1;  
    while(b){  
        if(b&1) ret=(ret*a)%MOD;  
        a=(a*a)%MOD;  
        b>>=1;  
    }  
    return ret;  
}  
LL fac[100005];  
LL Get_Fact(LL p){  
    fac[0]=1;  
    for(LL i=1;i<=p;i++)  
        fac[i]=(fac[i-1]*i)%p;  //预处理阶乘 
}  
LL Lucas(LL n,LL m,LL p){  
    LL ret=1;  
    while(n&&m){  
        LL a=n%p,b=m%p;  
        if(a

或者

#include 
#include 
#define ll unsigned long long int
const ll mod=1000000007;
using namespace std;
ll pow(ll a,ll b )//用迭代法求快速幂a^b%mod，即逆元
{
    ll ans=1,ret=a;
    while(b)
    {
        if(b&1) ans=ans*ret%mod;
        ret=ret*ret%mod;
        b>>=1;
    }
    return ans;
}
ll c(ll n,ll m)//用费马小定理求c（n,m)
{
    ll fac=1;
    for(ll i=1;i<=m;i++)
    {
        fac=fac*(n-m+i)%mod;
        fac=fac*pow(i,mod-2)%mod;//每一次都要求快速幂，这样最后就是m！的mod-2次方了，也可先求之前列出的费马小定理的公式的左半部分的阶乘，再求右半部分阶乘的快速幂
    }
    //cout<