Miller_Rabin质数测试

//Miller_Rabin的复杂度为S * log n
//注意快速幂中的取模n可能会乘法溢出
bool Miller_Rabin(ll n){//其中n为要测试的数字
	if(n <= 2){
		if(n == 2)return true;
		else return false;
	}
	if(n % 2 == 0)return false;
	ll u = n - 1;
	while(!(u & 1))u = u >> 1;
	int S = 1000;//S为测试次数
	ll x,y;
	for(int i = 1;i <= S;i++){
		int a = (rand() % (n - 1)) + 1;//随机获取一个2 ~ (n - 1)的数字
		x = qpow(a,u,n);//x = a ^ u % n,此处用快速幂
		while(u < n){
			y = qpow(x,2,n);//此处用快速幂
			if(y == 1 && x != 1 && x != n - 1)return false;
			x = y;
			u = u * 2;
		}
		if(x != 1)return false;
	}
	return true;
}