Lucas定理

总结:

C\bigl(\begin{smallmatrix} m\\n \end{smallmatrix}\bigr)\equiv C\bigl(\begin{smallmatrix} m/p\\n/p \end{smallmatrix}\bigr) * C\bigl(\begin{smallmatrix} m%p\\n%p \end{smallmatrix}\bigr)(mod \quad p)


为什么要用Lucas?

费马小定理求逆元要求n与mod互质,当n很大时,可能出现n=k*mod这种情况,因此不能直接通过求逆元求解组合数。


模版:

#include
using namespace std;
#define ll long long
int mod ;
ll ksm(ll a,ll b){
    ll kase = 1;
    while(b){
        if(b & 1){
            kase = kase * a % mod;
        }
        a = a * a % mod;
        b >>= 1;

    }
    return kase ;
}
ll ni(ll x){
    return ksm(x , mod - 2);
}
int main(){
    ios::sync_with_stdio(false);
    cin.tie(0);
    int T ;
    cin>>T;
    while(T--){
        int n , m ;
        cin>>n>>m>>mod;
        vectorfac(n + m + 1);
        fac[0] = 1;
        for(int i = 1; i <= n + m ; i++)fac[i] = fac[i - 1] * i % mod;
        auto get = [&](int a , int b){
            if(b > a)return (ll)0 ;
            ll res = fac[a] * ni(fac[b]) % mod * ni(fac[a - b]) % mod;
            return res;
        };
        ll ans = 1 ;

        m = n + m;
        while(n && m && ans){
            ans *= get(m % mod , n % mod);
            ans %= mod;
            n /= mod;
            m /= mod;
        }
        cout<

你可能感兴趣的:(学习,算法)