跳转至

阶乘取模

引入

本文讨论了某一模数下阶乘计算的相关结论,并提供一种时间复杂度线性相关于模数大小的计算方法,因而该方法主要适用于模数不太大(106)的情形。除了本文介绍的方法外,根据场景不同,还可以应用 多项式技术 进行快速计算。

根据 中国剩余定理,阶乘取模问题可以转化为模数为素数幂 pα 的情形。在处理这类问题时,常常需要对于素数 p 和正整数 n,将阶乘 n! 中的所有因子 p 都提取出来,进而得到分解:

n!=pνp(n!)(n!)p.

其中,νp(n!) 表示阶乘 n! 的素因数分解中 p 的幂次,(n!)p 表示在阶乘 n! 的结果中去除所有 p 的幂次得到的整数。本文将讨论 (n!)p 在素数(幂)模下的余数以及幂次 νp(n!) 的具体计算方法。

这种分解在解决阶乘同时出现在所求表达式的分子和分母的问题时尤为有用,比如 计算某一模数下的二项式系数。对于这类问题,分子和分母中 p 的幂次可以直接相减,而与 p 互素的部分 (n!)p 则可以利用 乘法逆元 计算。

本文还介绍了与上述问题相关的 Wilson 定理及其推广、Legendre 公式和 Kummer 定理等内容。

Wilson 定理

Wilson 定理给出了判断某个自然数是素数的一个充分必要条件。

Wilson 定理

对于自然数 n>1,当且仅当 n 是素数时,(n1)!1(modn)

证明

首先,证明对于素数 p(p1)!1(modp)。对于这一点,可以利用 同余方程原根 得到两种简洁的证明,此处略去不表。下面提供前置知识较少的一种证明方法:

p=2 时,命题显然成立。下面设 p3,继而要证明 Zp 中所有非零元素(即同余类)的积为 1。因为 Zp 中所有非零元素 a 都有逆元 a1,于是 Zp 中彼此互逆的元素乘积为 1。但是要注意 aa1 可能相等:a=a1,当且仅当 a21(modp),即

0a21(a+1)(a1),(modp)

从而,a1(modp)a1(modp)。这说明 Zp{0,1,1} 中所有元素的乘积为 1,进而 Zp 中所有非零元素的积为 1

反过来,对于合数 n 的情形,要证明 (n1)!1(modn)。利用反证法,不妨设 (n1)!1(modn),亦即存在整数 k 使得 (n1)!=kn1 成立。因为 n 是合数,必然存在素数 p<n 使得 n=pm,所以 (n1)!=kpm11(modp)。但是,乘积 (n1)! 中必然已经出现 p,故而一定有 (n1)!0(modp)。这一矛盾就说明了 (n1)!1(modn)

利用本文的记号,Wilson 定理可以写作 (p!)p1(modp)

推广

Wilson 定理可以推广到一般模数的情形。

定理(Gauss)

对于自然数 m>1,有

1k<m, kmk±1(modm).

而且,余数中的 ±1 取值为 1 当且仅当模 m原根存在,即 m=2,4,pα,2pα 时,其中 p 是奇素数且 α 是正整数。

证明

这个定理可以通过 n 整数乘法群 的结构简单地证明。此处给出思路相仿,但是较为初等的证明。

对于 m=2 的情形,有 1!=11(mod2)。对于其他存在原根的情形,设原根为 g,则所有满足小于 m 且与它互素的正整数 k 都可以唯一地表示为 gimodm 的形式,其中 0i<φ(m)φ(m)Euler 函数。直接验证可知,φ(m) 一定是偶数。因为 gigφ(m)i 互为乘法逆元,所以在乘积中将它们两两配对,就有

1k<m, kmki=0φ(m)1gi=gφ(m)/2i=1φ(m)/21gigφ(m)igφ(m)/2(modm).

因为 gφ(m)/2modm 是唯一的不等于 1modm 且乘法逆元就是它自身的元素,所以它就等于 1modm。这就说明了此时的余数等于 1

对于模 m 的原根不存在的情形,要证明余数等于 1。为此,可以首先做质因数分解 m=p1e1p2e2pses,然后应用 中国剩余定理 可知,只需要证明

1k<m, kmk1(modpjej)

对所有因子 pjej 都成立。中国剩余定理说明,每一个可能的余数组合 (r1,r2,,rs),其中,1rj<pjejpjrj,都唯一地对应着一个 1k<mkm 使得 krj(modpjej) 成立。所以,对于某个余数 rj,都恰好有 φ(m)/φ(pjej)k 使得 krj(modpjej) 成立。利用这一点,可以对乘积进行分组,就有

1k<m, kmk(1rj<pjej, rjpjrj)φ(m)/φ(pjej)(modpjej).

此处的指数 φ(m)/φ(pjej)=φ(m/pjej) 要成为奇数,必然要求 m/pjej=1,2,因为欧拉函数 φ(n) 对于 n3 都是偶数。如果 pj 是奇素数,因为模 m 的原根不存在,必然有 m/pjej1,2;如果 pjej=2,4,因为模 m 的原根不存在,必然有 m/pjej 含有某个奇素因子,故而大于 2:这两种情形指数 φ(m)/φ(pjej) 都是偶数。而上式中括号里的项已经证明是模 pjej1 的,所以这个幂模 pjej 的余数一定是 1。剩余的情形只有 pj=2ej>2 时,对于这个情形,可以直接证明 `

1rj<2ej, rj2rj1(mod2ej).

仿照前文的证明思路,可以将所有 1rj<2ej 的奇数 rj 两两配对而消去,那些无法配对的必然是方程 x21(mod2ej) 的解。该方程意味着 2ej(x1)(x+1)。令 x=2y+1,就必然有 2ej2y(y+1),而 yy+1 必然一奇一偶,所以 y=t2ej2y=t2ej21。故而,有 x=t2ej1±1t 是整数。模 2ej 的余数中,只有 ±12ej1±1 四个。因此,有

1rj<2ej, rj2rj(1)(2ej11)(2ej1+1)1(mod2ej).

这就完成了所有情形的证明。

在计算中,尤为重要的是模数为素数幂的情形:

推论

对于素数 p 和正整数 α,有

$$ \prod_{1\le k