51NOD1106 质数检测

基准时间限制：1 秒 空间限制：131072 KB 分值: 0 难度：基础题

 收藏

 关注

给出N个正整数，检测每个数是否为质数。如果是，输出"Yes"，否则输出"No"。

Input

第1行：一个数N，表示正整数的数量。(1 <= N <= 1000)
第2 - N + 1行：每行1个数(2 <= S[i] <= 10^9)

Output

输出共N行，每行为 Yes 或 No。

Input示例

5
2
3
4
5
6

Output示例

Yes
Yes
No
Yes
No

/*
*模板1
*Miller-Rabin算法是目前主流的基于概率的素数测试算法
*快速测试n是否满足素数的“必要”条件，出错概率极低 
*/ 
int modularExponent(int a, int b, int n) {
	int ret = 1;
	for (; b; b >>= 1, a = (int) ((long long) a * a % n)) {
		if (b & 1) {
			ret = (int) ((long long) ret * a % n);
		}
	}
	return ret;
}
bool millerRabin(int n,int a) {
	if (n == 1 || (n != 2 && !(n % 2)) || (n != 3 && !(n % 3)) || (n != 5 && !(n % 5)) || (n != 7 && !(n % 7))) {
		return false;
	}
	int r = 0, s = n - 1, j;
	if(!(n%a)) return false;
	while(!(s&1)){ s >>= 1; r++; }
	long long k = modularExponent(a, s, n);
	if(k == 1) return true;
	for(j = 0; j < r; j++, k = k * k % n)
	    if(k == n - 1) return true;
	return false;
}
bool miller_Rabin(int n)//是素数返回真值 
{
	int a[]={2,3,5,7},i;//能通过测试的最小素数为 3215031751(此数超int)
	for(i=0;i<4;i++){
		if(!millerRabin(n,a[i]))return false;
	}
	return true;
}