lucas定理及其拓展的推导

lucas定理及其拓展的推导

我的前一篇博客—— lucas定理 https://mp.csdn.net/mdeditor/100550317#主要是给出了lucas的结论和模板，不涉及推导。
本篇文章主要侧重lucas定理及其拓展的公式推导。先从$p^w$的一般情况进行推导，最后令$w=1$进一步推导出了lucas定理的结论。可能会比较枯燥。
想直接阅读结论的请转向这篇博客——拓展lucas结论及模板 https://blog.csdn.net/yuyilahanbao/article/details/100570075

---

$\left(\genfrac{}{}{0ex}{}{n}{m}\right)\mathrm{m}\mathrm{o}\mathrm{d}{p}^w$的求取

p是素数。主要是计算$\left(\genfrac{}{}{0ex}{}{n}{m}\right)\mathrm{m}\mathrm{o}\mathrm{d}{p}^w$的值。不想看推导过程的可以直接往下翻到下划线以后看结论。

对于$\left(\genfrac{}{}{0ex}{}{n}{m}\right)=\frac{n!}{m!(n-m)!}$.

对于$n!=1\times 2\times 3\times \dots \times n.$
我们将这些数进行分组。假设$k=\lfloor \frac{n}{p}\rfloor ,r=n\mathrm{m}\mathrm{o}\mathrm{d}p;u=\frac{n}{p^w},v=n\mathrm{m}\mathrm{o}\mathrm{d}{p}^w$

1. 把所有p的倍数抽出来，组成一组，它们的乘积记为$A$，易得$A=\prod _{i=1}^k(i\times p)=p^w\times k!$。

2. 剩下的数中,把在这个区间内$[u{p}^w+1,u{p}^w+v]$中的数分为一组，记乘积为$B$.由于每个数都不是$p$的倍数，则$B$与$p^w$互质,$B$关于$p^w$的逆元$B_{p^w}^{-1}$存在。

3. 剩下的每个数a,把$\lfloor \frac{a}{p^w}\rfloor$相同的数分在一组，并把$\lfloor \frac{a}{p^w}\rfloor$作为它们的组号。每一组的乘积分别记作$C_0,C_2,C_3,C_4,\dots ,C_{u-1}$。同理，显然$C_i$与$p^w$互质，逆元存在。另外，每一组的乘积关于$p^w$是同模的，因为每一组的元素都是$\left[p^{w}i+1,p^w(i+1)\right)$的所有元素中去掉了$p$的倍数。

于是$n!=AB\prod _{i=0}^{u-1}{C}_i=p^w\times k!\times B\prod _{i=0}^{u-1}{C}_i$.

回到原来组合数的问题，对于$n,m,n-m$依然有上面$k,r,u,v$的定义，只是分别用下标$1,2,3$区分。

于是$$\left(\genfrac{}{}{0ex}{}{n}{m}\right)=\frac{n!}{m!(n-m)!}=\frac{p^{k_1-k_2-k_3}\times k_1!\times B_1\times \prod _{i=0}^{u_1-1}{C}_i}{k_2!\times k_3!\times B_2\times B_3\times \prod _{i=0}^{u_2-1}{C}_i\times \prod _{i=0}^{u_3-1}{C}_i}$$

注意到$k_1-k_2-k_3\ge 0$,事实上，它只能取0或1.$B$系列，$C$系列都是与$p^w$互质的，而且$\frac{k_1!}{k_2!\times k_3!}$依旧是个整数，所以，可以直接把分母中$B$,$C$系列中的直接变成乘以逆元抹去。

假设$B_i\equiv b_i(\mathrm{m}\mathrm{o}\mathrm{d}{p}^w)$,$C_i\equiv c(\mathrm{m}\mathrm{o}\mathrm{d}{p}^w)$.$b,c\in \left[1,p^w\right)$.
$$\begin{gathered} \frac{p^{k_1-k_2-k_3}\times k_1!\times B_1\times \prod _{i=0}^{u_1-1}{C}_i}{k_2!\times k_3!\times B_2\times B_3\times \prod _{i=0}^{u_2-1}{C}_i\times \prod _{i=0}^{u_3-1}{C}_i} \\ \equiv p^{k_1-k_2-k_3}\times b_1\times {b_2}^{-1}\times {b_3}^{-1}\times c^{u_1}\times (c^{-1}{)}^{u_2}\times (c^{-1}{)}^{u_3}\times \frac{k_1!}{k_2!\times k_3!} \\ \equiv p^{k_1-k_2-k_3}\times b_1\times {b_2}^{-1}\times {b_3}^{-1}\times c^{u_1-u_2-u_3}\times \frac{k_1!}{k_2!\times k_3!}(\mathrm{m}\mathrm{o}\mathrm{d}{p}^w) \end{gathered}$$.

