威尔逊定理学习笔记

一、威尔逊定理内容及证明

内容

$p$ 质数当且仅当 $(p-1)!\equiv -1(\mathrm{m}\mathrm{o}\mathrm{d}p)$ 。

证明

（1）充分性

即 $p$ 是质数，则 $(p-1)!\equiv -1(\mathrm{m}\mathrm{o}\mathrm{d}p)$ 。

我们知道，在模奇素数 $p$ 的剩余系下，$1$ 到 $p-1$ 的逆元都存在。考虑到 $x{x}^{-1}\equiv 1(\mathrm{m}\mathrm{o}\mathrm{d}p)$ ，原式可化简为：
$$(p-1)!={\Pi }_{i^{-1}\ne i}i\times {\Pi }_{i^{-1}=i}i\equiv {\Pi }_{i^{-1}=i}i(\mathrm{m}\mathrm{o}\mathrm{d}p)$$
因此需要考虑哪些数的逆元是其本身，即满足二次剩余 $x^2\equiv 1(\mathrm{m}\mathrm{o}\mathrm{d}p)$ 。移向后因式分解可得
$$(x-1)(x+1)\equiv 0(\mathrm{m}\mathrm{o}\mathrm{d}p)$$

解得 $x=\pm 1$ 。

在 $1$ 到 $p-1$ 这 $p-1$ 个数中，除了 $1$ 和 $p-1$ 的逆元是本身外，对于 $\forall x\in [2,p-2],\exists y\in [2,p-2]$ 满足 $xy\equiv 1(\mathrm{m}\mathrm{o}\mathrm{d}p)$ 且 $y\ne x$ 。

综上，${\Pi }_{i^{-1}=i}i=1\times (p-1)\equiv -1(\mathrm{m}\mathrm{o}\mathrm{d}p)$ 。

（2）必要性

即 $(p-1)!\equiv -1(\mathrm{m}\mathrm{o}\mathrm{d}p)$，则 $p$ 是质数。直接证明有些困难，因此采用反证法，假设 $p$ 满足条件且为合数，推出矛盾“任何合数都不可能满足这个条件”。

① $p=1$

此时 $(p-1)!=0!=1\equiv 0(\mathrm{m}\mathrm{o}\mathrm{d}1)$ ，不满足上述条件。

② $p=4$

此时 $(p-1)!=3!=6\equiv 2(\mathrm{m}\mathrm{o}\mathrm{d}4)$ ，不满足上述条件。

③ $p>4$ 且 $p$ 是完全平方数

不妨设 $p=k^2$ ，显然 $k>2$ 。做差易得 $p>2k$ 。
那么有
$$(p-1)!=1\times \cdots \times k\times \cdots \times 2k\times \cdots \times (p-1)=2k^2\times (\dots )$$

为 $p$ 的整数倍，因此 $(p-1)!\equiv 0(\mathrm{m}\mathrm{o}\mathrm{d}p)$ 。

④ $p>4$ 且 $p$ 不是完全平方数

由于 $p$ 不是完全平方数且不是质数，不妨设 $p=a\times b,a<b$ 。
那么有
$$(p-1)!=1\times \cdots \times a\times \cdots \times b\times \cdots \times (p-1)=(a\times b)\times (\dots )$$

为 $p$ 的整数倍，因此 $(p-1)!\equiv 0(\mathrm{m}\mathrm{o}\mathrm{d}p)$ 。

二、例题和代码

(1) HDU2973 YAPTCHA

题目大意：
共 $T$ 组数据，求
$$\sum _{k=1}^n\lfloor \frac{(3k+6)!+1}{(3k+7)}-\lfloor \frac{(3k+6)!}{(3k+7)}\rfloor \rfloor$$

的值（$T,n⩽10^6$）。

题解

威尔逊定理模板题。
令 $3k+7=x$ 得：
$$原式=\sum _{k=1}^n\lfloor \frac{(x-1)!+1}{x}-\lfloor \frac{(x-1)!}{x}\rfloor \rfloor$$

当 $x$ 为质数时，根据威尔逊定理有 $(x-1)!\equiv -1(\mathrm{m}\mathrm{o}\mathrm{d}x)$ ，即 $(x-1)!+1\equiv 0(\mathrm{m}\mathrm{o}\mathrm{d}x)$ ，故 $x\mid (x-1)!+1$ ，因此
$$\frac{(x-1)!+1}{x}-\lfloor \frac{(x-1)!}{x}\rfloor =1$$

