组合数取模(卢卡斯定理)

问题描述：

给n,m,p，输出C(n,m)%p的值

思路：

卢卡斯定理

int lucas(int n,int m){//C(n,m)%mod
    if(m==0)return 1;
    return C(n%mod,m%mod)*lucas(n/mod,m/mod)%mod;
}

code:

递归

//https://loj.ac/problem/10228
#include<cstdio>
#include<iostream>
using namespace std;
#define int long long
int n,m,mod;
int ppow(int a,int b){
    int ans=1%mod;
    a%=mod;
    while(b){
        if(b&1){
            (ans*=a)%=mod;
        }
        (a*=a)%=mod;
        b>>=1;
    }
    return ans;
}
int inv(int a){//逆元
    return ppow(a,mod-2);
}
int C(int n,int m){
    if(m>n)return 0;
    int a=1,b=1;
    for(int i=n-m+1;i<=n;i++){
        (a*=i)%=mod;
    }
    for(int i=1;i<=m;i++){
        (b*=i)%=mod;
    }
    return a*inv(b)%mod;
}
int lucas(int n,int m){
    if(m==0)return 1;
    return C(n%mod,m%mod)*lucas(n/mod,m/mod)%mod;
}
signed main(){
    int T;
    cin>>T;
    while(T--){
        cin>>n>>m>>mod;
        cout<<lucas(n,m)<<endl;
    }
    return 0;
}