1. 当$r_1\ge r_2$时，$k_1=k_2+k_3$,故上面这个式子最后分式的部分$\frac{k_1!}{k_2!\times k_3!}=\left(\genfrac{}{}{0ex}{}{k_1}{k_2}\right)$.

2. 否则，$k_1=k_2+k_3+1$,最后的那个分式无法直接变成组合数，但是我们只需要分子分母同时乘以$k_1-k_2$，即可变成组合数。$\frac{k_1!}{k_2!\times k_3!}=(k_1-k_2)\times \left(\genfrac{}{}{0ex}{}{k_1}{k_2}\right)$

---

以上最后几行就是结论了，最后再次强调一下各个下标，字母代表的含义。

与n,m,n-m有关的量分别用下标1,2,3区分。

k,r是除以$p$的商与余数，u,v是除以模数$p^w$的商与余数。

b是最后剩下的v个数中不是p的倍数的数的乘积。

c是$\left[1,p^w\right]$中不是p的倍数的数的乘积。

---

从结论中的式子可以看到b,c我们只关注模$p^k$意义下的值，因此可以预先求出$[1..i]$中不是p的倍数的数的乘积$f(i)$(模$p^k$意义下的)。

lucas定理的推导

特殊地，当$w=1$时，$k$和$u$,$r$和$v$是相同的。
故$b\equiv r!(\mathrm{m}\mathrm{o}\mathrm{d}p),c\equiv (p-1)!(\mathrm{m}\mathrm{o}\mathrm{d}p)$，因此代入式子，可以得到lucas定理的结论。

1. 当$r_1\ge r_2$时，$k_1=k_2+k_3$,故上面最后的式子可以化为$ans\equiv r_1!\times (r_2!{)}^{-1}\times (r_3!{)}^{-1}\times \left(\genfrac{}{}{0ex}{}{k_1}{k_2}\right)(\mathrm{m}\mathrm{o}\mathrm{d}p)$

2. 否则，$k_1=k_2+k_3+1$,最后的那个分式无法直接变成组合数，但是我们只需要分子分母同时乘以$k_1-k_2$，即可变成组合数。$\frac{k_1!}{k_2!\times k_3!}=(k_1-k_2)\times \left(\genfrac{}{}{0ex}{}{k_1}{k_2}\right)$。而$ans\equiv p\times r_1!\times (r_2!{)}^{-1}\times (r_3!{)}^{-1}\times c\times (k_1-k_2)\times \left(\genfrac{}{}{0ex}{}{k_1}{k_2}\right)\equiv 0(\mathrm{m}\mathrm{o}\mathrm{d}p)$

$\left(\genfrac{}{}{0ex}{}{n}{m}\right)\mathrm{m}\mathrm{o}\mathrm{d}N$的求取

N是任意正整数。对$N$进行素数分解。$N=\prod _{i=1}^q{p}_i^{k_i}$.
对$\left(\genfrac{}{}{0ex}{}{n}{m}\right)\mathrm{m}\mathrm{o}\mathrm{d}{p}_i^{k_i}$问题，可以通过上一小节的拓展lucas求得，记答案是$c_i$.
于是得到了$q$个线性同余方程,即线性同余方程组$\left(\genfrac{}{}{0ex}{}{n}{m}\right)\equiv c_i(\mathrm{m}\mathrm{o}\mathrm{d}{p}_i^{k_i})(1\le i\le q)$.
对于线性同余方程组，并且注意到模数$p_i^{k_i}$两两互质，可以用中国剩余定理（也可以用拓欧）解出其通解$x=x_0+kt$。并且由于模数互质，$k=lcm(p_i^{k_i})=N(1\le i\le q)$.所以在$[0,N)$内只有一个特解$x_0$,而这个特解就是$\left(\genfrac{}{}{0ex}{}{n}{m}\right)\mathrm{m}\mathrm{o}\mathrm{d}N$.
至此，组合数对任意数的取模我们都予以解决了。