Miller-Rabin素性测试

根据费马小定理，若 $n$ 为质数，且 $a<n$，则有 $a^{n-1}\equiv 1(modn)$。

设 $n-1=d\times 2^r,gcd(d,2)=1$，此时若 $n$ 为质数，必有 $a^d\equiv 1(modn)$ 或 $a^{d\times 2^i}\equiv -1(modn),0\le i<r$ 成立。

但这只是 $n$ 为质数的必要不充分条件。

但若 $n$ 对于 $a$ 取 $2,3,5,7,13,29,37,89$ 均成立，则在 $int$ 范围内 $n$ 必为质数；

若对 $a$ 取前 $11$ 个质数均成立，则在 $longlong$ 范围内 $n$ 必为质数。

反之，若 $n$ 为偶数，则不一定均成立。

const int prime[8] = {2, 3, 5, 7, 13, 29, 37, 89};

inline int qpow(int a, int b, int p) {
	int ans = 1, x = a % p;
	while (b) {
		if (b & 1) ans = 1ll * ans * x % p;
		x = 1ll * x * x % p, b >>= 1;
	}
	return ans;
}

inline int miller_rabin(int x) {
	if (x <= 1) return 0;
	for (int i = 0; i < 8; ++i)
		if (x == prime[i]) return 1;
	for (int i = 0; i < 8; ++i) {
		int d = x - 1, r = 0;
		while (d % 2 == 0) d /= 2, ++r;
		int a = qpow(prime[i], d, x), flag = 0;
		if (a == 1) continue;
		for (int j = 0; j < r; ++j) {
			if (a == x - 1) {
				flag = 1;
				break;
			}
			a = 1ll * a * a % x;
		}
		if (!flag) return 0;
	}
	return 1;
}