Lucas 定理学习小结，组合数取模

图片中解释了Lucas定理，把n,m表示为质数p的p进制形式，然后各项系数的组合数的乘积就是C(n,m)%p的结果。

对于n,m比较大的时候，质数p不太大的时候可以使用Lucas计算。
具体的细节证明可以查阅维基，或者参考挑战上的介绍。

当然写熟练了Lucas也可以改成非递归的，也可以把C(n,m)合并到Lucas里面一块计算了。

小技巧，一般直接这样计算有可能时间不够，所以通常的做法是先打出阶乘表和逆元表，就可以大大简化代码，以及节约时间，比如这题，here

以下是代码，题目链接

#include<cmath>
#include<map>
#include<set>
#include<cstdlib>
#include<vector>
#include<cstdio>
using namespace std;
#define cl(a,b) memset(a,b,sizeof(a))
#define LL long long
#define pb push_back
#define gcd __gcd

#define For(i,j,k) for(int i=(j);i<k;i++)
#define lowbit(i) (i&(-i))
#define _(x) printf("%d\n",x)

const int maxn = 1e6+200;
const int inf  = 1 << 28;
LL p ;
LL qsm(LL a, LL b){//因为p是质数所以可以使用快速幂计算逆元，p不是质数的时候可以使用扩展欧几里得
    LL ans = 1;
    while(b){
        if(b&1)ans = ans*a%p;
        a = a*a%p;
        b>>=1;
    }
    return ans%p;
}
LL C(LL n,LL m){//计算C(n,m)%p
    if(m>n)return 0;
    LL ans = 1;
    for(int i=1;i<=m;i++){
        LL a = (n-m+i) % p;
        LL b = i % p;
        ans = ans * (a * qsm(b,p-2) % p ) % p;
    }
    return ans;
}
LL Lucas(LL n,LL m){//p进制拆分
    if(m==0)return 1;
    return C(n%p,m%p) * Lucas(n/p, m/p)%p;
}

int main(){
    int T;scanf("%d",&T);
    while(T--){
        LL n,m;
        scanf("%I64d%I64d%I64d",&n,&m,&p);
        printf("%I64d\n",Lucas(n,m)%p);
    }
    return 0;
}