卢卡斯定理

卢卡斯定理（Lucas定理）

链接中视频给出的代码是错误的，因为有可能出现$nmodp<mmodp$的情况，getC函数可能会越界。尽管这份代码不开O2确实是可以通过的，但是开了O2会无法通过。

卢卡斯定理解决大组合数模质数的问题。

卢卡斯定理指的是：
$\underset{i=0}{\stackrel{n}{\mathbb{I}}}\left(\begin{array}{c}{\overline{a_i}}_{(p)}\\ {\underleftarrow{\overline{b_i}}}_{(p)}\end{array}\right)\equiv \stackrel{n}{\prod _{i=0}}\left(\begin{array}{c}{a}_i\\ b_i\end{array}\right)(\mathrm{mod}p)(p\in \mathbb{P})$
其中$\mathbb{I}$表示循环记号，${\overline{a_n{a}_{n-1}...}}_{(p)}$表示将$a_n,a_{n-1},...$依照次序连成一个$p$进制数。

卢卡斯定理是一个把组合数分解$p$进制的过程，它表明，分解前后组合数模意义下相等。

卢卡斯定理还有一个常用版本：
$\left(\begin{array}{c}n\\ m\end{array}\right)\equiv \left(\begin{array}{c}n/p\\ m/p\end{array}\right)+\left(\begin{array}{c}nmodp\\ mmodp\end{array}\right)(\mathrm{mod}p)$
其中$/$表示整除。

证明一下：

（奥秘重重的证明过程）

设$p\in \mathbb{P}$

引理1

${(}_x^p)\equiv 0(\mathrm{mod}p)(0<x<p)$

这是因为，显然有${(}_x^p)\equiv \frac{p!}{x!(p-x)!}\equiv \frac{p}{x}\cdot \frac{(p-1)!}{(x-1)!((p-1)-(x-1))!}\equiv p\cdot x^{-1}\cdot {(}_{x-1}^{p-1})\equiv 0\cdot x^{-1}\cdot {(}_{x-1}^{p-1})\equiv 0$

在提出$p$的过程中，必须提出$x$，得到一个更小的组合数的原因是，必须要保证提出$p$之后仍然是整数。

引理2

$(1+x{)}^p\equiv 1+x^p(\mathrm{mod}p)$

结合引理1，做二项式展开易证。

结论

设$n=k_{1}p+r_1,m=k_{2}p+r_2(0\le r_1,r_2<p)$，则有：
$(1+x{)}^n\equiv \stackrel{n}{\sum _{m=0}}{(}_m^n)x^m\equiv \stackrel{k_{1}p+r_1}{\sum _{m=0}}\left({}_{k_{2}p+r_2}^{k_{1}p+r_1}\right)x^{k_{2}p+r}$

另有：
$(1+x{)}^n\equiv (1+x{)}^{k_{1}p}(1+x{)}^{r_1}$
$\equiv {\left((1+x{)}^p\right)}^{k_1}(1+x{)}^{r_1}$

应用引理2：
$\equiv (1+x^p{)}^{k_1}(1+x{)}^{r_1}$
$\equiv \stackrel{k_1}{\sum _{i=0}}\left({}_i^{k_1}\right)x^{i\cdot p}\times \stackrel{r_1}{\sum _{j=0}}\left({}_j^{r_1}\right)x^j$

对比系数：
$\equiv \stackrel{k_{1}p+r_1}{\sum _{m=0}}\left({}_{k_{2}p+r_2}^{k_{1}p+r_1}\right)x^{k_{2}p+r}$
$\equiv \stackrel{k_1}{\sum _{i=0}}\left({}_i^{k_1}\right)x^{i\cdot p}\times \stackrel{r_1}{\sum _{j=0}}\left({}_j^{r_1}\right)x^j$

可以看出，显然，模$p$下可以把一个数分解：
$\left(\begin{array}{c}n\\ m\end{array}\right)\equiv \left(\begin{array}{c}n/p\\ m/p\end{array}\right)+\left(\begin{array}{c}nmodp\\ mmodp\end{array}\right)(\mathrm{mod}p)$

QED.

时间复杂度O(logn)

组合数模合数较繁。

后记

于是皆大欢喜。