素数判定方法_Miller-Rabbin测试

原文：http://blog.csdn.net/sluqy671/article/details/41701655

首先了解一下伪素数的概念： 如果n是一个正整数，如果存在和n互素的正整数a满足a^(n-1)≡1(mod n)，我们说n是基于a的伪素数。

如果一个数是伪素数，它几乎就是素数。另一方面，如果一个数不是伪素数，它一定不是一个素数。那么在一定的条件下，如果我们选取了若干个基都发现n是伪素数， 那么n是素数的概率趋近于1。

Miller-Rabbin测试

现在我们只需要多次寻找不超过n-1基并检验是否有a^(n-1)≡1(mod n)， 如果一直有， 那么这个数就是一个素数， 否则可以立即判定这个数是个合数。

检验的时候使用快速幂进行优化， 可以使复杂度降至O(slogn)， 这里s为选取基的次数。 可以证明， 出错的概率不大于1/(2^s)， 因此s取到几十就差不多了。 再就是对于2要进行特判， 否则会对0求余， 那显然是不对的。

[html]  view plain  copy

1. #include <cstdio>  
2. #include <cmath>  
3. #include <cstdlib>  
4. #include <ctime>  
5.   
6. using namespace std;  
7.   
8. typedef long long LL;  
9. LL gcd(LL x, LL y)  
10. {  
11.     if (!y) return x;  
12.     return (y, x%y);  
13. }  
14. LL pow(LL a, LL x, LL mod)  
15. {  
16.     LL ans = 1;  
17.     while(x)  
18.     {  
19.         if (x & 1) (ans *= a) %= mod;  
20.         (a *= a) %= mod;  
21.         x >>= 1;  
22.     }  
23.     return ans;  
24. }  
25. bool MRT(LL x)  
26. {  
27.     if (x == 2) return true;  
28.     for (LL i = 1; i <= 30; ++i)  
29.     {  
30.         LL now = rand()%(x-2) + 2;  
31.         if (pow(now, x-1, x) != 1) return false;  
32.     }  
33.     return true;  
34. }  
35. int main()  
36. {  
37.     int n;  
38.     LL x;  
39.     scanf("%d", &n);  
40.     while(n--)  
41.     {  
42.         scanf("%I64d", &x);  
43.         if (MRT(x)) printf("YES\n");  
44.         else printf("NO\n");  
45.     }  
46.     return 0;  
47. }