MillarRabin模版随机测大素数

最新推荐文章于 2022-06-12 10:25:46 发布

原创

最新推荐文章于 2022-06-12 10:25:46 发布 · 1.2k 阅读

1 ·

CC 4.0 BY-SA版权

本文介绍如何利用Millar-Rabin算法模板高效地测试和验证长整型数值的素数属性，特别是针对大素数的判断。

支持long long

#include <algorithm>
typedef long long ll;
ll GCD(ll a, ll b) { return b ? GCD(b, a % b) : a; }
ll MultiMod(ll a, ll b, ll n) { // a * b % n
	ll res = 0; MillarRabin
	a %= n;
	while (b > 0) {
		if (b & 1) {
			res += a;
			if (res >= n) res -= n;
		}
		a <<= 1;
		if (a >= n) a -= n;
		b >>= 1;
	}
	return res;
}
ll QuickMod(ll a, ll b, ll n) { // a ^ b % n
	ll res = 1;
	a %= n;
	while (b > 0) {
		if (b & 1) res = MultiMod(res, a, n);
		a = MultiMod(a, a, n); b >>= 1;
	}
	return res;
}
bool MillarRabin(ll n) { // 判断是否素数
	if (n == 2 || n == 3 || n == 5 || n == 7 || n == 11) return true;
	if (n == 1 || !(n & 1) || !(n % 3) || !(n % 5) || !(n % 7) || !(n % 11)) return false;
	ll t = 0, m = n - 1, x, y;
	while (!(m & 1)) { m >>= 1; t++; }
	for (int i = 0; i < 10; i++) {
		ll a = rand() % (n - 2) + 2;
		x = QuickMod(a, m, n);
		for (ll j = 0; j < t; j++) {