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)
Lucas
發表評論
所有評論
還沒有人評論,想成為第一個評論的人麼? 請在上方評論欄輸入並且點擊發布.