【FZUOJ 2020】 组合

lucas的应用!!
中间组合数学开始竟然看不懂!!!! 约掉了一部分啊!
还用到 逆元 欧拉定理之类的(就一句话 没有名字这样高大上)

#include <cstdio>
#define LL long long
LL F(LL x,LL y,LL mo)
{
    LL ans=1;
    while(y)
    {
        if(y&1) ans=ans*x%mo;
        x=x*x%mo;
        y>>=1;
    }
    return ans;
}
LL C(LL n,LL m,LL p)
{
    LL x=1,y=1;
    for(LL i=1;i<=m;i++)
        x=x*(n-i+1)%p,y=y*i%p;
    return x*F(y,p-2,p)%p;
}
LL Lucas(LL n,LL m,LL p)
{
    if(m==0) return 1;
    return C(n%p,m%p,p)*Lucas(n/p,m/p,p)%p;
}
int main()
{
    int T;scanf("%d",&T);
    while(T--)
    {
        LL m,n,p;
        scanf("%I64d %I64d %I64d",&m,&n,&p);
        printf("%I64d\n",Lucas(m,n,p));
    }
    return 0;
}

你可能感兴趣的:(【FZUOJ 2020】 组合)