当 $x$ 不是质数时，在证明威尔逊定理的必要性时有 $\forall x>4,(x-1)!\equiv 0(\mathrm{m}\mathrm{o}\mathrm{d}x)$ ，因此 $x\mid (x-1)!$ ，因此
$$\lfloor \frac{(x-1)!+1}{x}-\lfloor \frac{(x-1)!}{x}\rfloor \rfloor =0$$

综上，所求的式子可以化简为
$$原式=\sum _{k=1}^n\left[x\text{ is prime}\right]=\sum _{k=1}^n[3k+7\text{ is prime}]$$

使用筛法筛出素数后统计前缀和即可，代码如下。

#include<bits/stdc++.h>
using namespace std;
const int maxn = 3e6 + 10;
bool vist[maxn];
int cnt;
int sum[maxn], prime[maxn];
int T, n;
void Euler(int n){
	vist[1] = true;
	for (int i = 2; i <= n; i ++){
		if (!vist[i]) prime[++ cnt] = i;
		for (int j = 1; j <= cnt && i * prime[j] <= n; j ++){
			vist[i * prime[j]] = true;
			if (prime[j] % i == 0) break;
		}
	}
	for (int i = 1; i < maxn; i ++) sum[i] = sum[i - 1] + (vist[3 * i + 7] == 0);
}
int main(){
	Euler(maxn - 1);
	cin >> T;
	while (T --){
		cin >> n;
		cout << sum[n] << endl;
	}
	return 0;
}

(2) HDU6608

题目大意：
共 $T$ 组数据，给定质数 $p$ ，求小于 $p$ 的最大质数 $q$ 并求下式的值（$T⩽10,10^9⩽p⩽10^{14}$）。
$$q!\mathrm{m}\mathrm{o}\mathrm{d}p$$

题解

由威尔逊定理
$$(p-1)!\equiv -1(\mathrm{m}\mathrm{o}\mathrm{d}p)$$

推导得
$$q!\times (q+1)\times (q+2)\times \cdots \times (p-1)\equiv -1(\mathrm{m}\mathrm{o}\mathrm{d}p)$$
$$q!\equiv -(q+1{)}^{-1}(q+2{)}^{-1}\times \cdots \times (p-1{)}^{-1}(\mathrm{m}\mathrm{o}\mathrm{d}p)$$
因此从 $p-1$ 开始倒序枚举，使用朴素的质数检验方法检验是否为质数并记录逆元的累乘，枚举到质数就输出并停止。不过观察到数据范围比较大，因此需要使用龟速乘，代码如下。

#include<bits/stdc++.h>
using namespace std;
const int maxn = 1e7 + 5;
bool vist[maxn];
int cnt;
int prime[1000005];

void Euler(int n){
	for (int i = 2; i <= n; i ++){
		if (!vist[i]) prime[++ cnt] = i;
		for (int j = 1; j <= cnt && i * prime[j] <= n; j ++){
			vist[i * prime[j]] = true;
			if (i % prime[j] == 0) break;
		}
	}
}

int T;
long long n;

bool is_prime(long long x){
	for (int i = 1; i <= cnt && prime[i] * prime[i] <= x; i ++)
		if (x % prime[i] == 0) return false;
	return true;
}

long long sMul(long long a, long long b, long long MOD){
	long long ans = 0;
	while (b){
		if (b & 1) ans = (ans + a) % MOD;
		a = (a + a) % MOD;
		b >>= 1;
	}
	return ans;
}

long long qPow(long long a, long long b, long long MOD){
	long long ans = 1;
	while (b){
		if (b & 1) ans = sMul(ans, a, MOD);
		a = sMul(a, a, MOD);
		b >>= 1;
	}
	return ans;
}

long long get_inv(long long a, long long p){
	return qPow(a, p - 2, p);
}

int main(){
	std::ios::sync_with_stdio(false);
	cin.tie(0);
	Euler(maxn - 1);
	cin >> T;
	while (T --){
		cin >> n;
		long long ans = n - 1;
		for (long long i = n - 1; ; i --){
			bool ret = is_prime(i);
			if (ret){
				cout << ans << '\n';
				break;
			}
			ans = sMul(ans, get_inv(i, n), n);
		}
	}
	return 0;
}