费马小定理求c(m,n)取模

ll fac[10000005];//阶乘数组
ll p(ll a,ll b){
    ll c=1;
    while(b)
    {
        if(b&1) c=c*a%mod;
        b>>=1;
        a=a*a%mod;
    }
    return c;
}
ll work(ll m,ll i)
{
    ll ans=((fac[m]%mod)*(p(fac[i]*fac[m-i]%mod,mod-2)%mod))%mod;//费马小定理求逆元
    return ans;
}