摘自队友devil的lucas模板

Posted 2020-07-10 very_czy

typedef long long LL;
const int MOD=1e9+7;
LL quick_mod(LL a,LL b)
{
LL ans=1;
a%=MOD;
while(b)
{
if(b&1)
{
ans=ans*a%MOD;
b--;
}
b>>=1;
a=a*a%MOD;
}
return ans;
}
LL C(LL n,LL m)
{
if(m>n) return 0;
LL ans=1;
for(int i=1; i<=m; i++)
{
LL a=(n+i-m)%MOD;
LL b=i%MOD;
ans=ans*(a*quick_mod(b,MOD-2)%MOD)%MOD;
}
return ans;
}
LL Lucas(LL n,LL m)
{
if(m==0) return 1;
return C(n%MOD,m%MOD)*Lucas(n/MOD,m/MOD)%MOD;
}