Lucas定理学习记录

序言

晚上的时间真的是不会利用……
和李某学习了很多积极向上的内容……

维基百科可以查到好~多~积极向上的内容……是学习（FuBai）的好工具！
前排膜拜百科
维基百科
百度百科
前排预警：这篇Blog没有任何公式！绝对没有公式！

Lucas Theorem

现给出公式：
${m\choose n}\equiv {m/p\choose n/p}\cdot {m\mod p\choose n\mod p}(\mod p)$
在$p$较小的情况下，可以通过递归在$O(m\log m)$的时间内求出组合数。

证明

设
$\begin{equation} \begin{aligned} m&=\sum\limits_{i=0}^km_ip^i\\ n&=\sum\limits_{i=0}^kn_ip^i\\ \end{aligned} \end{equation}$
在$m,n$化为$p$进制数下时，Lucas公式可变式为：
$\therefore {m \choose n}\equiv \prod\limits_{i=0}^{k}{m_i\choose n_i}(\mod p)$
由费马小定理可得：
$\begin{equation} \begin{aligned} x^p&\equiv x(\mod p)\\ \therefore x^p+1&\equiv x+1(\mod p)\\ \left(x+1\right)^p&\equiv x+1(\mod p)\\ \therefore \left(x+1\right)^p&\equiv x^p+1(\mod p)\\ \therefore \left(x+1\right)^{p^i}&\equiv x^{p^i}+1(\mod p)\\ \end{aligned} \end{equation}$
利用结论推导式子：
$\begin{equation} \begin{aligned} \sum\limits_{n=0}^{m}{m\choose n}x^n&=\left(x+1\right)^m\\ &=\prod\limits_{i=0}^{k}\left(\left(x+1\right)^{p^i}\right)^{m^i}\\ &\equiv \prod\limits_{i=0}^{k}\left(x^{p^i}+1\right)^{m^i}\\ &= \prod\limits_{i=0}^{k}\left(\sum\limits_{n_i=0}^{m_i}{m_i\choose n_i}x^{n_ip^i}\right)\\ &= \prod\limits_{i=0}^{k}\left(\sum\limits_{n_i=0}^{p-1}{m_i\choose n_i}x^{n_ip^i}\right)\\ \sum\limits_{n=0}^{m}{m\choose n}x^n&= \sum\limits_{n=0}^{m}\left(\prod\limits_{i=0}^{k}{m_i\choose n_i}\right)x^n(\mod p)\\ \therefore {m\choose n}&=\prod\limits_{i=0}^{k}{m_i\choose n_i}(\mod p)\\ \end{aligned} \end{equation}$
问题得证。

实际应用

给道例题吧：
P2183 [国家集训队]礼物
（就是怕你不看证明做直接用去做题，我才把例题放到后面。）
快速求大组合数还是很有用的……

局限性

当你想直接用的时候，你需要注意模数 p p <script type="math/tex" id="MathJax-Element-14">p</script>为质数。