Miller_Rabin整理笔记

[TOC]

问题

一个数到底是不是素数

别的

首先列一下我们可以求素数的东西根号暴力求 $O(nloglogn)$的埃氏筛 $O(n)$的欧拉筛还有我们要学习的Miller_Rabin算法对了，还有神奇的6倍法(也许叫这个吧)

bool pd(int x) {
    if(x==2||x==3) return 1;
    if(x==1||x%6!=1&&x%6!=5) return 0;
    for(int i=5;i*i<=x;i+=6)
        if(x%i==0||x%(i+2)==0) return 0;
    return 1;
}

正事

先去掉偶合数和2吧，现在剩下奇数了吧费马小定理对素数一定成立，但合数是不一定成立的，多试几次成立几率还是挺

大的，没啥说的 $Carmichael$数本身是合数,而且$1$到$n-1$都满足飞马小定理,suah as:561 着重说一下二次探测吧考虑$x^{2}≡1(mod n)$的根，如果n为奇素数，那跟就只有1和-1两个根(移项之后显然) 设$n-1=2^{r}*d$,d是奇数，如果存在$0≤k＜r，a^{2^{k}*d}≠1，-1(mod n)$ 但是$a^{2^{k+1}d}≡1（mod n）$,那就不满足费马小定理，那么n就是合数了凭借着二次探测和费马小定理，我们出错的几率就大大降低了但只测一次错误还是不够(应该是一次失败的几率是四分之一) 我们用前10个素数，在OI界就可以无忧了复杂度$O(klogn)$ 不知道百度百科搞得什么鬼,好多log

使用快速傅里叶变换能够将这个时间推进到$O(klognlog lognlog log logn) ="O"(klogn).$

代码

#include <bits/stdc++.h>
using namespace std;
typedef long long LL;
LL q_pow(LL a, LL k, LL p) {//calc a^k % p
	LL res = 1;
	for (; k ; k >>= 1, a = a * a % p)
		if (k & 1) res = a * res % p;
	return res;		
}
int prime[] = {2, 3, 5, 7, 11, 13, 17, 19, 23, 29};
bool detective(LL a, LL n) {
	int r = 0; LL d = n - 1; // n - 1 = 2 ^ r * d
	while (d % 2 == 0) d >>= 1, ++r;
	for (LL x = q_pow(a, d, n), y; r ; r--) {
		y = x * x % n;
		if (y == 1) {
			if (x == 1) return true;
			if (x == n - 1) return true;
			return false;
		}
		x = y;
	}
	return false;
}
bool miller_rabin(LL n) {
	for (int i = 0; i < 10; i++) {
		if (n == prime[i]) return true;
		if (n % prime[i] == 0) return false;
		if (!detective(prime[i], n)) return false;
	}
	return true;
}
int main() {
    LL wo,n,x;
    cin>>wo>>n;
    while(n--) {
	cin>>x;
        if (miller_rabin(x)) printf("Yes\n");
        else printf("No\n");
    }
}

来源：oschina

链接：https://my.oschina.net/u/4258874/blog/3655907

标签

Miller