Lucas 定理学习笔记

Lucas 定理

Lucas 定理用于求解大组合数取模的问题，其中模数必须为素数。正常的组合数运算可以通过递推公式求解（详见排列组合），但当问题规模很大，而模数是一个不大的质数的时候，就不能简单地通过递推求解来得到答案，需要用到 Lucas 定理。

求解方式

Lucas 定理内容如下：对于质数 $p$，有

$$\left(\genfrac{}{}{0ex}{}{n}{m}\right)\mathrm{m}\mathrm{o}\mathrm{d}p=\left(\genfrac{}{}{0ex}{}{\lfloor n/p\rfloor }{\lfloor m/p\rfloor }\right)\cdot \left(\genfrac{}{}{0ex}{}{n\mathrm{m}\mathrm{o}\mathrm{d}p}{m\mathrm{m}\mathrm{o}\mathrm{d}p}\right)\mathrm{m}\mathrm{o}\mathrm{d}p$$

观察上述表达式，可知 $n\mathrm{m}\mathrm{o}\mathrm{d}p$ 和 $m\mathrm{m}\mathrm{o}\mathrm{d}p$ 一定是小于 $p$ 的数，可以直接求解，$\left(\genfrac{}{}{0ex}{}{\lfloor n/p\rfloor }{\lfloor m/p\rfloor }\right)$ 可以继续用 Lucas 定理求解。这也就要求 $p$ 的范围不能够太大，一般在 $10^5$ 左右。边界条件：当 $m=0$ 的时候，返回 $1$。

时间复杂度为 $O(f(p)+g(n)\log n)$，其中 $f(n)$ 为预处理组合数的复杂度，$g(n)$ 为单次求组合数的复杂度。

ll qpow(ll x, ll n, ll mod) {
   
	ll ans = 1;
	while (n) {
   
		if (n & 1)ans = (ans * x) % mod;
		x = (x * x) % mod;
		n >>= 1;
	}
	return ans % mod;
}
ll C(ll a, ll b, ll p){
   //通过定理求组合数C(a, b)
    if (a < b) return 0;
    ll x = 1, y = 1;// x是分子，y是分母
    for (int i = a, j = 1; j <= b; i--, j++){
   
        x = x * i % p;
        y = y * j % p;
    }
    return x * qpow(y, p - 2, p) % p;
}
ll lucas(ll a, ll b, int p){
   //lucas定理求组合数
    if (a < p && b < p) return C(a, b, p);
    return C(a % p, b % p, p) * lucas(a / p, b / p, p) % p;
}

素数在阶乘中的幂次

Legengre 在 1808 年指出 $n!$ 中含有的素数 $p$ 的幂次为 ${\sum }_{j\ge 1}\lfloor n/p^j\rfloor$。

证明：将 $n!$ 记为 $1\times 2\times \cdots \times p\times \cdots \times 2p\times \cdots \times \lfloor n/p\rfloor p\times \cdots \times n$ 那么其中 $p$ 的倍数有 $p\times 2p\times \cdots \times \lfloor n/p\rfloor p=p^{\lfloor n/p\rfloor }\lfloor n/p\rfloor !$ 然后在 $\lfloor n/p\rfloor !$ 中继续寻找 $p$ 的倍数即可，这是一个递归的过程。为了方便记 $\nu (n!)={\sum }_{j\ge 1}\lfloor n/p^j\rfloor$。

另一种其他地方比较常见的公式，用到了 p 进制下各位数字和：

$v_p(n!)=\frac{n-S_p(n)}{p-1}$

与等比数列求和公式很相似。由于涉及各位数字和，利用数学归纳法可以轻松证明。

特别地，阶乘中 2 的幂次是：

$v_2(n!)=n-S_2(n)$

素数在组合数中的幂次

组合数对一个数取模的结果，往往构成分形结构，例如谢尔宾斯基三角形就可以通过组合数模 2 得到。

$v_p(C_m^n)=\frac{S_p(n)+S_p(m-n)-S_p(m)}{p-1}$

如果仔细分析，p 是否整除组合数其实和上下标在 p 进制下减法是否需要借位有关。这就有了 Kummer 定理。

Kummer 定理：p 在组合数 $C_m^n$ 中的幂次，恰好是 p 进制下 m 减掉 n 需要借位的次数。

特别地，组合数中 2 的幂次是：

$v_2(C_m^n)=S_2(n)+S_2(m-n)-S_2(m)$

Wilson 定理

对于素数 $p$ 有 $(p-1)!\equiv -1(\mathrm{m}\mathrm{o}\mathrm{d}p)$。

• 证明：我们知道在模奇素数 $p$ 意义下，$1,2,\dots ,p-1$ 都存在逆元且唯一，那么只需要将一个数与其逆元配对发现其乘积均为（同余意义下）$1$，但前提是这个数的逆元不等于自身。那么很显然 $(p-1)!\mathrm{m}\mathrm{o}\mathrm{d}p$ 就是逆元等于其自身的数的乘积，这两个数为 $\pm 1$。在 $p$ 为 $2$ 时单独讨论即可。对于整数 $n$，令 $(n!{)}_p$ 表示所有小于等于 $n$ 但不能被 $p$ 整除的正整数的乘积，即 $(n!{)}_p=n!/(\lfloor n/p\rfloor !p^{\lfloor n/p\rfloor })$。
• Wilson 定理指出 $(p!{)}_p=(p-1)!\equiv -1(\mathrm{m}\mathrm{o}\mathrm{d}p)$ 且可被推广至模素数 $p$ 的幂次。

推论 1

对于素数 $p$ 和正整数 $q$ 有 $(p^q!{)}_p\equiv \pm 1(\mathrm{m}\mathrm{o}\mathrm{d}{p}^q)$。

依然考虑配对一个数与其逆元，也就是考虑关于 $m$ 的同余方程 $m^2\equiv 1($