【ACM】求素数的多种方法

素数：质数（prime number）又称素数，有无限个。质数定义为在大于1的自然数中，除了1和它本身以外不再有其他因数的数称为质数。

下面给出求素数的几种不同的算法。
假设我们给定一位数N，现在判断这个数字是不是素数。

方法一：
从2开始遍历到N-1 。

bool isprime(int x)
{
        if(x == 1)
            return false;
        if(x>=2)
        {
             for(int i=2;i<x;i++)
             {
                  if(x % i == 0)
                      return false;
             }
             return true;
        }
}

这是一种非常慢的方法。

方法二：
对上述的方案进行稍微的优化

bool isprime(int x)
{
        if(x == 1)
            return false;
        if(x>=2)
        {
             for(int i=2;i<sqrt(x);i++)
             {
                  if(x % i == 0)
                      return false;
             }
             return true;
        }
}

因为一个数N = a*b，则只需找出a，b中的较小值就可以，而较小值的最大上限不会超过sqrt（N）。

方法三：
遍历从2到N的平方根下每一个素数，但是如果N很大，素数表的构建就会非常复杂。因此我们可以选择构建动态素数表，但是效率依然会很低。
1:试除法：试除sqrt(N)以内的素数
2:筛选法：（下面给出两个链接）
Sieve of Eratosthenes
Euler‘s Sieve

方法四：
费马小定理 + 快速幂运算

对于两个互质（两个正整数只有公约数1时）的素数n和p，必有n^(p-1) % p =1 。考虑到在运算高次幂运算时候可能会发生上溢，所以采用快速幂取模算法。

int Montgomery(int n,int p,int m)
{
        int ans= 1;
        n= n%m;
        while(p!=1)
        {
               if(p&1)
                   ans  = (ans*n)%m;
                n =(n*n)%m;
                p>>=1;
        }
        return ans;
}

以上是快速幂取模的算法

bool isprime(int x)
{
      if(x == 1)
          return false;
        int PrimeList[]={   2, 3, 5, 7, 11, 17, 19, 23, 29, 31, 41,
          43, 47, 53, 59, 67, 71, 73, 79, 83, 89, 97};
          const int len = sizeof(PrimeList)/sizeof(int);
          for(int i=0;i<len;i++)
          {
               if(Montgomery(PrimeList[i] , x-1 ,x) != 1)
               //如果不为一，则这两个数字不互质
               return false;
          }
          return true;
}

但是例如数字29341，等于13* 37*61，显然是一个合数，但是被判断为了素数，所以这种策略是一种不可靠的策略。这一类数字叫做卡米歇尔数。卡米歇尔数和素数非常相似，都满足了费马小定理，因此称之为伪素数。随着数越来越大，卡米歇尔数变得越来越少，1至10^17有585 355个卡米歇尔数。

方法五：
拉宾米勒测试（目前未知最快的判断素数的方法，没有之一）