卢卡斯定理

Dejavu___

已于 2023-02-24 16:35:39 修改

阅读量150

点赞数 1

分类专栏：算法中用到的数学文章标签：算法 c++

于 2023-02-22 21:21:47 首次发布

本文链接：https://blog.csdn.net/Dejavu___/article/details/129170032

版权

算法中用到的数学专栏收录该内容

2 篇文章 0 订阅

订阅专栏

文章目录

先上结论
应用场景
证明：

先上结论

求组合数时：
$C_{n}^{m}\equiv C_{n\quad mod\quad p}^{m\quad mod \quad p}*C_{n / p}^{m/p}(mod\quad p)$

应用场景

卢卡斯定理一般用于大组合数取模问题，当m，n很大且可能为p的倍数时，快速幂不适用，且递推法复杂度过大，这时应用卢卡斯定理就可以给运算带来极大的方便。

证明：

首先证明：
$1+x)^p=1+x^p$
按照二项式定理展开: $(1+x)^p\equiv C_{p}^{0}1+C_{p}^{1}x^1+C_{p}^{2}x^2+……C_{p}^{p}x^p\equiv1+x^p(mod\quad p)$
因为除 $C_{p}^{0}和C_{p}^{p}外其余$ 系数 $C_{p}^{1},C_{p}^{2},C_{p}^{3}……,C_{p}^{p-1}$ 都被p除掉

接着：
令n=ap+b,m=cp+d;
对于 $(1+x)^n\equiv C_{n}^{0}1+C_{n}^{1}x^1……+C_{n}^{n}x^n$ ……(1)
又 $(1+x)^n\equiv (1+x)^{a p+b}\equiv((1+x)^a)^p (1+x)^b$
$\equiv (1+x^p)^a(1+x)^b$

$\equiv \sum_0^aC_{a}^{i}x^{pi}\times \sum_0^aC_{b}^{i}x^i(mod\quad p)$ ……（2）

在(1)中 $X^m$ 的系数为 $C_{n}^{m}$
在(2)中 $X^m=X^{cp}X^d$ 的系数为 $C_{a}^{c}C_{b}^{d}$
故 $C_{n}^{m}\equiv C_{a}^{c}C_{b}^{d}(mod\quad p)$
即 $C_{n}^{m}\equiv C_{n\quad mod\quad p}^{m\quad mod \quad p}*C_{n / p}^{m/p}(mod\quad p)$ ，证毕.