扩展Lucas算法详解

原创

于 2025-04-30 22:06:25 发布 · 1.1k 阅读

标签

#算法

收录于

扩展Lucas算法

介绍：

这段内容主要是关于组合数取模的数学知识介绍：

当正整数 $n, m$ 很大，质数 $p$ 较小时，求 $C_{n}^{m}$ 对 $p$ 取模后的值可用Lucas定理；
若 $p$ 不是质数且 $m$ 较小，则可以用扩展Lucas定理,注意扩展Lucas定理与Lucas定理没有一定关心，你只要学会如下知识点就可以：1.质因数分解，2.逆元 ,3.CRT(中国剩余定理) 4.组合数的定义，5.组合数的定义。

第一步：中国剩余定理

设 $p_{1}^{r_{1}} p_{2}^{r_{2}} \cdots p_{k}^{r_{k}}$ ，其中 $p_{i}$ 为质数。我们可以先求出 $C_{n}^{m} \%p_{1}^{r_{1}}$ ， $C_{n}^{m} \%p_{2}^{r_{2}}$ ， $⋯\cdots$ ， $C_{n}^{m} \%p_{k}^{r_{k}}$ 的值 $a1，a2，⋯，aka_{1}，a_{2}，\cdots，a_{k}$ 。
我们把 $C_{n}^{m}$ 看作未知数x，可以得到以下方程组：

$x≡a1(modp1r1)x≡a2(modp2r2)x≡a3(modp3r3)⋯⋯x≡an(modpkrk)\begin{cases} x \equiv a_{1} \pmod {p_{1}^{r_{1}}}\\ x \equiv a_{2} \pmod {p_{2}^{r_{2}}}\\ x \equiv a_{3} \pmod {p_{3}^{r_{3}}}\\ \cdots\cdots\\ x \equiv a_{n} \pmod {p_{k}^{r_{k}}} \end{cases}$

利用中国剩余定理，我们可以求出 $a_i$ ，它是以 $pirip^{r_i}_i$ 为周期出现的无穷多个解。而在 $[0,piri)[0,p^{r_i}_i)$ 这个周期的解，就是 $C_{n}^{m} \ % p$ 后的值。
那么 $a1，a2⋯，aka_{1}，a_{2}\cdots，a_{k}$