卢卡斯(Lucas)定理

最新推荐文章于 2025-06-25 13:42:59 发布

jvruo_shabi

最新推荐文章于 2025-06-25 13:42:59 发布

阅读量1.3k

点赞数 5

CC 4.0 BY-SA版权

分类专栏：排列组合文章标签：数学

本文链接：https://blog.csdn.net/jvruo_shabi/article/details/109219110

这篇博客以易懂的方式介绍了Lucas定理，主要用于大组合数取膜，特别是当质数p很小的情况。文章指出Lucas定理可以通过递归实现，并引入了逆元的概念来处理模运算中的除法问题。还提到了费马小定理在求逆元中的应用，以及快速幂算法在计算中的作用。最后给出了组合数计算的代码实现，但警告Lucas定理仅适用于模数为质数且不超过1e6的情况。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

看了看网上关于 $L u c a s$ 的讲解，感觉很多充斥着大量难懂的公式，对于~~小学生~~初学者来说不太友好，于是决定写一篇没有任何算法数学基础的人都能看懂的 $L u c a s$ 。（当然本文中也会出现一些公式，不过没有 $∑\sum$ 这种令人难受的符号而且公式大多特别简单）

$L u c a s$ 定理主要用于大组合数取膜。

他的结论很简单：当且仅当 $p$ 为质数时， $pC_m^n\, mod\,p = C_{m\,mod\,p}^{n\, mod\, p} \times C_{m\div p}^{n\div p}\,mod\,p$

~~由于作者太弱并不会证明~~，如果对恶心的证明感兴趣可以左转百度。

在 $n, m$ 很大而 $p$ 又很小的时候就可以使用 $L u c a s$ 定理。

写个递归就可以实现 $L u c a s$ 定理啦（本文所有代码均默认已经添加了typedef long long ll）！

inline ll Lucas(ll m, ll n) {
   
   
	return n == 0 ? 1 : C(m % p, n % p) * Lucas(m / p, n / p) % p;
}

慢着，你以为这就完了？

$L u c a s$ 函数已经很优秀了，可是组合数如果按照杨辉三角计算…… $L u c a s$ 复杂度+n。

按公式计算呢？炸了。

这里科普一个叫逆元的东西，它是拿来干啥的呢？

我们知道：

$p(a+b)\,mod\,p=(a\,mod\,p+b\,mod\,p)\,mod\,p$
$p(a+b)\,mod\,p=(a\,mod\,p-b\,mod\,p + p)\,mod\,p$
$p(a\times b)\,mod\,p=(a\,mod\,p\times b\,mod\,p)\,mod\,p$