Lucas定理模板

最新推荐文章于 2022-10-20 19:24:36 发布

Megumin

最新推荐文章于 2022-10-20 19:24:36 发布

阅读量373

点赞数

CC 4.0 BY-SA版权

分类专栏： ACM模板

本文链接：https://blog.youkuaiyun.com/Megumin/article/details/51920294

ACM模板专栏收录该内容

1 篇文章

订阅专栏

输入p时，预处理所有n!的模p的余数，即调用init(p);

LL f[maxn];
void init(LL p)
{
    LL i;
    f[0]=1;
    for(i=1;i<=p;i++){
        f[i]=f[i-1]*i%p;
    }
}
LL powmod(LL a,LL b,LL p) {
    LL res=1;
    while(b!=0){
        if(b&1) res=(res*a)%p;
        a=(a*a)%p;
        b>>=1;
    }
    return
     res;
}

LL Lucas(LL n,LL m,LL p)
{
    LL ans=1;
    while(n&&m){
        LL nn=n%p,mm=m%p;
        if(nn<mm) return 0;
        ans=ans*f[nn]*powmod(f[mm]*f[nn-mm]%p,p-2,p)%p;
        n/=p;
        m/=p;
    }
    return ans;
}