O( nlogn )预处理， O( 1 ) 求组合数模板

最新推荐文章于 2022-03-08 20:36:06 发布

勇敢铁铁

最新推荐文章于 2022-03-08 20:36:06 发布

阅读量1k

点赞数

分类专栏：常用知识

本文链接：https://blog.csdn.net/qq_41105401/article/details/81711414

版权

常用知识专栏收录该内容

18 篇文章 0 订阅

订阅专栏

typedef long long LL;
const int maxn = 5e5;
const long long mod = 1e9 + 7;
LL A[maxn];
LL B[maxn];
void Init()
{
    A[0] = 1;
    for(int i=1; i<=maxn; i++){
        A[i] = (A[i-1] * i ) % mod;
    }
}

LL Ext_Gcd(LL a, LL b, LL &x, LL &y)
{
    if(b==0){
        x=1;
        y=0;
        return a;
    }
    LL d = Ext_Gcd(b, a%b, y, x);
    y-=a/b*x;
    return d;
}

LL Inv(LL a, LL n)
{
    LL x,y;
    LL d = Ext_Gcd(a,n,x,y);
    if(d == 1)
        return ((x%n)+n)%n;
    return -1;
}
LL get()
{
    for(int i=0;i<maxn;i++)
        B[i] = Inv(A[i],mod);
}
LL C(LL a, LL b, LL mod)
{
    if(a < b)
        return 0;
    return (A[a] * B[b] %mod)*B[a-b]%mod;
}