lucas定理、拓展lucas定理学习小结

lucas定理

正题

首先，这玩意就是下面这个式子：
$$C_m^n\%p=C_{m/p}^{n/p}\ast C_{m\%p}^{n\%p}\%p$$
当且仅当$p$为质数。

当然，还有一个形式：
首先把n写成一个p进制下的东东：$n=\sum a_i\ast p^i$
在把m写成一个p进制下的东东：$m=\sum b_j\ast p^j$
然后
$$C_m^n=\Pi C_{b_i}^{a_i}\%p$$

证明，由于则不是主要内容，所以在此不展开了。自己上百度。
或者这里有一张图：（转的）

其实是太懒了

拓展lucas定理

前言

这个东东是解决$C_m^n\%p$其中p不为质数的情况的。
其实和上面的lucas做法没有半点关系。

前置姿势

中国剩余定理、拓展gcd
没了。

中国剩余定理：（百度百科）

exgcd：（百度百科）

正文

直接求解$C_m^n\%P$我们需要两个条件，一个是能够快速求出$n！$，一个是能快速求出逆元。
然而当n和m都奇大无比时， 且p不为质数。那么就很头痛了。

然后拓展lucas定理（其实都可以不必叫定理了）就可以比较好地解决这个问题。
而这也不是万能的，lucas定理可以运用的前提是p这个合数分解质因数不能是一个巨大的质数乘上一个巨小的质数。

当然，这个可以当做一个判断是否用拓展lucas的一个条件。

讲这么多废话，那么现在开始。

首先，我们要求的东东长这样：
$$C_m^n\%P=\frac{m!}{n!(m-n)!}\%P$$
首先，我们考虑用拓展中国剩余定理来简化一下这个东东：
$$\{\begin{array}{c}x\equiv C_m^n(\%p_1^{q_1})\\ x\equiv C_m^n(\%p_2^{q_2})\\ x\equiv C_m^n(\%p_3^{q_3})\\ \dots \dots \end{array}$$

如果我们求出右式，那么意味着我们可以用中国剩余定理直接合并。

那么问题则转化为求解：
$$C_m^n\%p^q$$
那么我们现在的问题就有两个，一个是如何快速求出阶乘。其次是如何求出这个阶乘的逆元。

求逆元

就以求$n!$为例
由于$n!$中有许许多多的因子为$p$的数。
那么这就不与模数互质了，于是就不能直接套exgcd来求逆元。

怎么办？

考虑从$n!$把这些$p$的因子给提出来。
我们原式就变为：$\frac{n!}{p^x}\ast p^x$。其中$x$是一个未知量，且我们这个x要尽量大使得$\frac{n!}{p^x}$里没有$p$的因子。
我们可以知道，一共有$\lfloor \frac{n}{p}\rfloor$个这些p的因子，然后他们组成情况为：$${\Pi }_{i=1}^{\lfloor \frac{n}{p}\rfloor }i\ast p$$
化简得：$$(\lfloor \frac{n}{p}\rfloor )!\ast p^{\lfloor \frac{n}{p}\rfloor }$$
那么$n!$就被我们化成：
$$(\lfloor \frac{n}{p}\rfloor )!\ast p^{\lfloor \frac{n}{p}\rfloor }\ast {\Pi }_{i=1}^n[i\not\equiv 0\%p^q]i$$
然鹅$(\lfloor \frac{n}{p}\rfloor )!$还是有可能有$p$的因子，因此我们考虑设个$f(n)=\frac{n!}{p^x}$
于是原式可以通过递归的方法求解。其中我们的$p^{\lfloor \frac{n}{p}\rfloor }$可以在外面被$p^x$给化简掉，所以递归过程忽略。
那么递推式变为：

$$f(n)=f(\lfloor \frac{n}{p}\rfloor )\ast {\Pi }_{i=1}^n[i\not\equiv 0\%p^q]i$$

这样我们就可以求出$\frac{n!}{p^x}$
然后外面的$p^x$我们一样考虑通过递归来求出x的值。
设$g(n)=x$
那么递归过程中，每次消去$\lfloor \frac{n}{p}\rfloor$，且阶乘一样产生贡献。
因此递归式为：
$$g(n)=\lfloor \frac{n}{p}\rfloor +g(\lfloor \frac{n}{p}\rfloor )$$

至此，我们已经解决了求逆元的问题。

求阶乘
我们发现，上面的递推式中：${\Pi }_{i=1}^n[i\not\equiv 0\%p^q]i$长得不是很好看。
但是我们发现，由于是在模$p^q$意义下的东西，所以这玩意是有周期的。
怎么理解？
比如：$1、2、4、5、7、8、10、11、13、14、16、17$其中$p=3,q=2$
（去除了有因子3的数）
我们发现，$1\ast 2\ast 4\ast 5\ast 7\ast 8$和$10\ast 11\ast 13\ast 14\ast 16\ast 17$在模意义下是相同的。
那么我们就可以只算出一个周期的值，然后再把这些周期值快速幂乘起来。
然后可能剩下一小段，那么由于比较小，所以直接暴力乘起来即可。

$${\Pi }_{i=1}^n[i\not\equiv 0\%p^q]i=({\Pi }_{i=1}^{p^q}[i\not\equiv 0\%p^q]i{)}^{\lfloor \frac{n}{p}\rfloor }\ast ({\Pi }_{i=p^q\ast \lfloor \frac{n}{p}\rfloor }^n[i\not\equiv 0\%p^q]i)$$

代码

（看例题）

例题

luoguP4720 【模板】扩展卢卡斯
jzoj3214. 【SDOI2013】方程

参考资料

https://blog.csdn.net/qq_30115697/article/details/88942177
https://www.luogu.com.cn/problem/P4720 中的fading题解
各种百度百科