大组合数Lucas

https://blog.csdn.net/sr_19930829/article/details/39058487

 LL Lucas(LL n, LL m, int p){
        return m ? Lucas(n/p, m/p, p) * comb(n%p, m%p, p) % p : ;
 }

Saving Beans HDU3037

#include <iostream>
#include <string.h>
#include <cmath>
using namespace std;
typedef long long ll;
const int maxm = 1e5 + ;
ll fac[maxm];

void init(ll p) {
    fac[]=;
    for(int i = ; i <= p; i++)
        fac[i] = fac[i - ] * i % p;
}

ll qpow(ll x, ll n, ll mod) {
    ll res = ;
    while (n > ) {
        if (n & ) {
            res = res * x % mod;
        }
        x = x * x % mod;
        n >>= ;
    }
    return res;
}

ll lucas(ll n, ll m, ll p)
{
    ll res = ;
    while(n && m)
    {
        ll a = n % p;
        ll b = m % p;
        if(a < b) return ;
        res = (res * fac[a] * qpow(fac[b] * fac[a - b] % p, p - , p)) % p;
        n /= p;
        m /= p;
    }
    return res;
}

int t;
ll n, m, p;
int main() {
    scanf("%d", &t);
    while(t--) {
        scanf("%lld%lld%lld", &n, &m, &p);
        init(p);
        printf("%lld\n", lucas(n + m, m, p));
    }
    return ;
}