快速幂求组合数

利用乘法逆元

乘法逆元：
(a/b)% mod=a*(b^(mod-2)) % mod为素数

ll quick(ll a,ll b){ //计算a^b并且对mod取余
    ll ans = 1;
    while(b){
        if(b & 1) ans = ans * a % mod;
        a = a * a % mod;
        b >>= 1;
    }
    return ans % mod;
}
ll calc(ll n,ll m)
{//求组合数
	if(m > n) return 0;
	return (fac[n] * quick(fac[m], mod - 2) % mod * quick(fac[n - m], mod - 2) % mod) % mod;
}