Lucas定理

最新推荐文章于 2024-08-01 16:31:59 发布

愚末语

最新推荐文章于 2024-08-01 16:31:59 发布

阅读量106

点赞数

分类专栏：学习文章标签：算法

本文链接：https://blog.csdn.net/weixin_50547586/article/details/127701366

版权

学习专栏收录该内容

46 篇文章 0 订阅

订阅专栏

总结：

$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<bits/stdc++.h>
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;
        vector<ll>fac(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<<ans<<'\n';

    }

    return 0 ;
}