ll fac[N];
void N_()
{
    fac[0]=1;
    for(ll i=1;i<N;i++)
    fac[i]=fac[i-1]*i%mod;
}
ll quick_mod(ll a,ll b,ll mo){ll ans=1;for(;b;b>>=1,a=a*a%mo)if(b&1)ans=ans*a%mo;return ans; }
ll C(ll n, ll m)
{
    if(m > n) return 0;
   return fac[n]*quick_mod(fac[m]*fac[n-m],mod-2)%mod;
}
ll Lucas(ll n, ll m)
{
    if(!m) return 1;
    return C(n % mod, m % mod) * Lucas(n / mod, m / mod) % mod;
}
//C(n,m)%p=C(n/p,m/p)*C(n%p,m%p)