//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;
}