快速求排列组合 lucas定理

对于C(n, m) mod p。这里的n，m，p（p为素数）都很大的情况。

就不能再用C(n, m) = C(n - 1，m) + C(n - 1, m - 1)的公式递推了。

一般lucas定理的p不能大，在1e6以内，一下代码应该可以吧

typedef long long LL;
using namespace std;

LL exp_mod(LL a, LL b, LL p) {
    LL res = ;
    while(b != ) {
        if(b&) res = (res * a) % p;
        a = (a*a) % p;
        b >>= ;
    }
    return res;
}

LL Comb(LL a, LL b, LL p) {
    if(a < b)   return ;
    if(a == b)  return ;
    if(b > a - b)   b = a - b;

    LL ans = , ca = , cb = ;
    for(LL i = ; i < b; ++i) {
        ca = (ca * (a - i))%p;
        cb = (cb * (b - i))%p;
    }
    ans = (ca*exp_mod(cb, p - , p)) % p;
    return ans;
}

LL Lucas(int n, int m, int p) {
     LL ans = ;

     while(n&&m&&ans) {
        ans = (ans*Comb(n%p, m%p, p)) % p;
        n /= p;
        m /= p;
     }
     return ans;
}

上面是在线的算法，如果数字小还可以用以下写法

void pre()  //i表示n，j表示m
{
    C[][]=C[][]=;
    for(int i=;i<=;i++)
    {
        C[i][]=;
        for(int j=;j<=;j++)
          C[i][j]=(C[i-][j]+C[i-][j-])%mod;
    }
    return;
}