卢卡斯定理，费马小定理的应用

牛客多校第八场G：传送门

 

卢卡斯定理：C(n,m)%p=C(n/p,m/p)*C(n%p,m%p)%p，组合数的取模

费马小定理：当p为质数时候， a^(p-1)≡1(mod p)， a*a^(p-2)≡1(mod p)

费马小定理求解逆元：(a/b)%m 时，因b可能会过大，会出现爆精度的情况，所以需变除法为乘法：

设c是b的逆元，则有b*c≡1(mod m)；      b*b^(p-2)≡1(mod p)，故b^(p-2)为b的逆元

则(a/b)%m = (a/b)*1%m = (a/b)*b*c%m = a*c(mod m);

#include <iostream>
using namespace std;
#define p 1000000007
#define ll long long
 
ll qpow(ll a,ll b){    
    ll ans=1;
    a%=p;
    while(b){
        if(b&1)ans=ans*a%p;
        b>>=1;
        a=a*a%p
    }
    return ans;
}
ll C(ll n,ll m){
    if(n<m)return 0;
    ll ans=1;
    for(ll i=1;i<=m;i++)
        ans=ans*((n-m+i)%p*qpow(i,p-2)%p)%p;//费马小定理求逆元 
    return ans;
}
 
ll lacus(ll n,ll m){
    if(m==0)
        return 1;
    return (lacus(n/p,m/p)*C(n%p,m%p))%p;
}
int main(){
    ll n;
    cin>>n;
    cout<<(lacus(n-1,2)+lacus(n,4))%p;
    return 0;
}