【lucas定理】BZOJ2982 combination

最新推荐文章于 2021-09-09 22:50:27 发布

linkfqy

最新推荐文章于 2021-09-09 22:50:27 发布

阅读量196

点赞数 1

分类专栏： BZOJ lucas定理常见OJ题解专栏我的OI历程文章标签： lucas定理 BZOJ

本文链接：https://blog.csdn.net/linkfqy/article/details/78724152

版权

我的OI历程同时被 3 个专栏收录

259 篇文章 1 订阅

订阅专栏

常见OJ题解专栏

193 篇文章 4 订阅

订阅专栏

BZOJ

107 篇文章 0 订阅

订阅专栏

题面在这里

不说了，直接lucas

示例程序：

#include<cstdio>

const int maxn=10010,MOD=10007;
int tst,fac[maxn],inv[maxn];
void prepare(){
    fac[0]=inv[0]=inv[1]=1;
    for (int i=1;i<=MOD;i++) fac[i]=fac[i-1]*i%MOD;
    for (int i=2;i<=MOD;i++) inv[i]=-(MOD/i)*inv[MOD%i]%MOD;
    for (int i=1;i<=MOD;i++) (inv[i]*=inv[i-1])%=MOD;
}
int C(int x,int y){
    if (x>y) return 0;
    if (x<MOD&&y<MOD)
     return fac[y]*inv[x]%MOD*inv[y-x]%MOD;
    return C(x/MOD,y/MOD)*C(x%MOD,y%MOD)%MOD;
}
int main(){
    prepare();
    scanf("%d",&tst);
    while (tst--){
        int n,m;scanf("%d%d",&n,&m);
        printf("%d\n",(C(m,n)+MOD)%MOD);
    }
    return 0;
}