exLucas - 拓展Lucas

对普通Lucas不熟悉的可以看看这篇博客

拓展Lucas

在使用普通Lucas时，要求模数为质数，因为要用到模数$p$的逆元。拓展Lucas用于处理$p$不为质数的情况。

具体操作如下：

首先，我们对$p$进行质因数分解得到：$p={p_1}^{x1}\cdot {p_2}^{x2}\cdot ...\cdot {p_m}^{xm}$
我们欲求$C(n,m)modp=n!/m!/(n-m)!modp$，如果能分别求得$n!/m!/(n-m)!mod{p_i}^{xi}$，就能用$CRT$求得答案。那么我们将问题转化为求$n!/m!/(n-m)!$ 对${p_i}^{xi}$，下面简称为$p$和$x$。

求$n!/m!/(n-m)!mod{p}^x$，我们只需要会求$n!mod{p}^x$，另外两个阶乘分别求得。设 $i!mod{p}^x$ 为$F(i)$，可得到答案为 $F(n)\ast Inv(F(m))\ast Inv(F(n-m))$。下面我们考虑求解$F(n)$

$n!=1\ast 2\ast ...\ast n$，将其根据 $p$为长度分割得到 $n!=1\ast 2\ast ...\ast p\ast (p+1)\ast (p+2)\ast ...\ast 2p\ast (2p+1)\ast ...\ast n$。假设提出 $p$ 到 $2p-1$ 这一段，将其写作：$p\ast (p+1)\ast (p+2）\ast ...\ast (p+p-1)$，那么对于$p^t(0<=t<p)$，该项前的系数为$a_1\ast a_2\ast ...\ast a_{p-1-t}$。这里的$1-(p-1-t)$并不一定对应连续，因此我们要求一个$\sum$。设 $dp[i][j]$为前 $i$ 个 $a$中选 $j$ 个的和。易得$dp[i][j]=dp[i-1][j]+dp[i-1][j-1]+a[i]$，那么 $dp[t][p-t-1]$就为系数。对于$p^t$，$t>x$的项，我们可以直接模去，反之可以线性求解。而对于$2p-(3p-1)$，推到$kp-(k+1)p-1$，同次数$p^t$前的系数就为之前的$p-(2p-1)$的系数前乘上一个$k^t$。由此，我们可以求得一些段的答案。那么这么些“段”的定义是什么呢？我们发现，以$p^x$作为循环节分割，每个循环节的答案在模$p^x$下是等价的。如此，我们就学会了如何求$n!mod{p}^x$。

综上，我们就可以求$C(n,m)modp$了。

话是这么说，代码是真难打啊……