标签：frac Lucas cdot 定理 卢卡斯 mod equiv

卢卡斯定理/Lucas 定理

引入

求 \(C_{n+m}^n \mod p\)。

\(n,m,p \leq 10^5\)。

如果直接用阶乘求，可能在阶乘过程中出现了 \(p\)，而最后的结果没有出现 \(p\)，导致错误。

有两种解决方法：

1.求组合数时提前把 \(p\) 的质因子除掉。

2.Lucas 定理。

所以 Lucas 定理用于处理模数较小且为质数的情况下，求组合数的问题。

定理推导

先放结论：

\[C_a^b \equiv C_{\lfloor \frac{a}{p} \rfloor}^{\lfloor \frac{b}{p} \rfloor} \cdot C_{a \mod p}^{b \mod p} \mod p \]

先证明 \(C_p^i \equiv \frac{p}{i}C_{p-1}^{i-1} \equiv 0 \mod p\)。

\[C_p^i=\frac{p!}{i!(p-i)!}=\frac{p}{i} \cdot \frac{(p-1)!}{(i-1)!(p-i)!}=\frac{p}{i}C_{p-1}^{i-1} \]

得证。

考虑二项式定理，易得：

\[(1+x)^p \equiv C_p^0+C_{p}^1x+C_{p}^2x^2+\cdots+C_p^px^p \equiv 1+x^p \mod p \]

Ps：除去 \(C_p^p=1\) 以外，其他的项都被模为 \(0\)。

此时，令 \(a=lp+r,b=sp+j\)。

求证 \(C_a^b \equiv C_l^s\cdot C_r^j \mod p\)。

接着剥削二项式

\[(1+x)^a=(1+x)^{lp}(1+x)^r \]

展开 \((1+x)^{lp}\)

\[(1+x)^{lp} \equiv ((1+x)^p)^l \equiv (1+x^p)^l \mod p \]

\[\therefore (1+x)^a \equiv (1+x^p)^l(1+x)^r \mod p \]

通过上式，观察 \(x^b\) 项。

\[\because C_a^b x^b \equiv C_l^sx^{sp}\cdot C_r^jx^j \mod p \\ \therefore C_a^b x^b \equiv C_l^sC_r^jx^b \mod p \\ \therefore C_a^b \equiv C_{\lfloor \frac{a}{p} \rfloor}^{\lfloor \frac{b}{p} \rfloor} \cdot C_{a \mod p}^{b \mod p} \mod p \]

​

Ps：左边是直接二项式的 \(x^b\) 项，右边是二项式 \((1+x^p)^l\) 的 \(x^{sp}\) 和 \((1+x)^r\) 的 \(x^j\) 项。

标签：frac,Lucas,cdot,定理,卢卡斯,mod,equiv
From： https://www.cnblogs.com/binbinbjl/p/17822632.html