lucas定理模板

求c(n,m) %p 其中p为素数

#include <cstdio>
#include <algorithm>
using namespace std;
typedef long long LL
int T,n,m;
int fac[N],inv[N],invf[N];
inline int cmb(int a,int b)
{
    if(a<b)return 0;
    if(!b)return 1;
    return fac[a]*invf[b]%mod*invf[a-b]%mod;
}
int lucas(LL a,LL b)
{
    if(!b)return 1;
    return lucas(a/mod,b/mod)*cmb(a%mod,b%mod)%mod;
}

void init()
{
    fac[0]=fac[1]=1;
    invf[0]=invf[1]=1;
    inv[1]=1;
    for(int i=2;i<N;++i)
    {
        fac[i]=fac[i-1]*i%mod;
        inv[i]=inv[mod%i]*(mod-mod/i)%mod;
        invf[i]=invf[i-1]*inv[i]%mod;
    }
}