【总结】lucas 定理 + 中国剩余定理

先上学习资料：
Lucas’ Theorem
中国剩余定理 Chinese Remainder Theorem
组合数取模 - ACdreamer

lucas 定理用来计算组合数模素数。如果素数P可以先确定，则可以$O(P)$预处理，每次计算时间复杂度为 $O(logPlogN)$。不预处理的时间复杂度，$O(PlogN)$。

中国剩余定理用来解模方程组

其中 m1, m2 … mn 任意两数互质

结合两个定理，就可以求解组合数模任意数。通过将模分解成若干质数，用 Lucas定理 分别求解，再用 CRT 合并。

下面是一些习题

hdu 3944 DP?
先可以猜想两种路径，从终点斜向上走到边界，或竖直向上走到边界，则得到两种答案。
1) $\binom{n+1}{m} + n - m$
2) $\binom{n+1}{m+1} + m$
两者相减可以计算出，$2 * m < n$ 时，2）比较优，否则选 1）。

const int N = 1e4;
const int PN = 1300;
int *fc, Mod;
int fac[PN][N + 5]; // 阶乘表
int pnm, pi[PN]; // 素数总数，素数表

inline LL C(LL n, LL m) {
    if ( m > n ) return 0;
    if ( m == n || m == 0 ) return 1;
    return fc[n] * inv( fc[m] * fc[n-m] % Mod) % Mod;
}

LL lucas(LL n, LL m) {
    if ( m == 0 ) return 1;
    if ( n == 0 ) return 0;
    return C(n % Mod, m % Mod) * lucas(n / Mod, m / Mod) % Mod;
}

void init() {
    sieve(1e4);
    for (int i = 0; i < pnm; ++ i) fac[i][0] = 1;
    for (int i = 1; i <= N; ++ i) {
        for(int j = 0; j < pnm; ++ j)
            fac[j][i] = fac[j][i-1] * i % pi[j];
    }
}

int main() {
    init();
    int n, k, p, cas = 0;
    while ( scanf("%d%d%d", &n, &k, &p) != EOF ) {
        int i;
        for (i = 0; i < pnm; ++ i) if ( pi[i] == p ) break;
        fc = fac[i]; Mod = p;
        int ans;
        if ( k * 2 < n )
            ans = ( n - k + lucas( n + 1, k ) ) % p;
        else
            ans = ( k + lucas( n + 1, k + 1 ) ) % p;
        printf("Case #%d: %d\n", ++ cas, ans);
    }
    return 0;
}

hdu 5446 Unknown Treasure
lucas 结合 CRT 模板题，特点是 n, m, M(为分解的模)都达到了1e18
所以，在 CRT 合并的时候要采取避免 long long 溢出的姿势。

// x = ai (mod mi)
// M = m1 * m2 * ... * mk, Mi = M / mi
// x = sum { ai * inv(Mi) * Mi } + nM, 通解
LL RemindChina( const LL a[], const LL m[], int k ) {
    LL M = 1;
    for (int i = 0; i < k; ++ i) M *= m[i];
    LL ans = 0;
    for (int i = 0; i < k; ++ i) {
        LL Mi = M / m[i];
        Mod = m[i];
        LL _inv = inv(Mi);
        ans = ( ans + mul( mul(a[i], _inv, M), Mi, M ) ) % M;
    }
    return ans;
}

int main() {
    int t;
    LL n, m;
    cin >> t;
    while ( t -- ) {
        cin >> n >> m >> pnm;
        rep(i, 0, pnm-1) cin >> pi[i];
        rep(i, 0, pnm-1) fac[i][0] = 1;
        for (int i = 1; i <= N; ++ i) {
            for (int j = 0; j < pnm; ++ j) {
                fac[j][i] = fac[j][i-1] * i % pi[j];
            }
        }
        LL ai[16];
        rep(i, 0, pnm-1) {
            fc = fac[i];
            Mod = pi[i];
            ai[i] = lucas(n, m);
        }
        LL ans = RemindChina(ai, pi, pnm);
        cout << ans << endl;
    }
    return 0;
}