乘法逆元 求组合数取模 (模为质数)

为什么求组合数取模，可以用 乘法逆元来做。

所以求 组合数取模 就是求   % mod.

对于正整数  和 p，如果有 ，那么把这个同余方程中  的最小正整数解叫做 a 模 p 的逆元。

即如果  % p = 1， 那么x的最小正整数解就是 a 的逆元。

 可以得到    % M () % M;    //费马小定理

所以  % M =   (  % M )   % M.

typedef long long LL;
const LL maxn(1000005), mod(1e9 + 7);
LL Jc[maxn];
void calJc()    //求maxn以内的数的阶乘
{
    Jc[0] = Jc[1] = 1;
    for(LL i = 2; i < maxn; i++)
        Jc[i] = Jc[i - 1] * i % mod;
}
/*
//拓展欧几里得算法求逆元
void exgcd(LL a, LL b, LL &x, LL &y)    //拓展欧几里得算法
{
    if(!b) x = 1, y = 0;
    else
    {
        exgcd(b, a % b, y, x);
        y -= x * (a / b);
    }
}
LL niYuan(LL a, LL b)   //求a对b取模的逆元
{
    LL x, y;
    exgcd(a, b, x, y);
    return (x + b) % b;
}
*/
//费马小定理求逆元
LL pow(LL a, LL n, LL p)    //快速幂 a^n % p
{
    LL ans = 1;
    while(n)
    {
        if(n & 1) ans = ans * a % p;
        a = a * a % p;
        n >>= 1;
    }
    return ans;
}
LL niYuan(LL a, LL b)   //费马小定理求逆元
{
    return pow(a, b - 2, b);
}
LL C(LL a, LL b)    //计算C(a, b)
{
    return Jc[a] * niYuan(Jc[b], mod) % mod
        * niYuan(Jc[a - b], mod) % mod;
}

 以上即为逆元求取组合数的方法，无论使用拓展欧几里得还是费马小定理，一开始求取Jc数组是的复杂度是 O(n),拓展欧几里得算法和费马小定理的复杂度均为 O(log(mod))，如果要求取m次组合数，则总的时间复杂度为 O(mn log(mod)). 

转载自：https://www.zybuluo.com/ArrowLLL/note/713749