组合数模板

inv代表的逆元,逆元可以倒推出其他数的逆元
Comp(a,b) 就是C(a,b)
X 关于 mod 的逆元为 x^(mod - 2);

const int maxn=1e5+5;
const int mod=1e9+7;
ll inv[maxn],fac[maxn],ans[maxn],pos[maxn];
ll qpow(ll a,int b)
{
    ll ans=1;
    a%=mod;
    while(b)
    {
        if(b&1)
            ans=(ans*a)%mod;
        a=(a*a)%mod;
        b/=2;
    }
    return ans;
}
void pre()
{
    fac[0]=fac[1]=1;
    for(int i=2;i<maxn;++i) fac[i]=i*fac[i-1]%mod;
    inv[maxn-1]=qpow(fac[maxn-1],mod-2);
    for(int i=maxn-2;i>=0;i--) inv[i]=inv[i+1]*(i+1)%mod;
}
ll Comb(int n,int k)
{
    if (k < 0 || k > n) return 0;
    return fac[n]*inv[k]%mod *inv[n-k]%mod;
}
posted @ 2018-08-03 17:24  ffgcc  阅读(100)  评论(0编辑  收藏  举报