AcWing 886. 求组合数 II（乘法逆元）

最新推荐文章于 2022-07-08 17:04:55 发布

罚时大师月色

最新推荐文章于 2022-07-08 17:04:55 发布

阅读量287

点赞数

分类专栏：数论 acwing 零基础学习C++

本文链接：https://blog.csdn.net/qq_52358098/article/details/113917516

版权

组合数快速幂大数模运算乘法逆元费马小定理

关键词由CSDN通过智能技术生成

acwing 同时被 3 个专栏收录

128 篇文章 0 订阅

订阅专栏

数论

27 篇文章 1 订阅

订阅专栏

零基础学习C++

23 篇文章 2 订阅

订阅专栏

求组合数
在这里插入图片描述
这道题和上一篇很像，但是这道题的数据范围也变大了，那我们该怎么求呢，我们要求i(1-N)的 i! 和 i!^-1然后根据
C^b_a=a!/(b!*(a-b)!)这样求就会可以很方便的求组合数了，时间复杂度是N(log₂N）因为求快速幂的时间复杂度是log₂N级别的，输出结果是o(1)的
然后我们怎么求 i!，这个很简单，类似于递归的方式调用前一个数组，
但是难点是怎么求 i!^-1呢，我们还是要用到一个定理。
在这里插入图片描述
b是底数，m是模，如果b,m互质b^-1等于b^m-2%m这是类似于费马小定理的递推。
然后思路很好解决了，
我们求i的乘法逆元求i^-i
代码如下

#include<iostream>
#include<cmath>

using namespace std;

const int Mod=1e9+7;
const int N=1e5+1;

int fact[N],infact[N];

int quick_pow(int a,int k,int mod)
{
    int res=1;
    while(k)
    {
        if(k&1) res=res*1ll*a%mod;
        a=a*1ll*a%mod;
        k>>=1;
    }
    return res;
}


int main(void)
{
    infact[0]=fact[0]=1;
    for(int i=1;i<N;i++) fact[i]=(long long)fact[i-1]*i%Mod;
    for(int i=1;i<N;i++) infact[i]=(long long)infact[i-1]*quick_pow(i,Mod-2,Mod)%Mod;
    int n;
    scanf("%d",&n);
    while(n--)
    {
        int a,b;
        scanf("%d%d",&a,&b);
        printf("%d\n",(long long)fact[a]*infact[b]%Mod*infact[a-b]%Mod);
    }
}